R in action读书笔记（6）-第七章：基本统计分析（中）

2019-04-10 17:19:45 浏览数 (14)

7.2 频数表和列联表

> library(vcd)

> head(Arthritis)

ID Treatment Sex Age Improved

1 57 Treated Male 27 Some

2 46 Treated Male 29 None

3 77 Treated Male 30 None

4 17 Treated Male 32 Marked

5 36 Treated Male 46 Marked

6 23 Treated Male 58 Marked

7.2.1 生成频数表

函 数 描 述

table(var1, var2, …, varN) 使用 N 个类别型变量（因子）创建一个 N 维列联表

xtabs(formula, data) 根据一个公式和一个矩阵或数据框创建一个 N 维列联表

prop.table(table, margins) 依margins定义的边际列表将表中条目表示为分数形式

margin.table(table, margins) 依margins定义的边际列表计算表中条目的和

addmargins(table, margins) 将概述边margins（默认是求和结果）放入表中

ftable(table) 创建一个紧凑的“平铺”式列联表

  1. 一维列联表
代码语言：javascript复制
> mytable<-with(Arthritis,table(Improved))
代码语言：javascript复制
> mytable
代码语言：javascript复制
Improved
代码语言：javascript复制
  None   Some Marked
代码语言：javascript复制
    42     14     28

可以用prop.table()将这些频数转化为比例值：

代码语言：javascript复制
> prop.table(mytable)
代码语言：javascript复制
Improved
代码语言：javascript复制
     None      Some    Marked
代码语言：javascript复制
0.5000000 0.1666667 0.3333333

或使用prop.table()*100转化为百分比：

2. 二维列联表

对于二维列联表，table()函数的使用格式为：mytale<-table(A,B)

其中的A是行变量，B是列变量。除此之外，xtabs()函数还可使用公式风格的输入创建列联表，

格式为：mytable<-xtabs(~A B,data=mydata)

其中的mydata是一个矩阵或数据框。总的来说，要进行交叉分类的变量应出现在公式的右侧（即~符号的右方），以 作为分隔符。若某个变量写在公式的左侧，则其为一个频数向量（在数据已经被表格化时很有用）。

对于Arthritis数据，有：

代码语言：javascript复制
> mytable<-xtabs(~Treatment Improved,data=Arthritis)
代码语言：javascript复制
> mytable
代码语言：javascript复制
         Improved
代码语言：javascript复制
Treatment None Some Marked
代码语言：javascript复制
  Placebo   29    7      7
代码语言：javascript复制
  Treated   13    7     21

可以使用margin.table()和prop.table()函数分别生成边际频数和比例。行和与行比

例可以这样计算：

代码语言：javascript复制
> margin.table(mytable,1)
代码语言：javascript复制
Treatment
代码语言：javascript复制
Placebo Treated
代码语言：javascript复制
     43      41
代码语言：javascript复制
> prop.table(mytable,1)
代码语言：javascript复制
             Improved
代码语言：javascript复制
Treatment   None      Some    Marked
代码语言：javascript复制
  Placebo 0.6744186 0.1627907 0.1627907
代码语言：javascript复制
  Treated 0.3170732 0.1707317 0.5121951

列和与列比例可以这样计算：

代码语言：javascript复制
> margin.table(mytable,2)
代码语言：javascript复制
Improved
代码语言：javascript复制
  None   Some Marked
代码语言：javascript复制
    42     14     28
代码语言：javascript复制
> prop.table(mytable,2)
代码语言：javascript复制
         Improved
代码语言：javascript复制
Treatment      None      Some    Marked
代码语言：javascript复制
  Placebo 0.6904762 0.5000000 0.2500000
代码语言：javascript复制
  Treated 0.3095238 0.5000000 0.7500000

各单元格所占比例可用如下语句获取：

代码语言：javascript复制
> prop.table(mytable)
代码语言：javascript复制
         Improved
代码语言：javascript复制
Treatment       None       Some     Marked
代码语言：javascript复制
  Placebo 0.34523810 0.08333333 0.08333333
代码语言：javascript复制
  Treated 0.15476190 0.08333333 0.25000000

可以使用addmargins()函数为这些表格添加边际和

代码语言：javascript复制
> addmargins(mytable)
代码语言：javascript复制
         Improved
