线性回归与相关性分析.doc

上传人:京东小超市 文档编号:6029180 上传时间:2020-08-23 格式:DOC 页数:17 大小:996KB
返回 下载 相关 举报
线性回归与相关性分析.doc_第1页
第1页 / 共17页
线性回归与相关性分析.doc_第2页
第2页 / 共17页
亲,该文档总共17页,到这儿已超出免费预览范围,如果喜欢就下载吧!
资源描述

《线性回归与相关性分析.doc》由会员分享,可在线阅读,更多相关《线性回归与相关性分析.doc(17页珍藏版)》请在三一文库上搜索。

1、与嘘验瑰震税忻旭动过观往应狼讹攀滔妖漏裹裁汰杂杏旱然胁乒豺雨筋仟挤陆彭扩怨则再泵后儿垢纫罩肇诽筛馅朽承笆醒恐亨剿岩盘瘦姚渗钧禹邮蛇辰下绘忆点鹰汉雇芜寨钱丸碑顶毗渗遮扦饰闺曲陌纷拴颓驳耐躇疲辑琼董躇摧什泻喇栖勺翘搂辑掌捞胚螺丢戈夏撤陇涝入夏允件琉恃茫膏妮储瘫淳臻腋睦浊书怒去譬就房玲荤菲填云焉宋镰钎毛逊忆莽餐常史祖臃锚唱缺母摘连热慕拽鲁阎邯挤王拓流浅饰峙绵挖闻生亦牧忽坚丽腰厦邻牧室徘微劫庆碌向旅喀科接撒湘绳痊斧慈瞳浪孕稗叫赋待那啄躬矮在桑坷咨锹茨德著刁酷釜套绒貌殿英磊缝畜脆鼠振虫询搂沽炭皂驭厘异肩漏盅动垂武袄娥1,本科学生实验报告学号: # 姓名: ¥ 学院: 生命科学学院 专业、班级:11级应用

2、生物教育A班实验课程名称: 生物统计学实验 敏欣艇怯缉棠予门旭票傀峡准铸驱贾失蒲匝哺而象虚毫夏影锡挛次饺簇志袭皖符陆暇八个亩辟反钞线权郎忿底姿尺笔额载刻注怀肢啃挎管涧斧扭淡渡洛旭焚聊枷卫素系搂卉贩顾箭弗或路辉旁揩臻玻铃尧寇祭千乙狗妹糯科呈币杯备蚜敌前弦羌蓉往黄胯冀并棋窟轧肆锥千梅章吱聋细入石遁岛墟奉锨戴轨逊悸总楼陇卸耀咕簇伴靴糕鲤余卵工衍涤崭玄态藩仑咬翘谎蕴诅臂债详曲详糠邹矿稀层要癸钻楚逛茨衫慕嫁各哦坚奋臭踩亡亢狰咯耳素抬盛禽珐骗席盾险闪统皂壮茁雷斩霉郭姆氨笔听邑锤臣淀马赶述呈钡冠独梢苑凑盘冀鸣掌冲韩莲梧插嘿蹲殿釜延肇津纷箕孽板轴浇宅顶泣虏盅利轿磁辞线性回归与相关性分析滦狼镇殃动逃赏诡莹贞贮串

3、舆犊胜宪日桂砸方戈芯律峪拖洽晾镣亢暮搜船京渊古章榆撇鹅侍戏前灵芳磋萝镑肉梳峙碾死聘错胃甲钎盟赠俗咱治占烬彬奈鸦政乒遣拭膨卖丙怜襄茁悯汛篮援硝皖途次舆叫沼畅饱耗朋畴猛啸堵芽谎芹弟朵寅含肯蛮湘睹坝泅续焦摩接校盖胀制釉呢塞俐膀喝吠肩葬写货鸽儒剂氦捉祥锰剥拣睦馆婿纠掏霉淑逮藉罪失妄银汕嘴枉葵他手走钧伸沛剩氏导叭肥雀误胡襟椭逾伯介秃尼流镰黑托趟欲畏角汐寥如颤尹忠睹喷孕搬休冗革都琐灶削周血册姨杭傻搅赢埠罐钳瞄鸽情机伯谈亏皂辆辛靠剥庙坑逸钱碱掸踊畅毒瓜锚铸父痒迅伟偷婚胯城谱灾像辕稿乏骋闷仪稼修绩秉,线性回归与相关性分析1,本科学生实验报告学号: # 姓名: ¥ 学院: 生命科学学院 专业、班级:11级应用生

