R 效果图-固定效果

我正在使用MASS包对计数数据运行负二项回归模型。在我的1级（国家）中，大约有一半在所有观察结果中为0。因此，结果没有受试者内部的差异。 我想用R中的（effects）包绘制拟合值的效果。 然而，当使用固定效果模型时，仅在受试者内进行研究，我无法理解效果函数如何绘制拟合值。如果我只是尝试为固定效果模型绘制拟合值，它将平均化，并且无法显示效果（可能必须为一个“级别1”绘制）

我假设这是因为我必须为效果函数指定给定的.values或xlevels，以绘制特定级别1单元的效果，但我仍然无法使其工作

我可以很容易地为非固定效果模型绘制效果，但是如果有人能帮助我理解如何使效果函数在固定效果下工作，我将非常感激。谢谢

# Set Seed
set.seed(123)
# Create test data, z being the level 1 subjects
testdat <- data.frame(x=rnorm(n= 400), y=rep(c(2,0,0,0,1,2,5,0,0,0),40),    z=as.factor(rep(seq(1,10),40)))
#Load MASS and run fixed-effects
library(MASS)
nbmod <- glm.nb(y ~ x + z-1, data=testdat)
summary(nbmod)

Call:
glm.nb(formula = y ~ x + z - 1, data = testdat, init.theta = 42072.15381, 
    link = log)

Deviance Residuals: 
     Min        1Q    Median        3Q       Max  
0.00e+00  0.00e+00  0.00e+00  0.00e+00  2.98e-08  

Coefficients:
      Estimate Std. Error z value Pr(>|z|)    
x   -8.157e-17  5.308e-02   0.000        1    
z1   6.931e-01  1.125e-01   6.159 7.31e-10 ***
z2  -2.830e+01  1.342e+05   0.000        1    
z3  -2.830e+01  1.342e+05   0.000        1    
z4  -2.830e+01  1.342e+05   0.000        1    
z5  -1.563e-16  1.581e-01   0.000        1    
z6   6.931e-01  1.123e-01   6.170 6.83e-10 ***
z7   1.609e+00  7.072e-02  22.759  < 2e-16 ***
z8  -2.830e+01  1.342e+05   0.000        1    
z9  -2.830e+01  1.342e+05   0.000        1    
z10 -2.830e+01  1.342e+05   0.000        1    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for Negative Binomial(42072.15) family taken to be 1)
#Plot effect display for fixed effect
library(effects)
plot(effect(term = "x", mod = nbmod))
#Run and plot non-fixed effect model
nbmod2 <- glm.nb(y ~ x, data=testdat)
plot(effect(term = "x", mod = nbmod2))

#设置种子
种子集（123）
#创建测试数据，z为1级科目
testdat您需要查看您的模型（这就是我添加摘要输出的原因）。std错误在10^5范围内是病态拟合的标志。您不需要绘图建议。您需要统计教育。感谢您的洞察力。当兴趣在于单位内方差时，使用固定模型。固定效应“控制”对于单元（1级）内未观察到的影响。在这种情况下，neg.bin模型用于对城市内各个行政区的事件（即犯罪事件）进行建模。城市为1级和2级行政区。一些城市，如z2、z3等，在其行政区内没有任何犯罪，对平均数没有贡献。