#分段回归的代码 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,   ...
#有时候排序图的中物种或样方过多,造成标识重叠的情况,可以用vegan包中的orditorp将标识虚化 ## A cluttered ordination plot : library(vegan) data(BCI) mod - cca(BCI) plot(mod, dis=sp, type=t) # Now with orditorp and abbreviated species names cnam - make.cepnames(names( ...