问题描述
我想命名残差,因为我想将它与一阶差分模型进行比较。但是函数“plm”对名称进行排序,并且来自 plm 的残差 1 与来自 lm 的残差 1 不同。您可以推荐我哪种更优雅的比较两种残差的方法?
谢谢。
city <- c("B","D","E","C","A")
date <- rep(1,5)
value <- c(600,120,300,150,270)
count <- c(40,15,30,20,25)
city2 <- rep(city,each = 2)
date2 <- rep(1:2,5)
value2 <- c(600,650,350,200,270,300)
count2 <- c(40,36,10,27,25,22)
ex1 <- data.frame(city,date,value,count)
ex2 <- data.frame(city2,date2,value2,count2)
names(ex2) <- c("city","date","value","count")
fit.lm <- lm(count ~ value,data = ex1)
names(fit.lm$residuals) <- city
summary(fit.lm)
fit.plm <- plm(count ~ value,data = ex2,model = "fd"))
names(fit.plm$residuals) <- sort(city)
summary(fit.plm)
# attr(fit.plm$model,"index") #show order
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)