定量古生物学中的常用Q型多变量线性回归方法例析

提要Q型多变量线性回归方法是定量古生物学研究中经常使用的一组方法。不同于把研究重心放在变量或变量组之间关系的R型回归方法,Q型回归方法对样本进行整体分析,通过构建样本之间的相似性系数矩阵或者距离矩阵,并进行分解与变换,进而用可视化的方法展示样本的整体相似性关系。本文将介绍常用的Q型多变量线性回归方法,包括主坐标分析(Principal coordinate analysis)、Q型聚类分析(Q-Cluster analysis)、Q型因子分析(Q-Factor analysis),以及综合了R型和Q型的对应分析(Correspondence analysis),和它们在计算中的注意事项与适用场景。本文选取瓜迪亚纳河入海口大陆架地区底栖有孔虫分布特征研究作为实例,分别用上述方法进行分析,以解释它们在定量古生物学中的具体应用。     

我国实施医师多点执业的政策分析与思考

医师多点执业试点工作正在各地逐步推开。文章依据国家新医改政策,以医师多点执业为基点,对医师多点执业现状和对医疗卫生系统的运行管理可能带来的影响进行综合分析,结合医疗单位的实际,提出意见建议。     

多腺巨伪蝎(伪蝎目:苦伪蝎科)的雌性新发现（英文）

首次报道了采自中国贵州宽阔水国家级自然保护区一洞穴的多腺巨伪蝎Megachernes glandulosus Mahnert,2009雌性,并对其形态特征进行了详细描述。标本保存于河北大学博物馆(MHBU)。多腺巨伪蝎Megachernes glandulosus Mahnert,2009(图1~3)鉴别特征:该种触肢非常粗短,触肢膝节中央明显具1突起,腿节和膝节着生大量浓密的长羽状刚毛;背甲刚毛正常;基节、第Ⅳ~Ⅸ每半个腹板和触肢螯固定指基部内侧面具大量腺体刚毛,第三步足和第四步足腿节和膝节腹面着生有大量腺体微刚毛。     

太行山猕猴颅骨的形态学及多变量研究

对２１例成年太行山猕猴颅骨的主要形态特征进行了描述和测量,并进行了初步的统计学分析。结果表明,成年太行山猕猴颅骨在耳点间宽（Ｘ１３）、上面宽（Ｘ１９）、左眶宽（Ｘ２７）、右眶宽（Ｘ２８）、两眶宽（Ｘ３７）等变量间存在显的性别差异。与湘鄂、闽粤猕猴颅骨有关变量相比,有显的不同。３个猕猴种群的聚类分析结果显示,湘鄂、闽粤猕猴种群间的生物学距离相对较近,而太行山猕猴与前二种种群间的距离均较远。太行山     

单个协变量带有测量误差的多维结构回归模型

给出可交换条件下单个协变量的带有测量误差的多维结构回归模型,利用该模型研究总体平均处理效应的估计,给出当暴露组和对照组的协变量测量误差同分布时总体平均处理效应的拟极大似然估计及其性质.     

南丰蜜桔可溶性固形物近红外检测的多元线性回归变量筛选

文章采用反向区间偏最小二乘法结合连续投影算法,筛选南丰蜜桔近红外检测的多元线性回归变量。对南丰蜜桔近红外光谱进行多元散射校正后,利用反向间隔偏最小二乘法,从500～1750 nm中初选出7个光谱区间,用于多元线性回归变量筛选。利用通过遗传算法和连续投影算法筛选出的变量建立了多元线性回归模型。经比较发现,利用反向区间偏最小二乘法结合连续投影算法筛选出的变量建立的多元线性回归模型,预测结果最优,模型预测相关系数为0.937,模型预测均方根误差为0.613 oBrix。结果表明,反向区间偏最小二乘法结合连续投影算法,可以有效地筛选近红外光谱的多元线性回归变量,提高南丰蜜桔可溶性固形物模型的预测精度。     

基于哑变量和分位数回归的兴安落叶松更新幼树的树高-胸径模型

