怎样使用R软件运用Boostrap方法估计相关系数的标准误差？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

怎样使用R软件运用Boostrap方法估计相关系数的标准误差？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：boostrap boost 标准误差 相关系数 Trap 软件

用本x
N=length(x)
n=5000 重复5000次
mu=sapply(1:n,function(o){y=sample(x,N,replace=TRUE);mean(y)})
c(mean(mu),sd(mu))   #计算样本均值及均值标准差
SD=sapply(1:n,function(o){y=sample(x,N,replace=TRUE);sd(y)})
c(mean(SD),sd(SD))   #样本标准差点估计及其方差

样本X
N=length(x)
rou=sapply(1:5000,function(o){sample(x,N,replace=TRUE)})
sd(rou)

huyiustc 发表于 2014-4-27 19:37
样本X
N=length(x)
rou=sapply(1:5000,function(o){sample(x,N,replace=TRUE)})

不好意思第一次没看清题目

前面的不对，不好意思
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)

这个是对的
很抱歉啊

huyiustc 发表于 2014-4-27 19:52
前面的不对，不好意思
X为样本数据对N*2的矩阵
N=length(X)

谢谢你啦！