GEInfo (version 1.0)

TypeS: Construct Type_S prior information

Description

For G variables and G-E interactions, transform their prior information from counts(frequencies) into a set of significant variables (Type_S)

Usage

TypeS(
  G.count,
  GE.count,
  eta_G = 0.95,
  eta_GE = 0.95,
  varphi_G = NULL,
  varphi_GE = NULL
)

Arguments

G.count

A numeric vector, including the prior counts (frequencies) for G variables.

GE.count

A numeric matrix, including the prior counts (frequencies) for G-E interactions.

eta_G

A probability. The (eta_G)th quantile of G.count is used as a count (frequency) threshold (denoted by varphi_G) for G variables. Default is 0.95.

eta_GE

A probability. The (eta_GE)th quantile of GE.count is used as a data-dependent count (frequency) threshold (denoted by varphi_GE) for G-E interactions. Default is 0.95.

varphi_G

A user supplied count threshold for G variables. It is used to determine which G variables will be finally included in the Type_S prior information set. Default is NULL. Typical usage is to have the program calculate the (eta_G)th quantile of G.count as the threshold. Supplying a varphi_G value will override this.

varphi_GE

A user supplied threshold value used for G-E interactions. It is used to determine which G-E interactions will be finally included in the Type_S prior information set. Default is NULL. Typical usage is to have the program calculate the (eta_GE)th quantile of GE.count as the threshold. Supplying a varphi_GE value will override this.

Value

The outputs include the Type_S prior information sets for G variable and G-E interactions.

S_G

A numeric vector, denoting the Type_S set for G variables. For j in S_G, the jth G variable is suggested to be associated with the response.

S_GE

A numeric matrix, denoting the Type_S set for G-E interactions. For (l,k) in S_GE,the lth G variable and the kth E variable is suggested to have an interaction effect on the response.

Examples

# NOT RUN {
G.count<-c(100,300)
GE.count<-matrix(c(130,356,8,30,87,2),nrow=2)
TypeS(G.count,GE.count)
# }