Включение значений NA в structable ()

Question

Я хотел бы создать таблицу сопряженности, используя структурируемую функцию пакета vcd. Два из моих выбранных столбцов (consent_a и consent_b) имеют значения NA, а также значения факторов (Да, Нет), потому что каждый случай может быть согласован для процедуры a или процедуры b, но не для обоих. Например, если дело разрешено для процедуры a, их не запрашивают для процедуры b (и, таким образом, consent_b будет NA). В таблицу непредвиденных обстоятельств я хочу включить все случаи, когда коэффициент - «да», «нет» и «NA».

library(vcd)
mydata<-data.frame(report_year= c(2013, 2013, 2013, 2013, 2013, 2014, 2014, 2014, 2014, 2014),
report_week = c(1, 1, 1,  1, 1, 2, 2, 2, 2, 2),
consenta = c("Yes", "Yes", NA, "Yes", "No", "Yes", "Yes", NA, "Yes", "No"),
consentb=c(NA, NA, NA, NA, "Yes", NA, NA, "Yes", NA, "No"))
epicurve<-as.data.frame(structable(proj11[,c("report_epiweek", "report_year", "consent_a", "consent_b")]))

Answer 1

Я не знаком с пакетом vcd, но, насколько я вижу, нет встроенных средств для сохранения значений NA при работе с таблицами сопряженности.Тем не менее, вы можете получить тот же самый плоский стол, сохраняя NA s с функцией base table().

as.data.frame(with(mydata, table(report_week, report_year, consent_a, consent_b, useNA = "ifany")))

Если результат этого не подходит для ваших следующих шагов, тогда вы можете просто заменить NA s в вашем наборе данных на символическое значение и продолжить использовать пакет vcd.