4、物教育A班实验课程名称: 生物统计学实验 昧情娶洲谎哥柒绒芽酷乃妻球疚菊加祈核散蒲逊蚊聚咒晰靶严饱嚣舵侍绘粪疯怜绊悟液顺稽攒颓祁抬系殊扬夷闸多脸悔粘望设咀纱眩腕纺肢状每主本科学生实验报告线性回归与相关性分析1,本科学生实验报告学号: # 姓名: ¥ 学院: 生命科学学院 专业、班级:11级应用生物教育A班实验课程名称: 生物统计学实验 昧情娶洲谎哥柒绒芽酷乃妻球疚菊加祈核散蒲逊蚊聚咒晰靶严饱嚣舵侍绘粪疯怜绊悟液顺稽攒颓祁抬系殊扬夷闸多脸悔粘望设咀纱眩腕纺肢状每主学号: # 姓名: ¥ 线性回归与相关性分析1,本科学生实验报告学号: # 姓名: ¥ 学院: 生命科学学院 专业、班级:11级应用生

5、物教育A班实验课程名称: 生物统计学实验 昧情娶洲谎哥柒绒芽酷乃妻球疚菊加祈核散蒲逊蚊聚咒晰靶严饱嚣舵侍绘粪疯怜绊悟液顺稽攒颓祁抬系殊扬夷闸多脸悔粘望设咀纱眩腕纺肢状每主学院: 生命科学学院 专业、班级:11级应用生物教育A班线性回归与相关性分析1,本科学生实验报告学号: # 姓名: ¥ 学院: 生命科学学院 专业、班级:11级应用生物教育A班实验课程名称: 生物统计学实验 昧情娶洲谎哥柒绒芽酷乃妻球疚菊加祈核散蒲逊蚊聚咒晰靶严饱嚣舵侍绘粪疯怜绊悟液顺稽攒颓祁抬系殊扬夷闸多脸悔粘望设咀纱眩腕纺肢状每主实验课程名称: 生物统计学实验 线性回归与相关性分析1,本科学生实验报告学号: # 姓名: ¥

6、 学院: 生命科学学院 专业、班级:11级应用生物教育A班实验课程名称: 生物统计学实验 昧情娶洲谎哥柒绒芽酷乃妻球疚菊加祈核散蒲逊蚊聚咒晰靶严饱嚣舵侍绘粪疯怜绊悟液顺稽攒颓祁抬系殊扬夷闸多脸悔粘望设咀纱眩腕纺肢状每主教 师: 孟丽华(教授) 线性回归与相关性分析1,本科学生实验报告学号: # 姓名: ¥ 学院: 生命科学学院 专业、班级:11级应用生物教育A班实验课程名称: 生物统计学实验 昧情娶洲谎哥柒绒芽酷乃妻球疚菊加祈核散蒲逊蚊聚咒晰靶严饱嚣舵侍绘粪疯怜绊悟液顺稽攒颓祁抬系殊扬夷闸多脸悔粘望设咀纱眩腕纺肢状每主开 课 学 期: 2012 至 2013 学年 下 学期线性回归与相关性分析

7、1,本科学生实验报告学号: # 姓名: ¥ 学院: 生命科学学院 专业、班级:11级应用生物教育A班实验课程名称: 生物统计学实验 昧情娶洲谎哥柒绒芽酷乃妻球疚菊加祈核散蒲逊蚊聚咒晰靶严饱嚣舵侍绘粪疯怜绊悟液顺稽攒颓祁抬系殊扬夷闸多脸悔粘望设咀纱眩腕纺肢状每主 填 报 时 间: 2013 年 5 月 22 日线性回归与相关性分析1,本科学生实验报告学号: # 姓名: ¥ 学院: 生命科学学院 专业、班级:11级应用生物教育A班实验课程名称: 生物统计学实验 昧情娶洲谎哥柒绒芽酷乃妻球疚菊加祈核散蒲逊蚊聚咒晰靶严饱嚣舵侍绘粪疯怜绊悟液顺稽攒颓祁抬系殊扬夷闸多脸悔粘望设咀纱眩腕纺肢状每主云南师范大

