全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
1432 9
2013-12-30
chouyang<-function(x) {data1<-list(matrix(rep(0,3000),1000,3));
  for(i in 1:50)
  {xx<-x[sample(1:nrow(x),5),];hgfx<-lm(xx$V1~xx$V2+xx$V3);
   data1[i]<-hgfx$coefficients }
}
我的目的是要创建一个名为chouyang的函数,作用是从数据框x中抽取5行记录,然后对变量V1和V2,V3做回归分析,结果保存在hgfx中,然后只要回归分析的三个估计系数。这样做50次,每次结果保存在事先创建好的一个列表矩阵中,但是运行老不成功,请教大神,希望不吝赐教,谢谢。
二维码

扫码加我 拉你入群

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

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

全部回复
2013-12-30 22:13:49
试出来了,改动一个地方就行
data1[[1]][i,] <- hgfx$coefficients
其实楼主的问题可以用一句来完成
复制代码


二维码

扫码加我 拉你入群

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

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

2013-12-30 22:15:00
感觉易出的地方:
1、data1<-list(matrix(rep(0,3000),1000,3)) 不好
2、data1【i】有问题,且没有返回值
3、lm(xx$V1~xx$V2+xx$V3) 可能会出现名字不同,特别是 大小写
复制代码
chouyang(data) 你看看如何
二维码

扫码加我 拉你入群

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

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

2013-12-30 22:17:41
复制代码
稍微改了下  ,应该可以了吧,前提是你的x是数据框,且含有变量V1,V2,V3,否则得定义
二维码

扫码加我 拉你入群

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

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

2013-12-30 22:24:17
yywan0913 发表于 2013-12-30 22:17
稍微改了下  ,应该可以了吧,前提是你的x是数据框,且含有变量V1,V2,V3,否则得定义
复制代码
二维码

扫码加我 拉你入群

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

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

2013-12-30 22:45:33
我也要修改,
保留 (Intercept)  v2, v3

复制代码


二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

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