[R] Help to find correlation. Oops

Martin Henry H. Stevens HStevens at MUOhio.edu
Thu Feb 16 20:57:07 CET 2006


Hi,
Ignore the previous post-it contained a typo. Sorry! The following is  
ok.
Hank

See help files for read.csv for input.
See help files for cor.test
dat <- read.csv("whatever-your-file-is", other info...)

# Generate an upper triangular  matrix of Pearson product moment  
correlation coefficients
r.s <- matrix(NA, nr=dim(dat)[2], nc=dim(dat)[2])
rownames(r.s) <- names(dat)
colnames(r.s) <- names(dat)
for(i in 1:(dim(dat)[2]-1)) {
   for(j in i:dim(dat)[2]) r.s[i,j] <- {cor.test(dat[,i], dat[,j]) 
$estimate}}
r.s

# Generate an upper triangular  matrix of P values of the r's
p.s <- matrix(NA, nr=dim(dat)[2], nc=dim(dat)[2])
rownames(p.s) <- names(dat)
colnames(p.s) <- names(dat)
for(i in 1:(dim(dat)[2]-1)) {
   for(j in i:dim(dat)[2]) p.s[i,j] <- {cor.test(dat[,i], dat[,j]) 
$p.value}}
p.s

Please do read the help files.
Hank

On Feb 16, 2006, at 1:17 PM, Anusha Aiyaloo kannan ((RIT Student))  
wrote:

>
> Respected Sir,
>
> I am trying to import excel file into R, but I need to truncate  
> some columns from the original file.
>
>  How to delete unwanted columns when I import data from excel file.
>
> How to use cor.test for the data when I want the output rowwise.
>
> How to do grouping and use cor.test on that data
>
> I need some help regarding how to calculate the correlation.
> I don't know whether you understood my question, but I need help.
>
> Any help is appreciated.
>
> Thanks,
>
>
>
> 	[[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting- 
> guide.html

Dr. Martin Henry H. Stevens, Assistant Professor
338 Pearson Hall
Botany Department
Miami University
Oxford, OH 45056

Office: (513) 529-4206
Lab: (513) 529-4262
FAX: (513) 529-4243
http://www.cas.muohio.edu/~stevenmh/
http://www.muohio.edu/ecology/
http://www.muohio.edu/botany/
"E Pluribus Unum"




More information about the R-help mailing list