[BioC] error in runBioHMM() (package snapCGH)
Wolfgang Raffelsberger
wraff at titus.u-strasbg.fr
Mon May 4 10:31:41 CEST 2009
Dear list,
when analyzing Agilent 44k (human) CGH data* *using the function
runBioHMM() from the library snapCGH I've run into an error message
described below.
Curiously this error has never occurred before although I've been using
more or less the same code with previous data-sets where the code was
working well ! So I wonder if there's something special with may data,
but I can't find any particular anomalies in the data$M.observed (eg
I've checked that there are no NAs ...)
A summary of the data didn't indicate anything "completely" different to
the other arrays, at least nothing that gave the impression beyond those
differences due to the nature/biology of the biospies.
If you would like to get a copy of the real data, pls let me know, of
course the file-size would be way to large for attaching to this message.
As the function runBioHMM() doesn't contain a line containing "if
((phi[3, t - 1] + max(-1e+05, log(BFGS.trans.mat[ ..." I'm stuck to get
some clues of what may have caused the problem.
Note that other segmentation functions form the same package like
runDNAcopy() or runGALD() did work without problems !!
I'm running code like :
library(limma); library(snapCGH); library(DNAcopy); library(GLAD);
library(smoothseg)
RG1 <- read.maimages(file=mySamples,columns=extrCols1,other.columns=
extrOtherCols,source = "agilent")
MA1 <- normalizeBetweenArrays(RG1,method="Aquantile")
MA1$design <- rep(-1,ncol(MA1$M))
MA1_snpCGH <- processCGH(MA1, maxChromThreshold =
length(na.omit(unique(MA1$genes$Chr))) )
MA1_snpCGH$genes$Position <- round(MA1_snpCGH$genes$Position,6)
# transform to Mb scale
## checking (some of) my data
summary( MAagil_snpCGH$M.observed[,155:156] ) ## gives :
M.111_CN_18 M.111_CN_19
Min. :-4.856865 Min. :-2.919297
1st Qu.:-0.243359 1st Qu.:-0.235424
Median :-0.044505 Median :-0.039601
Mean : 0.007306 Mean : 0.001244
3rd Qu.: 0.163834 3rd Qu.: 0.152406
Max. : 7.773869 Max. : 4.161682
MA1_snpCGH$M <- MA1_snpCGH$M.observed
## and here I (consitently) get the error : (message truncated to show
the end)
sample is 153 Chromosomes: 1 2 3 4 5 6 7 8 9 10 11 12 13
14 15 16 17 18 19 20 21 22 23 24
sample is 154 Chromosomes: 1 2 3 4 5 6 7 8 9 10 11 12 13
14 15 16 17 18 19 20 21 22 23 24
sample is 155 Chromosomes: 1 2 3 4 5 6 7 8 9 10 11 12 13
14 15 16 17 18 19 20 21 22 23 24
sample is 156 Chromosomes: 1 2 3 4 5 6 7 8 9 10 11 Error
in if ((phi[3, t - 1] + max(-1e+05, log(BFGS.trans.mat[[t - 1]][3, :
missing value where TRUE/FALSE needed
In addition: There were 50 or more warnings (use warnings() to see the
first 50)
Error in res$nstates.list : $ operator is invalid for atomic vectors
## for completeness:
> sessionInfo()
R version 2.8.0 (2008-10-20)
x86_64-unknown-linux-gnu
locale:
C
attached base packages:
[1] splines grid tools stats graphics grDevices utils
[8] datasets methods base
other attached packages:
[1] smoothseg_0.0.1 OCplus_1.16.0 akima_0.5-1
[4] snapCGH_1.10.0 aCGH_1.16.0 sma_0.5.15
[7] multtest_1.22.0 cluster_1.11.11 GLAD_1.18.0
[10] DNAcopy_1.16.0 tilingArray_1.20.0 pixmap_0.4-9
[13] geneplotter_1.20.0 annotate_1.20.1 xtable_1.5-4
[16] AnnotationDbi_1.4.1 genefilter_1.22.0 survival_2.34-1
[19] vsn_3.8.0 lattice_0.17-15 strucchange_1.3-4
[22] sandwich_2.1-0 zoo_1.5-4 RColorBrewer_1.0-2
[25] affy_1.20.0 Biobase_2.2.1 limma_2.16.3
loaded via a namespace (and not attached):
[1] DBI_0.2-4 KernSmooth_2.22-22 RSQLite_0.7-1
[4] affyio_1.10.1 preprocessCore_1.4.0
Any suggestions how to resolve this problem would be very much appreciated !
Wolfgang Raffelsberger
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Wolfgang Raffelsberger, PhD
Laboratoire de BioInformatique et Génomique Intégratives
CNRS UMR7104, IGBMC,
1 rue Laurent Fries, 67404 Illkirch Strasbourg, France
Tel (+33) 388 65 3300 Fax (+33) 388 65 3276
wolfgang.raffelsberger (at) igbmc.fr
More information about the Bioconductor
mailing list