请教大神下列程序有何错误？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

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次，每次结果保存在事先创建好的一个列表矩阵中，但是运行老不成功，请教大神，希望不吝赐教，谢谢。

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：yang CHO Yan 程序

试出来了，改动一个地方就行
data1[[1]][i,] <- hgfx$coefficients
其实楼主的问题可以用一句来完成

1. sapply(1:50,function(y) coefficients(lm(V1~V2+V3,data=x[sample(1:nrow(x),5),])))

复制代码

感觉易出的地方：
1、data1<-list(matrix(rep(0,3000),1000,3)) 不好
2、data1【i】有问题，且没有返回值
3、lm(xx$V1~xx$V2+xx$V3) 可能会出现名字不同，特别是 大小写

1. chouyang<-function(x) {
   
2.   data1<-matrix(rep(0,300),100,3)
   
3.                        for(i in 1:50)
   
4.                        {xx<-x[sample(1:nrow(x),5),]
   
5.                         hgfx<-lm(xx$V1~xx$V2+xx$V3)
   
6.                         data1[i,]<-hgfx$coefficients                        
   
7.                        }
   
8.   data1
   
9. }
   
10. 
    
11. data <- data.frame(V1=1:1000,V2=sample(500,1000,replace=T),V3=sample(500,1000,replace=T))

复制代码

chouyang(data) 你看看如何

1. chouyang<-function(x) {data1<-matrix(rep(0,3000),1000,3);
   
2.   for(i in 1:50)
   
3.   {xx<-x[sample(1:nrow(x),5),];hgfx<-lm(xx$V1~xx$V2+xx$V3);
   
4.    data1[i,]<-hgfx$coefficients }

复制代码

稍微改了下  ，应该可以了吧，前提是你的x是数据框，且含有变量V1,V2,V3，否则得定义

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

1. chouyang<-function(x) {data1<-matrix(rep(0,3000),1000,3);
   
2.   for(i in 1:50)
   
3.   {xx<-x[sample(1:nrow(x),5),];hgfx<-lm(xx$V1~xx$V2+xx$V3);
   
4.    data1[i,]<-hgfx$coefficients }
   
5. data1[1:50,]
   
6. }

复制代码

我也要修改，
保留 (Intercept)  v2, v3

1. data <- data.frame(V1=1:1000,V2=sample(500,1000,replace=T),V3=sample(500,1000,replace=T))
   
2. 
   
3. chouyang <- function(x,n=50,m=5){
   
4.   hg=function(x){
   
5.     xx <- x[sample(1:nrow(x),m),]
   
6.     hgfx<-lm(V1~V2+V3,data=xx)
   
7.     hgfx$coefficients
   
8.   }
   
9. data1= replicate(n,expr=hg(x))
   
10. t(data1)
    
11. }
    
12. chouyang(data,50,5)

复制代码

jmpamao 发表于 2013-12-30 22:45
我也要修改，
保留 (Intercept)  v2, v3

修改要先站住沙发,啊哈哈

nuomin 发表于 2013-12-30 23:03
修改要先站住沙发,啊哈哈

更爱藤椅~

1. data <- data.frame(V1=1:1000,V2=sample(500,1000,replace=T),V3=sample(500,1000,replace=T))

复制代码

1. t(replicate(50, lm(V1~V2+V3,data=data[sample(1:nrow(data),5),])$coefficients))

复制代码

nuomin 发表于 2013-12-30 22:13
试出来了，改动一个地方就行
data1[[1]]

谢谢。问题找来了。

jmpamao 发表于 2013-12-30 22:15
感觉易出的地方：
1、data1

谢谢，我少了最后一个data1语句，没有输出语句。