如何执行mantel.haenszel卡方检验？

Question

我有一个数据集，其中的列是3种颜色的父鸟，分为三类（浅锈，中锈和深锈）。这些行是后代，分为这些相同的类。因此，一个光明的父亲生产了9个光亮的后代，14个中等的后代和6个黑暗的后代，等等。

在导入此列联表时，我执行了经典的卡方检验，以检查单元格之间的随机性。由于我的数据涉及鸟羽毛图案的编码颜色，因此与位于极端位置的像元相比，相邻的像元更容易被误码。例如，与将轻鸟称为深色鸟相比，我更可能将轻鸟称为中鸟。

因此，我想在我的表上执行mantel.haenszel测试，该测试应考虑这些基本关系，但它总是给我错误：

“ mantelhaen.test中的错误（dd $ Light，dd $ Medium，dd $ Dark，正确= TRUE ，：   每个阶层的样本量必须> 1“

我尝试了mantel.haenszel测试以及mantel.haenszel和add.cross.strata.test的组合，两者都给了我相同的错误。

mantelhaen.test（dd $ Light，dd $ Medium，dd $ Dark，                  正确= TRUE，精确=否，conf.level = 0.95）

以及

add.mantel.haenszel（dd，Z = NULL，保守= FALSE）

我希望能够或找到如何转换数据的答案，以便我可以对其进行测试，也欢迎提出其他更可靠的测试建议，我可以用它们来发现数据中的任何非随机效应。 / p>