R语言做固定效应模型

代码小白求教
需要用R studio做固定效应模型，想分析失业率人口数和犯罪率对宠物弃养的影响。
网上教程攻略也找了不少，可是实在基础太差看不懂，时间又比较紧
求大神给个代码，如果能讲得详细点最好不过了，刚注册的号论坛币不多，以后可以追加

自己看着网上的提示写的卡在了这一步
install.packages("ggplot2")
getwd()
read.csv("E:/NTHU/PROJECT/data1.csv",header=TRUE)
install.packages("plm")
library("plm")
mydata <- read.csv("E:/NTHU/PROJECT/data1.csv",header=TRUE)
mydata<-pdata.frame(mydata, index = c("city", "date"))
fixed<-pdata.frame(explained~unemployment+population+crime,data=mydata,index = c("country","date"), model="within")      
然后就提示

Error in pdata.frame(explained ~ unemployment + population + crime, data = mydata,  :   unused arguments (data = mydata, model = "within")

下面是老师给的提示

先将每县市的收容数数据除以该县市的最大收容量, 做为被解释变量.

解释变量 我目前想到的大概有各县的失业率(衡量景气好坏), 人口数(城市规模), 犯罪率(治安状况), 数据在总体经济数据库,  先利用这些变量做一下 fixed effect 模型, 看看结果如何. 但由于数据不算太多, 所以解释变量也不能太多

Panel Data Econometrics in R: The plm Package

https://cran.r-project.org/web/packages/plm/vignettes/plm.pdf


如果这可以帮助到你，请帮我评分，评分选项在我回复页面的右下方。

按照你的需求，建模那一步用到的主函数应该是`plm()`，而非导致错误的`pdata.frame()`，修改一下试试。

楼上说得对，同时数据本身已经是paneldata，所以不需要index了，即最后一步改为

1. fixed<-plm(explained~unemployment+population+crime,data=mydata, model="within")

复制代码

另外面板模型从计量的严谨性要进行hausman检验，

1.     random<-plm(explained~unemployment+population+crime,data=mydata,model="random")
   
2. phtest(fixed, random) #P值小于0.05选择固定效应模型

复制代码

风向南吹TP 发表于 2018-11-11 12:38
楼上说得对，同时数据本身已经是paneldata，所以不需要index了，即最后一步改为
另外面板模型从计量的严谨 ...

如果用去中心化方法做ols 多元回归该怎么写代码呢，我去写了每个变量减均值，但是跑不出来

最后一步好像有问题，应该是用plm（XXX~XXX+XXX），如果要双维度固定的话还要在里面加上effect="twoways"