LqG (version 0.1.0)

grsc.marg.MLqE: Group Screening based on marginal Maximum Lq-likelihood Estimation

Description

Group screening by ranking utility of each group. The group effect is defined based on the cumulation of the maximum Lq-likelihood estimate of the regression using only one predictor each time within the group.

Usage

grsc.marg.MLqE(
X,
Y,
n = dim(X)[1],
p = dim(X)[2],
q = 0.9,
m,
group,
eps = 1e-06,
d = n/log(n)
)

Arguments

X

A matrix of predictors.

Y

A vector of response.

n

A value of sample size

p

A value denoting the dimension of predictors

q

A value of distortion parameter of Lq function, default to 0.9.

m

A number of the predictor groups

group

A vector of consecutive integers describing the grouping of the coefficients (see example below).

eps

The iteration coverage criterion, default to 1e-06.

d

A value of the number of groups retained after screening, default to n/log(n)

Value

The grsc.marg.MLqE returns a list containing the following components:

beta.group

The vector of utility of each group, which is the criterion for the variable screening procedure.

group.screened

The vector of integers denoting the screened groups.

Details

grsc.marg.MLqE obtains the group effect of each group for subsequential group screening, based on the cumulative marginal MLqE coefficients within the group. It can work when both the correlation within groups and between groups are small. If group size equals to 1, individual screening is conducted.

Examples

# NOT RUN {
# This is an example of grsc.marg.MLqE with simulated data
data(LqG_SimuData)
X = LqG_SimuData$X
Y = LqG_SimuData$Y
n = dim(X)[1]
p = dim(X)[2]
m = 200
groups = rep(1:(p/5), each = 5)
result <- grsc.marg.MLqE(X = X,
                         Y = Y,
                         n = n,
                         p = p,
                         q = 0.9,
                         m = m,
                         group = groups,
                         eps = 1e-06,
                         d = 15)
result$beta.group
result$group.screened
# }