R语言中因子含3个水平但箱线图仅显示1个箱体的问题

问题情况 #

我正在按照R教程做嵌套ANOVA检验，先想用箱线图可视化数据，但x轴的location因子明明有3个水平，箱线图却只显示1个箱体。所有数据都已经用as.factor()转成因子类型，数据结构如下：

> str(ANOVADATArobin)
tibble [105 × 4] (S3: tbl_df/tbl/data.frame)
 $ location  : Factor w/ 3 levels "1","2","3": 1 1 1 1 1 1 1 1 1 1 ...
 $ site      : Factor w/ 12 levels "1","2","3","4",..: 1 1 1 2 2 2 3 3 3 4 ...
 $ repeat    : Factor w/ 3 levels "1","2","3": 1 2 3 1 2 3 1 2 3 1 ...
 $ flighttime: Factor w/ 73 levels "0","3","5","6",..: 73 73 31 57 73 56 73 65 73 73 ...

执行的绘图代码：

> boxplot(flighttime-location, xlab="location", ylab="flighttime")

现象：

• 初始代码生成的箱线图仅显示1个箱体
• 尝试添加x=factor(location)后仍未解决，仅出现第二条线
• 目标是生成按location分组、每个水平对应一个箱体的标准箱线图

问题根源 #

1. flighttime类型错误：把响应变量flighttime转成了因子类型，但箱线图需要数值型变量来计算中位数、四分位数等统计量，因子作为分类变量无法完成这个计算，只会被当作单一类别处理
2. 公式语法错误：基础R的boxplot分组语法应该用~（波浪线）表示“y随x变化”，而不是-（减号），flighttime-location实际是做数值减法，根本不是分组绘图

解决方案 #

步骤1：转换flighttime为数值型 #

注意要先转成字符再转数值，直接转因子会得到水平的编码值（比如数据里的73、31是因子的水平索引，不是实际飞行时间）：

ANOVADATArobin$flighttime <- as.numeric(as.character(ANOVADATArobin$flighttime))

步骤2：用正确语法绘制箱线图 #

使用~表示分组关系，并指定数据框：

boxplot(flighttime ~ location, data = ANOVADATArobin, xlab = "location", ylab = "flighttime")

可选：用ggplot2绘制更美观的图 #

如果需要更灵活的可视化效果，推荐用ggplot2包：

library(ggplot2)
ggplot(ANOVADATArobin, aes(x = location, y = flighttime)) +
  geom_boxplot() +
  labs(x = "location", y = "flighttime")

内容的提问来源于stack exchange，提问作者Kat Sandamas