We're trying to calculate the difference in residuals between two groups 
and generate mufScores (FIRMAGene) on those residual differences. 
Calculating the differences goes smoothly:

rsu_diff <- lapply(unlist(rsu, use.names=FALSE), byrow=FALSE, 

where cls is a grouping variable (e.g. cls <- c('A', 'A, 'B', 'B')

Next we apply the mufColumns function to the elements of the rsu_diff 

mufScores <- lapply(rsu_diff[w], FUN = function(u), c(mufColumns(u)))

w is an indicator variable that only keeps genes (this is gene level 
analysis) if the gene has more than 3 probes and less than some other 
number. If we set the upper limit of w to >500 we get the following error.

Error in mufC(x): Nan/NA/Inf in foreign function call

Having looked at the C code for the relevant function (mps.c) I can only 
see one line where that could cause this:


Can anyone shed any light on why we can't run mufColumns if we select genes 
with more than 500 probes over the gene?



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/

You received this message because you are subscribed to the Google Groups 
"aroma.affymetrix" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to aroma-affymetrix+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

Reply via email to