我有一个结果变量x
和三个解释变量a, b, c
,它们是分类变量。在我的示例a
中,有 8 个级别、b
4 个级别和c
35 个级别,但并非三个变量的所有组合都有观察值(这可能不重要)。
如果我在 Stata 中运行以下加法方差分析模型
anova x a b c
adjust, by(a b) gen(y)
然后我获得由变量和x
调整的变量的预测。调整命令在结果窗口中输出下表,并且它还生成一个带有调整预测的变量。a
b
y
| b
a | 2 4 8 16
----------+-----------------------------------
50 | .016655 .018487
75 | .008286 .011237
100 | .005937 .006677 .012467
150 | .001905 .004038 .009454
200 | .001774 .003107 .007592 .010081
400 | .004982 .006853 .009342
800 | .002126 .00521
1000 | .002732 .005221
----------------------------------------------
Key: Linear Prediction
我的问题是变量对于andy
的每个组合都有一个值,而上表只有 and 的每个组合的值。如何保存表格中的结果,以便可以使用这些结果?表中的值与表中的值有什么联系?a, b
c
a
b
y
提前致谢。
更新:我发现这个help adjust
:
估计命令中使用但未包含在 by() 变量列表或调整变量列表中的变量将保留其当前值,通过观察进行观察。这里 adjust 显示平均估计预测(或相应的概率或指数预测),用 by() 选项中的变量定义的每个组内替换这些未指定变量的平均值。
我的数据也是如此。例如 ifa=75
和b=2
, thenc
取值 12,13,14,15,16。y
对应的值c=14
(即平均值)正是表中显示的值。但是,如果这些值的平均值不是它所采用的值怎么办?