网状meta回归求教

Question

各位大神好！本人在用Stata16进行网状meta回归，参考的是15年的文献《基于频率学框架的网状Meta回归实现与应用》，在进行方差协方差矩阵的构建过程中，总是报错too many variables specified r(103);请教各位大神如何解决
代码如下：
foreach trt in A B C D E F {
  if " `trt' "=="A" continue
  gen y`trt' = log(r`trt'/(n`trt'-r`trt')) - log(rA/(nA-rA))
  gen S`trt'`trt' = 1/r`trt' + 1/(n`trt'-r`trt') + 1/rA + 1/(nA-rA)
  foreach trt2 in A B C D E F {
    if " `trt2' "=="A" continue
    if " `trt2' ">" `trt' " gen S`trt' `trt2' = 1/rA + 1/(nA-rA) if !mi(r`trt') & !mi(r`trt2')
}
}
format y* S* %6.2g
mat P7=0.5*(I(7)+J(7,7,1))

Answer 1

根据您提供的代码和错误信息，我们可以看到在构建 S 矩阵的过程中，变量的数量太多了，超过了 Stata 可以处理的限制，导致了错误。

解决这个问题的方法是将 S 矩阵的构建分成两部分，先只构建对角线上的元素，再分别构建上三角和下三角的元素。可以参考以下代码：

// 构建对角线上的元素
foreach trt in A B C D E F {
    if " `trt' "=="A" continue
    gen S`trt'`trt' = 1/r`trt' + 1/(n`trt'-r`trt') + 1/rA + 1/(nA-rA)
}

// 构建上三角和下三角的元素
foreach trt in A B C D E F {
    if " `trt' "=="A" continue
    foreach trt2 in A B C D E F {
        if " `trt2' "=="A" continue
        if " `trt2' ">" `trt' " {
            gen S`trt'`trt2' = 1/rA + 1/(nA-rA) if !mi(r`trt') & !mi(r`trt2')
        }
        else if " `trt2' "<" `trt' " {
            replace S`trt'`trt2' = S`trt2'`trt' if missing(S`trt'`trt2')
        }
    }
}
```bash

```

以上代码将对角线上的元素和上下三角的元素分别构建，可以避免一次构建太多变量的问题。如果还有其他问题，请随时提问。