《多元回归R案例》PPT课件.ppt

棉花的纤维强力 13级统计学 棉纤维的简介 棉纤维是我国纺织工业的主要原料 主要成份是纤维素 其元素组成为碳44 44 氢6 17 氧49 39 棉纤维的聚合度在6000 11000间 其组成与长度 线密度 成熟度有关 棉纤维长度是指纤维伸直时两端间的距离 是棉纤维的重要物理性质之一 棉纤维的长度主要由棉花品种 生长条件 初加工等因素决定 棉纤维长度与成纱质量和纺纱工艺关系密切 棉纤维长度长 整齐度好 短绒少 则成纱强力高 条干均匀 纱线表面光洁 毛羽少 棉纤维的线密度是指纤维的粗细程度 是棉纤维的重要品质指标之一 它与棉纤维的成熟程度 强力大小密切相关 且还是决定纺纱特数与成纱品质的主要因素之一 纤维较细 则成纱强力高 纱线条好 棉纤维的成熟度是指纤维细胞壁的加厚程度 即棉纤维生长成熟的程度正常成熟的棉纤维 截面粗 强度高 转曲多 弹性好 有丝光 纤维间抱合力大 成纱强力也高 棉纤维的强度是纤维具有纺纱性能和使用价值的必要条件之一 纤维强度高 则成纱强度也高 案例分析 某地区所产原棉的纤维强力Y与纤维的公制支数X1 指单位质量 g 的纤维所具有的长度 m 纤维成熟度X2有关 数据如下 R程序 library foreign w read csv 02211 csv 建立以Y为因变量 X1 X2为自变量做线性回归A lm Y X1 X2 data w summary A 输出结果 Call lm formula Y X1 X2 data w Residuals Min1QMedian3QMax 0 14092 0 05819 0 039050 021810 20788 Coefficients EstimateStd ErrortvaluePr t Intercept 6 595e 001 120e 005 8890 000105 X1 5 106e 049 623e 05 5 3060 000250 X22 155e 013 946e 010 5460 595816Signif codes 0 0 001 0 01 0 05 0 1 1残差标准差 决定系数及回归方程的显著性检验结果 Residualstandarderror 0 1225on11degreesoffreedomMultipleR squared 0 9163 AdjustedR squared 0 9011F statistic 60 24on2and11DF p value 1 186e 06自变量X1的回归系数在0 05的显著水平下是显著的 且该回归方程是显著的 方差分析 anova A AnalysisofVarianceTableResponse YDfSumSqMeanSqFvaluePr F X111 803401 80340120 17912 93e 07 X210 004480 004480 29840 5958Residuals110 165070 01501 Signif codes 0 0 001 0 01 0 05 0 1 1 标准化数据回归分析及散点图w as matrix w stdx scale w 2 4 center TRUE scale TRUE stdx as data frame stdx B lm Y X1 X2 data stdx summary B attach stdx par mfrow c 2 2 plot Y X1 abline lm Y X1 plot Y X2 abline lm Y X2 在实际问题中 不同变量的测量单位往往是不一样的 这就是所谓的量纲 而不同的量纲会引起各变量取值的分散程度和均值的差异较大 所以不能由偏回归系数的大小直接说明自变量对因变量线性影响的大小 回归方程为且方程是显著的 从该散点图可以看出 自变量 的分布趋势大致在一条直线上 且与因变量负相关 自变量 与因变量呈正相关 两个散点图 绘制残差图 标准化残差图和残差 图resid residuals B stdresid rstandard B C predict B par mfrow c 2 2 plot C resid plot C stdresid plot B 2 通过残差 图大概可以看出 以上的样本的标准残差都在 到 之间且图上的点近似的在一条直线附近 所以可认为样本数据来自正态分布的总体 所以知道该拟合的回归模型具有合理性与可靠性 w1 x data frame X1 0 48 X2 0 42 w1 y predict B w1 x interval prediction level 0 95 w1 y 若已知一组标准化后数据X1 0 48 X2 0 42 用拟合的线性模型预测Y的值及置信水平为95 的置信区间