[R] confused about CORREP cor.LRtest

From: Mark W Kimpel <mwkimpel_at_gmail.com>
Date: Fri, 07 Mar 2008 12:29:07 -0500

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


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)


R-help_at_r-project.org mailing list
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Received on Fri 07 Mar 2008 - 17:31:22 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Fri 14 Mar 2008 - 03:30:21 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive