Initial {hySAINT}

Creating initial parents

Description

This function gives initial parents.

Usage

Initial(X, y, EVAoutput, heredity = "Strong", r1, r2, numElite = 40)

Arguments

X

Input data. An optional data frame, or numeric matrix of dimension n observations by p main effects.

y

Response variable. A n-dimensional vector.

EVAoutput

The output from function EVA

heredity

Whether to enforce Strong, Weak, or No heredity. Default is "Strong".

r1

At most how many main effects do you want to include in your model?. For high-dimensional data, r1 cannot be larger than the number of screened main effects.

r2

At most how many interaction effects do you want to include in your model?

numElite

Number of elite parents. Default is 40.

Value

Initial parents. A numeric matrix with dimensions numElite by r1+r2.

See Also

EVA

Examples

set.seed(0)
interaction.ind <- t(combn(10,2))
X <- matrix(rnorm(100*10,1,0.1), 100, 10)
epl <- rnorm(100,0,0.01)
y <- 1+X[,1]+X[,2]+X[,3]+X[,1]*X[,2]+X[,1]*X[,3]+epl
EVAoutput <- EVA(X, y, r1 = 5, sigma = 0.01,
  interaction.ind = interaction.ind)
myParent <- Initial(X = X, y = y, EVAoutput, r1 = 5, r2 = 2)