我想创建一个函数,它接受 1 到 10 个向量并返回是或否,这意味着我的所有向量都代表一个“完整设计”。
这是我称之为“完整设计”的示例:
a <- c(1,1,1,1,2,2,2,2,3,3,3,3)
b <- c(1,2,1,2,1,2,1,2,1,2,1,2)
c <- c(1,1,2,2,1,1,2,2,1,1,2,2)
这是一个完整的设计,因为:对于所有向量的所有级别,(在同一位置)所有其他向量的所有级别的次数相同。
这是“不完整设计”的两个示例:(在以下两个示例中,a 和 b 匹配,但 c 不匹配 a 和 b)
示例 1:
a <- c(1,1,1,1,2,2,2,2,3,3,3,3)
b <- c(1,2,1,2,1,2,1,2,1,2,1,2)
c <- c(1,2,3,1,2,3,1,2,3,1,2,3)
示例 2:
a <- c(1,1,1,1,2,2,2,2,3,3,3,3)
b <- c(1,2,1,2,1,2,1,2,1,2,1,2)
c <- c(1,2,3,4,5,1,2,3,4,5,1,2)
希望我清楚。整个想法是我有一个数据集,这个数据集由 a、b、c、d、e 等因素解释......我想要一个函数来告诉我我是在测试一个完整的还是一个在其上运行 aov() 之前的设计不完整。
非常感谢 !
在我要创建的函数中要问的问题类似于:当等于给定级别时(例如 2)。我们查看 b[which(a==2)] 并检查 b 的所有级别是否包含在 b[which(a==2)] 中,并检查 b[which(a==2)] 的所有级别是否包含重复相同的次数。