Hi, I'll leave the details to FIRMA experts, but you are using really large values of argument 'ram'. It might be that you ran out of memory (the you got an error message). If you used cut'n'paste, instead of source(), to do the analysis, it might be that one of the fit() methods was preemptively finished. If so, not all units have been fitted. Rerun with ram=1 to see if you get a different result; the units already fitted will be skipped.
I'm also not sure if > 431210/(1190297*107) [1] 0.003385710 is an exceptionally large fraction. Note also that 431210/107 is exactly 4030; it could be that it is the exact same 4030 units that are NA in all samples. That could be explained by some units are not fitted or are Affymetrix control units. BTW, not all methods take argument 'ram'. Instead, use the global aroma settings for achieving the same, e.g. setOption(aromaSettings, "memory/ram", 50). More info at http://aroma-project.org/settings. This way your scripts are clean and can be ran "as-is" on other machines with less memory. Maybe this helps(?) /Henrik On Mon, Aug 2, 2010 at 7:24 PM, Adi Tarca <ata...@med.wayne.edu> wrote: > Hi all, > I have a batch of 107 mice exon arrays for which I computed FIRMA > scores and I got many NaN, Inf and 0 values which disable further > analysis based on log FIRMA values for some probesets. I was wondering > if this is a known issue or I am the only one to get these results. > > > Here is the code I use to get the FIRMA scores: > > library(aroma.affymetrix) > verbose <- Arguments$getVerbose(-8, timestamp=TRUE) > chipType <- "MoEx-1_0-st-v1" > cdf <- AffymetrixCdfFile$byChipType(chipType, > tags="fullR1,A20080718,MR") > cs <- AffymetrixCelSet$byName("mice2010", cdf=cdf) > bc <- RmaBackgroundCorrection(cs, tag="fullR1,A20080718,MR") > csBC <- process(bc,verbose=verbose,ram=500) > qn <- QuantileNormalization(csBC, typesToUpdate="pm") > csN <- process(qn, verbose=verbose,ram=500) > plmTr <- ExonRmaPlm(csN, mergeGroups=TRUE) > fit(plmTr, verbose=verbose,ram=500) > > firma <- FirmaModel(plmTr) > fit(firma, verbose=verbose,ram=500) > fs <- getFirmaScores(firma) > myres2=extractDataFrame(fs,addNames=FALSE) > myres=as.matrix(myres2[,-(1:3)]) > > > Here is the counts of NaN Inf and 0 values: > >> dim(myres) > [1] 1190297 107 >> sum(is.nan(myres)) > [1] 431210 >> sum(is.infinite(myres),na.rm=TRUE) > [1] 855 >> sum(myres==0,na.rm=TRUE) > [1] 214 > > > > R version 2.11.0 (2010-04-22) > x86_64-unknown-linux-gnu > > locale: > [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C > [3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8 > [5] LC_MONETARY=C LC_MESSAGES=en_US.UTF-8 > [7] LC_PAPER=en_US.UTF-8 LC_NAME=C > [9] LC_ADDRESS=C LC_TELEPHONE=C > [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C > > attached base packages: > [1] stats graphics grDevices utils datasets methods > base > > other attached packages: > [1] aroma.affymetrix_1.5.0 aroma.apd_0.1.7 > affxparser_1.20.0 > [4] R.huge_0.2.0 aroma.core_1.5.0 > aroma.light_1.16.0 > [7] matrixStats_0.2.1 R.rsp_0.3.6 > R.cache_0.3.0 > [10] R.filesets_0.8.1 digest_0.4.2 > R.utils_1.4.0 > [13] R.oo_1.7.2 R.methodsS3_1.2.0 > > > > > > > > > > > > > > > > -- > 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/