Блог

Я изо всех сил пытаюсь создать таблицу непредвиденных обстоятельств для следующей ситуации:

«Группы крови пациентов с язвенной болезнью и контрольной группой в Лондоне, Манчестере и Ньюкасле были зарегистрированы в исследовании «случай-контроль». Группы крови A и O были представлены в факторе B с двумя уровнями; упомянутые города были представлены как фактор C с тремя уровнями, L, M и N; и U представлял фактор с двумя уровнями, Контроль и Язва. В этом исследовании «случай-контроль» целесообразно рассматривать B как фактор ответа, а C и U как ковариационные факторы». Векторы задаются в виде

 B=c("A","A","A","A","A","A","O","O","O","O","O","O") C=c("L","L","M","M","N","N","L","L","M","M","N","N") U=c("C","U","C","U","C","U","C","U","C","U","C","U") count=c(4219,579,3775,246,5261,219,4578,911,4532,361,6598,396)  

Я попытался создать таблицу, используя код

 three=xtabs(~B C U,data=count) ftable(three)  

Но получите сообщение об ошибке

 Error in eval(predvars, data, env) :   numeric 'envir' arg not of length one  

Есть ли способ устранить эту проблему и создать точную таблицу?

Вы можете создать data.frame и использовать затем xtabs :

 xtabs(count~B C U, data.frame(B, C, U, count)) #, , U = C # # C #B L M N # A 4219 3775 5261 # O 4578 4532 6598 # #, , U = U # # C #B L M N # A 579 246 219 # O 911 361 396