grmsem.stpar {grmsem}

grmsem standardised parameters function

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)