时间戳

首页 > 解决方案 > R 相当于 proc genmod 中的 SAS 'estimate' 语句

r - R 相当于 proc genmod 中的 SAS 'estimate' 语句

问题描述

我正在尝试将我同事的一些 SAS 代码翻译成 R，以便我可以在我的 Rmarkdown 报告中使用它。我能够在 R 中拟合相同的模型，但不知道如何在estimate此处获取语句：

proc genmod data= smydata descending;
    class  id_var outcome outcome_name treatment1 treatment2;
    model outcome = outcome_name treatment1*outcome_name treatment2*outcome_name/ type3 dist=bin link=logit  noint covb;
    weight gbm_weights;
    repeated subject=id_var/ type=un corrw;

%let L= 0.1428 0.1428 0.1428 0.1428 0.1428 0.1428 0.1428;
estimate "6.2. H0: avg treatment1 effect =0"  treatment1*outcome_name  &L /e exp alpha = 0.0166;

这里treatment1和outcome是二进制的。outcome_name有7个级别。测试的目的是查看treatment1和之间的平均交互项outcome_name是否非零（或取幂时不同于 1）。

我知道如果没有数据，这有点难以做到，但我很感激任何建议。

谢谢！

标签： rsas

解决方案

推荐阅读