【求助】有关R做面板数据分析的问题！

11426

收藏 2010-07-02

各位大侠：

在下请教一个问题！！！

我在运行R里的以下程序的时候出现了警告信息，为什么呢？？？

程序：
fixef(ar1, effect="individual")

警告信息：
> fixef(ar1, effect="individual")
错误于lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) :
外接函数调用时不能有NA/NaN/Inf(arg1)

idividual的系数没有算出来，为什么呢？有没有什么建议啊？

还有我的自变量和因变量里面有一些零，在用plm函数的时候非要删掉这些零值吗？还是有什么其他方法？？？

因为用了form<-log(MOP)~log(MWP)来指定模型。

谢谢！！！

扫码加我拉你入群

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

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

全部回复

DM小菜鸟

2015-2-5 20:44:40

在R语言中缺失值通常以NA表示，用is.na判断是否缺失值
然后填补缺失值，比如——
library(mice)
imp=mice(sleep,seed=1234)
fit=with(imp,lm(Dream~Span+Gest))
pooled=pool(fit)
summary(pooled)

在R语言中实现方法是使用mice包中的mice函数，生成多个完整数据集存在imp中，再对imp进行线性回归，最后用pool函数对回归结果进行汇总。汇总结果的前面部分和普通回归结果相似，nmis表示了变量中的缺失数据个数，fmi表示fraction of missing information，即由缺失数据贡献的变异