Displaying 1 result from an estimated 1 matches for "genenames1000".
2005 Dec 26
0
problem with samr
Hello Everybody,
I am trying to perform SAM with the samr package.
I am using the following code:
sink ("R005")
library(siggenes)
library(samr)
library(nnet)
A <- as.matrix(read.table("D:\samrgenes1000.txt"))
B <- as.matrix(read.table("D:\genenames1000.txt"))
y1 <- c(rep(1,20),rep(2,6)) #there are 20 chips of one kind and 6 of the
other kind.
datasam = list(x=A,y=y1,genenames=B,logged2=TRUE)
testsamr <-samr(datasam,resp.type ="Two class unpaired",nperms=100)
del <- 2
samr.plot(testsamr,delta)
delta.table <- samr.comput...