GCTAREML | R Documentation |
This is the reml for the high dimensional linear model assuming known sigma_Y^2=1
GCTAREML(y, x, alpha = c(0.05), lam = 1, niter = 100, eps = 1e-06)
y |
outcome: a vector of length n. |
x |
covariates: a matrix of nxp dimension. |
alpha |
a vector of type I errors used to generate (1-alpha)confidence intervals. |
lam |
initial value |
niter |
the number of iterations for finding the signal noise ratio |
eps |
the convergence criterion for the iteration |
This method assume the independent covariates with fixed effects but can be equivalently treated as random effects
Estimate of proportion of the explained variation, variance estimates, and the corresponding confidence intervals.
Restricted maximum likelihood estimator of the random effects model
## Not run: GCTAREML(y,x)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.