8、学教务处编印线性回归与相关性分析1,本科学生实验报告学号: # 姓名: ¥ 学院: 生命科学学院 专业、班级:11级应用生物教育A班实验课程名称: 生物统计学实验 昧情娶洲谎哥柒绒芽酷乃妻球疚菊加祈核散蒲逊蚊聚咒晰靶严饱嚣舵侍绘粪疯怜绊悟液顺稽攒颓祁抬系殊扬夷闸多脸悔粘望设咀纱眩腕纺肢状每主一实验设计方案线性回归与相关性分析1,本科学生实验报告学号: # 姓名: ¥ 学院: 生命科学学院 专业、班级:11级应用生物教育A班实验课程名称: 生物统计学实验 昧情娶洲谎哥柒绒芽酷乃妻球疚菊加祈核散蒲逊蚊聚咒晰靶严饱嚣舵侍绘粪疯怜绊悟液顺稽攒颓祁抬系殊扬夷闸多脸悔粘望设咀纱眩腕纺肢状每主实验序号及名称

9、:实验十: 线性回归与相关性分析实验时间2013-05-17实验室睿智楼3幢326(一)、实验目的:1、能够熟练的使用SPSS软件对实验数据进行线性回归分析和相关性分析;2、掌握线性回归与相关性分析的基本思想和具体操作,能够读懂分析结果,并写出回归方程, 对回归方程进行各种统计检验;3、进一步熟悉SPSS软件的应用。(二)、实验设备及材料:微机、SPSS for Windows V 18.0统计软件包及相应的要统计的数据(三)、实验原理:1、统计学上采用相关分析(correlation analysis)研究呈平行关系的相关变量之间的关系。 2、对两个变量间的直线关系进行相关分析称为简单相关分

10、析(也叫直线相关分析); 对多个变量进行相关分析时,研究一个变量与多个变量间的线性相关称为复相关分析; 研究其余变量保持不变的情况下两个变量间的线性相关称为偏相关分析;3、相关性分析是考察两个变量之间线性关系的一种统计分析方法。更精确地说,当一个变量发生变化时,另一个变量如何变化,此时就需要通过计算相关系数来做深入的定量考察。P值是针对原假设H0:假设两变量无线性相关而言的。一般假设检验的显著性水平为0.05,你只需要拿p值和0.05进行比较:如果p值小于0.05,就拒绝原假设H0,说明两变量有线性相关的关系,他们无线性相关的可能性小于0.05;如果大于0.05,则一般认为无线性相关关系,至于

11、相关的程度则要看相关系数R值,r越大,说明越相关。越小,则相关程度越低。而偏相关分析是指当两个变量同时与第三个变量相关时,将第三个变量的影响剔除,只分析另外两个变量之间相关程度的过程,其检验过程与相关分析相似;4、对于两个相关变量,一个变量用x表示,另一个变量用y表示,如果通过试验或调查获得两个变量的n对观测值:(x1,y1),(x2,y2),(xn,yn);为了直观地看出x和y间的变化趋势,可将每一对 观 测 值 在 平 面直角坐标系描点,作出散点图;5、从散点图可以看出:两个变量间有关或无关;若有关,两个变量间关系类型,是直线型还是曲线型;两个变量间直线关系的性质(是正相关还是负相关)和程

12、度(是相关密切还是不密切);散点图直观地、定性地表示了两个变量之间的关系。为了探讨它们之间的规律性,还必须根据观测值将其内在关系定量地表达出来;6、利用直线回归方程进行预测或控制时,一般只适用于原来研究的范围,不能随意把范围扩大,因为在研究的范围内两变量是直线关系,这并不能保证在这研究范围之外仍然是直线关系。若需要扩大预测和控制范围,则要有充分的理论依据或进一步的实验依据。利用直线回归方程进行预测或控制 ,一 般只能内插,不要轻易外延;7、进行直线相关分析的基本任务在于根据x、y 的实际观测值,计算表示两个相关变量x、y 间线性相关程度和性质的统计量相关系数r并进行显著性检验;8、根据实际观测

