SAS作业

SAS作业2011级栽培生理班2011301110033唐红玲1SAS课后习题作业院系：植物科学技术学院学号：2011301110033姓名：唐红玲二O一一年十二月SAS作业2011级栽培生理班2011301110033唐红玲2习题2.1.1观测200头某品种仔猪的一月窝重（单位：kg）得到频数表如下：组下限组上限组中值频数7.523.515.51023.539.531.51939.555.547.53055.571.563.56171.587.579.54987.5103.595.524103.5119.5111.57试用SAS作样本观测值整理后的描述性统计，并画出样本观测值的频率分布直方图。解：dataex;inputxf@@;cards;15.51031.51947.53063.56179.54995.524111.57;procunivariate;varx;freqf;run;Variable:xQuantiles(Definition5)datahist01;inputx@@;cards;15.515.515.515.515.515.515.515.515.515.531.531.531.531.531.531.531.531.531.531.531.531.531.531.531.531.531.531.531.547.547.547.547.547.547.547.547.547.547.547.547.547.547.547.547.547.547.547.547.547.547.547.547.547.547.547.547.547.547.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.563.579.579.579.579.579.579.579.579.579.579.579.579.579.579.579.579.579.579.579.579.5SAS作业2011级栽培生理班2011301110033唐红玲379.579.579.579.579.579.579.579.579.579.579.579.579.579.579.579.579.579.579.579.579.579.579.579.579.579.579.579.579.595.595.595.595.595.595.595.595.595.595.595.595.595.595.595.595.595.595.595.595.595.595.595.595.5111.5111.5111.5111.5111.5111.5111.5;procgchart;vbarx/type=pctspace=0;run;习题2.4.1已知我国14岁女学生的平均体重为43.38kg，从该年龄的女学生中抽查10名运动员的体重，分别为39,36,43,43,40,46,45,45,42,41kg。试用SAS检验这些运动员的体重与上述平均体重是否有显著的差异（α=0.05）解：假设H0：μ=43.38，运动员体重与平均体重无显著差异HA：μ≠43.38，运动员体重与平均体重有显著差异dataex;inputx@@;y=x-43.38;cards;39364343404645454241;procmeansmeanstdtprt;vary;run;SAS结果：因为P｛|X||T|｝0.05,所以接受H0：μ=43.38，即运动员体重与平均体重无显著差异.SAS作业2011级栽培生理班2011301110033唐红玲4习题2.4.2某小麦品种经过4代选育，从第5代和第6代中分别抽出10株得到它们株高的观测值分别为66,65,66,68,62,65,63,66,68,62以及64,61,57,65,65,63,62,63,64,60.试用SAS检验这两代的株高是否有显著的差异（=0.05）解：dataex;doa=1to2;doi=1to10;inputx@@;output;end;end;cards;6665666862656366686264615765656362636460;procttestcochran;classa;varx;procprint;run;Variable:x假设H0：第5代和第6代的株高没有显著差异，即H0：μ1=μ2,，对HA：μ1≠μ2,显著水平α=0.05如果认为方差相等，则DF=18，Prob|T|为0.0193；0.010.01930.05所以放弃H0，即这两代的株高有显著差异；如果认为方差不相等，则：根据Satterthwaite检验法作近似的t检验得：DF=17.7，Prob|T|为0.0194；0.010.01940.05或根据Cochran检验法作近似的t检验得：DF=9，Prob|T|为0.0302；0.010.03020.05二者结果均表示：放弃H0综上可得：这两代的株高存在显著差异。ANMeanStdDevStdError11065.12.18330.690421062.42.50330.7916SAS作业2011级栽培生理班2011301110033唐红玲5习题2.5.1测定4种种植密度下金皇后玉米的千粒重（单位：g）如下：种植密度千粒重1234247238214210258244227204256246221200251236218210试用SAS对4种种植密度下金皇后玉米千粒重的观测值作正态性检验。解：假设H0：4种种植密度下金皇后玉米千粒重的观测值服从正态分布；H1：4种种植密度下金皇后玉米千粒重的观测值不服从正态分布dataex;doa=1to4;doi=1to4;inputx@@;output;end;end;cards;247258256251238244246236214227221218210204200210;procunivariatenormal;run;W检验的临界值ω0.05=0.859，P｛Wω0.05=0.859｝=0.05SAS结果表明P｛W0.867632｝=0.02500.05，因此接受H1，即4种种植密度下金皇后玉米千粒重的观测值不服从正态分布。SAS作业2011级栽培生理班2011301110033唐红玲6习题3.1.1试验3种猪饲料的饲养效果，得到9头猪的增重（单位：kg）如下：饲料月增重12351，40,43，4823,25，2623,28试用SAS作方差分析，估计各个总体的未知参数μi和μ。若有必要，试求出两两总体均值差的双侧0.95置信区间。解：设H0为各μi相等，也就是各个处理之间没有显著的差异。dataex;doa=1to3;inputn@@;doi=1ton;inputx@@;output;end;end;cards;451404348323252622328;procanova;classa;modelx=a;meansa/lsdcldiff;run;用SAS计算得到方差分析表：SourceDFSumofSquaresMeanSquareFValuePrFModel2934.722222467.36111131.100.0007A2934.722222467.36111131.100.0007Error690.16666715.027778Ctotal81024.888889SAS结果表明0.00070.01,即否定H0：各μi相等，推断各饲料饲养效果存在显著差异，且饲料1与2、1与3的饲养效果之间有显著差异，而饲料2与3的饲养效果无显著差异。SAS作业2011级栽培生理班2011301110033唐红玲7习题3.2.2某种药物3与药物1和药物2做有效率的对比试验，各做5组，每组100例，得到有效率的观测值如下：药物有效率（%）12323，21，29，28，2440，33，33，12，3462，45，27，69，29试用SAS先作Levene的F检验，再作平方根反正弦转换后的方差分析。解：dataex;doa=1to3;doi=1to5;inputx@@;output;end;end;cards;232129282440333312346245276929;procanova;classa;modelx=a;meansa/hovtest;run;DependentVariable:x从LeveneF检验可知，Pr|t|为0.02040.05应拒绝原假设，认为方差是不同质的。分析转换后的数据所用的SAS程序为：dataex;doa=1to3;doi=1to5;inputx@@;y=arsin(sqrt(x/100));output;end;end;cards;232129282440333312346245276929;procanova;classa;modely=a;meansa/dunnett(‘3’)hovtest;run;DependentVariable:ySAS作业2011级栽培生理班2011301110033唐红玲8由此得Pr|t|为0.05160.05接受原假设，认为方差是同质的，即变换后的方差是齐性的。Dunnett'stTestsfory根据上述结果，可知处理3与1有显著差异，处理2与1没有显著差异。各均值差的95%的区间估计为：Mu3-mu1(0.00825，0.44087)Mu2-mu1(-0.16223，0.27039)习题3.3.2有4个学生分析5个品种的稻米，得到含N量（单位：mg）的测量值如下：品种学生ABCDE甲2.42.53.23.42.0乙2.62.23.23.51.8丙2.12.73.53.81.8丁2.42.73.13.22.3试用SAS作方差分析，估计各个总体的未知参数μi。和μ.j及μ。若有必要，试求出两两总体均值差的双侧0.95置信区间。解：设H01为各μi。相等，也就是各学生分析的5个品种稻米的含氮量之间没有显著的差异，H02为μ.j相等，也就是品种之间的含氮量没有显著的差异。dataex;doa=1to4;dob=1to5;inputx@@;output;end;end;cards;2.42.53.23.42.02.62.23.23.51.8SAS作业2011级栽培生理班2011301110033唐红玲92.12.73.53.81.82.42.73.13.22.3;procanova;classab;modelx=ab;meansab/lsdcldiff;run;用SAS计算得到方差分析表为：DependentVariable:x因此接受H01，放弃H02，认为μi。相等，μ.j不相等，即各学生分析的5个品种稻米的含氮量之间没有有显著的差异，品种之间的含氮量有显著的差异。tTests(LSD)forx习题3.4.2某种化工过程在3种浓度、4种温度下成品的得率如下：温度浓度10℃24℃38℃52℃2%14,1011,1113,910,124%9,710,87,116,106%5,1113,1412,1314,10试先作得率数据的平方根反正弦变换后再用SAS作方差分析。解：设H01为各μi.相等，即各种浓度对应的观测值之间没有显著差异，H02为各μ。j相等，即各种温度对应的观测值之间没有显著的差异，H03为各γij都等于0，即各种浓度与温度搭SAS作业2011级栽培生理班2011301110033唐红玲10配试验的交互作用都等于0.dataex;doa=1to3;dob=1to4;doi=1to2;inputx@@;y=ars