MLqE.est | R Documentation |
The iterative algorithm for MLqE of coefficients of regression using each group of variables.
MLqE.est( X, Y, q = 0.9, eps = 1e-06 )
X |
The matrix of the predictor group. |
Y |
The vector of response. |
q |
The value of distortion parameter of Lq function, default to |
eps |
The iteration coverage criterion, default to |
The estimating equation of MLqE is a weighted version of that of the classical maximum likelihood estimation (MLE) where the distortion parameter q determines the similarity between the Lq function and the log function. When q = 1, MLqE is equivalent to MLE. The closer q is to 1, the more sensitive the MLqE is to outliers. As for the selection of q, there is presently no general method. However, MLqE is generally less sensitive to data contamination than MLE (to different degrees) when q is smaller than 1. Here, the default value of q is 0.9. Distortion parameter q can also be determined according to sample size n, choices of q_n with |1-q_n| between \frac{1}{n} and \frac{1}{√{n}} usually improves over the MLE.
The MLqE.est
returns a list containing the following components:
t |
The integer specifying the number of the total iterations in the algorithm. |
beta_hat |
The vector of estimated coefficients. |
sigma_hat |
The value of the estimated variance. |
OMEGA_hat |
The matrix of the estimated weight. |
# 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:( dim(X)[2] / 5), each = 5) Xb = X[ , which( groups == 1)] result = MLqE.est(Xb, Y, q = 0.9, eps = 1e-06) result$beta_hat result$sigma_hat result$OMEGA_hat result$t
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.