我正在使用 ISLR 包进行统计实践。我正在使用 OJ 数据集。我正在尝试为 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
我正在尝试找到几率比并最终应用 McNemar 检验。但我无法生成列联表。我可以手动完成，但需要在 R 中完成。
【问题讨论】：
标签： r statistics data-science
您正在尝试使用 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 购买都发生在两者都没有打折的情况下，但可能是销售促进了对打折品牌的购买。同样，统计测试会告诉您这可能只是随机机会还是不太可能。
【讨论】：