R语言学习系列2频率表和列联表

22.频率表和列联表分类变量包括无序分类变量、有序分类变量。对于分类变量的描述统计，主要是对分类变量各水平值分别进行频数和比例计算，再进步计算所需的一些相对频数指标。一、一维列联表频数统计分类变量的分析，首先要了解：各类别的样本数（频数），以及占总样本量的百分比；对有序分类变量，还需要了解：累积频数、累积百分比。使用table（）函数生成简单的频数统计表，基本格式为：table(var1,var2, •-,varN)将N个分类变量（因子）创建一个N维列联表。library(vcd)#使用数据集Arthritishead(Arthritis,3)574677SexAgeImprovedMale27574677SexAgeImprovedMale27SomeMale29NoneMale30NoneIDTreatmentTreatedTreatedTreated#统计频数#统计频数mytableNone SomeMarked42 14 28prop.table(mytable) #将频数转化为比例值None SomeMarked0.5000000 0.16666670.3333333prop.table(mytable)*100 #将频数转化为百分比NoneSomeMarked50.00000 16.6666733.33333集中/离散趋势观察原始频数，或者使用众数。对于分类变量，集中 /离散趋势是一体的。相对频数指标（1）比（Riatio）两个有关指标之比A/B,用来反映相对的大小关系，例如，月销售额/销售人数；（2）构成比用于描述事物内部各构成部分所占的比重，例如，百分比、累积百分比；（3）率（Rate）率是具有时间概念或速度、强度意义的指标，表示某个时期内某事件发生的频率或强度，例如速率、频率、费率、发病率等。二、二维列联表rxc二维列联表：(1)共n个样本；按两种属性A、B,属性A有r个水平值：Ai,…,Ar；属性B有c个水平值：Bi,…,Be.属性A=Ai，属性B=Bj的样本数为n,ni.=“属性A=Ai”的合计数，n.j=“属性B=Bj”的合计数。1二塞的FKf列联農■»*■含 it-"u…flh・・・"i-■；<■略1g*+*■ qA 计■1-n注：多分类变量对应高维列联表生成二维列联表使用函数table()，基本格式为：table(A,B)其中，A为行变量，B为列变量。或者使用函数xtabs()，基本格式为：xtabs(~A+B,data二mydata)其中，mydata为矩阵或数据框；要进行交叉分类的变量放在~右侧,以+作为分隔符，若某个变量写在~左侧，则为频数向量，常用于数据已经被表格化时。注：table函数默认忽略NA值，若要将NA值作为一个类别，需要使用参数useNA二"ifany"data=Arthritis)#生成data=Arthritis)#生成治疗方法与疗效的二维列联表（频数）mytableImprovedTreatmentNoneSomeMarkedTreatmentNoneSomeMarkedPlacebo29Treated1321prop.table(mytable)#Treated1321prop.table(mytable)#各单元格占的比例（所有行列总和=1）ImprovedTreatmentNoneSomeMarkedPlacebo0.345238100.08333333 0.08333333Treated0.154761900.08333333 0.25000000margin.table(mytable,1)#TreatmentNoneSomeMarkedPlacebo0.345238100.08333333 0.08333333Treated0.154761900.08333333 0.25000000margin.table(mytable,1)#对第1个变量（行变量）生成边际频数TreatmentPlaceboTreated414341prop.table(mytable,1)#prop.table(mytable,1)#对第1个变量（行变量）生成边际比例（每行和ImprovedTreatmentNoneSomeMarkedTreatmentNoneSomeMarkedPlacebo0.6744186 0.1627907 0.1627907Treated0.31707320.1707317 0.5121951#可见与安慰剂组相比，治疗组的治愈率从 16%提高到51%margin.table(mytable,2) #对第2个变量(列变量)生成边际频数ImprovedNoneSomeMarked421428421428prop.table（mytable,2）#对第2个变量（列变量）生成边际比例（每列和=1）ImprovedTreatmentNoneSomeMarkedPlacebo0.69047620.50000000.2500000Treated0.30952380.5000000 0.7500000TreatmentNoneSomeMarkedPlacebo0.69047620.50000000.2500000Treated0.30952380.5000000 0.7500000ImprovedTreatmentNoneSomeMarkedSumPlacebo297743Treated1372141Sum42142884#将边际求和结果添加到表格（频数）addmargins(mytable)#将边际求和结果添加到表格（比例）addmargins(prop.table(mytable))ImprovedTreatmentNoneSomeMarkedSumTreatmentNoneSomeMarkedSumPlacebo0.345238100.083333330.083333330.51190476Treated0.154761900.083333330.250000000.48809524Sum0.50000000 0.16666667 0.33333333 1.00000000addmargins(prop.table(mytable,1),2)#只添加列边际ImprovedTreatment None SomeMarked SumPlacebo0.67441860.16279070.16279071.0000000Treated0.31707320.17073170.51219511.0000000addmargins(prop.table(mytable.2),1)#只添加行边际ImprovedTreatmentNoneSome MarkedPlacebo0.69047620.50000000.2500000Treated0.30952380.50000000.7500000Sum1.00000001.00000001.00000002.gmodels包中的函数CrossTable()该函数仿照SAS中PROCFREQ或SPSS中CROSSTABS的形式生成二维列联表。它有很多选项，可以做许多事情：计算(行、列、单元格)的百分比；指定小数位数；进行卡方、Fisher和McNemar独立性检验；计算期望和(皮尔逊、标准化、调整的标准化)残差；将缺失值作为一种有效值；进行行和列标题的标注；生成SAS或SPSS风格的输出。library(gmodels)CrossTable(Arthritis$Treatment,Arthritis$lmproved)cellent?I N |Chi-squaretociributloftIM/RowTotalI\/ColTotalIN/rableTotalITotal0barrvat1005mTal>le:肌IArthritisllttprovedArthritl&ST'reatrftehlI Hone| Sottie|Mdrked| TotalITOC\o"1-5"\h\zPliccbn | 29 | 7 I 7 I 43 I| 2.61G | 0.004 | 3畀刃 |I 0.67^ I D.163 I 0.163 I 0r51? |0.690| 0.500I 0.250I0.345I O.OB3I 0.081Irr«LtedI 13 I 7 I 21 I 41 II 2.744 | 0.004 | J+935 II □.317 I 0,171 I 0.512 I0.480 II 0.310 I 0.500 I 0.750 I□.155I Q.O盟I IColurwiTotalI 42I 14 I 28I £4II 0.500| 0.16； | 0.133|三、多维列联表table()和xtabs()都可以基于三个或更多的分类变量(因子)生成多维列联表。margin.table()、prop.table()和addmargins()函数也可类似推广到多维。另外，ftable()函数可以用一种紧凑的“平铺式”输出多维列联mytable<-xtabs(~Treatment+Sex+lmproved,data=Arthritis)mytable,,Improved=NoneSexTreatmentFemaleMalePlacebo 19 10

