[R] confused about CORREP cor.LRtest

Mark W Kimpel mwkimpel at gmail.com
Fri Mar 7 18:29:07 CET 2008


After some struggling with the data format, non-standard in 
BioConductor, I have gotten cor.balance in package CORREP to work. My 
desire was to obtain maximum-likelihood p-values from the same data 
object using cor.LRtest, but it appears that this function wants 
something different, which I can't figure out from the documentation.

Briefly, my dataset consists of 36 samples from 12 conditions and I have 
  497 genes of interest to be correlated. The following works:

M <- cor.balance(stddata, m = 3, G=497)

The following does not:
M.p <- cor.LRtest(stddata, m1 = 3, m2 = 3)

Do I need to do something to stddata between example 1 and 2 or does "m" 
stand for something different in the two examples?

sessionInfo follows. Thanks, Mark

 > sessionInfo()
R version 2.7.0 Under development (unstable) (2008-03-05 r44683)
x86_64-unknown-linux-gnu

locale:
LC_CTYPE=en_US.UTF-8;LC_NUMERIC=C;LC_TIME=en_US.UTF-8;LC_COLLATE=en_US.UTF-8;LC_MONETARY=en_US.UTF-8;LC_MESSAGES=en_US.UTF-8;LC_PAPER=en_US.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=en_US.UTF-8;LC_IDENTIFICATION=C

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

other attached packages:
  [1] rat2302_2.0.1        Rgraphviz_1.17.13    graph_1.17.17
  [4] igraph_0.5           CORREP_1.5.0         e1071_1.5-17
  [7] class_7.2-41         affy_1.17.8          preprocessCore_1.1.5
[10] affyio_1.7.13        Biobase_1.99.1

loaded via a namespace (and not attached):
[1] cluster_1.11.10
-- 

Mark W. Kimpel MD  ** Neuroinformatics ** Dept. of Psychiatry
Indiana University School of Medicine

15032 Hunter Court, Westfield, IN  46074

(317) 490-5129 Work, & Mobile & VoiceMail
(317) 204-4202 Home (no voice mail please)

mwkimpel<at>gmail<dot>com



More information about the R-help mailing list