Re: [R] cross-correlation table with subscript or superscript to indicate significant differences

From: yoav baranan <ybaranan_at_hotmail.com>
Date: Thu, 05 May 2011 19:40:17 +0300

Here is an example for my earlier question. Say you have a 3x3 correlation matrix:corrs <- matrix(c(0.25,0.32,0.66,0.14,0.24,0.34,0.44,0.34,0.11), nrow=3, ncol=3, dimnames = list(c('varA','varB', 'varC'), c('varA','varB', 'varC')))And another matrix for the sample size of each correlation:sizes <- matrix(c(44,68,313,142,144,207,201,100,99), nrow=3, ncol=3, dimnames = list(c('varA','varB', 'varC'), c('varA','varB', 'varC')))corrs looks like this:      varA varB varCvarA 0.05 0.14 0.44varB 0.32 0.24 0.34varC 0.66 0.57 0.50sizes:     varA varB varCvarA   44  142  201varB   68  144  100varC  313  207   99i.e., the correlation between variables A and C was 0.66 with sample size of 313. (I got these tables from rcorr).What I want to do is to compare the correlations in each row (probably using r.test), and then create a correlation table with subscripts or superscripts indicating the significance "group" (again: correlations with different superscripts, in the same row, are signif!
 icantly different from each other). Something like this:       varA varB varCvarA 0.05b 0.14b 0.44avarB 0.32a 0.24a 0.34avarC 0.66a 0.57ab 0.50bOf course, I don't have a 3x3 table. I have about 20 tables of at least 7x7 each, so this is why I'm looking for methods to automate the process. Thanks,Yoav

> CC: r-help@r-project.org
> From: dwinsemius_at_comcast.net
> To: ybaranan_at_hotmail.com
> Subject: Re: [R] cross-correlation table with subscript or superscript to indicate significant differences
> Date: Thu, 5 May 2011 12:17:25 -0400
>
>
> On May 5, 2011, at 10:48 AM, yoav baranan wrote:
>
> >
> > Hi, I wonder whether the following is possible with R, and whether
> > anyone has done that and can share his/her code with me. I have a
> > correlation matrix, and I want to create a correlation table that I
> > can copy to Microsoft Word with a superscript above each
> > correlation, indicating significant differences in the same row.
> > That is, when correlations in the same row do not share superscript,
> > it means that they are significantly different from each other.
> > thanks,yoav
> > [[alternative HTML version deleted]]
>
> An example with data and the desired result might help focus the
> discussion.
>
> This shows how to set up an example showing how extract the row
> numbers from a correlation matrix with absolute values above 0.5 but
> less than 1 (to exclude the trivial cases).
>
> > set.seed(123)
> > X <- matrix(rnorm(100), 10)
> > apply(cor(X), 2, function(x) which(abs(x) > 0.5 & x < 1) )
> [[1]]
> [1] 2 4 8
>
> [[2]]
> [1] 1 3
>
> [[3]]
> [1] 2 6 9
>
> [[4]]
> [1] 1 7
>
> [[5]]
> integer(0)
>
> [[6]]
> [1] 3 10
>
> [[7]]
> [1] 4
>
> [[8]]
> [1] 1
>
> [[9]]
> [1] 3
>
> [[10]]
> [1] 6
>
> This would extract the rownames if they are letters[1:10]
>
> > lapply( apply(cor(X), 2, function(x) which(abs(x) > 0.5 & x < 1) ),
> function(x) rownames(X)[x])
> [[1]]
> [1] "b" "d" "h"
>
> [[2]]
> [1] "a" "c"
>
> [[3]]
> [1] "b" "f" "i"
>
> [[4]]
> [1] "a" "g"
>
> [[5]]
> character(0)
>
> [[6]]
> [1] "c" "j"
>
> [[7]]
> [1] "d"
>
> [[8]]
> [1] "a"
>
> [[9]]
> [1] "c"
>
> [[10]]
> [1] "f"
>
> Exactly how we are supposed to pass this to MS Word does not seem to
> be a proper question for this mailing list.
>
> --
>
> David Winsemius, MD
> West Hartford, CT
>
                                               

        [[alternative HTML version deleted]]



R-help_at_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. Received on Thu 05 May 2011 - 18:01:10 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 Thu 05 May 2011 - 18:10:06 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