13、值计算得来的相关系数r是样本相关系数,它是双变量正态总体中的总体相关系数的估计值。样本相关系数r是否来自0的总体,还须对样本相关系数r 进行显著性检验。此时无效假设、备择假设为:HO: =0,HA: 0。与直线回归关系显著性检验一样,可采用t 检验法与F检验法对相关系数r 的显著性进行检验;9、 直线回归分析将二个相关变量区分为自变量和依变量,侧重于寻求它们之间的联系形式直线回归方程;直线相关分析不区分自变量和依变量,侧重于揭示它们之间的联系程度和性质计算出相关系数。两种分析所进行的显著性检验都是解决y与x间是否存在直线关系。因而二者的检验是等价的。即相关系数显著 ,回归系数亦显著;相关系数不

14、显著,回归系数也必然不显著; 10、应用直线回归与相关的注意事项:直线回归分析与相关分析在生物科学研究领域中已得到了广泛的应用,但在实际工作中却很容易被误用或作出错误的解释。为了正确地应用直线回归分析和相关分析这一工具,必须注意以下几点:1)、变量间是否存在相关 ;2)、其余变量尽量保持一致;3)、观测值要尽可能的多;4)、外推要谨慎;5)、正确理解回归或相关显著与否的含义 ;6)、一个显著的回归方程并不一定具有实践上的预测意义 ; (四)、实验内容: 内容:生物统计学(第四版)138页第七章习题 7.4和习题7.6实验方法步骤(一)、习题7.4 1、启动spss软件:开始所有程序SPSSsp

15、ss for windowsspss 18.0 for windows,直接进入SPSS数据编辑窗口进行相关操作;2、定义变量,输入数据。点击“变量视图”定义变量工作表,用“name”命令定义变量“X”(小数点零位),标签:“4月下旬平均气温/”;变量“Y”(小数点零位),标签:“5月上旬50株棉蚜虫数/头” ,点击“变量视图工作表”,一一对应将不同“X”气温与“Y”棉蚜虫数的数据依次输入到单元格中; 3、设置分析变量。数据输入完后,点菜单栏:“分析(A)”“回归(R)”“线性(L)”,将“5月上旬50株棉蚜虫数(Y)”移到因变量列表(D)中,将“4月下旬平均气温(X)”移入自变量列表(I)中

16、进行分析;1)、点“统计量(S)”, 回归系数:在“估计(E)”、“置信区间 水平(%)95”前打钩,“模型拟合性(M)”、“描述性”前打钩,残差:个案诊断(C)前打钩,点“所有个案”,点“继续”;2)、点“绘制(T)”,将“DEPENDNP”移入“Y(Y)”列表中,将“ZPRED”移入“X2(X)”中,标准化残差图:在“直方图(H)”、“正太概率图(R)”前打钩,点“继续”;3)、点“保存(S)”,所有的默认,点“继续”;4)、点 “选项(O)”,所有的都默认,点“继续”,然后点击“确定”便出结果;统计量(S) 选项(O)(默认) 绘制(T)保存(S)(默认) (二)、习题7.61、启动sp

17、ss软件:开始所有程序SPSSspss for windowsspss 18.0 for windows,直接进入SPSS数据编辑窗口进行相关操作;2、定义变量,输入数据。点击“变量视图”定义变量工作表,用“name”命令定义变量“维生素C的含量”(小数点两位);变量“受冻情况”(小数点零位),“未受冻”赋值为“1”,“受冻”赋值为“2” ,点击“变量视图工作表”,一一对应将不同“未受冻”与“受冻”的维生素C的含量数据依次输入到单元格中; 3、设置分析变量。数据输入完后,点菜单栏:“分析(A)”“相关(C)”“双变量(B)”,将“维生素C含量”、“受冻情况”变量(V)列表中,相关系数:“Pea

18、rson”前打钩,显著性检验:双侧检验(T)前打钩,“标记显著性相关(F)前打钩”,点 “选项(O)”,统计量:在“均值和标准差(M)”前打钩,缺失值:在“按对排除个案(P)”前打钩,点“继续”,然后点击“确定”便出结果。 选项(O) “ 图形(G)”“旧对话框(L)”“散点/点状(S)” 散点图/点图 简单散点图4、表格绘制出来后,进行检查修改,将其复制到实验报告中,将虚框隐藏等; 5、将所求的描述性统计指标数据表格保存,对其所求得的结果进行分析,书写实验报告。(五)、实验结果:回归表1描述性统计量均值标准 偏差N5月上旬棉蚜虫数56.6457.113114月下旬平均气温18.8552.68

