看了一篇文章
“point and interval estimates of partial population attributable risks in cohort studies: examples and software"

里面举了一个例子：


根据这个计算出multivariate attributable risks

结果如下：



这个crude model我理解怎么做， 但我不知道她们这个multivariate model （多因素分析）结果是如何得出的

她们自己上传的macro如下：





但完全看不懂她们的%par是怎么调用的，求大神指点。。。。。

我自己也有一个数据集需要做这个归因危险度分析，希望对某种药物是否导致跌倒进行归因危险度分析，有暴露组和非暴露组的数据，但希望这个归因危险度能考虑到年龄，性别，居住地区和收入水平以及服用药物的种类这些因素（这些因素包括在multivariate model中间）

希望大神们指导一下，万分感谢！！

文章原文见附件