#分段回归的代码 x - c(1:10, 13:22) y - numeric(20) ## Create first segment y - 20:11 + rnorm(10, 0, 1.5) ## Create second segment y - seq(11, 15, len=10) + rnorm(10, 0, 1.5) ## Plot it library(segmented) lin.mod - lm(y~x) #psi starting values for the breakpoints to b ...
下面这些代码将展示由Francois Gillet编写的两个方差分析多重比较的函数 boxplert()和boxplerk(),来源Numerical Ecology with R (second Edition) #普通方差分析多重比较 boxplert-function(X, Y,   ...