代码语言：javascript复制
Treatment None Some Marked Sum
代码语言：javascript复制
  Placebo   29    7      7  43
代码语言：javascript复制
  Treated   13    7     21  41
代码语言：javascript复制
  Sum       42   14     28  84
代码语言：javascript复制
> addmargins(prop.table(mytable))
代码语言：javascript复制
         Improved
代码语言：javascript复制
Treatment       None       Some     Marked        Sum
代码语言：javascript复制
  Placebo 0.34523810 0.08333333 0.08333333 0.51190476
代码语言：javascript复制
  Treated 0.15476190 0.08333333 0.25000000 0.48809524
代码语言：javascript复制
  Sum     0.50000000 0.16666667 0.33333333 1.00000000

在使用addmargins()时，默认行为是为表中所有的变量创建边际和

代码语言：javascript复制
> addmargins(prop.table(mytable,1),2)#仅添加了各行的和
代码语言：javascript复制
         Improved
代码语言：javascript复制
Treatment      None      Some    Marked       Sum
代码语言：javascript复制
  Placebo 0.6744186 0.1627907 0.1627907 1.0000000
代码语言：javascript复制
  Treated 0.3170732 0.1707317 0.5121951 1.0000000

注意 table()函数默认忽略缺失值（NA）。要在频数统计中将NA视为一个有效的类别，请设定参数useNA="ifany"。.

使用gmodels包中的CrossTable()函数是创建二维列联表的第三种方法。CrossTable()

函数仿照SAS中PROC FREQ或SPSS中CROSSTABS的形式生成二维列联表

代码语言：javascript复制
> CrossTable(Arthritis$Treatment,Arthritis$Improved)

CrossTable()函数有很多选项，可以做许多事情：计算（行、列、单元格）的百分比；指

定小数位数；进行卡方、Fisher和McNemar独立性检验；计算期望和（皮尔逊、标准化、调整的

标准化）残差；将缺失值作为一种有效值；进行行和列标题的标注；生成SAS或SPSS风格的输出。

3.多维列联表

table()和xtabs()都可以基于三个或更多的类别型变量生成多维列联margin.table()、prop.table()和addmargins()函数可以自然地推广到高于二维的情况。另外，ftable()函数可以以一种紧凑而吸引人的方式输出多维列联表

代码语言：javascript复制
> mytable<-xtabs(~Treatment Sex Improved,data=Arthritis)
代码语言：javascript复制
, , Improved = None
代码语言：javascript复制
         Sex
代码语言：javascript复制
Treatment Female Male
代码语言：javascript复制
  Placebo     19   10
代码语言：javascript复制
  Treated      6    7
代码语言：javascript复制
, , Improved = Some
代码语言：javascript复制
         Sex
代码语言：javascript复制
Treatment Female Male
代码语言：javascript复制
  Placebo      7    0
代码语言：javascript复制
  Treated      5    2
代码语言：javascript复制
, , Improved = Marked
代码语言：javascript复制
         Sex
代码语言：javascript复制
Treatment Female Male
代码语言：javascript复制
  Placebo      6    1
代码语言：javascript复制
  Treated     16    5
代码语言：javascript复制
代码语言：javascript复制
> ftable(mytable)
代码语言：javascript复制
                 Improved None Some Marked
代码语言：javascript复制
Treatment Sex
代码语言：javascript复制
Placebo   Female            19    7      6
代码语言：javascript复制
          Male              10    0      1
代码语言：javascript复制
Treated   Female             6    5     16
代码语言：javascript复制
          Male               7    2      5

> margin.table(mytable,c(1,3))#治疗情况（Treatment） × 改善情况（Improved）的边际频数

代码语言：javascript复制
         Improved
代码语言：javascript复制
Treatment None Some Marked
代码语言：javascript复制
  Placebo   29    7      7
代码语言：javascript复制
  Treated   13    7     21

7.2.2独立性检验

1. 卡方独立性检验

可以使用chisq.test()函数对二维表的行变量和列变量进行卡方独立性检验

代码语言：javascript复制
> library(vcd)
代码语言：javascript复制
> mytable<-xtabs(~Treatment Improved,data=Arthritis)
代码语言：javascript复制
> chisq.test(mytable)
代码语言：javascript复制
       Pearson's Chi-squared test
