[BioC] affyPLM and exon array question

[Crispin Miller]

Hi Allen,
Does rma() work with your cdf? 

We've also produced one that works OK with rma() (see the 'exonmap'
package vignette for more details, including how to get it).  Don't know
if that helps?

Crispin


 

> -----Original Message-----
> From: bioconductor-bounces at stat.math.ethz.ch 
> [mailto:bioconductor-bounces at stat.math.ethz.ch] On Behalf Of Allen Day
> Sent: 01 May 2007 01:32
> To: bioconductor at stat.math.ethz.ch
> Subject: [BioC] affyPLM and exon array question
> 
> Hi,
> 
> I've been trying to get NUSE, RLE, and RMA values for 
> HuEx-1_0-st-v2 (Human "all exon") Affymetrix arrays.
> 
> So far I have successfully read the arrays into an affybatch object.
> This required creating the CDF environment, which I have 
> already done with makecdfenv.  I'll be submitting that for 
> inclusion shortly, but that's another topic.
> 
> After creating the AffyBatch, I try to use affyPLM to do an 
> RMA model fit.  R = 2.4.1, affyPLM = 1.12.0, affy = 1.12.2.  
> That's where there's trouble, and it appears to be caused by 
> the lack of mismatch probes on the array.  Here's code 
> illustrating the problem:
> 
> > library( 'affy' );
> > library( 'affyPLM' );
> > ab = read.affybatch( 
> filenames='/home/allenday/cel/0001.CEL' ); ab; # 
> > works, output omitted pm( ab ); # works, output omitted mm( ab ); # 
> > fails!
> Error in FUN(X[[1411190]], ...) : subscript out of bounds
> > plm = fitPLM( ab ); #same failure in fitPLM, caused by a 
> call to mm() 
> > on variable ab;
> Error in FUN(X[[1411190]], ...) : subscript out of bounds
> 
> I'm only proficient enough in R and C to track this down -- 
> I'm don't know R or Bioconductor well enough to know how to 
> fix it.  If I can get this going I will submit a new package 
> that provides just.nuse() and just.rle() functions.  Can 
> someone give me a pointer for how to make this work?
> 
> Thanks.
> 
> -Allen
> 
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: 
> http://news.gmane.org/gmane.science.biology.informatics.conductor
>
 
--------------------------------------------------------

 
This email is confidential and intended solely for the use o...{{dropped}}