[生物学]Algae海藻例子.ppt

首都经济贸易大学 32009122125 * library(leaps) lm2leaps-regsubsets(a1~factor(season)+factor(size)+factor(speed)+mxPH+mnO2+Cl+NO3+NH4+oPO4+PO4+Chla,data=clean.algae, nbest=1, nvmax=15,really.big=TRUE) ##方程中不含常数项共有15 个自变量，所以nvmax=15，而nbest 为相同个数的自变量的最佳组合排列的方程个数。系统默认为1，即取最好的，因为我要找出所有变量组合中最好的模型，所以为了计算简便，就取系统默认的，因为运算量较大，really.big=TRUE，系统默认really.big=FALSE）。 lm2leaps.sum-summary(lm2leaps) names(lm2leaps.sum) [1] which rsq rss adjr2 cp bic outmat obj # 建立仅有截距项的模型，即规模为0 的模型 lm2leapsdum-lm(a1~1, data=clean.algae) lm2leapsdum.sum-summary(lm2leapsdum) names(lm2leapsdum.sum) [1] call