R中提供了用于创建频数表和列联表的若干种方法。
可以使用table()函数生成简单的频数统计表。示例如下:
> mytable <- with(Arthritis, table(Improved))
> mytable
Improved
None Some Marked
42 14 28
可以用prop.table()将这些频数转化为比例值:
> prop.table(mytable)
Improved
None Some Marked
0.500 0.167 0.333
或使用prop.table()*100转化为百分比:
> prop.table(mytable)*100
Improved
None Some Marked
50.0 16.7 33.3
这里可以看到,有50%的研究参与者获得了一定程度或者显著的改善
对于二维列联表,table()函数的使用格式为:
mytable <- table(A, B)
其中的A是行变量,B是列变量。除此之外,xtabs()函数还可使用公式风格的输入创建列联表,
格式为:
mytable <- xtabs(~ A + B, data=mydata)
其中的mydata是一个矩阵或数据框。总的来说,要进行交叉分类的变量应出现在公式的右侧(即
~符号的右方),以+作为分隔符。若某个变量写在公式的左侧,则其为一个频数向量(在数据已