森林可燃物平均值调查的样本数确定方法技术

本发明专利技术公开了森林可燃物平均值调查的样本数确定方法，利用R软件制作了“正态分布单变量平均值调查所需样本数查询表”，用于森林可燃物外业调查现场的样本数计算，为不多不少地抽样调查并得到较为可靠的结果提供技术方法；查询表基于t值、预调查均值和方差，原理清晰；以森林可燃物样线数确定和可燃物平均直径调查为例说明查询表的使用，并以可燃物平均直径计算了可燃物表积比；采用R软件power函数对预调查和后续调查合并数据进行功效检验，检验平均值的可靠性；方法不仅适用于森林可燃物调查，也可用于其他平均值调查研究。

A method for determining the number of samples for the investigation of the average value of forest combustibles

The invention discloses a method to determine the average value of forest fuel sample survey, produced a univariate normal distribution survey average number of samples required query table by using R software for forest fuel industry survey site sample number for much less sample survey and get more reliable the results provide technical methods; query table based on t value, mean and variance, the pre survey principle is clear; the average diameter of investigation to forest fuel line and determine the number of combustible materials as an example to illustrate the use of the query table, and the average diameter of combustible fuel can be calculated surface area ratio; using R software of power function with the data pre investigation and follow-up investigation function test, reliability test mean value; method is not only applicable to the investigation of forest fuel, can also be used for other average investigation.

【技术实现步骤摘要】
森林可燃物平均值调查的样本数确定方法
本专利技术属于森林生态领域，涉及正态分布数据指定平均值精确度要求下的样本数计算，具体涉及森林可燃物平均值调查的样本数确定方法。
技术介绍
样本数大小和样本数的决定对森林生态调查有着重要影响。森林可燃物调查的样本数包括样线条数、某一径级可燃物直径测量个数、样方个数等等。森林可燃物调查通常要求样本平均值误差在总体真实平均值±20％范围内，增加样本数可以减少标准误，提高结果的准确性，但增加样本数意味着成本加大；如何以有限的成本，不多不少地抽样调查，最终得到较为可靠的结果是本专利技术的目的。植物生态调查一般预先在目标较多的地方少量调查5～6个样方，再在目标明显较少的地方调查几个样方，用预调查数据估算均值和方差，再以此计算样本数。正态分布单变量的样本数计算有多种方法，比如Elzinga(2001,MonitoringPlantandAnimalPopulations)用类型Ⅰ错误的标准正态系数、标准差、半幅置信区间来计算初始样本数，再通过复杂的查表将初始样本数转为最终样本数；Bonhan(2013,MeasurementsforTerrestrialVegetation)用t值、标准差、均值和样本数等式，反复替代不同t值和n至等式成立；SPSSSamplePower和R软件也有计算样本数的专门软件包。目前我国大多数森林可燃物调查仍然采用落后的样方收获法，沿样地对角线设置极小样方收获全部可燃物，由于这种方法费工费时，样方数通常偏少；改用样线法进行可燃物调查是大势所趋，但究竟需要调查多少条样线才能得到有相当把握的可燃物预测值？在我国气候条件下，常见树种各径级可燃物的平均直径、表积比是多少？这些问题在国内都无先例可借鉴，缺乏恰当技术方法和可燃物基本信息是我国森林防火研究的特点。本专利技术编制了样本数查询表，以可燃物样线数确定和可燃物平均直径调查为例说明查询表的使用方法，利用平均直径计算表积比，为我国林火预测和防控提供基础数据收集的基本方法。
技术实现思路
本专利技术的目的是针对森林可燃物调查中存在的技术问题，旨在提供一种森林可燃物平均值调查的样本数确定方法。为了实现上述目的，本专利技术的技术方案如下：一种森林可燃物平均值调查的样本数确定方法，其特征在于：包括样本数查询表制作原理、样本数查询表制作方法、可燃物调查所需样线数的计算、可燃物平均直径测量的样本数计算、可燃物表积比计算、单变量平均值的功效检验工序，主要操作步骤如下：(1)样本数查询表制作原理设定样本均值与总体真实均值μ之差为样本均值的k倍，样本数为n,得到去除方根号得到样本数计算公式：总体方差用预调查得到的样本方差来替代，森林可燃物调查中的k值一般可取0.2，因此通过预调查可以得到的值；但t值随预调查自由度(n-1)和所需置信水平(1-α)变化，因此，需要反复试验代入成对t、n值至公式成立；为了简化样本数计算，我们将不同置信水平的n和其对应的列成查询表，根据预调查的值和的关系，查表即可得到可燃物调查所需样本数；(2)样本数查询表制作方法利用R软件输入以下代码制作样本数查询表：n<-1:100；alpha<-0.20；a<-n/(qt(1-alpha/2,n-1))^2write.table(a,"clipboard",sep＝"",col.names＝NA)至电子表格粘贴得表1第一列和第二列，修改代码中的alpha值得到其他几列；表1为两尾检验所需样本数；如果进行的是单尾检验，则表中置信水平80％改为90％，90％改为95％，95％改为97.5％；表1：正态分布单变量平均值调查所需样本数查询表(3)可燃物调查所需样线数的计算设置1.8m样线调查直径0～0.6cm可燃物载量；预先调查20条样线，计算预调查样线截交0～0.6cm可燃物的平均点数方差s2和按调查要求精度查表1得到适合的样本数，继续完成剩余样本的调查；(4)可燃物平均直径测量的样本数计算沿1.8m样线用卡尺逐一量取与样线截交的0～0.6cm可燃物的直径，预先调查20个可燃物颗粒后，计算平均直径方差s2和按调查要求精度查表1得到适合的样本数；继续完成剩余样本的调查；(5)可燃物表积比计算可燃物平均直径是计算可燃物载量和表积比的基础，一次测量可重复使用；圆柱形可燃物表积比为4/d(m-1),d-为圆柱形某一径级可燃物的平均直径(m)；叶状可燃物表积比为2/t(m-1)，t-为叶状可燃物厚度(m)；(6)单变量平均值的功效检验功效检验又称为类型Ⅱ错误可能性检验；在森林可燃物外业调查完成后，在室内对合并的预调查数据和后续样本调查数据进行功效检验；森林可燃物调查数据的功效检验在R3.4.1中进行，代码为：power.t.test(delta＝？,sd＝？,n＝？,sig.level＝？,power＝？,type＝"one.sample")代码括号内有5个选择项(delta、sd、n、sig.level、power)，把需要计算的那一项去掉，其余4项在问号处代换具体数据，则可计算出去掉项的数值；当分析功效时，先算出样本数n、样本均值标准差sd，可接受误差范围将n、delta、sd、sig.level值代入以下代码回车即可得到功效power:power.t.test(delta＝？,sd＝？,n＝？,sig.level＝？,type＝"one.sample")。本专利技术具有的优点及有益效果如下：1、本专利技术利用R软件制作了“正态分布单变量平均值调查所需样本数查询表”，用于森林可燃物外业调查现场的样本数计算，为不多不少地抽样调查并得到较为可靠的结果提供技术方法。2、查询表基于t值、预调查均值和方差，原理清晰，既可用于双尾检验也可用于单尾检验；不仅适用于森林可燃物调查，也可用于其他平均值调查研究。3、查询表与原公式应用相比，免除了反复测试不同样本数值、配平公式的繁琐步骤。4、以森林可燃物调查样线数和可燃物平均直径调查为例说明查询表的使用，并以可燃物平均直径计算了可燃物表积比。5、采用R软件power函数对预调查和后续调查合并数据进行功效检验，检验平均值的可靠性。6、R软件在室内比查询表好用，但R与安卓系统不兼容，野外调查用查询表较为轻便。附图说明图1：为在不同置信水平的t值与样本数关系图。具体实施方式下面结合具体实施例对本专利技术做进一步说明。实施例1：可燃物调查所需样线数的计算设置1.8m样线调查直径0～0.6cm可燃物载量；预先调查20条样线，样线截交0～0.6cm可燃物的点数为：30,25,16,9,10,25,40,35,12,15,6,9,10,25,40,35,12,15,8,30；平均数等于20.35，方差131.5，为32.31，如果要求置信水平为90％，查表1得到适合的样本数为90条样线，继续完成剩余70条样线的调查。实施例2：可燃物平均直径测量的样本数计算沿1.8m样线用卡尺逐一量取与样线截交的0～0.6cm可燃物的直径，预先调查20个可燃物颗粒，直径分别为：0.3,0.4,0.6,0.3,0.2,0.6,0.6,0.5,0.3,0.2,0.3,0.2,0.1,0.2,0.3,0.3,0.2,0.6,0.5,0.4；平均直径0.3本文档来自技高网...

