calc: helper function to be called in 'update.alpha.coef'.

In matthiasspeidel/hmi: Hierarchical Multiple Imputation

Description

It (basically) calculates (Z'Z + sigma.e * sigma.b^-1)^-1 * z' resid.

It (basically) calculates (Z'Z + sigma.e * sigma.b^-1)^-1 * z' resid.

Usage

calc(x, sigma.y.new, sigma.alpha.new.inv)

calc(x, sigma.y.new, sigma.alpha.new.inv)

Arguments

x	A data.frame which is basically cbind(Z_obs, y_obs - X_obs * beta.new).
sigma.y.new	The numeric residual variance.
sigma.alpha.new.inv	The inverse random effects covariance matrix.
x	A data.frame which is basically cbind(Z_obs, y_obs - X_obs * beta.new).
sigma.y.new	The numeric residual variance.
sigma.alpha.new.inv	The inverse random effects covariance matrix.

Value

It returns (for this given cluster) a (cluster specific) random intercept and random slope.

It returns (for this given cluster) a (cluster specific) random intercept and random slope.

matthiasspeidel/hmi documentation built on Aug. 18, 2020, 4:37 p.m.