本篇学习数据分析， Excel 表格制作
Task： 创建一个 行 百分比 表格
row percentage table
先看一下 PlantGrowth 数据集
library(dplyr)
data("PlantGrowth")
view(PlantGrowth)
给数据集新加一列 weight_cat ，并用 case_when 自定义名称
PG = PlantGrowth |>
mutate(weight_cat = case_when(
weight <= 4.5 ~ "light",
weight > 4.5 & weight <= 5.5 ~ "medium",
weight > 5.5 ~ "heavy"
)
)
view(PG)
统计 group 和 weight_cat 的次数
table(PG$group, PG$weight_cat)
table_count1 = table(PG$group,PG$weight_cat)
view(table_count1)
addmargins
addmargins(table_count1, 1) # 在最后一行加
addmargins(table_count1, 2) # 在最后一列加
addmargins(table_count1) # 在最后一行和最后一列都加
创建 行 百分比 表格
即每行加起来为100%，但每列加起来不一定是100%
# row percentage table
table_count2 = addmargins(table_count1, 1)
table_count2
proportion 比例
prop.table(x,1) 按行计算百分比
prop.table(x,2) 按列计算百分比
*100 转换成百分制
最后一列加sum
table_count3 = addmargins(table_count2,2)
table_count3
table_per2 = addmargins(table_per1,2)
table_per2
可以看到每行 sum 都是100
paste0
如果直接把 count 和 percentage 拼接起来，结果有点奇怪，可以根据需求自定义修改
paste0(table_count3, "(", round(table_per2), "%)")
先用round去掉了小数点，再把百分比数值前边加上（，再把百分比数值右边加上%和）
matrix
matrix(paste0(table_count3, "(", round(table_per2), "%)"),dim(table_count3))
将输出变成矩阵的形式，可以指定输出的 shape，用dimention来指定shape更方便
数据有了，我们还缺少行名和列名，下面介绍一个快捷方式
dimnames
先看例子
dimnames(table_count3)
先输出行名，在输出列名
加到我们需要的表格上
matrix(paste0(table_count3, "(", round(table_per2), "%)"),dim(table_count3),dimnames = dimnames(table_count3))
以上，我们就得到了需要的表格，现在将最后一条代码规整化并命名新的表名
tab_combine = matrix(paste0(table_count3, "(", round(table_per2), "%)"),
dim(table_count3),
dimnames = dimnames(table_count3)
)
完整代码
library(dplyr)
data("PlantGrowth")
view(PlantGrowth)
PG = PlantGrowth |>
mutate(weight_cat = case_when(
weight <= 4.5 ~ "light",
weight > 4.5 & weight <= 5.5 ~ "medium",
weight > 5.5 ~ "heavy"
)
)
view(PG)
table_count1 = table(PG$group,PG$weight_cat)
table_count2 = addmargins(table_count1, 1)
table_per1 = prop.table(table_count2,1)*100
table_count3 = addmargins(table_count2,2)
table_per2 = addmargins(table_per1,2)
tab_combine = matrix(paste0(table_count3, "(", round(table_per2), "%)"),
dim(table_count3),
dimnames = dimnames(table_count3)
)
拓展任务：column percentage table & overall percentage table
column percentage table
tab_1 = table(PG$group,PG$weight_cat)
tab_2 = addmargins(tab_1,2)
tab_3 = addmargins(tab_2,1)
tab_per1 = prop.table(tab_2,2)*100
tab_per2 = addmargins(tab_per1,1)
tab_colmn_combine = matrix(paste0(tab_3, "(", round(tab_per2), "%)"),
dim(tab_3),
dimnames = dimnames(tab_3)
)
tab_colmn_combine
overall percentage table
tab_1 = table(PG$group,PG$weight_cat)
tab_2 = addmargins(tab_1)
tab_per1 = prop.table(tab_1)*100
tab_per2 = addmargins(tab_per1)
tab_colmn_combine = matrix(paste0(tab_2, "(", round(tab_per2), "%)"),
dim(tab_2),
dimnames = dimnames(tab_2)
)
tab_colmn_combine
