scores.wire: The weighted contrast

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

Description

This function caculates the weighted contrast W_j in order to find out the most significant dfferentially expressed genes.

Usage

scores.wire(obj, contrast = NULL, useZ = TRUE)

Arguments

obj	The return list of function emmixwire.
contrast	The vector of the specified contrast.
useZ	use the latent variable allocation Z_i (default, TRUE) or the posterior probability (FALSE).

Details

The number of classes of samples is either two or three.

Value

The vector of the statistic Wj

See Also

wire.init.fit scores.wire

Examples

data(hedenlc)
hedenlc<-hedenlc[seq_len(100),] #for speed
set.seed(123456)
obj<-emmixwire(hedenlc, g=5, ncov=3, nvcov=1, n1=7, n2=8, 
     debug=1, itmax=20, epsilon=1e-4) 
# calculate the weighted contrasts W_j
Wj  <- scores.wire(obj)

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