本研究以大兴安岭地区翠岗林场2018—2019年55块固定样地2054株兴安落叶松幼树为对象,采用四分位数法将林分密度指数(SDI)划分为4个等级,即等级Ⅰ(SDI₁<1863株·hm^-2)、等级Ⅱ(1863≤SDI₂<2155株·hm^-2)、等级Ⅲ(2155≤SDI₃<2459株·hm^-2)和等级Ⅳ(SDI₄≥2459株·hm^-2),并采用哑变量方法引入SDI构建兴安落叶松幼树树高-胸径的哑变量模型和分位数回归模型。结果表明：选取的5个代表性非线性树高曲线模型中,Richards模型的拟合效果最好,其R_a²、RMSE、MAE分别为0.7637、0.8250 m、0.5696 m;基于Richards模型构建的包含SDI的哑变量模型,其R_a²较基础模型提高了1.3%,而RMSE、MAE、AIC分别降低了2.1%、1.5%和11.2%;当分位点τ=0.5时,分位数回归模型的R_a²最大,RMSE、MAE、AIC最小,分别为0.7612、0.8294 m、0.5657 m、-767.19。相较于SDI₁,SDI₂～SDI₄林分内幼树的树高分别增加5.6%、5.6%和11.3%。因此,合理调控兴安落叶松林的林分密度有利于增加更新幼树的高生长。     

不同化学型樟树叶挥发性成分组成的多变量分析

为明确不同化学型樟树(Cinnamomum camphora)叶挥发性成分的异同,以5种不同樟树叶的25个样本为实验材料,采用静态顶空-气相色谱-质谱联用法(SHS-GC-MS)研究了其挥发性成分组成,并运用主成分分析(PCA)、聚类分析(CA)和判别分析(DA)等对不同化学型樟树叶挥发性成分组成进行了多变量分析。研究结果表明,不同化学型樟树叶的挥发性成分共有8种:β-芳樟醇、莰烯、β-愈创木烯、γ-松油烯、α-侧柏烯、2-乙基呋喃、α-石竹烯和大牛儿烯;不同化学型樟树叶的挥发性化合物种类和含量存在显著差异。检测出的平均化合物数量分别为52(异樟)、40(脑樟)、37(油樟)、34(芳樟)和33(龙脑樟)。主成分分析提取了3个主成分因子,累计贡献率达到86.40%。聚类分析和判别分析结果表明,其中的12种挥发性成分能够对5种化学型樟树进行较好的区分,准确率分别为96%和100%。实验结果说明,这12种挥发性成分可以较好地用于樟树化学型分类。     

社鼠种群年龄结构的回归分析

本文利用社鼠(Nivitverter confucianus)头骨度量的8项指标(单位：mm)：颅全长X1、颅基长X2、基底长X3、颧宽X4、眶间宽X5、齿隙长X6、上裂齿长X7、门齿孔长X8,用多元回归分析方法,建立了其种群年龄结构的判别函数为：Y=0．3577X2-0．4599X7 0．3364X8—8．4267将头骨度量X2、X7、X8带入回归方程进行计算,如INT(Y 0．5)=G1,则该个体属于G1年龄组。回归方程表明,社鼠的年龄仅与头骨的颅基长X2、上裂齿长X7、门齿孔长X8有关,而颧宽X4、眶间宽X5、齿隙长X6在各年龄组相对较稳定,可以作为分类的指标。     

利用变量分析探讨西双版纳州勐海等地鸟类生境变化

郑作新,郑宝赉(1961,1962),依万诺夫(1961)等在西双版纳的勐海、勐养、勐腊进行了鸟类资源调查。经二十多年的开发利用后,1983年作者在上述地区观察、采集、记录这些地区的留鸟种数,利用变量分析来对两次的留鸟种数进行比较分析,通过种数的差异性变化分布来反映三地开发后的鸟类生境变化。     

A Note on Maximum Likelihood Ratio Test of No Outliers in Regression Models

A derivation of the maximum likelihood ratio test for testing no outliers in regression models is given using the method of WETHERILL (1981, pp. 106–107) for estimating the regression parameters. This method is essentially similar to the one outlined in BARNETT and LEWIS (1978, p. 263), although by our detailed derivation it is easier to see that the maximum likelihood estimate of θ of model (3) under the hypothesis that the i^th observation in an outlier is the same as that obtained from model (1) when the i^th observation is removed.     

ELISA法检测HBeAg假性结果原因分析

目的:探讨ELISA法检测HBeAg假性结果原因方法:用ELISA法检测乙肝血清标志物,对HBsAg阳性而HBeAg阴性的标本以及HBeAg阳性的标本用ELISA法和电化学发光法复查。结果:136例HBsAg阳性而HBeAg阴性的标本经稀释复查后检出10例HBeAg阳性标本。23例溶血标本引起HBeAg假阳性。结论:钩状效应和标本溶血是引起HBeAg假阴性和假阳性的重要原因。必要时应加以复查,以减少HBeAg的错检和漏检。     

