因子分析和时间序列分析

第十二章因子分析因子分析是一种在许多变量中提取出隐藏的具有代表性的共性因子、构造因子模型的统计技术。因子分析模型设m个可能存在相关关系的观测变量z1,z2,......,zm（经过标准化后）含有p个独立的公共因子F1,F2,......,Fp（m斗），观测变量z.含有独特因子Ui（i=1...m），诸U间互不相关，且与Fj（j=1...p）也互不相关，每个z.可由p个公共因子和自身对应的独特因子U线性表出：Z=aF+aF+A+aF+cU1 111 122 1pp11Z=aF+aF+A+aF+cUV2 211 222 2pp22 （模型1）AAAAAAAAAAAAAAAZ=aF+aF+A+aF+cUPmm11 m22 mppmmfZ1:fF1:fc1U1]Z2M=(a.)-IJmxpF2M+C2U2MmF<p>IcmUm简记为Z=A-F+CU （模型2）（mx1） （mxp）（px1） （mx1）A称为因子负荷矩阵（即模型1中各方程的系数aq的矩阵），aq表示第i个变量z.在第j个公共因子Fj上的负荷，简称因子负荷。因子负荷反映了某一变量与某个因子的相关关系。每一个因子也可以表示各观测变量的线性组合：牛WjSW.2Z2+Wj3z3+....+wjmzmwj1:权重或因子得分系数，用于计算因子得分。因子分析步骤：一、 形成问题二、 基于原始数据构造相关矩阵CorrelationMatrixV1V2V3V4V5V6CorrelationV11.000-.053.873-.086-.858.004V2-.0531.000-.155.572.020.640V3.873-.1551.000-.248-.778-.018V4-.086.572-.2481.000-.007.640V5-.858.020-.778-.0071.000-.136V6.004.640-.018.640-.1361.000KMOandBartlett'sTestKaiser-Meyer-OlkinMeasureofSamplingAdequacy..660Bartlett'sTestofApprox.Chi-Square111.314Sphericitydf15Sig..000三、 确定因子分析方法主成分分析法(Principalcomponentsanalysis简称PCA)和主因子分析法(Principalfactoranalysis简称？?'人，也称commonfactoranalysis)o主成分解释了变量的总方差，主因子解释了协方差。当主要目的是要减少变量时，采用主成分分析；当为了寻找对协方差有贡献的潜在因子时，采用主因子分析。常用主成分分析法。四、 提取因子提取因子原则：按照能够解释方差的大小逐序提取因子。所有公因子(与其他变量所共有)能够解释某个变量方差的比例称为公因子方差(communality)，记作h2。公因子方差反映了各个因子对该变量的解释程度。某变量的h2越大，说明这些因子对该变量的解释程度越强，用这些因子来描述该变量就越有效。CommunalitiesInitialExtractionV11.000.926V21.000.723V31.000.894V41.000.739V51.000.878V61.000.790ExtractionMethod:PrincipalComponentAnalysis.因子的特征值揭示了各个因子能够解释总方差的多少，反映了因子的重要程度，可作为提取因子的依据。

TotalVarianceExplainedComponentInitialEigenvaluesTotal%ofVarian匚日Cumulative%12.73145.52045.52022.21B36.96982.48B3.4427.36089.8484.3415,68895,5365.1S33.044^8.50068.521E-021.420100.000Extra匚ticnMethod;PrincipalComponentAnalysis.五、确定因子的数量方法一：依据经验知识确定。方法二：选择特征值大于1的因子。方法三：scree检验，即将各因子的特征值用折线图表示出来，寻找平滑递减的特征值在图的右边停止不前的位置。六、旋转因子若确定两个因子，得出因子/成分矩阵图如下，展示了两个因子与6个变量之间的相关关系，即因子负荷，体现了因子对变量的解释作用。但是，由于公因子是从多个变量中提取出来的，很有可能出现很多变量甚至所有变量在一个因子上的负荷都较大的情形，从而难以将变量归类和解释因子。故此，将因子矩阵进行正交旋转，采用方差最大正交旋转法，使得每个因子只和少数变量之间显著相关，或者每个变量都只和个别因子（最好是一个）强相关。旋转不会影响公因子方差，但每一个因子能够解释的方差会发生变化，从而使得因子负荷发生变化。

ComponentMatrixaComponent12V1.928.253V2-.301.795V3.936.131V4-.342.789V5-.869-.351V6-.177.871ExtractionMethod:PrincipalComponentAnalysis.a.2componentsextracted.RotatedComponentMatri*Component12V1.962-2.66E-02V2-5.72E-02.848V3.934-.146V4-9.83E-02.854V5-.933-8.40E-02V68.337E-02.885ExtractionMethod:PrincipalComponentAnalysis.RotationMethod:VarimaxwithKaiserNormalization.a.Rotationconvergedin3iterations.六、 解释因子基于旋转后的因子矩阵进行解释。变量1、3、5和因子1相关，变量2、4、6和因子2相关。七、 计算因子得分可根据研究需要计算因子得分：即每一个样本在每一个因子上的得分。计算方法：对于每个样本，用其各个标准化变量值乘以相应的因子得分系数Wij，再加总。ComponentScoreCoefficientMatrixComponent12V1.358.011V2-.001.375V3.345-.043V4-.017.377V5-.350-.059V6.052.395ExtractionMethod:PrincipalComponentAnalysis.RotationMethod:VarimaxwithKaiserNormalization.ComponentScores.

ED1ffaw一SPSSDa.tfEii-lorDI4E.4LU硕”w 班/怔!£UllU.«EiII 山M囹二ti&j« □庄iri巧陶1:nx*b£ir ■nunber5■■JVSraci_irac2_ivarwar1136241M彻i-.296792213J目54-11&248-33S2133B77目131.27922•4日45■4625277661106785&12J3&2-I407566eB36目24JEI64Q3-202237753t3目3390B9-92T5T88G47141.314213=i：r|®369934233-101060-H1010282676-1301S51539D21111R43231.10561-J337S31：212231454■1MW.312941313726目13129229-&22S9佃d645361-J23S1358121515132264-132IBS-.918771616646334mi如-33446171753L：334.611E3-.6129718ie7344149635-.286591919243363-1.02367-.67CM020.'I.:，31.4--1011::■：!■：,212132353■1.18766一河04222254SA24@1431-013772-21■■:d4.皿244■4S47-.010081.88491252565214.83420-.33195hJjJXDfltflVictAVariableVIw/IaL _uSTSSFxica-aiM-iirvady 的k(B芒r'*心涅西Lscix1[^]t-34^b.也bactu., jime处. »]1 [赣八、模型评估即评估模型的拟合程度。因子分析的原假设：变量之间的相关关系可以由一个共同的因子来描述，即变量间的相关程度可以转化为变量与因子间的相关程度/可由变量与因子间的相关关系来反映变量间关系。把通过变量与因子的关系反映出的变量间关系称为再生相关系数，初始的变量间相关系数与再生相关系数之间的差异反映了模型的拟合效果。ReproducedCorrelations V1 V2 V3 V4 V5 V6 ReproducedCorrelationV1.926b-7.76E-02.902-.117-.8955.662E-02V2-7.76E-02.723b-.177.730-1.79E-02.746V3.902-.177.894b-.217-.859-5.13E-02V4-.117.730-.217.739b1.999E-02.748V5-.895-1.79E-02-.8591.999E-02.878b-.152V65.662E-02_ 746-5.13E-02_ 748 -152_ 790bResiduala V12.440E-02-2.92E-023.115E-023.770E-02-5.25E-02V22.440E-022.224E-02-.1583.763E-02-.105V3-2.92E-022.224E-02-3.13E-028.138E-023.327E-02V43.115E-02-.158-3.13E-02-2.66E-02-.107V53.770E-023.763E-028.138E-02-2.66E-021.574E-02 V6 -5.25E-02— -.1053.327E-02— -.1071.574E-02—ExtractionMethod:PrincipalComponentAnalysis.Residualsarecomputedbetweenobservedandreproducedcorrelations.Thereare5(33.0%)nonredundantresidualswithabsolutevaluesgreaterthan0.05.Reproducedcommunalities

第十三章时间序列分析时间序列分析是对于按时间顺序观测的系列数据进行分析的方法。例1:表1为中国海关出口商品总额的有关数据（亿美元），计算环比、定基、平均发展速度和增长速度。年份出口总额逐期增长量环比发展速度环比增长速度定基发展速度定基增长速度1990620.91———100.00—1991718.4397.52115.7115.71115.7115.711992849.40130.97118.2318.23136.8036.801993917.4468.04108.018.01147.7647.7619941210.06292.62131.9031.90194.8894.8819951487.80277.74122.9522.95239.62139.6219961510.6622.86101.541.54243.30143.3019971827.92317.26121.0021.00294.39194.3919981837.099.17100.500.50295.87195.8719991949.31112.22106.116.11313.94213.9420002492.03542.72127.8427.84401.35301.3520012661