MAIL_Split: MAIL_Split

In hwyneken/MAILPackage: Model-Averaged Inferential Learning

MAIL_Split

Description

MAIL_Split runs MAIL with data splitting

Usage

MAIL_Split(XMat, yVec)

Arguments

XMat	a n by p numeric matrix
yVec	a n by 1 numeric vector

Details

This is a specific use of the MAIL function. The function uses the following arguments with MAIL:

• splitOption = "Split"

• firstSOILWeightType = "BIC"

• smallestModelWeightType = "AIC"

• firstSOILPsi = 0.5

• smallestModelPsi = 0 sigma2EstFunc = "LPM_AIC_CV_50Split" verbose = FALSE

Value

resList a list

See Also

MAIL and MAIL_Full

Examples

#' two main ideas to explain with an example
 
 how to run MAIL_Split

 how to run diagnostics for MAIL_Split - is it a reliable method for the current data


#### how to run MAIL_Split

# first step - organize the data
# gene expression problem: response is the expression of gene probe X249.
# install package colonCA through Bioconductor (use these steps)
# https://bioconductor.org/packages/release/data/experiment/html/colonCA.html
require(colonCA)
data(colonCA)
dataSet <- t(colonCA@assayData$exprs)
colon_y <- dataSet[,249]
colon_x <- dataSet[,-1*249]
colon_x <- t(unique(t(colon_x))) # remove duplicate columns - original has p = 1999, reduce to p = 1990
colon_y <- scale(colon_y)
colon_x <- scale(colon_x)

colonMAILSplit <- MAIL_Split(XMat = colon_x,yVec = colon_y)

hwyneken/MAILPackage documentation built on July 27, 2022, 12:46 a.m.