有关R语言for循环问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

> p=8
> n=200
> tbeta=c(3,1.5,2,rep(0,p-3))
> Sig=matrix(rep(1,p*p),p)
> for (ii in 1:p) for (jj in 1:p) if (ii!=jj) Sig[ii,jj] = 0.5^(abs(ii-jj))
> x = mvrnorm(n, rep(0,p), Sig)
> y = x%*%tbeta + rnorm(n,0,3)
> cv.out<-cv.glmnet(x,y,nfolds=5)
> tbetahat<-coef(cv.out)[2:9]
> yhat = x%*%tbetahat + rnorm(n,0,3)
> a = y-yhat
> pse=sum(a^2)/n
> pse
[1] 17.4298
我想把这段代码重复50次，得到50个这个值，该怎么用for语句，尝试过很多次，总是不行

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：for循环 For R语言 matrix beta matrix

all_data=c()
for(k in 1:50)
{
p=8
n=200
tbeta=c(3,1.5,2,rep(0,p-3))
Sig=matrix(rep(1,p*p),p)
for (ii in 1:p) for (jj in 1:p) if (ii!=jj) Sig[ii,jj] = 0.5^(abs(ii-jj))
x = mvrnorm(n, rep(0,p), Sig)
y = x%*%tbeta + rnorm(n,0,3)
cv.out<-cv.glmnet(x,y,nfolds=5)
tbetahat<-coef(cv.out)[2:9]
yhat = x%*%tbetahat + rnorm(n,0,3)
a = y-yhat
pse=sum(a^2)/n
all_data=c(all_data,pes)

}
all_data

不知道是不是您想要的结果

cbq95229 发表于 2016-8-29 17:06
不知道是不是您想要的结果

是的，谢谢大神

王新军 发表于 2016-8-29 18:27
是的，谢谢大神

> s$coef[4,2:9]
                 b.1                  b.2                  b.3                  b.4                  b.5                  b.6                  b.7
   "Mean   :3.556  "   "Mean   :1.2632  "   "Mean   :1.6506  "  "Mean   : 0.1332  " "Mean   : 0.06361  " "Mean   : 0.03635  " "Mean   :-0.01616  "
                 b.8
"Mean   :-0.1406  "
我再问一下，我想把这里面的数据拿出来，不想要Mean:,该怎么弄

cbq95229 发表于 2016-8-29 17:06
不知道是不是您想要的结果

> s$coef[4,2:9]
                 b.1                  b.2                  b.3                  b.4                  b.5                  b.6                  b.7
   "Mean   :3.556  "   "Mean   :1.2632  "   "Mean   :1.6506  "  "Mean   : 0.1332  " "Mean   : 0.06361  " "Mean   : 0.03635  " "Mean   :-0.01616  "
                 b.8
"Mean   :-0.1406  "
我再问一下，我想把这里面的数据拿出来，不想要Mean:,该怎么弄

王新军 发表于 2016-8-29 18:44
> s$coef[4,2:9]
                 b.1                  b.2                  b.3                  b ...

try this:

x<-s$coef[4,2:9]
z<-sapply(x,function(y) as.numeric(unlist(strsplit(y,":"))[2]))
z

都很厉害，学习了！