代码语言：javascript复制
data:  mytable
代码语言：javascript复制
X-squared = 13.055, df = 2, p-value = 0.001463#治疗情况和改善情况不独立

2. Fisher精确检验

可以使用fisher.test()函数进行Fisher精确检验。Fisher精确检验的原假设是：边界固定

的列联表中行和列是相互独立的。其调用格式为fisher.test(mytable)，其中的mytable是

一个二维列联表

代码语言：javascript复制
> fisher.test(mytable)
代码语言：javascript复制
       Fisher's Exact Test for Count Data
代码语言：javascript复制
data:  mytable
代码语言：javascript复制
p-value = 0.001393
代码语言：javascript复制
alternative hypothesis: two.sided

这里的fisher.test()函数可以在任意行列数大于等于2的二维列联表上使用，但不能用于2×2的列联表。

3.Cochran-Mantel—Haenszel检验

mantelhaen.test()函数可用来进行Cochran—Mantel—Haenszel卡方检验，其原假设是，两

个名义变量在第三个变量的每一层中都是条件独立的。

代码语言：javascript复制
 > mantelhaen.test(mytable)
代码语言：javascript复制
       Cochran-Mantel-Haenszel test
代码语言：javascript复制
data:  mytable
代码语言：javascript复制
Cochran-Mantel-Haenszel M^2 = 14.6323, df = 2,
代码语言：javascript复制
p-value = 0.0006647

7.2.3 相关性的度量

如果可以拒绝原假设，那么你的兴趣就会自然而然地转向用以衡量相关性强弱的相关性度量。vcd包中的assocstats()函数可以用来计算二维列联表的phi系数、列联系数和Cramer’sV系数

代码语言：javascript复制
> mytable<-xtabs(~Treatment Improved,data=Arthritis)
代码语言：javascript复制
> assocstats(mytable)
代码语言：javascript复制
                    X^2 df  P(> X^2)
代码语言：javascript复制
Likelihood Ratio 13.530  2 0.0011536
代码语言：javascript复制
Pearson          13.055  2 0.0014626
代码语言：javascript复制
Phi-Coefficient   : 0.394
代码语言：javascript复制
Contingency Coeff.: 0.367
代码语言：javascript复制
Cramer's V        : 0.394

总体来说，较大的值意味着较强的相关性。vcd包也提供了一个kappa()函数，可以计算混

淆矩阵的Cohen’s kappa值以及加权的kappa值。（举例来说，混淆矩阵可以表示两位评判者对于一系列对象进行分类所得结果的一致程度。）

7.2.5将表转换为扁平格式

通过table2flat将表转换为扁平格式

代码语言：javascript复制
> table2flat<-function(mytable){
代码语言：javascript复制
  df<-as.data.frame(mytable)
代码语言：javascript复制
  rows<-dim(df)[1]
代码语言：javascript复制
  cols<-dim(df)[2]
代码语言：javascript复制
  x<-NULL
代码语言：javascript复制
  for(i in 1:rows){
代码语言：javascript复制
  for(j in 1:df$Freq[i]){
代码语言：javascript复制
  row<-df[i,c(1:(cols-1))]
代码语言：javascript复制
  x<-rbind(x,row)
代码语言：javascript复制
  }
代码语言：javascript复制
  }
代码语言：javascript复制
  row.names(x)<-c(1:dim(x)[1])
代码语言：javascript复制
  return(x)
代码语言：javascript复制
  }

使用table2flat()函数转换已发表的数据

代码语言：javascript复制
> treatment<-rep(c("Placebo","Treated"),times=3)
代码语言：javascript复制
> improved<-rep(c("None","Some","marked"),each=2)
代码语言：javascript复制
> Freq<-c(29,13,7,17,7,21)
代码语言：javascript复制
> mytable<-as.data.frame(cbind(treatment,improved,Freq))
代码语言：javascript复制
> mydata<-table2flat(mytable)
代码语言：javascript复制
> head(mydata)
代码语言：javascript复制
  treatment inmproved
代码语言：javascript复制
1   Placebo      None
代码语言：javascript复制
2   Placebo      None
代码语言：javascript复制
3   Placebo      None
代码语言：javascript复制
4   Placebo      None
代码语言：javascript复制
5   Treated      None
代码语言：javascript复制
6   Placebo      Some
编程算法 音视频点播加速 r语言

0 人点赞