Я использую пакет arsenal для использования функции tableBy.
У меня есть две переменные, по которым я хочу разделить: PROC и GROUP.
Мои данные выглядят так:
PATNO PROC LOS HC
1 20 LAP RNYGB 1 0
2 22 LAP RNYGB 1 0
3 23 LAP RNYGB 2 0
4 24 LAP RNYGB 2 0
5 26 LAP RNYGB 5 0
6 27 LAP RNYGB 1 0
Я хотел бы создать tableBy, который расслаивается как по процедуре, так и по HC, чтобы дать мне среднее значение LOS.Затем я хотел бы сравнить среднее значение LOS процедур непосредственно между группами HC (то есть mean(LOS) в LAP RNYGB существенно не отличается между HC=0 и HC=1, а mean(LOS) в ROBOTIC RNYGB значительно отличается между HC=0 и HC=1).
Я некоторое время возился с tableBy, но не могу понять.Это близко, поскольку я пришел (los_hc это подмножество los, где HC=0 и los_e, где HC=1):
los.hc.means <- tableby(PROC ~ LOS, data=los_hc)
los.e.means <- tableby(PROC ~ LOS, data=los_e)
los.means <- merge(los.hc.means,los.e.means)
summary(los.means, title='LOS by Procedure',text=TRUE)