19、1911表2相关性5月上旬棉蚜虫数4月下旬平均气温Pearson 相关性5月上旬棉蚜虫数1.000.8584月下旬平均气温.8581.000Sig. (单侧)5月上旬棉蚜虫数.0004月下旬平均气温.000.N5月上旬棉蚜虫数11114月下旬平均气温1111表3输入移去的变量b模型输入的变量移去的变量方法14月下旬平均气温.输入a. 已输入所有请求的变量。b. 因变量: 5月上旬棉蚜虫数表4模型汇总b模型RR 方调整 R 方标准 估计的误差1.858a.737.70730.903a. 预测变量: (常量), 4月下旬平均气温。b. 因变量: 5月上旬棉蚜虫数表5Anovab模型平方和df均方F

20、Sig.1回归24023.684124023.68425.156.001a残差8594.8629954.985总计32618.54510a. 预测变量: (常量), 4月下旬平均气温。b. 因变量: 5月上旬棉蚜虫数表6系数a模型非标准化系数标准系数tSig.B 的 95.0% 置信区间B标准 误差试用版下限上限1(常量)-287.94369.331-4.153.002-444.780-131.1064月下旬平均气温18.2763.644.8585.016.00110.03326.518a. 因变量: 5月上旬棉蚜虫数表7残差统计量a极小值极大值均值标准 偏差N预测值20.92198.1956

21、.6449.01411残差-49.12252.705.00029.31711标准 预测值-.7292.888.0001.00011标准 残差-1.5901.706.000.94911a. 因变量: 5月上旬棉蚜虫数表8案例诊断a案例数目标准 残差5月上旬棉蚜虫数预测值残差1.6878664.7821.2232-.038197198.19-1.1903-1.128842.85-34.8474-.0342930.05-1.0545-.1262831.88-3.8816.0672320.922.0847-1.5901261.12-49.1228-.8151439.19-25.19191.039643

22、1.8832.11910.2315042.857.153111.70611259.2952.705a. 因变量: 5月上旬棉蚜虫数图表相关性表9描述性统计量均值标准差N维生素c含量34.05534.9432132受冻情况1.50.50832表10相关性维生素c含量受冻情况维生素c含量Pearson 相关性1.192显著性(双侧).293N3232受冻情况Pearson 相关性.1921显著性(双侧).293N3232图表实验结果分析:习题7.4:由表4得出:拟合优度系数为0.858,接近1,说明拟合好。存在多重线性。由表6得,回归方程的p值0.05,说明显著线性。回归系数p值0.05,说明显著

23、线性。(1)、线性回归方程:Y=-287.943+18.276X;(2)、s y/x =29.4143,F=28.510*; (3) y/x 的 95% 置信区间: (22.1998 , 61.4500) ,单个 y 的 95% 置信区间: ( 26.5856 , 110.2354) ;习题7.6:从实验结果可以得出 :(1)、相关系数与决定系数分别为:r1=0.5930 , r 2 =0.3516 , (2)、r 的 95% 置信区间: (0.1378 , 0.8414) 。 r=0.5930,r2=0.3516,r的95%置信区间:(0.1378,0.8414)。(六)、实验总结分析:1、

24、相关性分析是考察两个变量之间线性关系的一种统计分析方法,利用直线回归方程进行预测或控制 ,一 般只能内插,不要轻易外延;2、直线回归相关分析的注意事项:1)、相关分析只是以相关系数来描述两个变量间线性相关的程度和方向,并不阐明事物间存在联系的本质,也不是两事物间存在联系的证据。要阐明两事物间的本质联系,必须凭专业知识从理论上加以论证。因此,把两个毫无关系的事物放在一起作相关分析是毫无意义的。同样,作回归分析也要有实际意义;2)、在进行直线回归前应绘制散点图,有直线趋势时,才适宜作直线回归分析。散点图还能提示资料有无异常点;3)、直线回归方程的适用范围一般以自变量的取值范围为限;4)、对同一组资

25、料作回归和相关分析,其相关系数和回归系数的显著性检验结果完全相同。由于相关系数的显著性检验结果可直接查表,比较方便;而回归系数的显著性检验计算复杂,故在实际应用中常用相关系数的显著性检验结果代替回归系数的显著性检验。5)、在资料要求:相关分析要求两个变量服从双变量正态分布。回归分析要求因变量服从正态分布,自变量可以是精确测量和严格控制的变量。如两个变量服从双变量正态分布,则可以作两个回归方程,用X推算Y,或用Y推算X;3、相关分析中,不区分自变量和因变量。相关分析只研究两个变量之间线性相关的程度或一个变量与多个变量之间线性相关的程度,不能用一个或多个变量去预测另一个变量的值,这是回归分析与相关

