- library('psych')
- cors <- seq(-.9,.9,.1)
- zs <- fisherz(cors)
- rs <- fisherz2r(zs)
- round(zs,2)
- n <- 30
- r <- seq(0,.9,.1)
- rc <- matrix(r.con(r,n),ncol=2)
- t <- r*sqrt(n-2)/sqrt(1-r^2)
- p <- (1-pt(t,n-2))/2
- r.rc <- data.frame(r=r,z=fisherz(r),lower=rc[,1],upper=rc[,2],t=t,p=p)
- round(r.rc,2)
-
-
- library('psych')
- n <- 30
- r <- seq(0,.9,.1)
- rc <- matrix(r.con(r,n),ncol=2)
- test <- r.test(n,r)
- r.rc <- data.frame(r=r,z=fisherz(r),lower=rc[,1],upper=rc[,2],t=test$t,p=test$p)
- round(r.rc,2)
- r.test(50,r)
- r.test(30,.4,.6) #test the difference between two independent correlations
- r.test(103,.4,.5,.1) #Steiger case A of dependent correlations
- r.test(n=103, r12=.4, r13=.5,r23=.1)
- #for complicated tests, it is probably better to specify correlations by name
- r.test(n=103,r12=.5,r34=.6,r13=.7,r23=.5,r14=.5,r24=.8) #steiger Case B
其中,函数r.con(r,n)是什么意思?
还有,r.test()函数也不太清楚,哪位大神解释一下?


雷达卡



京公网安备 11010802022788号







