【问题标题】:SPSS syntax of a quadratic term with interaction具有交互作用的二次项的 SPSS 语法
【发布时间】:2014-05-27 13:34:07
【问题描述】:

SPSS 中具有二次项和交互的回归的语法如何?在 R 中,代码为:

fit <- lm(c ~ a*b + a*I(b^2), dat)

fit <- lm(c ~ a*(b+I(b^2), dat)

感谢您的帮助。

【问题讨论】:

    标签: regression spss interaction quadratic


    【解决方案1】:

    使用REGRESSION,您需要在提交命令之前实际在SPSS数据文件中制作变量。因此,如果您的变量命名相同:

    COMPUTE ab = a*b. /*Interaction*/.
    COMPUTE bsq = b**2. /*squared term*/.
    COMPUTE absq = a*bsq. /*Interaction with squared term*/.
    

    然后这些可以放在回归方程的右侧。

    REGRESSION VARIABLES=a,b,bsq,absq,c
     /DEPENDENT=c
     /METHOD=ENTER a,b,bsq,absq.
    

    我以为你只能为交互作用做因子变量——但我错了,你也可以做连续变量(对不起!)。这是一个使用MIXED 的示例(如果使用REGRESSION,您仍然需要创建单独的变量)。

    INPUT PROGRAM.
    LOOP Case = 1 TO 200000.
    END CASE.
    END LOOP.
    END FILE.
    END INPUT PROGRAM.
    
    COMPUTE a = RV.BERNOULLI(0.5).
    COMPUTE b = RV.NORMAL(0,1).
    COMPUTE ab = a*b /*Interaction*/.
    COMPUTE bsq = b**2 /*squared term*/.
    COMPUTE absq = a*bsq /*Interaction with squared term*/.
    COMPUTE c = 0.5 + 0.2*a + 0.1*b -0.05*ab + .03*bsq -.001*absq + RV.NORMAL(0,1).
    
    VARIABLE LEVEL a (NOMINAL).
    RECODE a (0 = 2)(ELSE = COPY).
    MIXED c BY a WITH b bsq
      /FIXED = a b b*b a*b
      /PRINT SOLUTION.
    

    【讨论】:

    • 感谢您的回复。您是否有可能在 ENTER 中错过了 ab?您能给我建议混合模型的语法(SPSS:MIXED)吗?它看起来像:/MIXED c by a with b bsq absq /FIXED a b bsq ab absq?是否还有其他选择,例如:/MIXED c by a with b bsq \FIXED c by a with b bsq a*bsq?甚至:/MIXED a b /FIXED a b b*b a*b a*b*b?
    • 是的@giordano 我忘了在REGRESSION 命令中包含ab。是的,在MIXED(和其他GENLIN 命令)中,您可以指定因子 变量的交互,但不能指定连续变量。因此,鉴于您的示例,b*b 绝对不是有效的语法。 (只有BY 之后的变量才能估计这样的交互作用。)
    • 再次感谢。我纠正了你的答案,我希望这没问题。那么,前两个语法正确,第三个错误?也就是说,ab(分类numerical)是可能的但不是numericalnumerical(我的意图是平方b)?还有... b2 ab2 错了吗?抱歉这些问题,但是在网上很难找到SPSS编码语言的语法。
    • @giordano - 抱歉我错了,显然可以在使用 MIXED 时指定连续交互(我认为您也可以为 GENLIN 指定,尽管我还没有确定)。请参阅添加的示例。
    • @Andi。谢谢你的模拟。我还是很困惑。在新语法的倒数第二行:/FIXED = a b b*b a*b。 bb 与 bsq 相同吗?并且没有缺少 absq(或 abb 或 ab2)吗?我将 bb 解释为 b2。
    猜你喜欢
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 2022-09-09
    • 1970-01-01
    • 2021-09-20
    • 1970-01-01
    • 1970-01-01
    相关资源
    最近更新 更多