论文写作中存在的数理统计问题课件

1、报告人：环境科学学报编委会执行副主编、编辑部主任 2007-11-26重要假定n作者所处理的数据属于随机变量的特定样本。作者所处理的数据属于随机变量的特定样本。n作者已经掌握最基本的数理统计学常识，如概率、作者已经掌握最基本的数理统计学常识，如概率、假设检验、均值、方差、标准差、正态分布、相假设检验、均值、方差、标准差、正态分布、相关分析、回归分析、方差分析关分析、回归分析、方差分析。数理统计问题的重要性数理统计问题的重要性n在科学研究中，经常会涉及到对随机变量在科学研究中，经常会涉及到对随机变量大小大小、离散离散及及分布分布特征的描述以及对特征的描述以及对2 2个或多个随机变量之间的个或多个

2、随机变量之间的关系关系描述问题。描述问题。地学、环境科学研究也不例外地学、环境科学研究也不例外。n对随机变量及随机变量之间的关系进行定量描述的数学工具对随机变量及随机变量之间的关系进行定量描述的数学工具就是就是数理统计学数理统计学。n在科学研究中，能否正确使用各种数理统计方法关系到所得在科学研究中，能否正确使用各种数理统计方法关系到所得出结论的客观性和可信性。所以，出结论的客观性和可信性。所以，来稿中使用的数理统计方来稿中使用的数理统计方法是否正确法是否正确应是学术期刊编辑们极为重视的问题。应是学术期刊编辑们极为重视的问题。n目前，国内环境科学与技术类学术期刊对稿件中数理统计方目前，国内环境科

3、学与技术类学术期刊对稿件中数理统计方法问题的重视程度存在差异。法问题的重视程度存在差异。 1 统计软件的选择统计软件的选择 n统计分析通常涉及大量的数据，需要较大的计统计分析通常涉及大量的数据，需要较大的计算工作量。算工作量。n在进行统计分析时，尽管作者可以自行编写计在进行统计分析时，尽管作者可以自行编写计算程序，但在统计软件很普及的今天，这样做算程序，但在统计软件很普及的今天，这样做是毫无必要的。是毫无必要的。n出于对出于对工作效率工作效率以及对以及对算法的通用性、可比性算法的通用性、可比性的考虑，一些学术期刊要求作者采用专门的数的考虑，一些学术期刊要求作者采用专门的数理统计软件进行统计分析

4、。理统计软件进行统计分析。 1 统计软件的选择统计软件的选择n环境科学学报环境科学学报的编辑们在处理稿件时的编辑们在处理稿件时经常发现的问题是：经常发现的问题是：作者未使用专门的数作者未使用专门的数理统计软件，而采用理统计软件，而采用Excel这样的电子表格这样的电子表格软件进行数据统计分析。软件进行数据统计分析。n由于电子表格软件提供的统计分析功能十由于电子表格软件提供的统计分析功能十分有限，只能借助它进行较为简单的统计分有限，只能借助它进行较为简单的统计分析，故我们不主张作者采用这样的软件分析，故我们不主张作者采用这样的软件进行统计分析。进行统计分析。 1 统计软件的选择统计软件的选择n目

5、前，国际上已开发出的专门用于统计分析的商目前，国际上已开发出的专门用于统计分析的商业软件很多，比较著名有业软件很多，比较著名有SPSS(Statistical Package for Social Sciences)和和SAS(Statistical Analysis System)。此外，还有此外，还有BMDP和和STATISTICA等等。nSPSS是专门为社会科学领域的研究者设计的，但是专门为社会科学领域的研究者设计的，但此软件在自然科学领域也得到广泛应用。此软件在自然科学领域也得到广泛应用。nBMDP是专门为生物学和医学领域研究者编制的统是专门为生物学和医学领域研究者编制的统计软件。计软

6、件。 1 1 统计软件的选择统计软件的选择n目前，国际学术界有一条不成文的约定：凡目前，国际学术界有一条不成文的约定：凡是用是用SPSS和和SAS软件进行统计分析所获得的软件进行统计分析所获得的结果，在国际学术交流中不必说明具体算法。结果，在国际学术交流中不必说明具体算法。由此可见，由此可见，SPSS和和SAS软件已被各领域研究软件已被各领域研究者普遍认可。者普遍认可。n我们建议作者们在进行统计分析时尽量使用我们建议作者们在进行统计分析时尽量使用这这2个专门的统计软件。目前，有关这个专门的统计软件。目前，有关这2个软个软件的使用教程在书店中可很容易地买到。件的使用教程在书店中可很容易地买到。

7、2 2 均值的计算均值的计算 ：理论问题n均值（准确的称呼应为均值（准确的称呼应为“样本均值样本均值”）的统计学意义：反）的统计学意义：反映随机变量样本的大小特征。映随机变量样本的大小特征。n均值对应于随机变量总体的数学期望均值对应于随机变量总体的数学期望总体的数学期望客总体的数学期望客观上决定着样本的均值，反过来，通过计算样本的均值可观上决定着样本的均值，反过来，通过计算样本的均值可以描述总体的数学期望。以描述总体的数学期望。n在处理实验数据或采样数据时，经常会遇到对相同采样或在处理实验数据或采样数据时，经常会遇到对相同采样或相同实验条件下同一随机变量的多个不同取值进行统计处相同实验条件下同

