R:用于MI数据集的拆分工具和SD,拆分文件

时间:2018-07-19 09:40:10

标签: r pool imputation

我是这个论坛的新手,我刚刚在另一个领域提出了类似的问题。由于这不是问我问题的合适地点,因此我将其发布在这里。希望你能帮帮我!

我对R很陌生,并且一直在用R进行分析。我有一个多重插补数据集(mids对象),我想为其计算一些描述性数据。

更具体地说,我想根据特定变量(抑郁是/否)来划分均值和方差,并进行Mann-Whitney U检验以测试患有抑郁症的人的BMI是否比没有抑郁症的人更高

您知道如何在R中使用多个插补数据集吗?

由于使用以下语法对另一个问题的回答,我能够计算均值和方差(不使用拆分文件选项):

impF1lang <- complete(impF1, "long")
(means <- with(impF1lang, tapply(BMI, .imp, mean)))
mean(means)
(variances <- with(impF1lang, tapply(BMI, .imp, var)))
mean(variances)

我要用于拆分文件的变量的变量名是'PAIDhoog'。

这是我的数据集样本的外观(使用'complete')选项:

     .imp .id verrichtingen verpleegkanders Leeftijd Datuminvullen HbA1c  BMI
1       1   1             0               0       26    27-6-2016     69 26.7
2       1   2             0               0       69    22-7-2016     75 34.5
3       1   3             0               0       67    22-jul-16     62 24.3
4       1   4             6               0       38    22-7-2016     96 25.4
5       1   5             0               0       29    26-7-2016     63 19.1
6       1   6             0               0       50    16-7-2016     86 37.9
7       1   7             1               0       29    29-7-2016     44 29.1
3176    7 452             1               3       55    12-10-2016    50 36.3
3177    7 453             8               6       63    13-10-2016    36 29.3
3178    7 454             0               0       66    8-12-2016     86 28.3
     Type_Treat DurationDM          Socstatus PAID1 PAID2 PAID3 PAID4 PAID5
1       Insulin          5  wel achterstandsw     0     1     2     1     0
2       Insulin         17 geen achterstandsw     2     1     1     2     0
3       Insulin          1 geen achterstandsw     0     0     0     0     0
4       Insulin         10 geen achterstandsw     0     0     1     1     0
5       Insulin         25 geen achterstandsw     0     0     0     0     0
6          Both          9  wel achterstandsw     0     1     0     2     0
7          Both         33 geen achterstandsw     1     1     2     3     0
3176    Insulin         28 geen achterstandsw     2     2     3     3     3
3177    Insulin         12 geen achterstandsw     3     4     3     3     0
3178    Insulin         14 geen achterstandsw     0     0     2     3     0
     PAIDtot PHQ1 PHQ2 PHQ3 PHQ4 PHQ5 PHQ6 PHQ7 PHQ8 PHQ9 Hypoaantal Geslacht
1          4    0    0    0    1    1    0    0    0    0         20    vrouw
2          6    0    0    0    1    0    0    1    1    0         20      man
3          0    0    0    0    0    0    0    0    0    0          0      man
4          2    1    0    3    3    1    1    1    1    0         10    vrouw
5          0    1    1    3    0    0    0    3    0    0          0      man
6          3    2    0    3    1    1    1    0    0    0         10      man
7          7    1    1    3    3    0    0    2    0    0          0    vrouw
3176      13    0    1    1    1    3    0    1    1    0         10    vrouw
3177       9    0    2    3    1    0    3    1    3    2         20    vrouw
3178       5    0    0    3    1    0    0    0    0    0          0    vrouw
            Etnicit HAPOH Bedrijfsarts MW FysioErgo Diet Psychiat Psychol Dvk
1        Overigwest     0            0  0         0    5        0       0   5
2               Mar     0            1  4         0    2        0       0   2
3        Overigwest     2            3  0         5    0        0       0   1
4    Overignietwest     0            0  0         0    0        1       2  11
5        Overigwest     0            0  0         0    0        0       0   4
6              Turk     0            0  0         0    1        0       0   2
7        Overigwest     0            0  0         0    3        0       0   4
3176     Overigwest     2            0  0        26    1        0       6   4
3177 Overignietwest    36            0  4         6    3        0       0   6
3178     Overigwest     0            0  6        48    0        0       0   2
     VPtot Internist Specialist ICUopname Opnamegewoon SEH    Comorb DMtype
1        5         2          3         0            0   0      geen    DM1
2        2         3          8         0            0   0      geen    DM2
3        1         2          3         0            0   0      geen    DM1
4       11         6         25         0            1   0      geen    DM2
5        4         2          6         0            0   0      geen    DM1
6        2         2          0         0            0   0      geen    DM2
7        4         2          3         0            0   0      geen    DM2
3176     7         5          6         0            1   0         1    DM2
3177    12         2         22         0           13   2         1    DM1
3178     2         4          5         0            7   3 2 of meer    DM2
     PAIDtotaal PHQtotaal PAIDhoog PHQhoog kostenopnames kosteninternist
