数学建模-统计.pptVIP

数学建模 概率统计专题 专题板块系列 概率统计专题 part1：方差分析与协方差分析 单因素方差分析--计算 data ex；do a=1 to 3；input n @@； do i=1 to n； input x @@； Output；end；end； Cards； 8 21 29 24 22 25 30 27 26 10 20 25 25 23 29 31 24 26 20 21 6 24 22 28 25 21 26 ； proc anova; class a;model x=a； means a/duncan cldiff；run; 双因素方差分析-考虑交互作用-理论 可设双因素试验的一个因素为A，共有A1、A2、…、Ar等r个水平，另一个因素为B，共有B1、B2、…、Bs等s个水平。这两个因素的水平互相搭配各安排m次试验，其中A因素的A水平与B因素的B水平搭配安排试验所得到的样本为X，相应的观测值为x。 协方差分析 如果在单因素、双因素或多因素试验中有无法控制的因素x影响试验的结果Y，且x可以测量、x与Y之间又有显著的线性回归时，常常利用线性回归来矫正Y的观测值、消去x的差异对Y的影响。 例如，研究施肥对苹果树产量的影响，由于苹果树的长势不齐，必须消去长势对产量的影响。又如，研究饲料对动物增重的影响，由于动物的初重不同，必须消去初重对增重的影响。 data ex; do a=1 to 3;do i=1 to 8; input x y @ @;output ;end;end; cards; 47 54 58 66 53 63 46 51 49 56 56 66 54 61 44 50 52 54 53 53 64 67 58 62 59 62 61 63 63 64 66 69 44 52 48 58 46 54 50 61 59 70 57 64 58 69 53 66 ; proc glm;class a;model y=x a/solution; lsmeans a/stderr pdiff;run; 双因素协方差分析-不考虑交互作用 双因素协方差分析-考虑交互作用 data ex; do a=1 to 4; do b=1 to 2; do i=1 to 2;input x y @@;output;end;end;end; cards; 14.6 97.8 12.1 94.2 19.5 113.2 18.8 110.1 13.6 100.3 12.9 98.5 18.5 119.4 18.2 114.7 12.8 99.2 10.7 89.6 18.2 122.2 16.9 105.3 12.0 102.1 12.4 103.8 16.4 117.2 17.2 117.9 proc glm; class a b;model y=x a b a*b/solution; lsmeans a b/stderr pdiff;run; 一元非线性回归 data ex;input x y @@; x1=1/x;lx=log(x);ly=log(y); cards; 1 1.85 2 1.37 3 1.02 4 0.75 4 0.56 6 0.41 6 0.31 8 0.23 8 0.17 ; proc gplot; plot y*x; symbol i=spline v=star; proc reg;model y=x1; proc reg;model ly=lx; proc reg;model ly=x;run; 计算剩余平方和Q data ex;input x y @@; x1=1/x;lx=log(x);ly=log(y); y1=0.1159+1.9291*x1;q1+(y-y1)**2; y2=exp(0.9638-1.1292*lx);q2+(y-y2)**2; y3=exp(0.9230-0.3221*x);q3+(y-y3)**2; cards; 1 1.85 2 1.37 3 1.02 4 0.75 4 0.56 6 0.41 6 0.31 8 0.23 8 0.17 ; proc print;var q1-q3;run; 多元线性回归 多元线性回归模型的一般形式 多元回归分析数据格式 多元线性回归分析的一般步骤 多元回归分析的几个公式 实例：湖北省油菜投入与产出的统计分析 1投入指标 （1）土地（S）。土地用播种面积来表示。农作物播种面积是指当年从事农业 （2）劳动（L）。劳动用劳动用工数（成年劳动力一人劳动一天为一个工）来表示。劳动用工中包含着直接和间接生产用工。 （3）资本（K）。资本