What is the name of your custom CDF, and what does print(cdf) report after loading it with cdf <- AffymetrixCdfFile$byName(...)?
/Henrik On Wed, May 4, 2011 at 3:56 PM, Fong <fongchunc...@gmail.com> wrote: > Hi all, > > I've created a customCDF following the procedure listed at: > http://www.aroma-project.org/node/40 > > The customCDF I made was based on the Brainarray group CDF (ense14 for > Human Exon arrays to be exact). I wanted to create a CDF which I > could run FIRMA on since the customCDF provided by the BrainArray > doesn't appear map the probeset (group ids) onto the appropriate gene > (unit IDs) which is necessary to run FIRMA. I create one that works > in the flat2CDF() script provided by aroma.affymetrix and then when I > used it on aroma affymetrix I got the following error: > > 20110504 15:46:17| Setting up ExonChipEffectFile...done > 20110504 15:46:17| Retrieving chip-effect #19 of 20 (8_4N)...done > 20110504 15:46:17| Retrieving chip-effect #20 of 20 (9_5T)... > 20110504 15:46:17| Setting up ExonChipEffectFile... > 20110504 15:46:17| Pathname: plmData/ > gardina_colon_cancer,ense14,QN,RMA,merged/HuEx-1_0-st- > v2/9_5T,chipEffects.CEL > 20110504 15:46:17| Setting up ExonChipEffectFile...done > 20110504 15:46:17| Retrieving chip-effect #20 of 20 (9_5T)...done > 20110504 15:46:17| Retrieving chip-effects from data set...done > 20110504 15:46:17| Getting chip-effect set from data set...done > 20110504 15:46:17| Updating ExonChipEffectSet... > 20110504 15:46:17| Scanning for and applying sample annotation > files... > 20110504 15:46:17| Defining 0 files... > > 20110504 15:46:17| Defining 0 files...done > 20110504 15:46:17| No sample annotation files found. > 20110504 15:46:17| Scanning for and applying sample annotation > files...done > 20110504 15:46:17| Updating ExonChipEffectSet...done > 20110504 15:46:17| Identifying non-fitted units in chip-effect > file... > 20110504 15:46:17| Pathname: plmData/ > gardina_colon_cancer,ense14,QN,RMA,merged/HuEx-1_0-st- > v2/9_5T,chipEffects.CEL > 20110504 15:46:17| Found indices cached on file > 20110504 15:46:17| Reading data for these 442179 cells... > Error in readCel(getPathname(this), indices = idxs, readIntensities = > FALSE, : > Argument 'indices' is out of range [1,308580]. > Calls: fit ... findUnitsTodo -> findUnitsTodo.ChipEffectFile -> > readCel > 20110504 15:46:17| Reading data for these 442179 cells...done > 20110504 15:46:17| Identifying non-fitted units in chip-effect > file...done > 20110504 15:46:17| Identifying non-estimated units...done > 20110504 15:46:17|Fitting model of class ExonRmaPlm...done > > Has anyone run into this error before? Any ideas on what I am doing > wrong? > > Thanks, > > Fong > > -- > When reporting problems on aroma.affymetrix, make sure 1) to run the latest > version of the package, 2) to report the output of sessionInfo() and > traceback(), and 3) to post a complete code example. > > > You received this message because you are subscribed to the Google Groups > "aroma.affymetrix" group with website http://www.aroma-project.org/. > To post to this group, send email to aroma-affymetrix@googlegroups.com > To unsubscribe and other options, go to http://www.aroma-project.org/forum/ > -- When reporting problems on aroma.affymetrix, make sure 1) to run the latest version of the package, 2) to report the output of sessionInfo() and traceback(), and 3) to post a complete code example. You received this message because you are subscribed to the Google Groups "aroma.affymetrix" group with website http://www.aroma-project.org/. To post to this group, send email to aroma-affymetrix@googlegroups.com To unsubscribe and other options, go to http://www.aroma-project.org/forum/