pid <- 1:13 sex <- as.factor(c(rep(1,6),rep(2,7))) levels(sex) <- c("男","女") disease <- as.factor(c(1,1,1,2,2,2,1,1,1,1,2,2,2)) levels(disease) <- c("有","無") print(ctab <- table(sex,disease)) mosaicplot(ctab,main="2×2クロス集計表のモザイクプロット例")