在r中，如何根据地层中的观测值数量从3维列联表中删除地层？

Question

我正在用数据框中的3个变量构建3d联想表。让我们假设我正在通过

table(x,y,z)

其中z是我进行分层的变量。我希望去掉该层中观察值数量为1的任何(,,z(i))。

我该怎么做呢？首先，我很难弄清楚如何计算观测值，我认为我可以使用子集来削减我的偶然表。

Answer 1

假设您的数据包含在名为data的数据框对象中，此代码应通过一次观测删除图层中的所有数据。

data <- data[-which(data$z %in% which(table(data$z)==1)),]

编辑

这似乎现在起作用了。我不确定这是否会在一般情况下工作，但它在这种情况下是有效的。

data <- read.csv(file='~/Downloads/juveniles2forMax.csv')

data <- data.frame(
  Urban = data$Urban,
  RecidivismPlacement = data$RecidivismPlacement,
  timeinjj = data$timeinjj
)

removeStrata <- function(data, z) {
    data[-which(data[,z] %in% as.numeric(attr(which(table(data[,z])==1),"names"))),]
}

removeStrata(data=data, z='timeinjj')