“如何修复”;未定义的列已选定";用于R中的网络元分析?

“如何修复”;未定义的列已选定";用于R中的网络元分析?,r,networking,R,Networking,我正在用两个软件包gemtc和rjags对R进行网络元分析。然而,当我打字时 Model <- mtc.model (network, linearmodel=’fixed’). 阅读手册怎么样??mtc.model。它明确规定如下: 所需列[响应者,样本大小] 所以你的响应者变量应该是响应者,你的样本大小变量应该是样本大小 接下来,您的绘图(网络)将帮助您确定无法进行某些比较。在您的数据中,有两组试验进行了比较。治疗18和21未与任何其他治疗进行比较。因此,您只能对21和18进行元分

我正在用两个软件包gemtc和rjags对R进行网络元分析。然而,当我打字时

Model <- mtc.model (network, linearmodel=’fixed’). 

阅读手册怎么样?
?mtc.model
。它明确规定如下:

所需列[响应者,样本大小]

所以你的
响应者
变量应该是
响应者
,你的
样本大小
变量应该是
样本大小

接下来,您的
绘图(网络)
将帮助您确定无法进行某些比较。在您的数据中,有两组试验进行了比较。治疗18和21未与任何其他治疗进行比较。因此,您只能对21和18进行元分析,或对其余部分进行网络元分析

network <- mtc.network(data.ab=SAE[!SAE$treatment %in% c(21, 18), ])
model.fe <- mtc.model(network, linearModel="fixed")
网络
    study treatment responder sample.size
1      1         3         0          76
2      1        30         2          72
3      2         3        99        1389
4      2        23       132        1383
5      3         1         6         352
6      3        30         2         178
7      4         2         6         106
8      4        30         3          95
9      5         3        49         393
10     5        25        18         198
11     6         1        20          65
12     6        22        10          26
13     7         1         1          76
14     7        30         3          76
15     8         3         7         441
16     8        26         1         220
17     9         2         1          47
18     9        30         0          41
19    10         3        10         156
20    10        30         9         150
21    11         1         4          85
22    11        25         5          85
23    11        30         4          84
24    12         3         6         152
25    12        30         5         160
26    13        18         4         158
27    13        21         8         158
28    14         1         3         110
29    14        30         2         111
30    15         3         3          83
31    15        30         1          92
32    16         1         3         124
33    16        22         6         123
34    16        30         4         125
35    17         3       236        1553
36    17        23       254        1546
37    18         6         5         398
38    18         7         6         403
39    19         1        64         588
40    19        22        73         584
network <- mtc.network(data.ab=SAE[!SAE$treatment %in% c(21, 18), ])
model.fe <- mtc.model(network, linearModel="fixed")