# 無相関の場合の順位相関係数の分布
# 2006 年 6 月 11 日
repp=100000
xx=c(7,8,6,4,3,5,2,1)
yy=c(5,3,4,2,8,1,6,7)
sl=cor(xx,yy)
sl
sample=8
corr=rep(0,repp)
for(i in 1:repp){
x<-rnorm(sample,0,1)
y<-rnorm(sample,0,1)
corr[i]=cor(rank(x),rank(y))
}
mean(corr)
sd(corr)
hist(corr)
sum(corr<=sl)/repp