1 / 6
文档名称:

用R软件进行一元线性回归实验报告 (2).doc

格式:doc   大小:34KB   页数:6页
下载后只包含 1 个 DOC 格式的文档,没有任何的图纸或源代码,查看文件列表

如果您已付费下载过本站文档,您可以点这里二次下载

分享

预览

用R软件进行一元线性回归实验报告 (2).doc

上传人:莫比乌斯 2023/3/13 文件大小:34 KB

下载得到文件列表

用R软件进行一元线性回归实验报告 (2).doc

文档介绍

文档介绍:该【用R软件进行一元线性回归实验报告 (2) 】是由【莫比乌斯】上传分享,文档一共【6】页,该文档可以免费在线阅读,需要了解更多关于【用R软件进行一元线性回归实验报告 (2) 】的内容,可以使用淘豆网的站内搜索功能,选择自己适合的文档,以下文字是截取该文章内的部分文字,如需要获得完整电子版,请下载此文档到您的设备,方便您编辑和打印。数理统计上机报告
上机实验题目:用R软件进行一元线性回归
上机实验目的:
进一步理解假设实验的基本思想,学会使用实验检验和进行统计推断。
学会利用R软件进行假设实验的方法。
一元线性回归基本理论、方法:
基本理论:假设预测目标因变量为Y,影响它变化的一个自变量为X,因变量随自变量的增(减)方向的变化。一元线性回归分析就是要依据一定数量的观察样本(Xi,Yi),i=1,2…,n,找出回归直线方程Y=a+b*X  
方法:对应于每一个Xi,根据回归直线方程可以计算出一个因变量估计值Yi。回归方程估计值Yi与实际观察值Yj之间的误差记作e-i=Yi-Yi。显然,n个误差的总和越小,说明回归拟合的直线越能反映两变量间的平均变化线性关系。据此,回归分析要使拟合所得直线的平均平方离差达到最小,据此,回归分析要使拟合所得直线的平均平方离差达到最小,简称最小二乘法将求出的a和b代入式(1)就得到回归直线Yi=a+bXi。那么,只要给定Xi值,就可以用作因变量Yi的预测值。
(一)
实验实例和数据资料:
有甲、乙两个实验员,对同一实验的同一指标进行测定,两人测定的结果如下:
实验号
1
2
3
4
5
6
7
8


















试问:甲、乙两人的测定有无显著差异?取显著水平α=.
上机实验步骤:
设置假设:H0:u1-u-2=0:H1:u1-u-2<0
确定自由度为n1+n2-2=14;显著性水平a=
计算样本均值样本标准差和合并方差统计量的观测值
alpha<-;
n1<-8;
n2<-8;
x<-c(,,,,,,,);
y<-c(,,,,,,,);
var1<-var(x);
xbar<-mean(x);
var2<-var(y);
ybar<-mean(y);
Sw2<-((n1-1)*var1+(n2-1)*var2)/(n1+n2-2)
t<-(xbar-ybar)/(sqrt(Sw2)*sqrt(1/n1+1/n2));
tvalue<-qt(alpha,n1+n2-2);
计算临界值:tvalue<-qt(alpha,n1+n2-2)
比较临界值和统计量的观测值,并作出统计推断
实例计算结果及分析:
alpha<-;
>n1<-8;
>n2<-8;
>x<-c(,,,,,,,);
>y<-c(,,,,,,,);
>var1<-var(x);
>xbar<-mean(x);
>var2<-var(y);
>ybar<-mean(y);
>Sw2<-((n1-1)*var1+(n2-1)*var2)/(n1+n2-2)
>t<-(xbar-ybar)/(sqrt(Sw2)*sqrt(1/n1+1/n2));
>var1
[1]
>xbar
[1]
>var2
[1]
>ybar
[1]
Sw2
[1]
>t
[1]-
tvalue
[1]-
分析:t=->tvalue=-,所以接受假设H1即甲乙两人的测定无显著性差异。
(二)
实验实例和数据资料:
%,且其服从正态分布,现对一批该批号的玻璃纸测得100个数据如下:
(x%横向延伸率)