8、一随机变量的多个不同取值进行统计处理的问题。理的问题。n为找到代表这些观测值总体大小特征的代表值（统计量，为找到代表这些观测值总体大小特征的代表值（统计量，该统计量根据样本数据算出），多数作者会不假思索地直该统计量根据样本数据算出），多数作者会不假思索地直接给出算术平均值和标准差。显然，这种做法是不严谨接给出算术平均值和标准差。显然，这种做法是不严谨的的不一定总是正确的。不一定总是正确的。2 均值的计算：技术问题均值的计算：技术问题n在数理统计学中，作为描述随机变量样本的在数理统计学中，作为描述随机变量样本的总体大小特征的统计量有算术平均值、几何总体大小特征的统计量有算术平均值、几何平均值和中

9、位数等多个。平均值和中位数等多个。n何时用算术平均值？何时用几何平均值？以何时用算术平均值？何时用几何平均值？以及何时用中位数？及何时用中位数？这不能由研究者根据主观这不能由研究者根据主观意愿随意确定，而要根据随机变量的分布特意愿随意确定，而要根据随机变量的分布特征确定征确定。 2 均值的计算：技术问题均值的计算：技术问题n反映随机变量总体大小特征的统计量是数学期望，而在随机反映随机变量总体大小特征的统计量是数学期望，而在随机变量的分布服从正态分布时，其数学期望就可以用样本的算变量的分布服从正态分布时，其数学期望就可以用样本的算术平均值描述。此时，可用样本的术平均值描述。此时，可用样本的算术平

10、均值算术平均值描述随机变量描述随机变量的大小特征。的大小特征。n如果所研究的随机变量不服从正态分布，则算术平均值不能如果所研究的随机变量不服从正态分布，则算术平均值不能准确反映该变量的大小特征。在这种情况下，可通过假设检准确反映该变量的大小特征。在这种情况下，可通过假设检验来判断随机变量是否服从对数正态分布。如果服从对数正验来判断随机变量是否服从对数正态分布。如果服从对数正态分布，则几何平均值就是数学期望的值。此时，就可以计态分布，则几何平均值就是数学期望的值。此时，就可以计算变量的算变量的几何平均值几何平均值。n如果随机变量既不服从正态分布也不服从对数正态分布，则如果随机变量既不服从正态分布

11、也不服从对数正态分布，则按现有的数理统计学知识，尚无合适的统计量描述该变量的按现有的数理统计学知识，尚无合适的统计量描述该变量的大小特征。此时，可用大小特征。此时，可用中位数中位数来描述变量的大小特征。来描述变量的大小特征。 3 3 相关分析相关分析：相关系数的选择：相关系数的选择 n在相关分析中，作者们常犯的错误是：简单地计算在相关分析中，作者们常犯的错误是：简单地计算Pearson 积矩相关系数，而且既不给出正态分布检验结果，积矩相关系数，而且既不给出正态分布检验结果，也往往不明确指出所计算的相关系数就是也往往不明确指出所计算的相关系数就是Pearson 积矩相积矩相关系数。关系数。n在数

12、理统计学中，除有针对数值变量设计的在数理统计学中，除有针对数值变量设计的Pearson 积矩积矩相关系数（对应于相关系数（对应于 “参数方法参数方法”）外，还有针对顺序变量）外，还有针对顺序变量（即（即“秩变量秩变量”）设计的）设计的Spearman秩相关系数和秩相关系数和Kendall秩相关系数（对应于秩相关系数（对应于 “非参数方法非参数方法”）等。）等。nPearson 积矩相关系数可用于描述积矩相关系数可用于描述2个随机变量的线性相关个随机变量的线性相关程度，程度，Spearman或或Kendall秩相关系数用来判断两个随机秩相关系数用来判断两个随机变量在二维和多维空间中是否具有某种共

13、变趋势。变量在二维和多维空间中是否具有某种共变趋势。 3 相关分析：相关系数的选择相关分析：相关系数的选择n在相关分析中，计算各种相关系数是有前提条件在相关分析中，计算各种相关系数是有前提条件的。的。n在相关分析中，对于秩变量，一般别无选择，只在相关分析中，对于秩变量，一般别无选择，只能计算能计算Spearman或或Kendall秩相关系数。秩相关系数。n对于数值变量，只要条件许可，应尽量使用对于数值变量，只要条件许可，应尽量使用检验检验功效最高功效最高的参数方法，即计算用的参数方法，即计算用Pearson 积矩相积矩相关系数。只有计算关系数。只有计算Pearson 积矩相关系数的前提积矩相关

14、系数的前提不存在时，才考虑退而求其次，计算专门为秩变不存在时，才考虑退而求其次，计算专门为秩变量设计的量设计的Spearman或或Kendall秩相关系数（秩相关系数（尽管尽管这样做会导致检验功效的降低这样做会导致检验功效的降低）。）。 3 相关分析：相关系数的选择相关分析：相关系数的选择n对于数值变量，相关系数选择的依据是变量是否服从正态对于数值变量，相关系数选择的依据是变量是否服从正态分布，或变换后的数据是否服从正态分布。分布，或变换后的数据是否服从正态分布。n对于二元相关分析，如果对于二元相关分析，如果2个随机变量服从二元正态分布假个随机变量服从二元正态分布假设，则应该用设，则应该用Pe

15、arson 积矩相关系数描述这积矩相关系数描述这2个随机变量个随机变量间的相关关系。间的相关关系。n如果样本数据不服从二元正态分布，则可尝试进行数据变如果样本数据不服从二元正态分布，则可尝试进行数据变换，看变换后的数据是否符合正态分布？如果是，则可以换，看变换后的数据是否符合正态分布？如果是，则可以针对变换后的数据计算针对变换后的数据计算Pearson 积矩相关系数；否则，就积矩相关系数；否则，就不能计算不能计算Pearson 积矩相关系数，而应改用检验功效较低积矩相关系数，而应改用检验功效较低的的Spearman或或Kendall秩相关系数（此时，如果强行计秩相关系数（此时，如果强行计算算P

