r - Rに2つの変数を持つ3x3分割表を作成する方法

Question

例：

x <- c( 1, NA, 0, 1)
y <- c(NA, NA, 0, 1)
table(x,y, useNA="always") # --->
#       y
# x      0 1 <NA>
#   0    1 0    0
#   1    0 1    1
#  <NA>  0 0    1

私の質問は：

a <- c(NA, NA, NA, NA)
b <- c(1, 1, 1, 1)
table(a, b, useNA="always") ## --> It is 1X2 matrix. 
#       b
# a      1 <NA>
#   <NA> 4    0

上記の例と同じ列名、行名、ディメンションを持つ3X3テーブルを取得したいと思います。次に、テーブルにchisq.testを適用します。ご回答ありがとうございます！

Answer 1

aこれは、とbを同じレベルの因子に変換することで実現できます。これが機能するのは、因子ベクトルが、実際にはそれらのサブセットのみを含んでいる場合でも、要素が取る可能性のあるすべての可能な値（別名レベル）を追跡するためです。

a <- c(NA, NA, NA, NA)
b <- c(1, 1, 1, 1)

levs <- c(0, 1)

table(a = factor(a, levels = levs), 
      b = factor(b, levels = levs), 
      useNA = "always")
#       b
# a      0 1 <NA>
#   0    0 0    0
#   1    0 0    0
#   <NA> 0 4    0