sureEst: SURE model parameter estimation

In qtlmt: Tools for Mapping Multiple Complex Traits

Description

Estimate parameters in a SURE model.

Usage

sureEst(y, x, v, sigma, iter=10000, tol=1e-12)

Arguments

y	a n by p matrix, whose columns are dependent variables.
x	a n by m matrix, whose columns are predictor variables to select from.
v	a list; v[[j]] indicates which x's in the model for y[,j].
sigma	initial residual variance-covarance matrix (if given).
iter	maximum number of iterations in a numerical process to estimate model parameters.
tol	convergence tolerance.

Value

a list with the following components:

loglik	log-likelihood of the model
b	estimates of model coefficients
sigma	estimated residual variance-covariance
fitted.values	fitted mean values

Examples

data(etrait)
x<- as.matrix(mdat-1/2)
y<- as.matrix(traits)[,1:3]
v<- list(c(1,25,50),numeric(0),3)
## Not run: 
o<- sureEst(y, x, v=v, iter=250, tol=1e-12)

## End(Not run)

qtlmt documentation built on May 2, 2019, 2:23 p.m.