频数
7
8
11
9
9
12
17
14
5
3
2
0
2
0
1
上机实验步骤:
(1)设置假设:H0:u=65,H1:u<65.
(2)确定自由度为n=100-1=99;显著性水平a=
(3)输入数据x<-c(,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,)
(4)用R软件计算临界值
(5)比较临界值和统计量的观测值,并作出推断
实例计算结果及分析:
计算过程如下:
alpha<-;
n<-100;
x<-c(,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,)
sd1<-sd(x);
xbar<-mean(x);
t<-(xbar-)/(sd1/sqrt(n));
tvalue<-qt(alpha,n-1);
sd1
[1]
xbar
[1]
t
[1]-
tvalue
[1]-
分析推断:因为t=-<tvalue=-。即该批玻璃纸的横向延伸率不符合要求
(三)
实验实例和数据资料:
为了检验一种杂交作物的两种新处理方案,在同一地区随机选择16块地段在各实验地段,按两种方案处理作物,这8块地段的单位面积产量(单位:公斤)是:
一号方案产量:86875693849375798178799068658790;
二号方案产量:80795891778274665859647876808255;
假设两种方案的产量都服从正态分布,分别为N(u1,a^2),N(u2,a^2),a^2未知,求均值差u1-u2的置信区间;
实例计算结果及分析:
利用R软件求解过程如下:
>alpha<-;
>x<-c(86,87,56,93,84,93,75,79,81,78,79,90,68,65,87,90);
>y<-c(80,79,58,91,77,82,74,66,58,59,64,78,76,80,82,55);
>n1<-length(x);
>n2<-length(y);
>xbar=mean(x);
>ybar=mean(y);
>sw<-sqrt((n1-1)*sqrt(var(x))+(n2-1)*sqrt(var(y)))/(n1+n2-2);
>q<-qt(1-alpha/2,(n1+n2-2));
>left<-xbar-ybar-q*sw*sqrt(1/n1+1/n2);
>right<-xbar-ybar+q*sw*sqrt(1/n1+1/n2);
>n1
[1]16
>n2
[1]16
>left
[1]
>right
[1]
所以置信区间【,,】

最近更新

养鸡大棚薄膜安装工艺 23页

粤教粤科版小学科学四年级下册全册练习(附答案.. 48页

给初中毕业孩子的一封信(精选9篇) 12页

绿色家园果蔬超市商业计划书 26页

考研之教育心理学摸拟题及答案解析55 24页

自考护理教育导论复习资料新编 44页

学校述职报告2022年完整版汇编5篇 25页

学校运动会开幕式领导讲话稿(全文共1811字) 4页

二色补血草的化学成分及质量标准研究的中期报.. 1页

先声药业抗体工艺 29页

2024年高中毕业生登记表自我鉴定00字6篇 8页

辽海版六年级下册劳动教育《主题二、我学会包.. 5页

部编人教版四年级上册语文第六单元试卷(含答案.. 12页

部编版四年级语文上册 第八单元语文素养评估.. 7页

二元氟铝酸盐微结构的量子化学从头计算与高温.. 2页

安全文明施工、环境保护措施及雨季施工措施 28页

安全文化是企业平安的基础 27页

除数是两位数的除法单元测试题及答案 4页

安全教育的一般规定 27页

安全教育培训小组责任 23页

安全教育与培训安全隐患排查 28页

早期失智老年人护理ppt 23页

股骨转子间合并转子下骨折手术 27页

婴幼儿常见疾病及护理ppt 21页

胰腺炎的病情观察及护理PPT 20页

潜在茶庄团购方案策划相关7篇 12页

股东及其权利与义务 25页

大班音乐我是快乐的小蜗牛教案 31页

新生儿水痘护理病案讨论ppt 19页

肝衰竭的血液净化治疗从理论到实践课件 27页