16、earson 积矩相关系数有可能会得出完全错误的结论）。积矩相关系数有可能会得出完全错误的结论）。4 相关分析与回归分析的区别相关分析与回归分析的区别 n相关分析和回归分析是极为常用的相关分析和回归分析是极为常用的2种数理种数理统计方法，在环境科学及其它科学研究领域统计方法，在环境科学及其它科学研究领域有着广泛的用途。然而，由于这有着广泛的用途。然而，由于这2种数理统种数理统计方法在计算方面存在很多相似之处，且在计方法在计算方面存在很多相似之处，且在一些数理统计教科书中没有系统阐明这一些数理统计教科书中没有系统阐明这2种种数理统计方法的内在差别，从而使一些研究数理统计方法的内在差别，从而使一些

17、研究者不能严格区分相关分析与回归分析者不能严格区分相关分析与回归分析 。4 相关分析与回归分析的区别相关分析与回归分析的区别n最常见的错误是最常见的错误是:用回归分析的结果解释相用回归分析的结果解释相关性问题。例如，作者将关性问题。例如，作者将“回归直线（曲线）回归直线（曲线）图图”称为称为“相关性图相关性图”或或“相关关系图相关关系图”；将回归直线的将回归直线的R2(拟合度，或称拟合度，或称“可决系可决系数数”)错误地称为错误地称为“相关系数相关系数”或或“相关系相关系数的平方数的平方”；根据回归分析的结果宣称；根据回归分析的结果宣称2个个变量之间存在正的或负的相关关系。变量之间存在正的或负

18、的相关关系。 4 相关分析与回归分析的区别相关分析与回归分析的区别n相关分析与回归分析均为研究相关分析与回归分析均为研究2个或多个个或多个变量间关联性的方法，但变量间关联性的方法，但2种数理统计方种数理统计方法存在本质的差别，即它们用于不同的法存在本质的差别，即它们用于不同的研究目的。研究目的。n相关分析的目的在于检验两个随机变量相关分析的目的在于检验两个随机变量的共变趋势（即共同变化的程度），回的共变趋势（即共同变化的程度），回归分析的目的则在于试图用自变量来预归分析的目的则在于试图用自变量来预测因变量的值。测因变量的值。 4 相关分析与回归分析的区别相关分析与回归分析的区别n在相关分析中，

19、两个变量必须同时都是在相关分析中，两个变量必须同时都是随机变量，如果其中的一个变量不是随随机变量，如果其中的一个变量不是随机变量，就不能进行相关分析。这是相机变量，就不能进行相关分析。这是相关分析方法本身所决定的。关分析方法本身所决定的。 4 相关分析与回归分析的区别相关分析与回归分析的区别n对于回归分析，其中的因变量肯定为随对于回归分析，其中的因变量肯定为随机变量（这是回归分析方法本身所决定机变量（这是回归分析方法本身所决定的），而自变量则可以是普通变量（有的），而自变量则可以是普通变量（有确定的取值）也可以是随机变量。确定的取值）也可以是随机变量。 4 相关分析与回归分析的区别相关分析与回

20、归分析的区别n如果自变量是普通变量，即模型如果自变量是普通变量，即模型回归分析，回归分析，采用的回归方法就是最为常用的最小二乘法。采用的回归方法就是最为常用的最小二乘法。n如果自变量是随机变量，如果自变量是随机变量，即模型即模型回归分析，回归分析，所采用的回归方法与计算者的目的有关。所采用的回归方法与计算者的目的有关。n在以预测为目的的情况下，仍采用在以预测为目的的情况下，仍采用“最小二乘法最小二乘法”（但精度下降（但精度下降最小二乘法是专为模型最小二乘法是专为模型 设计的，设计的，未考虑自变量的随机误差）；未考虑自变量的随机误差）；n在以估值为目的（如计算可决系数、回归系数等）在以估值为目的

21、（如计算可决系数、回归系数等）的情况下，应使用相对严谨的方法（如的情况下，应使用相对严谨的方法（如“主轴法主轴法”、“约化主轴法约化主轴法”或或“BartlettBartlett法法” ” ）。）。4 相关分析与回归分析的区别相关分析与回归分析的区别n显然，对于回归分析，如果是模型显然，对于回归分析，如果是模型回归分析，鉴于两个回归分析，鉴于两个随机变量客观上存在随机变量客观上存在“相关性相关性”问题，只是由于回归分析问题，只是由于回归分析方法本身不能提供针对自变量和因变量之间相关关系的准方法本身不能提供针对自变量和因变量之间相关关系的准确的检验手段，因此，若以预测为目的，最好不提确的检验手段

22、，因此，若以预测为目的，最好不提“相关相关性性”问题；问题；若以探索两者的若以探索两者的“共变趋势共变趋势”为目的为目的，应该改，应该改用相关分析。用相关分析。n如果是模型如果是模型回归分析，就根本不可能回答变量的回归分析，就根本不可能回答变量的“相关相关性性”问题，问题，因为普通变量与随机变量之间不存在因为普通变量与随机变量之间不存在“相关性相关性”这一概念这一概念（问题在于，大多数的回归分析都是模型（问题在于，大多数的回归分析都是模型回归回归分析！）。此时，即使作者想描述分析！）。此时，即使作者想描述2个变量间的个变量间的“共变趋共变趋势势”而改用相关分析，也会因相关分析的前提不存在而使而

