祸劫 发表于 2011-10-18 21:18

求sas小程序！急！

本人现在急求sas的几个小程序，请各位高手帮忙！多谢啦！！
主要想要下面几种：
1 过原点的线性加平台，求拐点
2 过原点的二次加平台，求拐点
3 对数加平台，求拐点

下面是一个线性加平台的程序，但是不过原点，贴上来，希望能有个参考。我想要拟合的数据也写在里面了。

DATA A ;
INPUT aa ry;
cards;
11.8        58.72210953
20.1        88.0967142
20.1        81.24529722
20.1        84.29360695
20.1        81.10726644
15.6        98.75298664
31.6        68.86474304
69.9        92.33568267
8.9        84.69947471
43.2        90.84231146
39.1        109.1894109
86.6        93.08064686
27.5        86.64394553
65.1        76.04877217
144.8        91.16522565
60        78.80126767
52        94.33156385
32.7        95.63820283
68.7        94.13253669
101.4        89.24645622
84.25        101.196025
133.1        85.42372881
75.6        85.5045225
139.5        96.59318637
56.4        107.7097506
112        94.01759531
45        97.63031155
76.6        85.12898331
14        92.92788563
264        84.72287086
11.8        27.16763006
8.9        66.1675435
8.4        69.86062718
48.2        95.39295393
85        90.46265185
56.66666667        96.16099071
21.5        89.88411761
167        87.70897566
130.5        97.58812616
167.5        99.29691725
27.5        94.3345668
27.8        83.32801701
33.7        84.03827751
121.3        62.1226565
199.8        79.47899302
270.5        72.63651182
206        87.08424454
107.5        72.13753714
134.9        88.94919802
179        107.8421955
3.4        77.79831933
12.5        69.97430957
12.5        41.72274563
12.5        42.59395726
21.5        96.77574591
18        80.52768162
17.5        91.15776212









;
TITLE1 'Yield Response Data';
TITLE2 'LRP Model';
PROC NLIN DATA=A BEST=3;
PARMS A=27 to 112 by 5 B=0 to 30 by 1 K 3.4 to 271 by 2;
IF AA<K THEN DO;
MODEL RY = A + B*AA;
DER.A = 1;
DER.B = AA;
DER.K = 0;
END;
ELSE DO;
MODEL RY = A + B*K;
DER.A = 1;
DER.B = K;
DER.K = B;
END;
FILE PRINT;
IF _OBS_ = 1 AND _MODEL_= 0 THEN DO;
PLATEAU = A + B*K;
PUT PLATEAU=;
END;
PLATEAU = A + B*K;
ID PLATEAU;
OUTPUT OUT=NEW P=PRY PARMS=A B K SSE=SSE;
RUN;
PROC PLOT;
PLOT RY*AA='+' PRY*AA='*'/OVERLAY;
RUN;
/*   For you information:
     */
/*   This part of the program outputs Rsquares, significance,
plateau,  */
/*   and soil-test value at the joining point into a nice dataset.
    */
PROC MEANS NOPRINT;
VAR RY SSE A B K PLATEAU;
OUTPUT OUT=NEW2 N=TDF
                MEAN=RY SSE A B K PLATEAU
                CSS=CSST;
DATA NEW3; SET NEW2;
INTERCPT=A; SLOPE=B; JOINT=K;
RSQ=(CSST-SSE)/CSST;
EDF=TDF-3;
SSR=CSST-SSE;
MSR=SSR/2;
MSE=SSE/EDF;
F=MSR/MSE;
PROBF=1-(PROBF(F,2,EDF));
KEEP INTERCPT SLOPE JOINT PLATEAU RSQ F PROBF;
PROC PRINT;
RUN;

神秘了一场雨 发表于 2011-12-21 19:12

{:3_41:}{:3_41:}

自生为役~ 发表于 2012-1-14 21:59

我要体力啊啊

lynli 发表于 2012-2-2 20:21

{:3_41:}{:3_41:}{:3_41:}

老牦 发表于 2012-2-19 22:18

。。。。。。。。。。。。。。。。。。。。。。。。。。。

071029 发表于 2012-9-14 13:48

额，不清楚
页: [1]
查看完整版本: 求sas小程序！急！