Question: paCalls with Clariom Chips
gravatar for kmvaldez2
2.4 years ago by
United States / Bethesda, MD / National Institutes of Health
kmvaldez20 wrote:

Good Afternoon,

I am trying to do present absent calls on my Clariom chips with oligo's paCalls function, but each time I receive a subscript out of bounds error.  I've tried the Human S HT and Mouse S Clariom chips.  Does anyone know of a work around?






pd<-read.AnnotatedDataFrame('Phenotype.txt', row.names = 'SampleName', header = T, sep = '\t')

celfiles <- read.celfiles(filenames = c('sample1.CEL', 'sample2.CEL', 'sample3.CEL', 'sample4.CEL', 'sample5.CEL', 'sample6.CEL'), phenoData=pd)

 PApset <- paCalls(celfiles,"MAS5")



Getting probe level data... OK.
Error in mms[o, , drop = FALSE] : subscript out of bounds


microarray oligo clariom • 400 views
ADD COMMENTlink modified 2.4 years ago by James W. MacDonald52k • written 2.4 years ago by kmvaldez20
Answer: paCalls with Clariom Chips
gravatar for James W. MacDonald
2.4 years ago by
United States
James W. MacDonald52k wrote:

You can't do a MAS5 style P/A call with an array that has no MM probes.

ADD COMMENTlink written 2.4 years ago by James W. MacDonald52k

As a followup, you should just get an error saying that paCalls doesn't work with this particular array. The ExpressionFeatureSet was originally intended for the old style 3'-biased arrays, which mostly did have MM probes (although even back then Affy was phasing them out).

When I added the Clariom S arrays to pdInfoBuilder, I used this class, because in most respects the Clariom S arrays are very similar to the old 3'-biased arrays. I'll add a more informative error message.

ADD REPLYlink written 2.4 years ago by James W. MacDonald52k
Please log in to add an answer.


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