1             4         2     laag    laag             0             259
2             6         3     laag    laag             0             240
3             0         0     laag    laag             0             160
4             2        11     laag    laag           443             480
5             0         8     laag    laag             0             160
6             3         8     laag    laag             0             160
7             7        10     laag    laag             0             160
3176         13         8     hoog    laag           443             400
3177          9        15     hoog    hoog          5759             160
3178          5         4     laag    laag          3101             320
     kostenspecialist kostenhuisarts kostenMW kostenfysioergo kostendvk
1                 240             NA       NA              NA       240
2                 640             NA       NA              NA        96
3                 240             NA       NA              NA        48
4                2000             NA       NA              NA       528
5                 480             NA       NA              NA       192
6                   0             NA       NA              NA        96
7                 240             NA       NA              NA       192
3176              480             66        0             858       192
3177             1760           1188      260             198       288
3178              400              0      390            1584        96
     kostendietist totaalkosten jaarHAPOH jaarbedrijfsarts jaarMW jaarfysioergo
1              240          946        NA               NA     NA            NA
2               96          129        NA               NA     NA            NA
3                0        10711        NA               NA     NA            NA
4                0         4994        NA               NA     NA            NA
5                0         1073        NA               NA     NA            NA
6               48          545        NA               NA     NA            NA
7              144          946        NA               NA     NA            NA
3176            48         2967         8                0      0           104
3177           144        10275       144                0     16            24
3178             0         6668         0                0     24           192
     totaalverbruikjaar kostenHAjaar kostenMWjaar kostenjaarfysioergo
1                     9           NA           NA                  NA
2                     2           NA           NA                  NA
3                   174           NA           NA                  NA
4                    45           NA           NA                  NA
5                    11           NA           NA                  NA
6                     4           NA           NA                  NA
7                     8           NA           NA                  NA
3176                134          264            0                3432
3177                229         4752         1040                 792
3178                237            0         1560                6336
     kostenopnameICU kostenpsycholoog kostenpsychiater kostenvpanders
1                  0                0                0              0
2                  0                0                0              0
3                  0                0                0              0
4                  0              188               94              0
5                  0                0                0              0
6                  0                0                0              0
7                  0                0                0              0
3176               0              564                0            144
3177               0                0                0            288
3178               0                0                0              0
     kostenverrichtingen totaalutilization kostenseh totaalkostennieuw
1                      0                15         0               979
2                      0                30         0              1072
3                      0                81         0               448
4                    876                60         0              4609
5                      0                13         0               832
6                      0                 6         0               304
7                    146                14         0               882
3176                 146               139         0              6113
3177                1168               246       518             16669
3178                   0               237       777             12590
     contactprimarycare contactsecondarycare totaalhealthcareutilization
1                     1                   15                          NA
2                   184                   15                          NA
3                   108                    6                          NA
4                     1                   52                          NA
5                     1                   12                          NA
6                     1                    5                          NA
7                     1                   13                          NA
3176                112                   27                         139
3177                184                   62                         246
3178                216                   21                         237
     kostenprimarycare kostensecondarycare totaalkostenhealthcare etndich
1                 0.01                 979                     NA       1
2              1056.00                1072                     NA       2
3              1056.00                 448                     NA       1
4                 0.01                4609                     NA       2
5                 0.01                 832                     NA       1
6                 0.01                 304                     NA       2
7                 0.01                 882                     NA       1
3176           3696.00                2417                   6113       1
3177           6584.00               10085                  16669       2
3178           7896.00                4694                  12590       1
     comorbdich etnicitnieuw paidhoogimp phqhoogimp totkosthealthcare
1             1       Overig           0          0            979.01
2             1          Mar           0          0           2128.00
3             1       Overig           0          0           1504.00
4             1       Overig           0         NA           4609.01
5             1       Overig           0          0            832.01
6             1         Turk           0          0            304.01
7             1       Overig           0          0            882.01
3176          2       Overig           1          0           6113.00
3177          2       Overig           1         NA          16669.00
3178          2       Overig           0          0          12590.00
     tothealthcareutil util paidnieuw paiddich
1                   16   20         4        0
2                  199   21         6        0
3                  114   14         0        0
4                   53   63         2        0
5                   13   16         0        0
6                    6    7         3        0
7                   14   17         7        0
3176               139   62        13        1
3177               246  120        13        1
3178               237   77         5        0

我希望这个问题足够具体,如果需要,我当然可以提供一些额外的信息。

编辑:NA问题已解决。

0 个答案:

没有答案