make.dmmobj {dmm}R Documentation

Construct an object of class dmm from user-supplied data

Description

Construct an object of class dmm containing all attributes needed to run the gresponse() function. The user must supply a phenotypic covariance matrix, and a genetic covariance matrix for each of the components needed by gresponse() Warning; this function does not currently support class specific genetic parameters.

Usage

make.dmmobj(p = NULL, components = c("VarG(Ia)"), ...)

Arguments

p

A phenotypic covariance matrix. Dimnames for rows and columns must be set to the trait names

components

A character vector specifying names of all of the genetic variance components for which a genetic covariance matrix is to be supplied. If there are genetic covariance components ( eg "CovG(Ia,Ma)"), these must be listed after all the genetic variance components

...

A variable number of genetic covariance matrices, one for each of the names listed in the components argument. Each matrix must have dimnames the same as the phenotypic covariance matrix

Details

Only the minimal attributes for use by the gresponse() function are constructed. The remainder are set to NULL

Value

An object of class dmm. Only the following attributes contain data

b

A dummy set of fixed effects

siga

A matrix of genetic variance (and covariance) components set out one component per row and the traitpairs in columns

variance.components

A matrix of genetic variance (and covariance) components with the phenotypic covariance matrix appended

phenotypic.variance

A matrix of phenotypic covariances set out one trait per row and per column

Note

If matrices are supplied which are not positive definite, function make.dmmobj() will alter them to the nearest positive definite matrix using routine neadPD() from package Matrix. No message is given. Check the value returned to see if matrices have been adjusted.

Author(s)

Neville Jackson

See Also

Functions gresponse() and dmm().

Examples

library(dmm)
p <- matrix(c(3,2,2,4),2,2)
dimnames(p) <- list(c("T1","T2"),c("T1","T2"))
gia <- matrix(c(2,1,1,3),2,2)
dimnames(gia) <- dimnames(p)
myobj <- make.dmmobj(p,components=c("VarG(Ia)"),gia)
myresp <- gresponse(myobj,psd=list(dp=c(0.5,0.1)))
print(myresp)
#cleanup
rm(p)
rm(gia)
rm(myobj)
rm(myresp)


[Package dmm version 2.1-9 Index]