Problem using RankProd
1
0
Entering edit mode
icxa0409 • 0
@icxa0409-8933
Last seen 4.9 years ago
United States

I am trying use the rankprod package but i am unable to do so. Its like initially I had .txt files and I think RankProd does not takes .txt format. Please tell me if it does and how shall I proceed? Also if I take data in .cel format then I do not have the sample-description.csv file. So what shall I do?

rankprod • 668 views
ADD COMMENT
1
Entering edit mode
@james-w-macdonald-5106
Last seen 2 hours ago
United States

If you have celfiles, then you should use either the affy or oligo or xps package to normalize and summarize your data. You then need a 'classlabel' vector of 0s and 1s, indicating which group a given sample comes from. RankProd won't accept an ExpressionSet, so you have to extract it. As an example, say you have six samples, three of one type and three of another. Also assume these are Affy Gene 2.0 ST arrays. This is how I would analyze the data.

library(oligo)
library(RankProd)
dat <- read.celfiles(filenames = list.celfiles())
eset <- rma(dat)
cl <- rep(0:1, each  = 3) ## assuming that the first three and last three are different
results <- RP(exprs(eset), cl)
ADD COMMENT

Login before adding your answer.

Traffic: 342 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