多元统计应用.pdfVIP

与小麦室交流 统计学应用 汇报人：陈庭木 连云港市农业科学院 报内容 1. 主成分分析与因子分析的应用 2. 最优分类方法的应用 3. 偏相关与通径分析的应用 4. 生长方程的最优拟合 主成分分析与因子分析的应用 主主成成分分分分析析法法的的基基本本思思想想 主成分分析 （PrincipalComponentAnalysis）是利 用 维的思想，将多个变量转化为少数几个综合变量 （即主成分），其中每个主成分都是原始变量的线性 组合，各主成分之间互不相关，从而这些主成分能够 反映始变量的绝大部分信息，且所含的信息互不重叠。 分分析析相相关关系系数数矩矩阵阵特特征征值值及及相相应应特特征征向向量量 采用Ja obi算法 aa (())AAssDDoouubbllee方方阵阵EEiinnggeennVVaalluueess(())AAssDDoouubbllee特特征征根根,, ‘‘’’ ‘‘vv (())AAssDDoouubbllee特特征征向向量量,, BByyVVaall oouunnttAAssLLoonngg最最大大迭迭代代次次数数,, ’’ ‘‘’’BByyVVaallnn AAssLLoonngg方方阵阵阶阶数数 Sub Cy li Ja obi(a()AsDouble,EingenValues()AsDouble,v()AsDouble,ByVal ountAsLong,ByValnAsLong) Const toleran e= 0.0000000001 Dimmax ount, self ount, i,p,qAsLong Dim Sumsq,avp(),avq()AsDouble ReDim avp(n,n),avq(n,n) Dim su essAsBoolean:su ess=False max ount= ount self ount= 0 Forp = 1Ton 特征向量初始化 For q= 1Ton v(p,q)= 0 Next v(p,p)= 1 Next Do selfcount=selfcount+1 Sumsq=0 Forp= 1Ton主对角线元素平方 Sumsq=Sumsq+a(p,p)*a(p,p) Next MsgBox主对角线元素平方 = Sumsq success=True Forp= 1Ton- 1中心，JACOBI计算 Forq=p+1Ton IfAbs(a(p,q)/Sumsq)toleranceThen success=False jacobirotationa,p,q,v,n EndIf Next Next LoopWhileNot (successOrselfcountmaxcount) IfsuccessThen Forp= 1Ton EingenValues(p)=a(p,p) Next EndIf EndSub 刘刘艳艳博博士士耐耐盐盐试试验验数数据据分分析析 根根据据以以下下特特征征向向量量与与下下页页计计算算程程序序计计算算主主成成分分得得分分 计计算算主主成成分分得得分分方方法法 Function countCompoScore(v()AsDouble,x()AsDouble,ByValnxAs Integer,_ B