家畜试验中的一种相关扰动回归模型

本文根据数量遗传学原理,从家畜群体中同一家系内的个体之间存在的遗传相关性出发,建立了扰动项相关的回归模型,运用回归理论讨论一类相关样本的统计推断及其应用。     

不育男性精浆锌与其它精液参数的回归分析

对２２３例男性不育患者的精浆锌进行了测定,并用回归分析的方法研究了精浆锌与其它精液参数的关系．结果显示：精浆锌与精液ｐＨ、精子存活率之间存在显著的线性回归（Ｐ＜０．０１）,而精浆锌与精子密度、前向运动精子百分率、精子畸型率具有显著的曲线回归（Ｐ＜０．０１）．     

Estimation of a Conditional Mean in a Linear Regression Model after a Preliminary Test on Regression Coefficient

Consider the two linear regression models of Y_ij on X_ij, namely Y_ij = β_io + β_ij, X_ij + E_ij = 1, 2,…, n_i, i = 1, 2, where E_ij are assumed to be normally distributed with zero mean and common unknown variance σ². The problem of estimating the conditional mean of Y₁ for a given value of X₁ is considered when it is a priori suspected that β₁₀ = β₂₀ and β₁₁ = β₂₁. The preliminary test estimator is proposed. The exact expressions for the bias and the mean square error of the estimator are derived. The relative efficiency of the new estimator to the usual least square estimator based on the first regression alone is computed and is used to determine the appropriate value of the significance level of the preliminary test β₁₀ = β₂₀ and β₁₁ = β₂₁.     

肾细胞癌淋巴结转移预测指标的Logistic回归分析

目的：研究肾细胞癌淋巴结转移的危险因素,并建立Logistic回归模型。方法：2002年2月-2010年10月我院手术治疗的肾细胞癌163例,对其临床病理资料进行单因素和多因素的Logistic回归分析。结果：淋巴结转移的发生率为20．9％（34／163）。单因素分析显示：肿瘤大小、临床分期、Fuhrman核分级和贫血与肾细胞癌淋巴结转移的风险有关（P〈0．05）;多因素分析显示：肿瘤大小、临床分期和Fuhrman核分级是RCC淋巴结转移独立的风险因素。结论：肾细胞癌淋巴结转移的风险与肿瘤大小、临床分期和Fuhrman核分级有关,Logistic回归模型对于判断预后、指导术后治疗及随访方案的制订具有重要作用。     

Forecast of Sugarcane Yield Using Eye Estimate Alongwith Plant Characters

The present study proposes to include farmers' appraisal (eye estimate) as one of the variables in the models alongwith plant characters. This will reduce the number of plant character in the model, thus reducing the cost on data collection. This has been demonstrated through regression and probability models.     

人工养殖大鲵全长与体重关系的回归分析

大鲵(Andrias davidianus)隶属两栖纲,有尾目,隐鳃鲵科,大鲵属,为我国特有的珍稀濒危两栖动物。目前在我国一些主要历史分布区掀起了一股大鲵养殖高潮,但是对当前养殖效果缺乏成熟的评价手段,因此本研究通过对陕西省汉中市和安康市两大鲵养殖场养殖的1530尾大鲵的体重和全长进行实际测量,利用SPSS分析软件对数据进行了回归分析。结果表明,体重与全长之间存在极显著的相关关系,体重与全长之间的关系主要表现为幂函数关系,其关系式为:Y=0.010X2.867。本研究的顺利完成,为今后评价大鲵的人工养殖效果提供了科学的方法。     

Ames试验假阳性的判断方法

目的:确立一种判断Ames试验假阳性的可靠方法。方法:在Ames平皿掺入试验中发现经受试物处理的TA97和TA1535的菌落数相比溶媒对照组明显增加,为了证实其是因为受试物致突变性导致的回变菌落数增加还是由受试物毒性导致的菌落数增加,对分别经受试物和阳性对照物处理的Ames菌落进行增菌培养后接种于无组氨酸的培养基上,观察比较细菌的生长情况。结果:经受试物处理的菌株不能生长在无组氨酸的培养基中,而经阳性对照物处理的菌株则可以生长在无组氨酸的培养基上,说明经受试物处理的菌株没有发生突变。结论:此方法能可靠地验证菌株是否为突变菌株,由本研究可以发现经受试物处理的菌株没有发生突变,Ames平皿掺入试验中所观察到的菌落数增加是由于受试物毒性造成的假阳性。