我对R非常陌生,所以如果有一种方法可以使用现有代码/包的轻微变化,我会道歉。
我使用predict命令为129个国家创建了变量(学生注册)的年度预测,然后我将它们绑定。我这样做是因为我预测使用多元回归。
以下是我正在做的事情(如果有帮助的话)
`fm1=lm(log(y+1)~Var.Ind)
XNew=data.frame(Var.Ind)
(rse<-summary(fit)$sigma(fm1)* df.residual(fm1))/2
rse<-summary(fm1)$sigma
yhat1=exp(predict(fm1,XNew)+rse*rse/2)-1
pos2014=which(Var.Ind[,1]==c(2014))
Var.Ind.2015=model.matrix(~as.matrix(Imp.Data4[pos2014,-2])-1)
head(Var.Ind.2015)
Var.Ind.2015=data.frame(Var.Ind.2015)
Var.Ind.2015.Ord=as.data.frame(Var.Ind.2015[order(Var.Ind.2015[,3],Var.Ind.2015[,1]), ])
head(Var.Ind.2015.Ord)
X.New.New=data.frame(cbind(model.matrix(~as.matrix(Var.Ind.2015.Ord))))
head(X.New.New)
ColNames.N=ColNames[-2]
colnames(X.New.New)=c("Int",ColNames.N,"Lag1","Lag2")
head(X.New.New)
Beta.Coef=matrix(as.numeric(fm1$coefficients),ncol=1)
Beta.Coef
Pred2015=as.data.frame(cbind(X.New.New[,3],exp(as.matrix(X.New.New)%*%Beta.Coef+rse*rse/2)-1))
dim(Pred2015)
colnames(Pred2015)=c("country","Yhat")
*以后的几年直到2030年)
cbind(Pred2015, Pred2016, Pred2017, Pred2018, Pred2019)`
我需要弄清楚是否有办法理解这些结果:
a)如何将预测结果导出为excel b)或者,如果我可以使用R将这些结果放入表中。
此外,这些结果不会出现在全局环境中,仅出现在程序的结果部分中,这就是为什么我不询问如何导出数据,而是询问这些特定结果。
如前所述,我的编码知识仅限于我1周的R经验(我通常使用STATA)。
任何帮助将不胜感激!