怎么合并两个矩阵

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

x=c(1:10)
y=c(101:110)
x1=sanmple(x,10)##抽样
y1=sample(y,10)##抽样
fit1=lm(y1~x1)##建立OLS
x2=sample(x,10)##抽样
y2=sample(y,10)##抽样
fit2=lm(y1~x1)#建立OLS
问题：将fit1和fit2预测值进行去均值？难点：由于每次是重抽，不能直接将fit1和fit2直接取均值。所以需要将fit1和fit2预测值按照x进行标识列出，然后再取均值。
谢谢

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Sample AMPL fit PLE MPL

x=c(1:10)
y=c(101:110)
x1=sample(x,10)##抽样
y1=sample(y,10)##抽样
fit1=lm(y1~x1)##建立OLS
#读取列名称
coef1<-coef(summary(fit1))
write.csv(coef1,"coef1.csv")
coef1<-read.csv("coef1.csv",header=T)
#修改以示区别
names<-names(coef1)
names(coef1)<-c("Variable",paste(names[2:5],"1",sep=""))

x2=sample(x,10)##抽样
y2=sample(y,10)##抽样
fit2=lm(y1~x1)#建立OLS
##读取列名称
coef2<-coef(summary(fit2))
write.csv(coef2,"coef2.csv")
coef2<-read.csv("coef2.csv",header=T)
#修改以示区别
names<-names(coef2)
names(coef2)<-c("Variable",paste(names[2:5],"2",sep=""))

#合并两个估计参数表
coef<-merge(coef1,coef2,by="Variable")
#对估计参数取均值
coef<-within(coef,
             {Estimate_Mean<-(Estimate1+Estimate2)/2
             })
write.csv(coef,"coef.csv")

楼主：看看我的程序以及数据包。 Yinuo9257.rar (5.83 KB)

2013-10-21 10:28:28 上传

nieqiang110 发表于 2013-10-21 10:30
楼主：看看我的程序以及数据包。

谢谢楼主。不过预测是预测值，不是参数。

哈哈哈哈。具体一点，你需要把结果中的什么求均值

nieqiang110 发表于 2013-10-21 10:37
哈哈哈哈。具体一点，你需要把结果中的什么求均值

每一个OLS模型都对y有预测值，现在对两个y预测值求均值。

那是你的OLS的知识决定的，预测值我看过一些，你直接找本回归分析的书，把程序抄下来。按照我上面的方法，就可以求均值了。

nieqiang110 发表于 2013-10-21 10:42
那是你的OLS的知识决定的，预测值我看过一些，你直接找本回归分析的书，把程序抄下来。按照我上面的方法，就 ...

谢谢。

If I understand correctly you would like to predict your original x with the model built on your resampled data. If that is the case, simply using the "predict" function. You do not need even need to worry about the sequence of your predictions. The only thing you need to keep in mind is that the predict function requires your independent variable(s) as a data.frame, which your "x" is not:

1. (predict(fit1, data.frame(x = x)) + predict(fit2, data.frame(x = x)))/2

复制代码

zhangyangsmith 发表于 2013-10-25 04:33
If I understand correctly you would like to predict your original x with the model built on your res ...

谢谢！