1。3。2 麦长管蚜选择性实验的数据分析

由于对同一重复的4次不同时间的调查数据存在相关性，故将调查次数作为随机效应变量，用广义线性混合效应模型（poisson分布）分析生物质炭处理（固定效应变量）与各处理上有翅蚜的数目间关系。用R软件的lme4功能包拟合模型，用似然比测验比较模型的简约性。用lsmeans功能包比较逐对处理之间的差异。数据分析用R软件[29]。

2 结果与分析

2。1 生物质炭处理对麦长管蚜体型大小的影响

各处理情况下，麦长管蚜具体体型大小情况见表2。通过单因素方差分析，我们发现麦长管蚜体型大小在施用生物质炭的处理与未施用生物质炭处理间存在显著差异（F = 7。8869，p < 0。05），施用生物质炭的处理，平均大小显著长于未施用生物质炭的处理0。0791mm（8。33%，如图1）。