Re: [R] test for difference between non-independent correlations

From: Jonathan Baron <baron_at_psych.upenn.edu>
Date: Tue 27 Jul 2004 - 23:05:33 EST

On 07/27/04 13:35, Avril Coghlan wrote:
>Hello,
>
> I am wondering whether there is a way to
>test whether two non-independent correlation
>coefficients are significantly different, in R?

About 20 years ago I asked the same question (without the "R"), found Steiger's paper (below) and wrote this BASIC program to do it. (If you try to type it in each time you're bound to make a mistake. This has been checked repeatedly.) It is trivial to write this as an R script, of course.

I and my colleagues have been using this for about 20 years, and it seems to be OK as judged by comparing it to the jackknife. (Remember, when we did this, the "bootstrap" was a newfangled idea that we didn't quite trust.) My hunch is that something better has come along, and I'd be interested to see how others respond. Hence I haven't bothered to translate this into R. Or maybe everyone who has this problem just uses the bootstrap now. That is certainly easy enough. (I will eventually add this to our "Notes on R for psychology...," since it is a very common problem in psychology research.)

10 print "looks for difference between r12 and r13, given r23"
12 print "(input r12>1 to quit)"
20 input "r12 "; r12

22 if r12>1 then end
30 input "r13 "; r13
40 input "r23 "; r23
50 input "N "; n

52 let rd=1-r12*r12-r13*r13-r23*r23+2*r12*r13*r23
54 let rb=(r12+r13)/2
56 let cube=(1-r23)*(1-r23)*(1-r23)
60 let t2=(r12-r13)*sqrt((n-1)*(1+r23)/(2*rd*(n-1)/(n-3)+rb*rb*cube))
70 print "t (";n-3;") = ";t2

80 print "Steiger, J.H. (1980). Tests for comparing elements of a 82 print "correlation matrix. Psychological Bulletin, 87, 245-251." 100 goto 12
-- 
Jonathan Baron, Professor of Psychology, University of Pennsylvania
Home page: http://www.sas.upenn.edu/~baron
R search page: http://finzi.psych.upenn.edu/

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.htmla>
Received on Tue Jul 27 23:15:06 2004

This archive was generated by hypermail 2.1.8 : Wed 03 Nov 2004 - 22:55:17 EST