1、查看数据，以R自带的数据集HairEyeColor为例。ftable(HairEyeColor)，表示把分类数据按列联表的形式显示。从数据中可以看到该数据集有三个分类变量，即：Hair、Eye、Sex。
3、同样，利用vcd包中的函数mosaic绘制马赛克图。添加选项shade = TRUE将根据拟合模型的皮尔逊残差值对图形上色。添加选项legend = TRUE将展示残差的图例。代码如下：library(vcd)mosaic(HairEyeColor,shade=TRUE,legend=TRUE,color=TRUE)
5、检验性别与头发和眼睛颜色联合作用的独立性。结果表明：男性在棕色头发和眼睛的人群中很少，在棕色头发和蓝色眼睛的人群中占了过多的比例。fm <- loglin(HairEyeColor, list(1:2, 3))pchisq(fm$pearson, fm$df, lower.tail = FALSE)
