Entering edit mode
hi gurus
i work with a matrix of data and there are no NAs or Infinite or
missing values
i tried to work with this package but i get all sorts of errors:
1) the initializeHP() function, called as HP <- initializeHP(D, cond)
when a value in the matrix created with the makeMyD() function is ==
1.000000 or -1.00000, the initializeHP() function will somehow create
NAs and quit giving me the following error
Error in mvnX(data = data, prior = prior):
NA/NaN/Inf in foreign function call (arg 1)
the condition variable was created as
cond <- c(rep(1, 3), rep(2, 3))
and the input matrix for the initializeHP function was created as
D <- makeMyD(myData, cond)
2) the HP variable (matrix) contains no NAs/Inf but when i call the
ebCoexpressOneStep(D, cond, pat, HP) function it gives me an error:
Error in if (mix[1] < 0.8) { : missing value where TRUE/FALSE needed
and i find out that the function creates a matrix of ONLY NaN..
pat was created with:
pat <- ebPatterns(c("1,1","1,2")) as requested by the documentation
so....any help on how to solve these issues?
did anyone use this package successfully in the past?
thanks
Seb
R version 2.15.1 (2012-06-22)
Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit)
locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] EBcoexpress_1.2.0 minqa_1.2.1 Rcpp_0.10.1
mclust_4.0
[5] EBarrays_2.22.0 lattice_0.20-10 Biobase_2.18.0
BiocGenerics_0.4.0
loaded via a namespace (and not attached):
[1] cluster_1.14.3 grid_2.15.1 tools_2.15.1