Description Usage Arguments Details Value See Also Examples
This function caculates the weighted contrast W_j in order to find out the most significant dfferentially expressed genes.
1 | scores.wire(obj,contrast=NULL)
|
obj |
The return list of function emmixwire |
contrast |
The vector of the specified contrast |
The number of classes of samples is either two or three.
The vector of the statistic Wj
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 | ## Not run:
dat <- read.table("GSE36703_37628_col.txt",header=FALSE,sep='\t')
rownames(dat) <- 1:nrow(dat)
###normalize the rows
x <- DoRows(dat)
set.seed(12345)
ret <-emmixwire(x,g=3,ncov=3,nvcov=1,n1=5,n2=6,n3=3,
debug=1,itmax=1000,epsilon=1e-5)
###calculate the W_j
wj <- scores.wire(ret,contrast=c(0.5,0.5,-1))
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.