设计性试验.docx
- 文档编号:12163304
- 上传时间:2023-06-04
- 格式:DOCX
- 页数:14
- 大小:144.04KB
设计性试验.docx
《设计性试验.docx》由会员分享,可在线阅读,更多相关《设计性试验.docx(14页珍藏版)》请在冰点文库上搜索。
设计性试验
课程名称:
数据分析与sas实验
实验编号
及实验名称
设计性试验
系别
姓名
学号
班级
实验地点
实验日期
实验时数
6
指导教师
同组其他成员
成绩
一、实验目的及要求
目的:
1.使学生掌握服从正态分布总体的均值,方差的区间估计与假设检验以及非参数检验
2.掌握单因素,两因素及正交实验的方差分析方法。
利用REG过程进行多元线性回归。
要求:
实验形式
二、实验环境及相关情况(包含使用软件、实验设备、主要仪器及材料等)
装有sas系统的电脑一台
三、实验内容及步骤(包含简要的实验步骤流程)
实验内容:
1、用sas软件进行服从正态分布的单个总体的均值和方差进行参数估计与区间估计,与假设检验以及非参数检验。
2、单因素,两因素及正交实验的方差分析方法。
用REG过程进行多元线性回归
步骤:
1.
datazt;
retainseed0;
mul=0;
mu2=2;
sigmal=1;
sigma2=4;
doi=1to1000;
normal1=mul+sigmal*rannor(seed);
normal2=mu2+sigma2*rannor(seed);
output;
end;
dropseedimu1sigma1mu2sigma2;
run;
procunivariatedata=ztnormal;
varnormal1normal2;
histogramnormal1normal2;
probplotnormal1normal2;/*正态性假设检验*/
run;
procttestdata=zth0=0alpla=0.01;/*总体均值的假设检验*/
varnormal1;
run;
procttestdata=zth0=0;/*成配对样本的假设检验*/
pairnormal1*normal2;
run;
procmeansdata=lwh;
varpriceopen;
run;
dataszsh;
setlwh;
select(region);
when('sz')outputsz;
when('sh')outputsh;
end;
dropsex;
run;
procunivariatedata=sznormal;
varprice;
histogramprice;
probplotprice;
run;
procunivariatedata=shnormal;
varprice;
histogramprice;
probplotprice;
run;
procttestcochrandata=lwh;
classregion;
varprice;
run;
procnpar1waydata=lwh;/*两独立样本的假设检验(非参数)*/
classregion;
varprice;
run;
procmeansvardata=zt;/*总体方差的假设检验*/
varnormal2;
outputout=testvar=varex;
run;
dataa(drop=_type_);
settest;
k=_freq_-1;
chisq=k*varex/16;
p=1-probchi(chisq,k);
cil=cinv(0.025,k);
ci2=cinv(0.975,k);
procprintdata=anoobs;
run;
2.
datalwh;%导入要分析的数据
inputgdpsr@@;
cards;
5294.71212.335934.51366.9571711642.868964.42004.8210202.22122.0111962.52199.3514928.32357.2416909.22664.918547.92937.121617.83149.4826638.13438.37
34634.44348.9546759.45218.158478.16242.267884.67407.9974462.68651.1478345.29875.9582067.4611444.0889442.213395.2395933.316386.04
;
run;
procgplotdata=lwh;%作散点图、直观认识gdp*sr的回归关系
plotgdp*sr;
symboli=joinv=dotcolor=red;
run;
datazs;
setlwh;
a=1/gdp;b=1/sr;c=log(gdp);d=log(sr);e=exp(-sr);f=sr**2;
g=sr**3;h=sr**4;
i=sr**5;j=sr**6;
run;
proccorrdata=zs;%相关分析
vargdpsrfghij;
run;
procregdata=zs;%回归分析
modelgdp=srfghij;
run;
procregdata=zs;
modelgdp=ij/noint;
run;
procregdata=zs;
modelgdp=sr;
modelgdp=sr/noint;
modela=b/noint;
modela=b;
modelc=d;
modelc=sr;
modelc=b;
modelgdp=d;
modelgdp=ghij/noint;
run;
datalwh;%导入方差分析的数据
inputpainlevel$1.0codeine1.0acupuncture1.0relief2.1@@;
datalines;
a1100a2105a1206a2212b1103b2106b1207b2213c1104c2108c1208c2216d1104d2107d1209d2215e1106e2110e1215e2219
f1109f2114f1216f2223g1110g2118g1217g2221h1112h2117h1216h2224
;
run;
procanova;%单因素方差分析
classpainlevel;
modelrelief=painlevel;
meanspainlevel;
meanspainlevel/t;
run;
procanovadata=lwh;%双因素的交互作用检验
classpainlevelcodeineacupunctuer;
modelrelief=painlevelcodeine|acupuncture;
run;
procanovadata=lwh;%无交互作用的双因素方差分析
classpainlevelcodeineacupuncture;
modelrelief=painlevelcodeineacupuncture;
meanspainlevelcodeineacupuncture;
run;
四、实验结果(包括程序或图表、结论陈述、数据记录及分析等,可附页)
1.
normal1的正态性检验结果
Normal1的直方图
Normal1的QQ图
由于以上数据pvalue>0.05总成立,且Normal1的直方图是对称的,QQ图也是一条直线,故Normal1是服从正态分布的。
Normal1均值的假设检验
因为0.2617>0.05,所以normal1的均值u=0是被接受的
Normal1均值的假设检验
从以上数据<.0001显示,变量normal1和normal2的均值相等是不被接受的
Priceopen的means过程
上海的price正态性检验结果
由此可知,上海的price不服从正态分布
深圳的price正态性检验结果
由此可知,深圳的price也不服从正态分布
深圳的price正态性检验结果
深圳和上海市场股价的npar1way过程
以上pvalue小于0.001,从检验的数量结果显示深圳和上海市场的股价是有显著性的区别的
Normal2的方差假设检验
P>0.05显示接受normal2的方差等于16.79
2.
Gdp与sr的散点图
相关分析的结果
显然这些函数的相关性是相当显著的
回归分析结果
拟合效果非常不错,但是单个变量有些是不显著的
单因素方差分析结果
显示因素painlevel对relief是有显著性的影响
单因素方差分析中因素排名
单因素方差分析中因素分类
带交互作用的双因素方差分析
0.0923>0.05显示他们的交互作用是不显著的
五、实验总结(包括心得体会、问题回答及实验改进意见,可附页)
1.学会运用sas软件画散点图和直方图
2.对sas的作图命令还不熟练,作图过程缺乏细心
3.分析数据需对数据进行重新的整理,应该选择有代表性的变量进行分析
4.对运行出来的结果还不能完全读懂,可以从网上翻查资料
六、教师评语
1.完成所有规定的实验内容,实验步骤正确,结果正确;
2.完成绝大部份规定的实验内容,实验步骤正确,结果正确;
3.完成绝大部份规定的实验内容,实验步骤基本正确,结果基本正确;
4.基本完成规定的实验内容,实验步骤基本正确,所完成的结果基本正确;
5.未能很好地完成规定的实验内容或实验步骤不正确或结果不正确。
评定等级为:
优秀、良好、中等、及格、不及格
教师签名:
年月日
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 设计 试验
![提示](https://static.bingdoc.com/images/bang_tan.gif)