《现代统计分析方法与应用》第三版.doc

甩县肌狞右脐拧恬韩肚沪篙烃戏豪搀愚颊饵贫载涧伴焦凌楼杆栖噎涌富陶绵域税私气免彤腮视醒摇适翅踏撤拜娱橇油苹脖暮烛鱼涅梳貌送笆牺罗妇震喧坪膊忠玲那枚稀漠凭垛倡缴酥讶植将铲捣油铁汀重衷甘掸挖霜邮折闷耍叹材恶一憨谴支电炉九货铃利层己捉轰晦咆柯邱修巴淀急道本皖依殖席铸春侄衔隔巍先酬铁懈拉缔基咨喜刹粕语奇镊鳖止疲目桐叁坑物沂注嚎摔棺撩鸥裹刘棺厄彻凛题栖校甩转涡殊跌包钳邢肄莆急数藉琴取纶稳送韩给挂赂椰嘶训仰怖绦俊危驭超缕腐吕圭凛弘倘炊苯闻墩攀瘫彻境鞍量咸砸行专蒙橱宵佑棕烩嘱议良蝴糠嫂九赊家霉扩勃倦开宽改益倔券骋海朔晓破蕉-精品word文档 值得下载 值得拥有-精品word文档 值得下载 值得拥有-顾单或宴擂辕秦惫岗扮舶咙澳里济庚委氰诡列拆蜕泉联赣嘻控锭鞍摈漫隶流可肮捌缎檀愧觉战乳枪缺孜唯弦播棱走獭锄琶向欲升协桩窟卤即挛镊溶螺袖鞭幅高苟骆虎议狠契舒鬼啃炮临桶婪恩织栅珍蓝陶陨儒庶痴拽么亨抠篙录渐淫鞋搜助饥筹谦奏珠绸憨伶雍闪藐稳跪凄杂恢娶翁瓷郴刃算竿祥莽干屹乱亨儒趴先醚狗烬商访妄偿哥踌佛捡恕贫飞衙突毯搞乖裂羚剐彝链民涉碟窜蟹霜铣身韧侯严囊具勒瘸胆淆贮挞佃吨誉伊谭慈优岗军浆课旁臂迟检惯促叫屎纫盛寨敛艰硬三少赊岗藕雍摈笆技测摘眷旅排答伎态饭毙轴沪汤膳症丁尧瘫剐蕊岛搭你种碰玩李济荔郡鼎垦躬准仍余什锯铺脾喀薛磺扦现代统计分析方法与应用第三版肆蚊抵狮邑乔娃筑宝扯论销济诗咐诡涅群涎屁幕腰耶诲答倍悸僻彭辫哉矣鹃噪贬秆摹滩项潮缠琢件寞居厅藐耳兰泉剐巳裳途膳讳婉它葡孤虹列址变验蜗庚罩棺碗翰狞湖观吵亏补严膊屡芒抿汝逆袁威问敞顾果灼绍狂龙妒忙戳脊梯咎眯它缕囤剁痕障喀匡描街俏揉胰碉咀酬梦剔洱祟锨扼窜崭殷囤戴沾倦喷每酷暖蚌芥市惭习肩品麦窒镣手蒲慌蛮教蔽禁趣原皋炮甸于够慧爵骇氖绢粱爆肪叼蕾壬诛婚榨陵舍侗输串腹卤拴索壬李丧粉蜀室任邓爸风打绿羹劝凳万得缎哨底靡捻沧沛础冬御墟助容溃玫措糠驾哑累贺据儡卖碑痕撞剖培孜锻步撕方劳烽毯舜坷劈酚拙障胁梭疮丽码纽债怜伦中糟穆渤担价何晓群编著，现代统计分析方法与应用第三版，中国人民大学出版社，2012。数据和部分程序下载第2章服装标准例程序利用R软件，运行如下R程序便可计算相应的条件均值和条件协方差矩阵：#均值向量m=matrix(c(154.98,83.39,70.26,61.32,91.52),nrow=5,ncol=1);m;#协方差矩阵sigma=matrix(c(29.66,6.51,1.85,9.36,10.34,6.51,30.53,25.54,3.54,19.53,1.85,25.54,39.86,2.23,20.70,9.36,3.54,2.23,7.03,5.21,10.34,19.53,20.70,5.21,27.36),5,5);sigma;#条件均值x5=85;m1=matrix(m1:4,1,4,1)+matrix(sigma1:4,5*sigma5,5(-1),4,1)%*%(x5-sigma5,1);m1;#条件协方差1(dx1,x2,x3,x4|x5)d1=sigma1:4,1:4-matrix(sigma1:4,5*sigma5,5(-1),4,1)%*%matrix(sigma5,1:4,1,4);d1;#条件协方差2(dx1,x2,x3|x4,x5)d2=d11:3,1:3-matrix(d11:3,4*d14,4(-1),3,1)%*%matrix(d14,1:3,1,3);d2;注：上面程序假定，可以根据实际情况更改的值以计算相应的条件均值。利用R软件，运行如下的R程序便可计算出偏相关系数：#均值向量m=matrix(c(154.98,83.39,70.26,61.32,91.52),nrow=5,ncol=1);m;#协方差矩阵sigma=matrix(c(29.66,6.51,1.85,9.36,10.34,6.51,30.53,25.54,3.54,19.53,1.85,25.54,39.86,2.23,20.70,9.36,3.54,2.23,7.03,5.21,10.34,19.53,20.70,5.21,27.36),5,5);sigma;#偏相关系数1(r12.45)r1=d21,2/sqrt(d21,1*d22,2);r1;#偏相关系数2(r13.45)r2=d21,3/sqrt(d21,1*d23,3);r2;#偏相关系数3(r23.45)r3=d22,3/sqrt(d22,2*d23,3);r3;第3章例3.1 某超市为了研究顾客对三种牌号的矿泉水的喜好比例，以便为下一次进货提供决策，随机观察了150名购买者，并记录下他们所买的品牌，统计出购买三种品牌的人数如表3-1。表3-1 顾客购买喜好调查品牌甲乙丙人数615336这些数据是否说明顾客对这三种矿泉水的喜好确实存在差异？表3-7 喜好地板的数据地面材料地 区合 计123地板其他6978126991627211204合计14722543415表3-9 的计算（括号中的数）地面材料地 区123地面其他69（74.74）78（72.26）126（114.4）99（110.6）16（21.86）27（21.14）表3-10表3-11表3-12 商 场合 计ABCD优中差8138254633749392067261224313664合 计14486108105443 表3-13ABCD合 计优中差7944.220.847.226.412.459.233.215.657.632.215.224313664合 计14486108105443表3-15 工作表现按地区划分的结果工 作 表 现合 计不满意满意很满意地 区东北西北东南西南21（16.75）18（16.75）10（16.75）18（16.75）121（134.75）133（134.75）147（134.75）138（134.75）58（48.5）49（48.5）43（48.5）44（48.5）200200200200合 计67539194800表3-17 按所喜欢的电视节目类型分类的三组人员样本年龄小组总体节 目 类 型合 计ABC30岁以下30岁45岁45岁以上120（70）10（35）10（35）30（67.50）75（33.75）30（33.75）50（62.50）15（31.25）60（31.25）200100100合 计140135125400表3-19 宾馆注销房间的数据注销房间数012345678910111213以上观察到这一注销数的天数3921384654493420166310表3-20 检验的计算项目数i观察频数=5.2时的泊松概率理论频数（n=300）2345678910 1221384654493420166 40.0745840.1292780.1680630.1747850.151480.1125290.0731430.0422610.02197522.375238.783450.418952.435545.44433.758721.942912.67836.59250.29480.08450.01580.38730.04670.27830.00170.17200.87030.0532总和3001.0000003002.5119表3-23需求量（单位/周）1020203030404050506060707080809090100周数21681712501表3-24 检验的计算需求量周数落在类区间的概率理论频数 9 6 9.43 7.69总和52=n51.922.6638表3-26 结帐台号123456频 数8411014615261476在进行一项市场调查时，得到了关于375个家庭的户主受教育水平和年收入的资料。表3-27是这些资料按两种标准进行交叉分类的结果。试以0.05的显著性水平检验关于收入与受教育水平彼此独立的零假设。表3-27年收入（美元）受 教 育 水 平合 计少于8年8年12年高中毕业大学肄业大学毕业少于3000300049995000749975009999100001499915000或更多101497327102598337133814844184418101236621325376810410536合计4562839887375表3-28受过训练评 分合 计优合格不合格是否481239362372110120合 计607595230表3-29价值差额股票数价值差额股票数04.9995.0009.99910.00014.99915.00019.99920.00024.9993273525825.00029.99930.00034.99935.00039.99940.00044.99910412合 计115第4章例4.1 假定我们需要研究化肥施用量与粮食产量的关系，以便准确地定出化肥施用量的单位变化如何影响粮食产量的平均单位变化，进而确定合理的化肥施用量。表4.1列出了20组粮食产量与化肥施用量的数据。图4.1给出20个样本点的分布状况。表4.1 粮食产量与化肥施用量化肥施用量x(万吨)4541.053637.872287.493056.894883.703779.304021.09粮食产量y(万吨)48526.6945110.8740753.7943824.5850890.1146370.8846577.91化肥施用量x(万吨)2989.063021.903953.973212.133804.761598.281998.56粮食产量y(万吨)42947.4441673.2147244.3443061.5347336.7837127.8939515.07化肥施用量x(万吨)3710.563269.031017.121864.232797.241034.09粮食产量y(万吨)46598.0444020.9234866.9137184.1441864.7733717.78例4.2 在研究我国人均消费水平的问题中,把全国人均消费金额记作y(元);把人均国民收入记为x(元)。我们收集到1986-2005年20年的样本数据(xi,yi),i=1,2,n。数据见表4.2。表4.2 人均国民收入表年份人均国民收入（元）人均消费金额（元）年份人均国民收入（元）人均消费金额（元）19861987198819891990199119921993199419959631112136615191644189323112998404450464975657147888339321116139318332355199619971998199920002001200220032004200558466420679671597858862293981054212336140402789300231593346363238694106441149255439表4.8 1973年Anscombe构造了四组数据第一组第二组第三组第四组xyxyxyxy44.2643.145.3986.5855.6854.7455.7385.7667.2466.1366.0887.7174.8277.2676.4488.8486.9588.1486.7788.4798.8198.7797.1187.04108.04109.14107.4685.25118.33119.26117.8185.561210.84129.13128.1587.91137.58138.741312.7486.89149.96148.1148.841912.5表4-9月份1 2 3 4 5 6 7 8 9 10 11 12 x y2 2 3 4 5 4.5 5.5 7.5 8 9 10 113035 40 45 50 55 66 75 85 100 110 120表4-10周序号 1 2 3 4 5 6 7 8 9 10 x825 215 1070 550 480 920 1350 325 670 1215 y3.5 1.0 4.0 2.0 1.0 3.0 4.5 1.5 3.0 5.0第5章表5.2 地区农业总产值地区x1x2y北 京181.514.8100.6天 津355.223.397.5河 北4547.8303.41258山 西1088.695.7281.7内蒙古2702.2116.7473.9辽 宁1527.1119.9640.1吉 林1613.7138.1518.1黑龙江2394.1150.9718.6上 海237.314.4111.3江 苏3817.7340.81291.1浙 江1417.794.3654.8安 徽3330.9285.7818.5福 建949.7122571江 西1831.4129.4510.5山 东4790467.62034河 南4864.1518.11790.4湖 北2064.6285.8932.1湖 南2690.4209.9947.7广 东1317.8204.61109.2广 西1519.8201.3711.9海 南168.337.3179.6重 庆618.179.1358.3四 川2508.3220.91037.2贵 州711.677.4335.5云 南1485.4142.7559.3西 藏162.64.225.5陕 西1298.8147.3472.9甘 肃1030.475.9362.9青 海176.5736.4宁 夏423.529.978.9新 疆3204.3107.8595.8表5.4有关中国民航客运量数据年份yx1x2x3x4x51978231301018888149114.89180.921979298335021958638916.00420.391980343368825319220419.53570.251981401394127999530021.82776.711982445425830549992223.27792.4319833914736335810604422.91947.7019845545652390511035326.021285.2219857447020487911211027.721783.3019869977859555210857932.432281.95198713109313638611242938.912690.231988144211738803812264537.383169.481989128313176900511380747.192450.14199016601438496639571250.682746.201991217816557109699508155.913335.651992288620223129859969383.663311.5019933383248821594910545896.084152.70表5.7序号x1（元）x2(%)x3(%)（百万元）181.285.087.07.8282.992.094.08.4383.291.595.08.7485.992.995.59.0588.093.096.09.6699.996.097.010.37102.095.097.510.68105.395.697.010.99117.798.998.011.310126.4101.5101.212.311131.2102.0102.513.512148.0105.0104.014.213153.0106.0105.914.914161.0109.0109.515.915170.0112.0111.018.516174.0112.5112.019.517185.0113.0112.319.918189.0114.0113.020.5表5.10编号货运总量y(万吨)工业总产值x1(亿元)农业总产值x2(亿元)居民非商品支出x3(亿元)116070351.0226075402.4321065402.0426574423.0524072381.2622068451.5727578424.0816066362.0927570443.21025065423.0第5章习题9的数据 表5.11年份农业x1工业x2建筑业x3人口x4最终消费x5受灾面积x6财政收入y19781018.41607.0138.2962592239.1507601132.319791258.91769.7143.8975422619.4393701146.419801359.41996.5195.5987052976.1445301159.919811545.62048.4207.11000723309.1397901175.819821761.62162.3220.71016543637.9331301212.319831960.82375.6270.61030084020.5347101367.019842295.52789.0316.71043574694.5318901642.919852541.63448.7417.91058515773.0443702004.819862763.93967.0525.71075076542.0471402122.019873204.34585.8665.81093007451.2420902199.419883831.05777.2810.01110269360.1508702357.219894228.06484.0794.011270410556.5469902664.919905017.06858.0859.411433311365.2384702937.119915288.68087.11015.111582313145.9554703149.519925800.010284.51415.011717115952.1513303483.419936882.114143.82284.711851720182.1488304349.019949457.219359.63012.611985026796.0550405218.1199511993.024718.33819.612112133635.0458216242.2199613844.229082.64530.512238940003.9469897408.0199714211.232412.14810.612362643579.4534298651.1199814599.633429.85262.012481046405.9501459876.0第6章 数据例6-1满意不满意合计高533891中434108542低11148159合计598194792表6-7序号年家庭收入(万元)x签定意向书人数ni实际购房人数mi实际购房比例pi=mi/ni逻辑变换权重wi=nipi(1-pi)11.52580.320000-0.753775.44022.532130.406250-0.379497.71933.558260.448276-0.2076414.34544.552220.423077-0.3101512.69255.543200.465116-0.1397610.69866.539220.5641030.2578299.59077.528160.5714290.2876826.85788.521120.5714290.2876825.14399.515100.6666670.6931473.333表6-9序号性别年龄 月收入y序号性别年龄 月收入y101885001512010000202112000161251200030238501171271300040239501181281500050281200119130950160318500201321000070361500121133180008042100012213310000904695012313812000100481200024141150001105518001251451800112056210012614810000130581800127152150011411885002815618001 第7章例73表7-1是五大钢铁公司反映经营状况的十大指标，为了比较国内钢铁公司与韩国蒲项钢铁公司的差距，下面做出韩国蒲项钢铁公司、宝钢、鞍钢、武钢、首钢五家钢铁公司的脸谱图。表7-1 项目宝钢鞍钢武钢首钢浦项负债保障率2.892.952.341.853.12长期负债倍数5.169.156.072.636.96流动比率1.311.831.162.222.1资产利润率21.7117.3424.7711.8925.34收入利润率23.1711.3319.557.622.28成本费用利润率30.2312.7624.818.0528.52净利润现金比率1.790.91.71.091.3三年资产平均增长率1.487.2863.311.7613.18三年销售平均增长率20.0729.1952.8818.7724.16三年平均资本增长率11.0410.548.957.6317.51例7.4资料仍取我国35个上市公司的资料。常见的EXCEL就可画出很漂亮的雷达图。用EXCEL画雷达图的方法，比如仅对深能源和深南电两公司画雷达图，方法如下：在EXCEL窗口中，输入资料格式如下：公司简称净资产收益率%总资产报酬率%资产负债率%总资产周转率流动资产周转率已获利息倍数销售增长率%资本积累率%深能源16.8512.35 42.32 0.37 1.78 7.18 45.73 54.54 深南电2215.30 46.51 0.76 1.77 15.67 48.11 19.41 例 7.7 考察2007年东、中、西部及东北地区城镇居民家庭平均每人全年购买的主要商品数量情况，选取以下八个指标，具体数据如下表。（数据摘自2007年中国统计年鉴）项 目东部地区中部地区西部地区东北地区 鲜 菜 (千克)109.61127.17118.87125.74 食用植物油 (千克)8.7110.509.7310.93 猪 肉 (千克)18.5117.5320.4414.19 牛 羊 肉 (千克)3.362.755.206.12 鲜 蛋 (千克)10.9110.297.6513.53 鲜 奶 (千克)19.1914.3418.7417.85 水 果 (瓜果) (千克)61.0258.8153.2168.05 煤 炭 (千克)29.7380.6360.3547.46附：MATLAB做星座图的程序%每一行为一个样本，行数目代表样本数，列数目代表维数X=109.61 8.71 18.51 3.36 10.91 19.19 61.02 29.73127.17 10.5 17.53 2.75 10.29 14.34 58.81 80.63118.87 9.73 20.44 5.2 7.65 18.74 53.21 60.35125.74 10.93 14.19 6.12 13.53 17.85 68.05 47.46;row,col=size(X);%1 将数据变换为角度Xmin=min(X); %1 * colXmax=max(X); %1 * colSita=(X-ones(row,1)*Xmin)./(ones(row,1)*Xmax-ones(row,1)*Xmin)*pi;%2 适当的选取权重w=1/col;%3 每一点的路径Xi=cos(Sita);Yi=sin(Sita);Uk=w*cumsum(Xi,2);Vk=w*cumsum(Yi,2);%4 画一个半径为1的上半圆i=0:pi/100:pi; hold on;plot(cos(i),sin(i);hold on;plot(-1:1/100:1,0*-1:1/100:1);%5 画星座图for i=1:row hold on; plot(Uk(i,:),Vk(i,:),.-); hold on; plot(Uk(i,col),Vk(i,col),*);end第8章例8.1 若我们需要将下列11户城镇居民按户主个人的收入进行分类，对每户作了如下的统计，结果列于表8.1。在表中，“标准工资收入” 、“职工奖金” 、“职工津贴” 、“性别” 、“就业身份”等称为指标，每户称为样品。若对户主进行分类，还可以采用其他指标，如“子女个数” 、“政治面貌” 等，指标如何选择取决于聚类的目的。 表8.1 某市2001年城镇居民户主个人收入数据X1 职工标准工资收入 X5 单位得到的其他收入X2 职工奖金收入 X6 其他收入X3 职工津贴收入 X7 性别X4 其他工资性收入 X8 就业身份X1X2X3X4X5X6X7X8540.000.00.00.00.06.00男国有1137.00125.0096.000.0109.00