全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
5529 5
2014-04-21
怎样使用R软件运用Boostrap方法估计相关系数的标准误差?
二维码

扫码加我 拉你入群

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

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

全部回复
2014-4-27 19:35:52
用本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))   #样本标准差点估计及其方差
二维码

扫码加我 拉你入群

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

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

2014-4-27 19:37:51
样本X
N=length(x)
rou=sapply(1:5000,function(o){sample(x,N,replace=TRUE)})
sd(rou)
二维码

扫码加我 拉你入群

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

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

2014-4-27 19:38:25
huyiustc 发表于 2014-4-27 19:37
样本X
N=length(x)
rou=sapply(1:5000,function(o){sample(x,N,replace=TRUE)})
不好意思第一次没看清题目
二维码

扫码加我 拉你入群

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

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

2014-4-27 19:52:20
前面的不对,不好意思
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)

这个是对的
很抱歉啊
二维码

扫码加我 拉你入群

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

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

2015-10-23 09:22:25
huyiustc 发表于 2014-4-27 19:52
前面的不对,不好意思
X为样本数据对N*2的矩阵
N=length(X)
谢谢你啦!
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群