grmsem.biher: grmsem bivariate heritability estimation function.
In grmsem: Genetic-Relationship-Matrix Structural Equation Modelling (GRMSEM)

Description Usage Arguments Details Value Examples

This function estimates the bivariate heritability.

1	grmsem.biher(ph, grmsem.var.out = NULL)

`ph`	Phenotype file as R dataframe (columns: >=2 phenotypes, rows: ni individuals in the same order as G). No default.
`grmsem.var.out`	A grmsem.var object with unstandardised parameters (factor loadings). Default NULL.

The grmsem.biher function estimates the bivariate heritability (DS, Cholesky, IP and IPC models) from the observed phenotype data and a grmsem.var object. All standard errors are derived with the Delta method.

grmsem.biher returns a list object consisting of the following matrices:

`VPO`	observed phenotypic variance/covariance matrix
`VA`	estimated genetic variance
`BIHER`	estimated bivariate heritability (off-diagonals): VA / VPO
`BIHER.se`	standard error of estimated bivariate heritability
`BIHER.Z`	Z (wald) of estimated bivariate heritability
`BIHER.p`	p (Wald) of estimated bivariate heritability

#(runtime should be less than one minute)

out <- grmsem.fit(ph.small, G.small, LogL = TRUE, estSE = TRUE)
var.out <- grmsem.var(out)
grmsem.biher(ph.small, var.out)