如何用sas中的proc混合计算回归系数？

Question

这是我的数据.数据的结构如下:id x1 x2 x3 y.

我曾经proc mixed分析它,但现在想确定回归系数,我不知道该怎么做.我只是一个初学者sas.从结果我发现x1,x2,x3和x1x2x3是显著的影响,但是如何确定系数alpha, beta, gamma, delta, theta:

y = theta + alpha*x1 + beta*x2 + gamma*x3 + delta*x1*x2*x3

这是我的代码:

ods graphics on;
proc mixed data=test;
  class x1 x2 x3;
  model y = x1 | x2 | x3 / solution residual;
  random id;
run;
ods graphics off;

编辑1:这是表格的一部分Solutions for Fixed Effects:

由于x1有两个级别,因此表中有两行.x1通过总结这两个值-109.07来获得效果:对于第一行和0第二行,还是应该做其他事情？请注意,这是2^k设计.x1应该将效果计算为ywhen x1为高(20)和低(10)时的平均值之差的一半.

Answer 1

该solution选项应该生成您的估计。您需要将其包含在model和random报表中。您应该看到两个表，固定效应的解决方案和随机效应的解决方案，其中保存了估计值。

proc mixed data=test;
class x1 x2 x3;
model y = x1 | x2 | x3 / solution residual;
random id / s;
run;

文档中的随机系数示例与您的问题很接近。