Я пытаюсь сопоставить контрольную группу и группу больных по возрасту (дискретно), LDH, образованию, CAPS и полу (категориально) в R. Есть ли простой способ сделать это? У меня 59 пациентов в моей больной группе и 336 в моем контроле.
dput(droplevels(df_cluster10_lh[1:10, ]))
Ouput:
structure(list(SUBJID = c("TRT_0001_recon", "TRT_0002_recon",
"TRT_0003_recon", "TRT_0004_recon", "TRT_0006_recon", "TRT_0008_recon",
"TRT_0010_recon", "TRT_0011_recon", "TRT_0012_recon", "TRT_0021_recon"
), Binge_Group_15 = c("Late", "Late", "Late", "Late", "Late",
"Early", "Early", "Late", "Early", "Late"), Age = c(23, 24, 42,
45, 29, 62, 44, 23, 45, 49), LDH_TOTKG = c(1897.741079, 1872.442119,
179.8455581, 4830.892644, 3189.005661, 897.6581303, 5815.605113,
3382.794477, 2768.912015, 696.3706047), CAPS_C = c(56, 82, 0,
11, 12, 7, 22, 85, 7, 4), BG_GEND = c("Male", "Male", "Female",
"Male", "Male", "Male", "Male", "Male", "Female", "Male"), BG_EDUC = c(12, 12, 14, 12, 12, 17, 12, 13, 13, 12)), row.names = c(NA, -10L), class = c("tbl_df",
"tbl", "data.frame"))
Я использовал R's match, но результат не выглядит соответствующим. Вот команда
library(MatchIt)
m.out <-matchit(Binge_Group_15 ~ Age + CAPS_C + BG_GEND + LDH_TOTKG, data = binge)
s.out <- summary(m.out, standardize= TRUE)
binge.match <- match.data(m.out)
head(m.out$match.matrix)
binge[12,]
binge[380,]
i <- rownames(binge.match) %in% c(12,380)
binge.match[i,]
m.out$match.matrix
Вывод, который я получил, был SUBJID Binge_Group_15 LDH_TOTAL Возраст LDH_TOTKG PTSD_C LifetimePTSD CAPS_C Binge_History Расстояние BG_GEND:
12 TRT_0021_recon 1 Мужской 0 681 0 068 0 068 0 068 0 0 6 0 6 0 6 0 6 0 6 0 6 0 6 0 6 0 6 0 6 0 6 0 6 0 6 0 6 0 6 0 6 0 0 0 6,0 0 0*
380 TRT_0529_recon 0 18336,6 33 3159,3119 1 1 64 1 Мужской 0,8304429
вес 12 = 1, 380 = 1