R语言零膨胀模型的随机效应怎么做，用什么函数？

Python016

R语言零膨胀模型的随机效应怎么做，用什么函数？,第1张

面板数据确定采用固定效应还是随机效应需要做hausman test(豪斯曼检验)。过程是，先对面板数据做随机性检验，在结果窗口的PROC菜单下选择hausman test就可以了，检验的原假设是应该采用随机效应，备则假设是固定效应。

F检验是判断方程是混合方程、变截距方程还是变系数方程，需要求解三个方程的残差平方和。豪斯曼检验是确定常数项的固定或者随机效应的。两个是不同的问题。

具体方法是：用EVIEWS先对回归方程做混合模型求解，在结果中有一项Sum squared resid(在结果的下面，R平方值的旁边)，这个就是残差平方和，这个值就是S3；然后在用变截距模型求解，得出S3，最后是变系数模型，得出S1。有了这三个值，F值自己手算就可以了。

我是用的pscl包，zeroinfl()函数

零膨胀负二项模型（ZINB）

mod <- zeroinfl(ReportedNumber~ A+B+C+D+E | F+G+H+I, data = zinb, dist = "negbin", EM = TRUE)

ZINB模型由点模型和零膨胀模型两部分结合而成，ABCDE是点模型内变量， 影响因变量发生次数的多少，FGHI是零膨胀模型内变量，决定因变量是否能够发生（为0还是非0）。

http://www.ats.ucla.edu/stat/r/dae/zinbreg.htm 这个网站里讲的很清楚

r语言数据分析是查看数据的结构、类型，数据处理。根据查询相关资料信息显示：R语言是一个开源、跨平台的科学计算和统计分析软件包，具有丰富多样、强大的的统计功能和数据分析功能，数据可视化可以绘制直方图、箱型图、小提琴图等展示分数的分布情况可以通过散点图和线性拟合来展示分数和年龄之间的关系。