grmsem.stpar: grmsem standardised parameters function

In grmsem: Genetic-Relationship-Matrix Structural Equation Modelling (GRMSEM)

Description

This function estimates standardised parameters for a grmsem.fit object.

Usage

grmsem.stpar(grmsem.out = NULL)

Arguments

grmsem.out

grmsem.fit object as provided by the grmsem.fit function. Default NULL.

Details

grmsem.stpar standardises grmsem.fit estimates so that derived or estimated A and E variance components will add up to phenotypic unit variance. The SEs of standardised parameters are derived using the Jacobian matrix.

Value

grmsem.stpar returns a list object consisting of:

model.in	list of input parameters
stand.model.out	dataframe of fitted grmsem model with standardised parameters and SEs
stVCOV	standardised variance/covariance matrix
k	number of phenotypes
n	total number of observations across all phenotypes
model	type of grmsem model
ph.nms	vector of phenotype names

model.in list of input parameters:

part	a - genetic, e - residual parameters
label	parameter label
value	starting values
freepar	free parameters

stand.model.out data.frame of fitted grmsem model with standardised parameters and SEs:

label	parameter label
estimates	standardised estimated parameters
se	SE
Z	Z (Wald)
p	p (Wald)

Examples

#(runtime should be less than one minute)

out <- grmsem.fit(ph.small, G.small, LogL = TRUE, estSE = TRUE)
stout <- grmsem.stpar(out)
print(stout)

grmsem documentation built on Jan. 29, 2021, 5:07 p.m.