23、改用相关分析，也会因相关分析的前提不存在而使分析结果毫无意义。分析结果毫无意义。4 相关分析与回归分析的区别相关分析与回归分析的区别n需要特别指出的是，回归分析中的需要特别指出的是，回归分析中的R2在数学上恰好是在数学上恰好是Pearson积矩相关系数积矩相关系数r的平方。因此，这极易使作者们错的平方。因此，这极易使作者们错误地理解误地理解R2的含义，认为的含义，认为R2就是就是 “相关系数相关系数”或或“相关系相关系数的平方数的平方”。问题在于，对于自变量是普通变量（即其取值。问题在于，对于自变量是普通变量（即其取值有确定性的变量）、因变量为随机变量的模型有确定性的变量）、因变量为随机变量的

24、模型回归分析，回归分析，2个变量之间的个变量之间的“相关性相关性”概念根本不存在，又何谈概念根本不存在，又何谈“相关相关系数系数”呢？呢？n更值得注意的是，一些早期的教科书作者不是用更值得注意的是，一些早期的教科书作者不是用R2来描述回来描述回归效果（拟合程度，拟合度）的，而是用归效果（拟合程度，拟合度）的，而是用Pearson积矩相关积矩相关系数来描述。这就更容易误导读者。系数来描述。这就更容易误导读者。5 重要的数理统计学常识n1）假设检验n假设检验的基本思想假设检验的基本思想n统计推断统计推断：是根据样本数据推断总体特征的一种方法。：是根据样本数据推断总体特征的一种方法。n假设检验：假设

25、检验：是进行是进行统计推断统计推断的途径之一（另一种途径是参数估计，的途径之一（另一种途径是参数估计，如点估计和区间估计）。如点估计和区间估计）。n假设检验的基本思路是假设检验的基本思路是：首先，对总体参数值提出假设（原假设）；：首先，对总体参数值提出假设（原假设）；然后，利用样本数据提供的信息来验证所提出的假设是否成立（统然后，利用样本数据提供的信息来验证所提出的假设是否成立（统计推断）计推断）-如果样本数据提供的信息不能证明上述假设成立，则应如果样本数据提供的信息不能证明上述假设成立，则应拒绝该假设；如果样本数据提供的信息不能证明上述假设不成立，拒绝该假设；如果样本数据提供的信息不能证明上

26、述假设不成立，则不应拒绝该假设。则不应拒绝该假设。n接受或拒绝原假设的依据接受或拒绝原假设的依据：小概率事件不可能发生。显然，这样做：小概率事件不可能发生。显然，这样做是有风险的（小概率事件真的发生了）。是有风险的（小概率事件真的发生了）。n假设检验中的关键问题假设检验中的关键问题：1）在原假设成立的情况下，如何计算样）在原假设成立的情况下，如何计算样本值或某一极端值发生的概率？本值或某一极端值发生的概率？2）如何界定小概率事件？）如何界定小概率事件？5 重要的数理统计学常识n1）假设检验n假设检验的基本步骤n1）提出原假设（或称）提出原假设（或称“零假设零假设”，H0）；）；n2）选择检验统

27、计量；）选择检验统计量；n3）根据样本数据计算检验统计量观测值的发）根据样本数据计算检验统计量观测值的发生概率（相伴概率，生概率（相伴概率，p）；）；n4）根据给定的小概率事件界定标准（显著性）根据给定的小概率事件界定标准（显著性水平，如水平，如0.05，0.01）做出统计推断。）做出统计推断。假设检验的基本步骤n为什么要设计并计算检验统计量？为什么要设计并计算检验统计量？n在假设检验中，样本值（或更极端的取值）发生的概率在假设检验中，样本值（或更极端的取值）发生的概率不能直接通过样本数据计算，而是通过计算不能直接通过样本数据计算，而是通过计算检验统计量检验统计量观测值观测值的发生概率而间接得

28、到的。的发生概率而间接得到的。n所设计的检验统计量一般服从或近似服从某种已知的理所设计的检验统计量一般服从或近似服从某种已知的理论分布（如论分布（如t-分布、分布、F-分布、卡方分布），易于估算其分布、卡方分布），易于估算其取值概率。取值概率。n对于不同的假设检验和不同的总体，会有不同的选择检对于不同的假设检验和不同的总体，会有不同的选择检验统计量的理论和方法验统计量的理论和方法。假设检验的基本步骤n计算检验统计量观测值的发生概率计算检验统计量观测值的发生概率n在假定原假设成立的前提下，利用样本数据计算检验统计量观测值发生的概率（即p值，又称“相伴概率”指该检验统计量在某个特定的极端区域在原假

29、设成立时的概率）。该概率值间接地给出了在原假设成立的条件下样本值（或更极端值）发生的概率。假设检验的基本步骤n进行统计推断进行统计推断n依据预先确定的 “显著性水平” （即值），如0.01或0.05，决定是否拒绝原假设。n如果p值小于值，即认为原假设成立时检验统计量观测值的发生是小概率事件，则拒绝原假设。否则，就接受原假设。显著性水平：概念与意义显著性水平：概念与意义n在假设检验中，在假设检验中，显著性水平（显著性水平（Significant level，用用表示）的确定是假设检验中至关重要的问题。表示）的确定是假设检验中至关重要的问题。n显著性水平是在原假设成立时检验统计量的制落在显著性水平

30、是在原假设成立时检验统计量的制落在某个极端区域的概率值。因此，如果取某个极端区域的概率值。因此，如果取= 0.05，如果计算出的如果计算出的p值小于值小于 ，则可认为原假设是一个，则可认为原假设是一个不可能发生的小概率事件。当然，如果真的发生了，不可能发生的小概率事件。当然，如果真的发生了，则犯错误的可能性为则犯错误的可能性为5%。显然，显著性水平反映显然，显著性水平反映了拒绝某一原假设时所犯错误的可能性，或者说，了拒绝某一原假设时所犯错误的可能性，或者说， 是指拒绝了事实上正确的原假设的概率。是指拒绝了事实上正确的原假设的概率。显著性水平：通常的取值显著性水平：通常的取值n值一般在进行假设检

