R twang 包的 mnps 函数显示绝对标准差异图的警告错误
mnps function of R twang package shows warning errors for Absolute Standardized Difference plot
我有一个包含 3465 个观察值和 5 个变量的数据集（命名为：数据集）
变量
$ 治疗 --> 因子（3 级：“1”、“2”、“3”）
$ 年龄 --> 整数（例如 24,54 等）
$死亡率 --> 整数（0=否，1=是）
$ LDL --> 整数（例如 200,120,143 等）
$ 吸烟 --> 整数（0=否，1=是）
我使用 twang 包中的 mnps 函数根据 treat
对我的观察进行加权
mnps 代码
mnps.dataset <- mnps(treat ~ Age + LDL + Smoking, data = dataset, estimand = "ATE", stop.method = c("es.mean", "ks.mean"), n.trees=3000)
```
＃＃＃ 阴谋 ＃＃＃
然后我用绘图函数来评估加权后数据的平衡性。
图 1 和图 2 没有问题（每次迭代处理的平衡，以及处理的倾向得分分布。）
剧情代码
plot(mnps.dataset, plots = 1)
plot(mnps.dataset, plots = 2)
但是当我想 运行 绘制 3 以获得绝对标准化差异时，它出现了一些警告和空白图的错误：
警告错误
plot(mnps.dataset, plots = 3)
警告消息：
1: In order(as.numeric(x)) : NAs introduced by coercion
2: In split(as.numeric(x) [ord], gg[ord]) : NAs introduced by coercion
3: ...(same as error 1)
4: ... (same as error 2)
5: ... (same as error 1)
6: ... (same as error 2)
7: ... (same as error 1)
8: ... (same as error 2)
9: In panel.superpose(x, y, type = type, groups = groups, ... : NAs introduced by coercion
10: ...
11: ...
12: ...
___________________
注意：“...”表示重复之前的错误
Image of empty plot generated by mnps plot=3
这是 twang 中的错误（我能够在自己的数据集上复制它）。您对此无能为力。
您可以改用 cobalt 包，它有助于平衡检查并与 twang 中的对象兼容。与类型3情节传达相同信息的情节是爱情情节，可以使用cobalt中的love.plot()调用。
如果您 运行 love.plot(mnps.dataset, which.treat = .all) 您可以看到所有成对比较的标准化均值差异。有关如何使用它的更多信息，请参阅 multi-category treatments 上的 cobalt 小插图部分。
我有一个包含 3465 个观察值和 5 个变量的数据集（命名为：数据集）
变量
$ 治疗 --> 因子（3 级：“1”、“2”、“3”）
$ 年龄 --> 整数（例如 24,54 等）
$死亡率 --> 整数（0=否，1=是）
$ LDL --> 整数（例如 200,120,143 等）
$ 吸烟 --> 整数（0=否，1=是）
我使用 twang 包中的 mnps 函数根据 treat
mnps 代码
mnps.dataset <- mnps(treat ~ Age + LDL + Smoking, data = dataset, estimand = "ATE", stop.method = c("es.mean", "ks.mean"), n.trees=3000)
``` ＃＃＃ 阴谋 ＃＃＃
然后我用绘图函数来评估加权后数据的平衡性。
图 1 和图 2 没有问题（每次迭代处理的平衡，以及处理的倾向得分分布。）
剧情代码
plot(mnps.dataset, plots = 1)
plot(mnps.dataset, plots = 2)
但是当我想 运行 绘制 3 以获得绝对标准化差异时，它出现了一些警告和空白图的错误：
警告错误
plot(mnps.dataset, plots = 3)
警告消息：
1: In order(as.numeric(x)) : NAs introduced by coercion
2: In split(as.numeric(x) [ord], gg[ord]) : NAs introduced by coercion
3: ...(same as error 1)
4: ... (same as error 2)
5: ... (same as error 1)
6: ... (same as error 2)
7: ... (same as error 1)
8: ... (same as error 2)
9: In panel.superpose(x, y, type = type, groups = groups, ... : NAs introduced by coercion
10: ...
11: ...
12: ...
___________________
注意：“...”表示重复之前的错误
Image of empty plot generated by mnps plot=3
这是 twang 中的错误（我能够在自己的数据集上复制它）。您对此无能为力。
您可以改用 cobalt 包，它有助于平衡检查并与 twang 中的对象兼容。与类型3情节传达相同信息的情节是爱情情节，可以使用cobalt中的love.plot()调用。
如果您 运行 love.plot(mnps.dataset, which.treat = .all) 您可以看到所有成对比较的标准化均值差异。有关如何使用它的更多信息，请参阅 multi-category treatments 上的 cobalt 小插图部分。