eq8.wire: Calculate the lambda values

In andrewthomasjones/EMMIXcontrasts: Contrasts in Mixed Effects for EMMIX Model with Random Effects

Description

This function calculates the lamda values in equation 2.8.

Usage

eq8.wire(m, g, nb, X, W, U, V, sigma.e, sigma.b, sigma.c, nh, contrast)

Arguments

m	The number of columns in the data
g	The number of components in the mixture model
nb	The number of columns in the design matrix X
X	The design matrix X
W	The design matrix W
U	The design matrix U
V	The design matrix V
sigma.e	The covariance parameters of error item e
sigma.b	The covariance parameters of random item b
sigma.c	The covariance parameters of random item c
nh	A vector with each element as the number of genes in each cluster
contrast	The contrast vector, for example, c(1/2, 1/2, -1) or c(1, 0, -1), etc.

Details

The default contrast for two classes is c(1, -1), and three classes c(1/2, 1/2, -1). And it is called by function scores.wire and wj2.permuted.

Value

A vector of lamda values.

See Also

wj2.permuted scores.wire

andrewthomasjones/EMMIXcontrasts documentation built on June 26, 2019, 5:46 a.m.