我是这个论坛的新手,我刚刚在另一个领域提出了类似的问题。由于这不是问我问题的合适地点,因此我将其发布在这里。希望你能帮帮我!
我对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问题已解决。