我有一个包含因子和数值列的R数据框架,我想通过各种因子分组来计算数值列的某些汇总。特别是,我希望能够一次通过多个因素组合进行总结,并在单个对象中返回结果。如果我想要得到每种毛线和张力组合的平均折断次数,我知道我可以使用aggregate(breaks ~ wool + tension, data = warpbreaks, mean),它会给我类似的结果L 36.388895 A NA 31.037047 A L 4
这是我的密码;y <- runif(100)r <- rep(1,times=100) # random effectr[51:100] <- 2 # replace 1 into 2, making two groups
gam_fit <- gam(z~ s(x) + s(y) + s(r,bs="re"), d