31、验前由研究者根据实际的需值一般在进行假设检验前由研究者根据实际的需要确定。要确定。n常用的取值是常用的取值是0.05或或0.01。对于前者，相当于在。对于前者，相当于在原假设事实上正确的情况下，研究者接受这一假原假设事实上正确的情况下，研究者接受这一假设的可能性为设的可能性为95%；对于后者，则研究者接受事；对于后者，则研究者接受事实上正确的原假设的可能性为实上正确的原假设的可能性为99%。n显然，降低显然，降低值可以减少拒绝原假设的可能性。因值可以减少拒绝原假设的可能性。因此，在报告统计分析结果时，必须给出此，在报告统计分析结果时，必须给出值。值。 显著性水平：进行统计推断显著性水平：进行统

32、计推断n在进行假设检验时，各种统计软件均会给出在进行假设检验时，各种统计软件均会给出检验统检验统计量观测值计量观测值以及原假设成立时该检验统计量取值的以及原假设成立时该检验统计量取值的相伴概率相伴概率（即（即检验统计量检验统计量某特定取值及更极端可能某特定取值及更极端可能值出现的概率，用值出现的概率，用p p表示）。表示）。np p值是否小于事先确定的值是否小于事先确定的值，是接受或拒绝原假值，是接受或拒绝原假设的依据。设的依据。n如果如果p p值小于事先已确定的值小于事先已确定的值，就意味着检验统值，就意味着检验统计量取值的可能性很小，进而可推断原假设成立的计量取值的可能性很小，进而可推断原

33、假设成立的可能性很小，因而可以拒绝原假设。相反，如果可能性很小，因而可以拒绝原假设。相反，如果p p值大于事先已确定的值大于事先已确定的值，就不能拒绝原假设。值，就不能拒绝原假设。 统计推断：过去的回忆统计推断：过去的回忆n在计算机技术十分发达，以及专业统计软件功能十分强大的今天，计算在计算机技术十分发达，以及专业统计软件功能十分强大的今天，计算检验统计量及其相伴概率是一件十分容易的事情。检验统计量及其相伴概率是一件十分容易的事情。n然而，在然而，在20世纪世纪90年代以前，只有服从标准正态分布的检验统计量，年代以前，只有服从标准正态分布的检验统计量，人们可以直接查阅事先准备好的人们可以直接查

34、阅事先准备好的标准正态分布函数表标准正态分布函数表，从中获得特定计，从中获得特定计算结果的相伴概率。而对于的服从算结果的相伴概率。而对于的服从t-分布、分布、F-分布、卡方分布或其它特分布、卡方分布或其它特殊的理论分布的检验统计量（大多数的假设检验是这样），人们无法直殊的理论分布的检验统计量（大多数的假设检验是这样），人们无法直接计算相伴概率。人们通常查阅各类假设检验的临界值表进行统计推断。接计算相伴概率。人们通常查阅各类假设检验的临界值表进行统计推断。这些表格以自由度和很少的几个相伴概率（通常为这些表格以自由度和很少的几个相伴概率（通常为0.1、0.05和和0.01）为自变量，以检验统计量的

35、临界值为函数排列。为自变量，以检验统计量的临界值为函数排列。n在进行统计推断时，人们使用上述临界值表根据事先确定的显著性水平，在进行统计推断时，人们使用上述临界值表根据事先确定的显著性水平，查阅对应于某一自由度和特定相伴概率的检验统计量的临界值，然后将查阅对应于某一自由度和特定相伴概率的检验统计量的临界值，然后将所计算出的检验统计量与该临界值相比较。如果检验统计量的计算值大所计算出的检验统计量与该临界值相比较。如果检验统计量的计算值大于临界值，即实际的相伴概率小于事先规定的显著性水平，便可拒绝原于临界值，即实际的相伴概率小于事先规定的显著性水平，便可拒绝原假设。否则，可接受原假设。假设。否则，

36、可接受原假设。显著性水平：举例显著性水平：举例n在根据显著性水平进行统计推断时，应注意原假设的性质。在根据显著性水平进行统计推断时，应注意原假设的性质。n以二元相关分析为例，相关分析中的原假设是以二元相关分析为例，相关分析中的原假设是“相关系数为相关系数为零零”（即（即2个随机变量间不存在显著的相关关系）。如果计个随机变量间不存在显著的相关关系）。如果计算出的检验统计量的相伴概率（算出的检验统计量的相伴概率（p值）低于事先给定值）低于事先给定值值（如（如0.05），就可以认为），就可以认为“相关系数为零相关系数为零”的可能性很低，的可能性很低， 既既2个随机变量之间存在显著的相关关系。个随机变

37、量之间存在显著的相关关系。n在正态分布检验时，原假设是在正态分布检验时，原假设是“样本数据来自服从正态分布样本数据来自服从正态分布的总体的总体”。此时，如果计算出的检验统计量的相伴概率（。此时，如果计算出的检验统计量的相伴概率（p值）低于事先给定值）低于事先给定值（如值（如0.05），则表明数据不服从正态），则表明数据不服从正态分布。只有分布。只有p值高于值高于值时，数据才服从正态分布。值时，数据才服从正态分布。这与相这与相关分析的假设检验不同。关分析的假设检验不同。显著性水平显著性水平n作者在描述相关分析结果时常有的失误是：仅给出相关系数作者在描述相关分析结果时常有的失误是：仅给出相关系数的

