在对线性模型应用对比度之前，解决组间差异不相等的问题？（r）
我的目标：我有一个顺序因子变量（5个级别），我想对其应用对比度来测试线性趋势。然而，因子组具有方差异质性 我所做的：根据推荐，我使用在对线性模型应用对比度之前，解决组间差异不相等的问题？（r）,r,regression,robustness,R,Regression,Robustness,我的目标：我有一个顺序因子变量（5个级别），我想对其应用对比度来测试线性趋势。然而，因子组具有方差异质性 我所做的：根据推荐，我使用robustpckg中的lmRob（）创建了一个稳健的线性模型，然后应用对比度 # assign the codes for a linear contrast of 5 groups, save as object contrast5 <- contr.poly(5) # set contrast property of sf1 to contain th
robustlmRob（）# assign the codes for a linear contrast of 5 groups, save as object
contrast5 <- contr.poly(5)
# set contrast property of sf1 to contain the weights
contrasts(SCI$sf1) <- contrast5
# fit and save a robust model (exhaustive instead of subsampling)
robmod.sf1 <- lmRob(ICECAP_A ~ sf1, data = SCI, nrep = Exhaustive)
summary.lmRob(robmod.sf1)
#为5组线性对比度分配代码，另存为对象
对比度5vcovHC
是处理异方差的正确函数。HC代表异方差一致估计量。这不会降低模型效应估计中的异常值，但会以不同的方式计算CIs和p值，以适应此类外围观测的影响lmRob
不降低外围值，不处理异方差
请参阅此处的更多信息：
vcovHC
是处理异方差的正确函数。HC代表异方差一致估计量。这不会降低模型效应估计中的异常值，但会以不同的方式计算CIs和p值，以适应此类外围观测的影响lmRob
不降低外围值，不处理异方差
请参阅此处的更多信息：
日志转换您的响应变量有帮助吗？日志转换您的响应变量有帮助吗？谢谢您，AdamO！我交叉发布了我的问题（我是新来的，不确定溢出和交换的汇款方式）。。。在使用vcovHC
后，我有一个新的担忧，即显著（出乎意料地）不同的p值，以及一个关于报告平均值/std DEV的问题。你能帮我一下吗？更新：谢谢你，阿达莫！我交叉发布了我的问题（我是新来的，不确定溢出和交换的汇款方式）。。。在使用vcovHC
后，我有一个新的担忧，即显著（出乎意料地）不同的p值，以及一个关于报告平均值/std DEV的问题。你能帮我一下吗？更新：