Treated 67,,Improved=SomeSexTreatmentFemaleMalePlacebo 70Treated 52,,Improved=MarkedSexTreatmentFemaleMalePlacebo 61Treated 165ftable(mytable)ImprovedNoneSomeMarkedTreatment SexPlacebo Female19 76Male1001Treated Female6 516Male7 25margin.table(mytable,1)#边际频数TreatmentPlaceboTreated43 41margin.table(mytable,2)SexFemale Male59 25margin.table(mytable,3)ImprovedNone SomeMarked42 14 28margin.table(mytable,c(1,3))#相当于1:治疗方法x3:疗效的二维列联表ImprovedTreatmentNoneSomeMarkedPlacebo2977Treated13721ftable(prop.table(mytable,c(1,2)))#紧凑形式显示1:治疗方法x2:性别的各类疗效情况(比例)ImprovedNoneSomeMarkedTreatmentSexPlaceboFemale0.593750000.218750000.18750000Male0.909090910.000000000.09090909TreatedFemale0.222222220.185185190.59259259Male0.500000000.142857140.35714286

ftable(addmargins(prop.table(mytable,#继续给第3变量添加边际c(1,Some2)),3))ImprovedNoneMarkedSumTreatmentSexPlaceboFemale0.593750000.218750000.187500001.00000000Male0.909090910.000000000.090909091.00000000TreatedFemale0.222222220.185185190.592592591.00000000Male0.500000000.142857140.357142861.00000000ftable(addmargins(prop.table(mytable.c(1,2)),3))*100#继续改成百分比ImprovedNoneSomeMarkedSumTreatmentSexPlaceboFemale59.37500021.87500018.750000100.000000Male90.9090910.0000009.090909100.000000TreatedFemale22.22222218.51851959.259259100.000000Male50.00000014.28571435.714286100.000000多维列联表的图形展示一一马赛克图单个分类变量，可以使用柱状图或者饼图；若存在两个分类变量，可以使用三维柱状图；两个以上的分类变量可以使用马赛克图。马赛克图中，嵌套矩形面积正比于单元格频率，其中该频率即多维列联表中的频率。颜色和/或阴影可表示拟合模型的残差值。使用mosaic()函数，基本格式为：mosaic(table)其中，table为数组形式的列联表；或者用mosaic(formula,data=)其中，formula为形如~var1+var2+…的公式，data为数据框或表格；可选项shade二TRUE将根据拟合模型的皮尔逊残差值对图形上色；legend二TRUE将展示残差的图例library(vcd)ftable(Titanic)SurvivedNoYesClassSexAge1stMaleChild05Adult11857FemaleChild01Adult41402ndMaleChild011Adult15414FemaleChild013Adult13803rdMaleChild3513

Adult 38775TOC\o"1