38、值，而不给出显著性水平。这就无法判断的值，而不给出显著性水平。这就无法判断2个随机变量间个随机变量间的相关性是否显著。的相关性是否显著。n有时作者不是根据显著性水平判断相关关系是否显著，而是有时作者不是根据显著性水平判断相关关系是否显著，而是根据相关系数的大小来推断（相关系数越近根据相关系数的大小来推断（相关系数越近1，则相关关系，则相关关系越显著）。问题是，相关系数本身是一个基于样本数据计算越显著）。问题是，相关系数本身是一个基于样本数据计算出的观测值，其本身的可靠性尚需检验。出的观测值，其本身的可靠性尚需检验。n此外，作者在论文中常常用此外，作者在论文中常常用“显著相关显著相关”和和“极显

39、著相关极显著相关”来描述相关分析结果，即认为来描述相关分析结果，即认为p值小于值小于0.05就是显著相关关就是显著相关关系（或显著相关），小于系（或显著相关），小于0.01就是极显著相关关系（或极就是极显著相关关系（或极显著相关）。显著相关）。 统计推断的注意事项统计推断的注意事项n在假设检验中，只有在假设检验中，只有 “显著显著”和和 “不显著不显著”，没，没有有“极显著极显著”这样的断语。只要计算出的检验统这样的断语。只要计算出的检验统计量的相伴概率（计量的相伴概率（p值）低于事先确定的值）低于事先确定的值，就值，就可以认为检验结果可以认为检验结果“显著显著”（相关分析的原假设（相关分析的

40、原假设是是“相关系数为零相关系数为零”，故此处的，故此处的“显著显著”实际意实际意味着味着“相关系数不为零相关系数不为零”，或说，或说“2个随机变量间个随机变量间有显著的相关关系有显著的相关关系”）；同样，只要计算出的检）；同样，只要计算出的检验统计量的相伴概率（验统计量的相伴概率（p值）高于事先确定的值）高于事先确定的值，值，就可以认为检验结果就可以认为检验结果“不显著不显著”。 统计推断的注意事项统计推断的注意事项n在进行相关分析时，不能同时使用在进行相关分析时，不能同时使用0.05和和0.01这这2个显著性水平来决定是否拒个显著性水平来决定是否拒绝原假设，只能使用其中的绝原假设，只能使用

41、其中的1个。个。 有关相关分析的断语有关相关分析的断语n1）显著和不显著：描述相关关系是否存在。n2）相关性强或不强：在存在相关关系的前提下，这种相关关系的强或弱。可以认为，相关系数越接近1，则相关性越强。n声明：第声明：第1 1）条是公认的数理统计常识，但第）条是公认的数理统计常识，但第2 2）条是个人理）条是个人理解，仅供参考。本文不对第解，仅供参考。本文不对第2 2）条承担责任。）条承担责任。5 重要的数理统计学常识n1）假设检验n统计推断：单侧检验与双侧检验n对于假设检验，其检验统计量的异常取值有2个方向，即概率分布曲线的左侧（对应于过小的值）和右侧（对应于过大的值）。检验统计量的极端

42、取值n检验统计量在左侧和右侧均有可能取值检验统计量的取值空间单侧检验与双侧检验n一般情况下，概率分布函数曲线两侧尾端的小概率事件都要考虑（即一般情况下，概率分布函数曲线两侧尾端的小概率事件都要考虑（即双侧检验）。如果事先有把握确定其中的一侧不可能取值，则仅需对双侧检验）。如果事先有把握确定其中的一侧不可能取值，则仅需对另一侧的小概率事件进行检验即可（单侧检验）。另一侧的小概率事件进行检验即可（单侧检验）。n在用在用 “ “查表法查表法”进行统计推断时，基于单侧小概率事件检验的临界值进行统计推断时，基于单侧小概率事件检验的临界值表称表称“单尾表单尾表”，基于双侧小概率事件检验的临界值表称，基于双

43、侧小概率事件检验的临界值表称“双尾表双尾表”。除除t-t-分布临界值表是双尾表外，大多数的检验临界值表均为单尾表分布临界值表是双尾表外，大多数的检验临界值表均为单尾表。n在显著性水平一定的情况下（例如在显著性水平一定的情况下（例如 =0.05 =0.05），对于单尾表，单侧检），对于单尾表，单侧检验时仍使用验时仍使用进行统计推断，双侧检验则用进行统计推断，双侧检验则用 /2 /2进行统计推断；对于进行统计推断；对于双尾表，单侧检验时改用双尾表，单侧检验时改用2 2进行统计推断，双侧检验则用进行统计推断，双侧检验则用 进行统进行统计推断。计推断。n在统计软件（如在统计软件（如SPSS或或SAS统

44、计软件）给出的计算结果中，已标注出统计软件）给出的计算结果中，已标注出所计算的相伴概率是单侧还是双侧，对应于上述的单尾表和双尾表。所计算的相伴概率是单侧还是双侧，对应于上述的单尾表和双尾表。单侧检验与双侧检验n以下是以下是SPSS 中的单样本中的单样本t检验输出结果：检验输出结果：nOne-Sample Test（原假设：储户1次平均存取的现金与2000元无显著差异）nTest Value=2000（均值比较的参比值）nt=1.240(检验统计量的观测值)ndf=312(自由度，样本量N=313)nSig.(2-tailed)=0.216（双侧相伴概率p ）nMean Difference=4

45、73.78（均值的标准误差）n95% Confidence Interval of the Difference（总体均值与原假设值之差的95%的置信区间）:-278.131225.69（有95%的把握可认为：储户1次平均存取的金额为1721.873225.69元）n上述检验属 “均值比较”，是双侧检验（大于或小于2000元都算拒绝原假设），计算的相伴概率也是双侧的。因此，可直接用p与比较。取=0.05,则因p大于，故不能拒绝原假设（不是小概率事件）。统计推断结果：根据313个储户调查数据，每个储户一次平均存取金额大体为2000元。单侧检验与双侧检验n在统计软件中，可通过选择Test of S