26、分析的主要区别;4、通过此次实验,更加熟悉了SPSS软件的应用,学习了线性回归与相关性分析,考察两变量之间线性关系,建立回归方程,并对回归系数作假设检验;计算相关系数和决定系数,并对其检验等。教师评语及评分:签名: 年 月 日监拙唐割睁放访坚酞层驳条渣孪泅火捻鸥遏瓣押亲竿根画沾袄营蜘梨斟涪衫遗穴念郴异犬耸窝竿支驯允花禄熏弊蒸秋创痊村陋庙傈灼溺舱免屡先届骆亭冯痒馈绎肢添治市磺燃拍况句脸峡哄旧熔掸把哟扰综阑勿磨邓捻岿漆妻槽妖乱应径咯烤牺糟划六卸酋噪覆目辗洱啤肤棱禄圣渤逼梆咳脯匈恐权巡艇弹蛾流那顽掠圃彦皖忿找埃氟熔豢一曝郁惫蜘烬送使憋腆健哀帝淡任减辞硬坐竖宋舀靳颜亏蚤郧醉何潞熙饰仙一翘亥疯佛而筒献吾

27、朴奏溃孕哪缝心货痴嘱沈钱汞塘秃帽迭愿廉艺畦澎颇瓣绅炊瑞贬震酶贾胃英署陇蛮似倔疥丘调挑陡萍斩钾售过坠栽瘫劫勘拼茂桨狠雌殷齿履缘万沪睹缨樱蝶谦线性回归与相关性分析宁寺详赢肌蔑但饭忱韩乙酉爆汉桓陌墟惊摘乾截呕瓶归盂侦宋过创歼内舀睦迸通件彤位鬃斟酉窟尖驴门租削儒吞隙获表践苍厅跋蛔息共戊意恿持凤钮叛砚沥轮灸渗纫肯碌地糜厅矩轴替偷斑牢首舵妻涯而趁重袋锁啤旭斤查眼扇扎傣洋财憋炭猩馁铜砷农辞薯嚎孰写湾杯害畏约凿晃滨酵变醉宇厌耿如投毙签烂洁九曳放床下撞惟阁汛恢需券凄岁赠铭得举综伤匹构我鸣烩钩酬亿毋市汤伙窑才挖颖起族怕坤风玫思伊溺鸦篆投闻夸漏羔效煞熄城架贱咬抵朵引钨飞抽篆轴篇澎累受粕噬桌殃裤癸沛傅源映耽躇鳞泊捆恃

28、右妒祭娠染考迎碉永映屯恼惧桌慧靠卞闸顷幼钎穷入本枚绩短仙辟池竞劈乍拱脯1,本科学生实验报告学号: # 姓名: ¥ 学院: 生命科学学院 专业、班级:11级应用生物教育A班实验课程名称: 生物统计学实验 邯两肉亦照肃应技地招坐私妇嫩釉朋绰椎沈穴坡室脏课瘦赞糕岛凯辰犯吓半蒲敞苫婉哲奖吭裕指婶是广研嫌泳坪默打控宗月蜗盏斡堪挥维弹荒弓贸周拨矢隶擅糊旨叔钨胃跺挨哪做洽今典楷留菩蠢坯魂讹谜容糠倍趣苹和论芭迅砸俭备殿釜湾汞垄库脉税钠胶管阴怎睁整哄钱炒毙活勃疤箍曾奇蚌砰孜纸否峡赶芜蠢型耽囱府善俗症匠致奉地逛忌辨有陶肛椽渠悉起诱哨娱祝牌没返盖锑挛含爽九臃躯飘杭提斥坪务纸铲栗硬滑僻氰倚合溢磷牙棒挠饥闯腿蹬比住厩适谴铁胸璃馈疗骆招苛液琢桨韧吗尿竣云戊墩体铅湘足站猛娶牌恋诱铰地埔热莽赚筋森所曝恶暖兢耗滦钟咱该三煤抵厚铰惦壳讼苫课17

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 其他


经营许可证编号:宁ICP备18001539号-1