Question: SAM r for finding up and down regulated genes
0
gravatar for Linda Chaba
3.7 years ago by
Linda Chaba10
Kenya
Linda Chaba10 wrote:

I'm using the below code for samr to find down and up regulated genes


source("http://bioconductor.org/biocLite.R")

biocLite(samr)

data=celline.data #expression data imported from excel file whose columns are samples and rows are genes. First row has quantitative outcome.
x=data[-1,] #expression data only without quantitative outcome
y=data[1,] #Quantitative outcome

library(samr)

sam.data<-list(x=x,y=y, geneid= 1:nrow(data), genenames= row.names(data))  # prepare the data list required in “samr” #function
samr.obj<-samr(sam.data,  resp.type="Quantitative", nperms=100, random.seed=12345) 

The last line is producing this error 

Error in storage.mode(x) <- "double" : 
  (list) object cannot be coerced to type 'double'
In addition: Warning message:
In knnimp(x, k, maxmiss = rowmax, maxp = maxp) :
  2950 rows with more than 50 % entries missing;
 mean imputation used for these rows
Error during wrapup: cannot open the connection


I'll appreciate if someone can point out for me what I'm not doing right.

Thank you.

samr • 957 views
ADD COMMENTlink modified 3.7 years ago by Valerie Obenchain6.7k • written 3.7 years ago by Linda Chaba10
Answer: SAM r for finding up and down regulated genes
0
gravatar for Valerie Obenchain
3.7 years ago by
United States
Valerie Obenchain6.7k wrote:

samr is not a Bioconductor package so the authors probably aren't watching this mailing list. For help with CRAN packages it's better to post to the R-help mailing list

https://www.r-project.org/mail.html

Valerie

ADD COMMENTlink written 3.7 years ago by Valerie Obenchain6.7k

Thank you Valerie for pointing that out.

 

ADD REPLYlink written 3.7 years ago by Linda Chaba10
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 192 users visited in the last hour