eq8.wire: Calculate the lamda values

In EMMIXcontrasts2: Contrasts in Mixed Effects for EMMIX Model with Random Effects 2

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.

See Also

wj2.permuted

EMMIXcontrasts2 documentation built on May 30, 2017, 3 a.m.