使用reshape2转换airquality数据时大量值为NA的问题排查与解决

你这次的问题主要是踩了两个小坑：

• 变量名拼接时的错位问题：你用airquality$Month和airquality$Day去拼接aql$variable，但melt后的aql行数是原数据集的4倍（对应Ozone、Solar.R、Wind、Temp这4个观测变量），而原airquality只有153行。R会自动循环重复原数据集的月份和日期值，导致每个变量对应的日期月份完全错位，后续dcast根本找不到对应的匹配值，自然就出现大量NA了。
• dcast语法错误：你写的Month~variable的内容是无效语法，正确的写法应该是直接指定Month~variable，用来定义行和列的对应关系。

直接用下面的代码就能解决问题：

# 加载包和数据集
library(reshape2)
data(airquality)

# 第一步： melt数据集，保留Month和Day作为id变量
aql <- melt(airquality, id.vars = c("Month", "Day"))

# 关键修正：用aql自身的Month和Day拼接变量名，保证对应关系正确
aql$variable <- paste0(aql$Month, aql$Day, aql$variable)
# （可选）如果想让列名可读性更强，可以加分隔符：
# aql$variable <- paste0(aql$Month, "_", aql$Day, "_", aql$variable)

# 去掉Day列
aql <- aql[, -2]

# 正确执行dcast，将Month设为行，拼接后的variable设为列
aqw <- dcast(aql, Month ~ variable)

运行后aqw里的数值就会正确填充，不会再出现大量NA了。

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