【技术保护点】
一种森林可燃物平均值调查的样本数确定方法，其特征在于：包括样本数查询表制作原理、样本数查询表制作方法、可燃物调查所需样线数的计算、可燃物平均直径测量的样本数计算、可燃物表积比计算、单变量平均值的功效检验工序，主要操作步骤如下：(1)样本数查询表制作原理假设样本均值

【技术特征摘要】
1.一种森林可燃物平均值调查的样本数确定方法，其特征在于：包括样本数查询表制作原理、样本数查询表制作方法、可燃物调查所需样线数的计算、可燃物平均直径测量的样本数计算、可燃物表积比计算、单变量平均值的功效检验工序，主要操作步骤如下：(1)样本数查询表制作原理假设样本均值与总体真实均值μ之差为样本均值的k倍，样本数为n,那么去除方根号得到样本数计算公式：通过预调查得到的值；将不同置信水平的n和其对应的列成查询表，根据预调查的值和的关系，查表可以得到可燃物调查所需样本数；(2)样本数查询表制作方法利用R软件输入以下代码制作样本数查询表：n<-1:100；alpha<-0.20；a<-n/(qt(1-alpha/2,n-1))^2write.table(a,"clipboard",sep＝"\t",col.names＝NA)至电子表格粘贴得表1第一列和第二列，修改代码中的alpha值得其他几列；表1为两尾检验所需样本数；表1：正态分布单变量平均值调查所需样本数查询表(3)可燃物调查所需样线数的计算设置1.8m样线调查直径0～0.6cm可燃物载量；预先调查20条样线，计算预调查样线截交0～0.6cm可燃物的平均点数方差s2和按调查要求精度查表1得到适合的样本数，继续完成剩余样本的调查；(4)可燃物平均直径测量的样本数计算沿1.8m样线用卡尺逐一量取与样线截交的0～0.6cm可燃物的直径，预先调查20个可燃物颗粒后，计算平均直径方差s2和按调查要求精度查表1得到适合的样本数；继续完成剩余样本的调查；(...

【专利技术属性】
技术研发人员：黄小荣，彭玉华，侯远瑞，
申请(专利权)人：广西壮族自治区林业科学研究院，
类型：发明
国别省市：广西,45