【问题标题】:SAS CONTRAST: "Weighting" for linear combination of parameter estimates in PROC GLMSAS CONTRAST:PROC GLM 中参数估计线性组合的“加权”
【发布时间】:2014-02-13 02:50:53
【问题描述】:

背景:我有一个分类变量,X,有四个级别,我适合作为单独的虚拟变量。因此,共有三个虚拟变量表示 x=1、x=2、x=3(x=0 是基线)。

问题/问题:我想测试我的模型参数线性组合的重要性,例如:2*B1+2*B2+B3=0。

在 Stata 中,第一个问题可以在使用以下方法拟合模型后轻松完成:

test 2*B1 + 2*B2 + B3 = 0

现在,如果我想在 SAS 中使用 CONTRAST 语句为 PROC GLM 执行此操作,我知道我的“权重”(因为没有更好的术语)总和必须为 0。例如,如果,在不相关的示例,我想测试以下四个连续变量:C1 + C2 = C3 + C4,我的对比语句如下所示:

CONTRAST 'Contrast1' C1 0.5 C2 0.5 C3 -0.5 C4 -0.5

在这种情况下,每个变量应该如何加权是显而易见的。但是,当我想将上面模型中给出的系数 (2*B1 + 2*B2 + B3 = 0) 与这些权重相结合时,我不清楚如何在 CONTRAST 语句中对函数进行加权,特别是对于虚拟变量编码的分类变量,如问题中最初所述。

【问题讨论】:

    标签: testing sas contrast categorical-data


    【解决方案1】:

    使用 PROC REG。

    proc reg data=mydata;
    model y = b1 b2 b3;
    test 2*b1+2*b2+b3=0;
    run;
    quit;
    

    【讨论】:

    • 谢谢。我很难跟踪 SAS 的众多回归过程中的每一个都支持哪些语句。我没有想到 PROC REG,因为我之前曾尝试对参数的线性组合进行简单的 计算(而不是测试)(即使用 SAS 作为“计算器”和存储的 beta 估计值),就像 Stata 的 lincom 命令,而 PROC REG 不支持。这实际上也是我需要能够做到的,所以我会将它添加到我的问题中。
    • 没关系 - 我认为它更适合作为自己的问题。现在将发布新问题。
    • 附加的,有些相关的,问题:stackoverflow.com/questions/21758486/…
    • 为您的其他问题添加了答案。简短的回答是 PROC SCORE。
    • 谢谢。我已经评论了该问题的发布。
    猜你喜欢
    • 1970-01-01
    • 2013-04-23
    • 2017-05-23
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    相关资源
    最近更新 更多