Fwd: error with Two class paired in samr function
0
0
Entering edit mode
@john-antonydas-gaspar-3144
Last seen 9.6 years ago
Dear friends, Kindly help me out to trouble shoot. I am interested to find out Differentially Expressed genes betwenn control group(triplicates) and treated group (triplicates). I am trying with the followin coding using samr but ended up with error. ----- library(samr) Data1<-read.table("normalized.txt",sep="t",header=TRUE,dec = ".",as.is =TRUE,na.strings = "NA", colClasses = NA,check.names = FALSE,strip.white = FALSE, blank.lines.skip = TRUE, allowEscapes = FALSE, flush = FALSE,encoding = "unknown") dim(Data1) [1] 27046 7 > de=Data1[,2:7] > de.gnames=Data1[,1] > de.cl=c(-1,-1,-1,1,1,1) > d=list(x=de,y=de.cl,geneid=as.character(de.gnames),genenames=de.gnames ,logged2=FALSE) > samr.obj<-samr(d,resp.type="Two class paired", nperms=100, return.x=FALSE, knn.neighbors = 50, random.seed=NULL) Error in o1[j] <- (1:ncol(x))[y == -o[j]] : replacement has length zero In addition: Warning message: In o1[j] <- (1:ncol(x))[y == -o[j]] : number of items to replace is not a multiple of replacement length ---------------------------------------- Herewith I have attached the normalixed data set also. I would greately appreciate the help With kind regards, antony -- John Antonydas Gaspar, Phd Student AG: Prof.A.Sachinidis Institute of Neurophysiology University of Cologne Robert-Koch-Str. 39 50931 Cologne/Germany Tel: 004922125918042 Handy: 004917683142627 ----- End forwarded message ----- ----- End forwarded message ----- -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: normalized.txt URL: <https: stat.ethz.ch="" pipermail="" bioconductor="" attachments="" 20090810="" d5fa5f9f="" attachment.txt="">
• 849 views
ADD COMMENT

Login before adding your answer.

Traffic: 891 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6