thank you john.
however, I am finding it difficult to automate on a matrix.
Pardon my ignorance in R computing:
I do not know how to automate on a matrix.
If I do the following it works:
x = cor.test(d6[1,],d6[2,])
x
Pearson's product-moment correlation
data: d6[1, ] and d6[2, ]
t = 10.5196, df = 10, p-value = 9.973e-07
alternative hypothesis: true correlation is not equal to 0
95 percent confidence interval:
0.8520623 0.9883592
sample estimates:
cor
0.9576655
If I want to run it on all rows, I do not know how to do it.
I tried following,
lapply(d6,cor.test)
Error in cor.test.default(X[[1L]], ...) :
element 1 is empty;
the part of the args list of 'length' being evaluated was:
(y)
sapply(d6,cor.test)
Error in cor.test.default(X[[1L]], ...) :
element 1 is empty;
the part of the args list of 'length' being evaluated was:
(y)
for(i in 1:14659){
+ k = i+1
+ cor.test(d6[i,],d6[k,])
+ x = cor.test(d6[i,],d6[k,])
+ return(x)}
Error: no function to return from, jumping to top level
I appreciate your help.
thank you.
Adrian
On Sun, Sep 20, 2009 at 6:17 PM, Adrian Johnson
<oriolebaltim...@gmail.com> wrote:
thank you john.
however, I am finding it difficult to automate on a matrix.
Pardon my ignorance in R computing:
I do not know how to automate on a matrix.
If I do the following it works:
x = cor.test(d6[1,],d6[2,])
x
Pearson's product-moment correlation
data: d6[1, ] and d6[2, ]
t = 10.5196, df = 10, p-value = 9.973e-07
alternative hypothesis: true correlation is not equal to 0
95 percent confidence interval:
0.8520623 0.9883592
sample estimates:
cor
0.9576655
If I want to run it on all rows, I do not know how to do it.
I tried following,
lapply(d6,cor.test)
Error in cor.test.default(X[[1L]], ...) :
element 1 is empty;
the part of the args list of 'length' being evaluated was:
(y)
sapply(d6,cor.test)
Error in cor.test.default(X[[1L]], ...) :
element 1 is empty;
the part of the args list of 'length' being evaluated was:
(y)
for(i in 1:14659){
+ k = i+1
+ cor.test(d6[i,],d6[k,])
+ x = cor.test(d6[i,],d6[k,])
+ return(x)}
Error: no function to return from, jumping to top level
I appreciate your help.
thank you.
Adrian
On Sun, Sep 20, 2009 at 5:13 PM, John Kane <jrkrid...@yahoo.ca> wrote:
?cor
?cor.test
--- On Sun, 9/20/09, Adrian Johnson <oriolebaltim...@gmail.com> wrote:
From: Adrian Johnson <oriolebaltim...@gmail.com>
Subject: [R] correlation help
To: r-help@r-project.org
Received: Sunday, September 20, 2009, 5:00 PM
Dear group,
I have a matrix like the following:
Name Sample1
sample2 sample3 sample4 .....
sample(n)
nm1 10.5
13.5
30 31
nm2 8
11
34
29
nm3 9
10.3
27.8 35
nm(j)
I want to be able to calculate correlation between
all pairs of names.
For example (nm1,nm2), (nm1,nm3), (nm1,nmj), (nm2,nm3),
(nm2,nmj)....
Then I want to calculate the significance of correlation
using t-score
or p-value.
I can calculate correlation coeffecient in excel but not
significance
in both excel and R.
I want to be able to do it in R, I appreciate your help.
thank you.
Ad.
______________________________________________
R-help@r-project.org
mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
>>> PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained,
reproducible code.
__________________________________________________________________
Be smarter than spam. See how smart SpamGuard is at giving junk
email the boot with the All-new Yahoo! Mail. Click on Options in
Mail and switch to New Mail today or register for free at
http://mail.yahoo.ca
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
