>From what I can see there is a bug in the function
findUnitsTodo.ChipEffectFile.
The following code fragment is supposed to (I think) get one index per
unit, to be used to access the standard deviation for that index, and
test if that value is <= 0.
idxs <-
lapplyInChunks(units,
Hi,
the factor of two has to do with the number of parameters
estimated/stored, i.e. storing theta, vs (thetaA, thetaB). I cannot
reproduce this, but I one guess is that you have tried both
combineAlleles=TRUE and combineAlleles=FALSE and somehow there exist
cached indices that are reused from o
Hi,
I don't know what the problem is. I tried your recipe on the HapMap
250K_Nsp and everything went fine:
library("aroma.affymetrix")
log <- Arguments$getVerbose(-8, timestamp=TRUE)
cdf <- AffymetrixCdfFile$byChipType("Mapping250K_Nsp")
cs <- AffymetrixCelSet$byName("HapMap270,500K,CEU,testSet"
