县域土壤有效磷测试数据数理统计分析应用.docVIP

县域土壤有效磷测试数据数理统计分析应用摘要：将县域测土配方施肥土壤有效磷测试数据视为大样本正态总体抽样分布，应用概率密度曲线函数变量X落入（μ-3σ, μ+3σ）概率占99.73%，据此标准对数据进行多轮次筛选，提高样本的代表性。结合专业背景知识，数据统计分析，作频数分布图表分析，研究分析县域土壤有效磷状况，提出施肥指导建议。 关键词：有效磷；正态分布；筛选；频数分布图表 中图分类号：S153.6 文献标识码：A 文章编号： 1674-0432（2011）-10-0083-2 土壤有效磷在土壤养分中属比较稳定的，适宜作为有效养分恒量监控指标，根据土壤有效磷测试结果和养分丰缺指标进行分级，指导实际生产的配方施肥有重要意文。为此应用原平市2008-2011年测土配方施肥项目土壤有效磷测试数据6347个做样本源，根据生物统计样本与总体的关系，科学研究所观测的数据只是随机变量X全部可能的一部分，统计学要做的事就是以“部分”所具有的统计特征来推断“全部”所具有的统计特征。根据统计学原理，对于实际问题中的非正态总体抽样，只要样本容量N充分大，根据大数定律，能将统计数的抽样分布视作正态总体抽样分布，据此以6347个土壤有效磷样本进行描述性统计分析，结果见土壤有效磷统计分析结果表（表1）第2列。 为了克服采样、化验等环节不可避免带来的异常数据，根据正态总体抽样分布概率密度曲线函数中变量X落入（μ-3σ,μ+3σ）概率占99.73%(μ为总体平均数)，据此对数据进行筛选，其原理参见标准正态分布概率密度曲线函数（图1）。 为方便计算，足够大样本取σ=S（σ为总体标准差，S为样本标准差），将6347个数据进行描述性统计分析，计算样本平均值x=13.63406mg/?，计算样本标准差S=10.29324mg/?，则可计算出x±3S的区间值，对其中大于x-3S和小于x+3S的数值进行剔除，再次对构成的新样本进行描述性统计，重新得到新的x和s，再依照x±3S标准对数据再次进行剔除，如此进行，直到所有数据值都在x±3S范围内，即从统计学角度保证所用数据值中有99.73%的概率保证分布在x±3S的区间内，提高了样本的代表性。在数据的剔除与结果分析中专业背景知识是不可或缺的，可更好地助于数据筛选。本次数据筛选经12轮次筛选，从6347个数据中共剔除439个，保留了5908个数据作为这次统计的样本基数。5908个数据描述统计分析有下列特征，结果见土壤有效磷统计分析结果表（表1）第3列。 土壤有效磷平均值11.5mg/?,中位值11.0mg/?，标准差6.58mg/?，最小值1.0mg/?，最大值31.2 mg/?，极差30.2 mg/?，有效统计数5908个。峰度0.298526，偏度0.93598（标准正态分布偏度为0）。 与未经筛选前的数据相比，结果见土壤有效磷统计分析结果表（表1）第4列。有效统计数减少439个，有效磷平均值降低了2.1mg/?，极差缩小36.0mg/?，最大值降低了36.0mg/?，变异系数（%）降低了18.3，置信度(95.0%)提高了0.08544，偏度降低了0.93069，曲线分布更加接近正态分布曲线。对数据筛选利用有了更科学的方法作保证，与此结果相对应的是土壤有效磷统计数更能反映当地土壤养分的真实水平，可以更好地用于指导实践。 为了了解土壤有效磷的详细分布，虽然对筛选出的数据用EXCEL表可以进行多角度数据分析，但我们采用了做数据的频数分布图表的方法，本法能更直观地看问题，解决问题。为了方便实际应用，频数分布按分10组进行（实际由于频数分布计算的方法最后形成分了11组）。选择组距30.2/10=3.02，第一组组中点值1.0，上限值2.51，依此类推。统计次数分布如下（表2）： 根据以上土壤有效磷次数分布表做频数分布直方图（图2），进行图表结合说明问题。 图2 从频数分布直方图和表可明显看出全市土壤有效磷次数分布的整体趋势情况，以2.51-5.53-8.55-11.6-14.6四段次数都在800次以上，经统计共占到总次数的约60%。说明土壤有效磷含量在2.51 mg /?-14.6 mg /?之间的相对较多。土壤有效磷平均值11.5 mg /?以上的约占到26%（9.24+6.16+4.35+3.25+2.47+1.22≈26），约1/4的土地。而绝大部分土壤有效磷在平均值11.5mg/?以下，约占74%，换言之约3/4的土地有效磷在平均值以下，全市整体有效磷水平偏低，全市整体施磷水平偏低，需补施磷。针对原平的主栽作物玉米，如果以土壤有效磷≤10mg/?为磷偏低，10-15mg/?为磷适宜，≥15mg/?为磷偏高作为丰缺界线指标，≤10mg/?磷偏低土