Question: Problem using RankProd
0
gravatar for icxa0409
4.2 years ago by
icxa04090
United States
icxa04090 wrote:

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 • 556 views
ADD COMMENTlink modified 4.2 years ago by James W. MacDonald52k • written 4.2 years ago by icxa04090
Answer: Problem using RankProd
1
gravatar for James W. MacDonald
4.2 years ago by
United States
James W. MacDonald52k wrote:

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 COMMENTlink written 4.2 years ago by James W. MacDonald52k
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: 286 users visited in the last hour