求sas小程序！急！

祸劫

本人现在急求sas的几个小程序，请各位高手帮忙！多谢啦！！
1 D* `$ h, f6 n; L8 `  a) W主要想要下面几种：
1 过原点的线性加平台，求拐点
7 G9 t3 n( K' D4 T2 过原点的二次加平台，求拐点
3 对数加平台，求拐点

下面是一个线性加平台的程序，但是不过原点，贴上来，希望能有个参考。我想要拟合的数据也写在里面了。
1 Q6 \: o% I9 q1 X4 b, }
! g" W' M0 q5 ]6 f# v: F! M  l: `DATA A ;
INPUT aa ry;
' I$ F' ]9 p3 X/ h7 _cards;
11.8        58.72210953
5 n  A  Z/ ]3 i5 q1 q20.1        88.0967142
! N$ a$ b4 h5 L20.1        81.24529722
9 C. v$ g0 [& m% q, F7 s$ ^20.1        84.29360695
20.1        81.10726644
15.6        98.75298664
31.6        68.86474304
7 m  D- O( y0 {69.9        92.33568267
8.9        84.69947471
43.2        90.84231146
39.1        109.1894109
% [1 ^' g  C+ \' q86.6        93.08064686
! j4 n; `" D3 y27.5        86.64394553
65.1        76.04877217
' T7 S& H2 t6 V/ e  q; W7 g8 k144.8        91.16522565
! i' n7 Q" j) X8 g4 R  I8 O0 A7 K60        78.80126767
5 p8 D% x/ J$ J* I" j52        94.33156385
32.7        95.63820283
& Z4 Y: }& V% ?! X& P) M. ?68.7        94.13253669
" V, E( c1 V; U' i2 j7 _101.4        89.24645622
( ?2 C7 I, [- g8 E$ n% N& ]) K84.25        101.196025
/ W* W2 ]5 U5 H133.1        85.42372881
75.6        85.5045225
5 I$ W. H- t& M3 W139.5        96.59318637
1 g* s" f' ?8 l3 g56.4        107.7097506
+ p+ E' E) R- N6 B4 L  o9 o6 ]112        94.01759531
) I" a4 Y, q5 `0 i4 o& L45        97.63031155
76.6        85.12898331
14        92.92788563
264        84.72287086
% D0 `  Y, D- N2 j. h" G11.8        27.16763006
8 u, n6 s8 a& U. r8 |3 V8.9        66.1675435
2 I1 ]+ a! p" q2 C1 T* j& f$ w8.4        69.86062718
48.2        95.39295393
# s# E! i% C8 N8 A3 Q: N( r85        90.46265185
56.66666667        96.16099071
21.5        89.88411761
167        87.70897566
2 x: o8 u+ a. I1 {5 b130.5        97.58812616
# `& n5 _. W( @/ `1 d. _167.5        99.29691725
7 m9 C$ p  ?, Y, F& o# t2 N* c' C27.5        94.3345668
27.8        83.32801701
) L- [% r. {/ G/ e# `& ?33.7        84.03827751
121.3        62.1226565
199.8        79.47899302
+ H3 W" H1 v3 Y( y1 ^% I& _" m270.5        72.63651182
206        87.08424454
107.5        72.13753714
3 m  h. P4 v2 W. n8 a3 t4 {: H134.9        88.94919802
! \  _  a# E" L- X! o8 S179        107.8421955
3.4        77.79831933
; Z  v% ^5 |. q) v- y' R" \( l12.5        69.97430957
12.5        41.72274563
12.5        42.59395726
21.5        96.77574591
18        80.52768162
: n6 }  H1 D1 N$ K; [- u7 C, Y17.5        91.15776212
0 I) \& w" i5 t6 o' A
( ~$ f" s# P* S8 v  k  o

. C1 R& |6 D4 d$ y% E, l) |# C

: h, n5 ~# s( P6 S- T4 w  Q- V
$ n& Q1 P$ z7 T' C2 L2 m* n! u2 i

' v, V2 y! G! y& e2 C& z
! h2 N4 H" y! S. E;
TITLE1 'Yield Response Data';
! Z& H9 d& d. @9 p+ CTITLE2 'LRP Model';
PROC NLIN DATA=A BEST=3;
! {4 }* F% G) ]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;
2 d. d( Z' k& @& H+ JDER.B = AA;
DER.K = 0;
END;
* l- o3 h* n9 tELSE DO;
MODEL RY = A + B*K;
DER.A = 1;
" V+ c+ ~) }5 o7 v  tDER.B = K;
DER.K = B;
END;
5 V& S$ P% T$ zFILE PRINT;
IF _OBS_ = 1 AND _MODEL_= 0 THEN DO;
9 a9 P- u. M9 S8 K4 p$ G# WPLATEAU = A + B*K;
PUT PLATEAU=;
5 l& H& x& p* dEND;
# t. M( r- _" Q, v0 G1 l4 V: M6 PPLATEAU = A + B*K;
1 E# T% Z# N1 U0 b- }& B! F& bID PLATEAU;
OUTPUT OUT=NEW P=PRY PARMS=A B K SSE=SSE;
RUN;
8 `) f/ b+ a: m2 c: ~5 Q) l/ uPROC PLOT;
4 r: V5 L- W* m% k* G$ {0 DPLOT RY*AA='+' PRY*AA='*'/OVERLAY;
, ?+ d% M, i1 c7 z6 t- ?RUN;
8 S# i/ D  r6 J* s+ n$ d/*   For you information:
& W& q( s: Z2 `" @) }     */
/*   This part of the program outputs Rsquares, significance,
$ K! A9 R3 ~9 c! [8 P, Vplateau,  */
/*   and soil-test value at the joining point into a nice dataset.
    */
. M1 F; {1 K  N0 I% G/ YPROC MEANS NOPRINT;
VAR RY SSE A B K PLATEAU;
OUTPUT OUT=NEW2 N=TDF
  w) i  m4 e- M/ s                MEAN=RY SSE A B K PLATEAU
- |2 W5 {! k, @- x) R* D" p                CSS=CSST;
DATA NEW3; SET NEW2;
, x) l3 m5 f$ u' i3 d! N1 OINTERCPT=A; SLOPE=B; JOINT=K;
: v) a8 e1 L2 w; v4 X3 A" gRSQ=(CSST-SSE)/CSST;
6 s7 F" r/ Y4 xEDF=TDF-3;
: G8 T' _  g/ D/ J3 [SSR=CSST-SSE;
. D2 ?. a" D' r+ ]. {$ _0 n0 pMSR=SSR/2;
MSE=SSE/EDF;
F=MSR/MSE;
PROBF=1-(PROBF(F,2,EDF));
/ H5 e5 x6 T) rKEEP INTERCPT SLOPE JOINT PLATEAU RSQ F PROBF;
PROC PRINT;
& B; v( W  a( Y2 MRUN;

神秘了一场雨

自生为役~

我要体力啊啊

lynli

老牦

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

071029

额，不清楚