不知道这样合不合你的原意
> aaov <- function(y)
+ { airquality$x<-airquality[[y]]
+ aov(x~Month,data=airquality)
+ }
>
> aaov("Wind")
Call:
aov(formula = x ~ Month, data = airquality)
Terms:
Month Residuals
Sum of Squares 59.9702 1826.5836
Deg. of Freedom 1 151
Residual standard error: 3.478014
Estimated effects may be unbalanced


雷达卡


京公网安备 11010802022788号