46、ignificance选项来控制所输出的相伴概率是单尾（1 tailed）概率还是双尾（2 tailed ）概率。5 重要的数理统计学常识n2）正态分布检验n目的：检验样本是否来自正态分布的总体n原假设：样本来自正态分布的总体n分布检验只能使用非参数方法（只有分布形式已知时才能使用参数方法）。n不同的统计软件给出了不同的检验方法。正态分布检验n在SAS中，提供了Shapiro-Wilk(适用于样本量小于50的情形)检验法。此检验无单尾、双尾之分。n在SPSS中提供了卡方检验（Chi-Square Test）和单样本的 Kolmogorov-Smirnov(柯尔莫哥洛夫-斯米尔诺夫，简称K-S)

47、检验。后者比前者精确一些，建议采用。正态分布检验n单样本的 Kolmogorov-Smirnov(柯尔莫哥洛夫-斯米尔诺夫，简称K-S)检验属于双侧检验，计算检验统计量（Z）的双尾概率。5 重要的数理统计学常识n3)均值比较na)将样本均值与某一特定值相比：t-检验（参数检验）n原假设：总体均值与特定值无显著差异n前提：样本来自正态分布的总体n双侧检验：是否等于。n单侧检验：已知不可能大于（或不可能小于），检验是否等于。nb)比较2个独立样本均值： t-检验（参数检验）n原假设：2个样本所代表的2个总体的均值无显著差异n用于对2个来自正态分布总体的样本的大小进行比较，且2个样本相互独立（无相关

48、关系）。n改检验有单侧和双侧之分。5 重要的数理统计学常识n3)均值比较nC) 比较2个独立总体大小的非参数检验n适用于对2个顺序变量的大小进行比较或对2个不服从正态分布的数值变量的大小进行比较n“Mann-Whitney U” 检验:适合样本量较大的样本。n “Wilcoxon秩和”检验：与“Mann-Whitney U” 检验在本质上完全等价。nKolmogorov单侧检验：适用于样本量较小的样本。5 重要的数理统计学常识n3)均值比较nd)比较多个来自正态分布总体的样本均值的检验方法：单因子方差分析（single-factor anova）。n对于将因子作为固定处理（而不是随机变量）的情

49、形，即模型1单因子方差分析，实际上可以看作比较2个总体均值的t-检验的直接推广。该方法属于参数检验。n有关假定：多个样本相互独立、样本均服从正态分布、方差同质性（各个样本的方差大小没有显著差异）等。n原假设：各样本的均值间无显著差异，即某影响因子的不同取值（等级）对各样本的大小没有影响。5 重要的数理统计学常识n3)均值比较nd)比较多个来自非正态分布总体的样本均值的检验方法：nKruskal-Wallis检验：该方法基于顺序变量设计，用于检验3个以上独立样本是否来自大小相同的总体，是应用最广泛的非参数检验方法。n推广的中位数检验：用于检验3个以上的独立样本是否来自中位数无显著差异的样本。该方

50、法检验功效低，不推荐采用。原假设：各独立样本所代表的总体的中位数无显著差异。nFriedman秩方差分析：用于检验3个以上相关样本是否来自大小相同的总体。案例1n2.4 Cd、Pb之间的交互作用之间的交互作用n如表如表4所示，三种花卉植物各部位对重金属所示，三种花卉植物各部位对重金属Cd、Pb的积累量与培养溶的积累量与培养溶液中所投加的液中所投加的Cd、Pb量之间，可以很恰当地被各量之间，可以很恰当地被各多元回归方程多元回归方程表示出表示出来，它们之间呈来，它们之间呈极显著相关关系极显著相关关系（P0.05p0.05）; ;n有相关性，显著（有相关性，显著（0.05p0.010.05p0.01

51、）; ;n有相关性，极显著（有相关性，极显著（p0.01p0.200）、珠江口（）、珠江口（p0.091）和澳门水域）和澳门水域（p0.110）呈正态分布（）呈正态分布（=0.05）.因此对珠江、珠江口和澳因此对珠江、珠江口和澳门水域进行门水域进行Pearson相关分析，对东江、西江和南海北部海域进相关分析，对东江、西江和南海北部海域进行行Kendall相关分析相关分析.从表从表2可以看出，除澳门水域外其它研究区域，可以看出，除澳门水域外其它研究区域，BDE209与与PBDEs相关性不显著（相关性不显著（r0.047）,这是这是由于由于BDE209与其它与其它PBDEs同系物分别来自不同的溴代

52、阻燃剂；同系物分别来自不同的溴代阻燃剂；但澳门水域沉积物中的但澳门水域沉积物中的BDE209与与PBDEs相关性显著（相关性显著（r=0.955，p=0）（图）（图5），表明澳门水域），表明澳门水域BDE209和其它其它和其它其它PBDEs同系物同系物具有相同的输入途径，正如上述，它们主要都是通过水体中颗粒具有相同的输入途径，正如上述，它们主要都是通过水体中颗粒物输入的，它们之间较高的相关性是物输入的，它们之间较高的相关性是PBDEs在水体颗粒物中再分在水体颗粒物中再分配的结果，这也证实了澳门水域是珠三角水体环境中配的结果，这也证实了澳门水域是珠三角水体环境中PBDEs的的“汇汇”. 案例3图

