以a=0.05为显著性水平,求出各个理化因子的置信区间如下表:
表3-3 每个主要理化因子置信区间统计表 总磷 磷酸盐磷 总氮 硝态氮 亚硝态氮 铵态氮
因为1号池塘为发生轻微水华,而由上表得每个理化因子的置信区间上限值分别为14.93045404,0.201131467,11.20798359,2.604117291,8.382504421,32.4057669,下限值分别为4.99921263,-0.018104801,3.208309738,0.202096042,-0.324324421,9.832539765。因此,我们预测发生重度水华时的主要理化因子的含量将大于置信区间的上限值,而不发生水华的主要理化因子范围小于置信区间的下限值。 4.4 问题四的解答 建立鱼类生长模型
针对问题四,由附件4里体长和体重的数据,通过matlab里的cftool工具对鳙鱼数据进行回归。具体处理结果如下: Linear model Poly3:
f(x) = p1*x^3 + p2*x^2 + p3*x + p4 Coefficients (with 95% confidence bounds):
p1 = 0.0144 (0.01062, 0.01817) p2 = -2.367 (-2.827, -1.906) p3 = 156 (140.5, 171.6) p4 = 891.1 (766.8, 1015)
Goodness of fit: SSE: 9.385e+07 R-square: 0.9154
Adjusted R-square: 0.915 RMSE: 428.6
得出方差SSE的数值是9.385e+07,这个值很小,回归可决系数R=0.9154,说明拟合效果很好,此次鳙鱼的生长和体重关系拟合可行。故
y=0.0144*x^3 -2.367*x^2 + 156*x + 891.1 (4)
10
置信区间1 4.99921263 -0.018104801 3.208309738 0.202096042 -0.324324421 9.832539765
置信区间2 14.93045404 0.201131467 11.20798359 2.604117291 8.382504421 32.4057669
用matlab画出函数图像如下:
图4-1 鳙鱼体长与体重函数关系拟合图
同理对鲢鱼的数据进行拟合,数据处理结果如下:
Linear model Poly3:
f(x) = p1*x^3 + p2*x^2 + p3*x + p4 Coefficients (with 95% confidence bounds):
p1 = 0.006887 (0.006181, 0.007592) p2 = -1.058 (-1.139, -0.9772) p3 = 80.2 (77.58, 82.82) p4 = 106.3 (85.02, 127.6)
Goodness of fit: SSE: 1.461e+06 R-square: 0.9954
Adjusted R-square: 0.9954 RMSE: 53.21
得出方差SSE的数值是1.461e+06,这个值很小,回归可决系数R=0.9954,说明拟合效果很好,此次鲢鱼的生长和体重关系拟合可行。故
y=0.006887*x^3 -1.058*x^2 + 80.2*x + 106.3 (5)