我正在努力為以下情況創建列聯表:
“在一項病例對照研究中記錄了倫敦、曼徹斯特和紐卡斯爾的消化性潰瘍患者和對照患者的血型。血型 A 和 O 代表具有兩個水平的因子 B;提到的城市代表具有三個水平的因子 C。水平,L、M 和 N;U 代表具有兩個水平的因子,Control 和 Ulcer。在此病例對照研究中,將 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
有沒有辦法解決這個問題并創建準確的表格?
uj5u.com熱心網友回復:
您可以創建一個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
轉載請註明出處,本文鏈接:https://www.uj5u.com/gongcheng/354153.html
標籤:r
上一篇:將每行的值寫入每行x天前
下一篇:將表轉換為資料框R
