Dear list,

I try to run gcrma using the affinity computation without taking into account MM probes, and following commands lead to segfault.

> library(simpleaffy)

> rawData <- ReadAffy() # Read in PrimeView CELs (no MM probes)

> eset <- gcrma(rawData, affinity.info=compute.affinities2(cdfName(rawData)))

Adjusting for non-specific binding.
 *** caught segfault ***
address 0x2116d0850, cause 'memory not mapped'

Traceback:
 1: runmed(x[O], 21)
 2: bg.parameters.ns(mms[index.affinities], amm, apm)
 3: bg.adjust.fullmodel(pms[, i], mms[, i], ncs = ncs[, i], apm = pm.affinities,     amm = mm.affinities, anc = anc, index.affinities, k = k,     rho = rho, fast = fast)
 4: gcrma.engine(pms = pm(object), mms = mm(object), ncs, pm.affinities = pm(affinity.info),     mm.affinities = mm(affinity.info), anc, type = type, k = k,     stretch = stretch, correction = correction, GSB.adjust = GSB.adjust,     rho = rho, verbose = verbose, fast = fast)
 5: bg.adjust.gcrma(object, affinity.info = affinity.info, affinity.source = affinity.source,     NCprobe = NCprobe, type = type, k = k, stretch = stretch,     correction = correction, GSB.adjust = GSB.adjust, rho = rho,     optical.correct = optical.correct, verbose = verbose, fast = fast)
 6: gcrma(rawData, affinity.info = compute.affinities2(cdfName(rawData)))

> sessionInfo()
R version 3.1.1 (2014-07-10)
Platform: x86_64-unknown-linux-gnu (64-bit)

locale:
 [1] LC_CTYPE=zh_CN.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=zh_CN.UTF-8        LC_COLLATE=zh_CN.UTF-8    
 [5] LC_MONETARY=zh_CN.UTF-8    LC_MESSAGES=zh_CN.UTF-8   
 [7] LC_PAPER=zh_CN.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=zh_CN.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] parallel  stats     graphics  grDevices utils     datasets  methods  
[8] base     

other attached packages:
[1] simpleaffy_2.42.0   gcrma_2.38.0        genefilter_1.48.1  
[4] affy_1.44.0         Biobase_2.26.0      BiocGenerics_0.12.0

loaded via a namespace (and not attached):
 [1] affyio_1.34.0         annotate_1.44.0       AnnotationDbi_1.28.0 
 [4] BiocInstaller_1.16.0  Biostrings_2.34.0     DBI_0.3.1            
 [7] GenomeInfoDb_1.2.2    IRanges_2.0.0         preprocessCore_1.28.0
[10] RSQLite_1.0.0         S4Vectors_0.4.0       splines_3.1.1        
[13] stats4_3.1.1          survival_2.37-7       XML_3.98-1.1         
[16] xtable_1.7-4          XVector_0.6.0         zlibbioc_1.12.0      

For PM-only arrays you have to specify the NCprobes that you want to use. The primeview array doesn't really have any negative control probes, so this can be a bit problematic. There are the ERCC probes, and if you didn't use the ERCC spike in cDNA, then you could hypothetically use those.

ind <- grep("ERCC", featureNames(rawData), value = TRUE)
eset <- gcrma(rawData, affinity.info=compute.affinities2(cdfName(rawData)), NCprobe = ind, type = "affinities")