在R中创建列联表的问题

Question

我正在使用ISLR包进行我的统计练习。我使用的是OJ数据集。我正在尝试为weach of population的Purchase列和specialPrice列创建一个联想表。

我正在尝试找出如果有特价，CH被出售的可能性。

到目前为止，这是我的代码。

library(ISLR)
CH <- table(OJ[OJ$Purchase == 'CH', "SpecialCH"])
MM <- table(OJ[OJ$Purchase == 'MM', "SpecialMM"])
table (MM, CH)

我得到的输出有点奇怪。

     CH
MM    121 532
  101   1   0
  316   0   1

我试着找出赔率比，并最终应用麦克内马尔检验。但是我无法生成列联表。我可以手工完成，但需要用R来完成。

Answer 1

您正在尝试处理3个变量，但列联表仅使用2个变量。我建议使用xtabs，因为公式方法省去了一些输入，而且它在标记表格方面做得更好：

xtabs(~SpecialMM+SpecialCH, OJ)  # Only 4 weeks are both on special
#          SpecialCH
# SpecialMM   0   1
#         0 743 154
#         1 169   4 
xtabs(~Purchase+SpecialCH, OJ)   # When CH is on special ca 75% CH
#         SpecialCH
# Purchase   0   1
#       CH 532 121
#       MM 380  37
# xtabs(~Purchase+SpecialMM, OJ)  # When MM is on special ca 58% MM
#         SpecialMM
# Purchase   0   1
#       CH 581  72
#       MM 316 101

第一个表提出了这个问题。是一个品牌与另一个品牌相关联的特价商品。有1070个OJ的购买代表。CH特效158次，MM特效173次。但只有4次是两个品牌的特价。此表表明MM和CH不是同时处于特殊状态。您可以使用卡方或其他测试来查看这是否与随机分配的特色菜有显着偏差。

第二个和第三个表查看OJ的购买情况，以查看一个品牌在打折时是否比另一个品牌更有可能被购买。请注意，大多数OJ购买都发生在两个都不打折的时候，但可能是销售促进了对打折品牌的购买。再一次，统计测试会告诉你这是随机的机会还是不太可能的机会。