53、2b表明，1/qN对1/D有很好的线性关系，相关系数R=0.9922（R2=0.9845）。 案例4nHA对有机农药甲基对硫磷、西维因、克百威的吸附等温线见图3，用线性吸附方程拟合甲基对硫磷、西维因和克百威的吸附等温线，拟合结果见表3。线性吸附方程为：n Q = KdCe + A （1）n式中Q为吸附量（mg/kg）；Ce为平衡浓度（mg/l）；Kd为线性吸附平衡常数，A为线性方程待定常数。用HA总有机碳标化有机农药的吸附系数Kd得有机碳标化吸附系数Koc，有机农药在HA上的Koc值见表3。由图3、表3可知，甲基对硫磷、西维因和克百威在HA上的吸附等温线较好的符合线性吸附方程，相关系数在0.8

54、7480.9940之间，但是克百威的相关系数要小于甲基对硫磷和西维因；从整体上看有机农药在HA上的Kd大小顺序为：水解处理HA原始HA肟化处理HA氧化处理HA。 案例6在下表中，作者将回归方程的可决系数误称为“相关系数”。案例7n早期的研究表明有机污染物通过分配作用吸附到土壤/沉积物有机质上，其吸附量与有机碳含量和有机污染物的辛醇-水分配系数成正比3。从甲基对硫磷、西维因和克百威分配系数Kd与改性HA有机碳含量的关系可知（见图4）， Kd与HA的有机碳含量成正比，但相关性不高分别为：0.7429、0.8870和0.6900，这表明有机农药在HA上的吸附行为不是由HA的有机碳含量唯一确定，还受到

55、其他因素的影响。图5为有机农药在处理前后HA上n的有机碳标化吸附系数Koc对数（lgKoc）与三种有机农药辛醇-水分配系数Kow对数（lgKow）之间的关系曲线，lgKoc与lgKow呈现较好的线性关系，相关系数分别为：0.8573、0.8367、0.8420和0.9408，可见用辛醇-水分配系数来预测有机污染物在土壤/沉积物上的吸附具有一定的合理性3。 案例8 0500100015002000250035404550甲基对硫磷010020030040050035404550西维因0153045607535404550克百威O / %Kd050010001500200025000.060.07

56、0.080.09甲基对硫磷01002003004005000.060.070.080.09西维因015304560750.060.070.080.09克百威H/C / %Kd图6 Kd与腐殖酸O元素含量和H/C比的相关性Fig 6 Correlation of Kd and O content and H/C rate of the humic acids 案例9图4为取每天19：00的DO值与叶绿素值做的趋势图，通过分析它们数据得出它们的相关性为0.8899，在一定程度上能反映藻类的变化趋势。可以作为藻类增长趋势的预报指标。 案例10n作图得到一条直线见图6，二级动力学速率方程可很好的描述C

57、u2、Cd2在生物膜上的吸附(RCu=0.9989，RCd0.9978)。 案例11案例11（续）n 由表由表1 1可知可知,0,05cm5cm土壤层中土壤层中, ,活动区土壤微生物生物量碳和缓活动区土壤微生物生物量碳和缓冲区土壤微生物生物量碳分别比背景区土壤微生物生物量碳冲区土壤微生物生物量碳分别比背景区土壤微生物生物量碳降低了降低了65.96%65.96%和和20.05%20.05%，而活动区土壤微生物生物量碳比缓，而活动区土壤微生物生物量碳比缓冲区土壤微生物生物量碳降低了冲区土壤微生物生物量碳降低了57.42%57.42%，并且并且3 3个试验区的差个试验区的差异均达到显著水平（异均达到

58、显著水平（P0.05P0.05）. 5. 515cm15cm土壤层中土壤层中, ,活动区土活动区土壤微生物生物量碳比缓冲区土壤微生物生物量碳降低了壤微生物生物量碳比缓冲区土壤微生物生物量碳降低了43.14%43.14%，而缓冲区土壤微生物生物量碳比背景区土壤微生物，而缓冲区土壤微生物生物量碳比背景区土壤微生物生物量碳降低了生物量碳降低了13.85%13.85%，3 3个试验区的差异也均达到显著水平个试验区的差异也均达到显著水平（P0.05P0.05）. 15. 1525cm25cm土壤层中土壤层中, ,活动区土壤微生物生物量碳活动区土壤微生物生物量碳比缓冲区土壤微生物生物量碳降低了比缓冲区土壤

59、微生物生物量碳降低了18.58%18.58%，而缓冲区土壤，而缓冲区土壤微生物生物量碳只比背景区土壤微生物生物量碳降低了微生物生物量碳只比背景区土壤微生物生物量碳降低了11.06%11.06%，但，但3 3个试验区的差异均达到显著水平（个试验区的差异均达到显著水平（P0.05P0.05）. . 案例11（续）n由表由表2可知，在可知，在05cm土壤层和土壤层和515cm土土壤层壤层,旅游踩踏对土壤微生物生物量氮的影响旅游踩踏对土壤微生物生物量氮的影响与对土壤微生物生物量碳的影响是相似的与对土壤微生物生物量碳的影响是相似的.但但在在1525cm土壤层，活动区土壤微生物生物土壤层，活动区土壤微生物

60、生物量氮比背景区土壤微生物生物量氮低，并且达量氮比背景区土壤微生物生物量氮低，并且达到显著水平（到显著水平（P0.05）；缓冲区土壤微生物）；缓冲区土壤微生物生物量氮与活动区土壤微生物生物量氮的差异生物量氮与活动区土壤微生物生物量氮的差异也达到显著水平（也达到显著水平（P0.05） .案例12n2.10 相关性分析n所有相关数据分析，通过SPSS10.0软件分析完成,采用t测验法检验相关系数的显著性。 案例12（续）表2：喷洒菌株与TSNA与硝酸盐、亚硝酸盐的相关性及显著性分析Table2: The correlation and significant analysis of sprayin