前面的不对,不好意思
X为样本数据对N*2的矩阵
N=length(X)
rou=sapply(1:5000,function(o){y=X[sample(1:N,N,replace=TRUE),];cor(y)})
sd(rou)
给个例子
rm(list=ls())
x=rnorm(1000)
y=rnorm(1000)
X=cbind(x,y)
N=length(x)
rou=sapply(1:5000,function(o){v=X[sample(1:1000,N,replace=TRUE),];cor(v)})
sd(rou)
这个是对的
很抱歉啊