协方差分析非线性回归分析.ppt

八、协方差分析 协方差分析的原理：方差分析+回归分析 协方差分析的应用：校正非处理因素UN; 协方差分析GLM过程的调用 PROC GLM; CLASS 效应名; MODEL 因变量=效应名 自变量/SOLUTION; MEANS 效应名 /DUNCAN; LSMEANS效应名 /STDERR PDIFF TDIFF; INPUT T$ X Y@@; CARDS; CK 1.50 12.40 CK 1.85 12.00 CK 1.35 10.80 CK 1.45 10.00 CK 1.40 11.00 CK 1.45 11.80 CK 1.50 12.50 CK 1.55 13.40 CK 1.40 11.20 CK 1.50 11.60 CK 1.60 12.60 CK 1.70 12.50 1 1.35 10.20 1 1.20 9.40 1 1.45 12.20 1 1.20 10.30 1 1.40 11.30 1 1.30 11.40 1 1.15 12.80 1 1.30 10.90 1 1.35 11.60 1 1.15 8.50 1 1.35 12.20 1 1.20 9.30 2 1.15 10.00 2 1.10 10.60 2 1.10 10.40 2 1.05 9.20 2 1.40 13.00 2 1.45 13.50 2 1.30 13.00 2 1.70 14.80 2 1.40 12.30 2 1.45 13.20 2 1.25 12.00 2 1.30 12.80 3 1.20 12.40 3 1.00 9.80 3 1.15 11.60 3 1.10 10.60 3 1.00 9.20 3 1.45 13.90 3 1.35 12.80 3 1.15 9.30 3 1.10 9.60 3 1.20 12.40 3 1.05 11.20 3 1.10 11.00 ; PROC GLM; CLASS T; MODEL Y= T X/SOLUTION; MEANS T/DUNCAN; LSMEANS T/STDERR PDIFF TDIFF; RUN; 实行内容 一元协方差分析 多因素协方差分析 多元协方差分析 九、非线性回归分析 非线性模型---线性模型 非线性模型的建立 模型参数的估计 牛产犊月份-产奶量 t ：牛产犊月份 Y: 305天产奶量 奶牛泌乳曲线 t: 产奶时间 yt: 产奶量 a*b*c0 a:反映奶牛的产奶潜力 b:反映奶牛达到产奶高峰的速率 c:反映奶牛达到产奶高峰后下降的速率 动物生长曲线 yt:时间t时的体重 a:反映体重极限 k:反映接近这一极限的速率 b:反映达到最大生长率的时间 产蛋曲线 y(t):时间t时的产蛋量 t:产蛋时间 a:反映群体最大的产蛋潜力 b:反映群体产蛋能力的衰减速度 c:反映了开产日龄的变异情况 d:反映开产日龄 非线性模型参数的估计方法 线性化后的最小二乘估计法 — REG 过程 或者 GLM过程 直接求解法 — NLIN 过程 线性化后的最小二乘估计法— GLM过程 data ; input t y @@; t2= t*t; cards; 1 3833 2 3812 3 3769 4 3566 5 3482 6 3373 7 3477 9 3466 10 3809 11 3817 12 3885 ; PROC REG; MODEL Y=T T2; RUN; PROC GLM; MODEL Y=T T2/SOLUTION; RUN; REG过程的参数估计结果 直接求解法 — NLIN 过程 PROC NLN 选项； MODEL 因变量=模型表达式； PARMS 参数=初始值； DER．参数=偏导表达式； 奶牛泌乳曲线 t: 产奶时间 yt: 产奶量 a*b*c0 a:反映奶牛的产奶潜力 b:反映奶牛达到产奶高峰的速率 c:反映奶牛达到产奶高峰后下降的速率 DATA； INPUT T Y@@; CARDS; 1 1200 2 1260 3 1310 4 1420 5 1320 6 1280 7 1200 8 1140 9 930 10 840 ; PROC NLIN METHOD=DUN; MODEL Y=a*t**b*EXP(-c*t); PARMS a=1000 b= 0.4 c=0.2; RUN; DUD Initialization Dependent Variable Y DUD Initialization