求sas小程序！急！

祸劫

本人现在急求sas的几个小程序，请各位高手帮忙！多谢啦！！
主要想要下面几种：
1 过原点的线性加平台，求拐点
2 过原点的二次加平台，求拐点
7 @1 B1 V: `5 i3 对数加平台，求拐点
/ g" s* W4 {2 B: {3 k5 l
% B1 ], I$ l% @% @下面是一个线性加平台的程序，但是不过原点，贴上来，希望能有个参考。我想要拟合的数据也写在里面了。
  S0 |/ Q& j' k" i
! r7 ?! z. w5 H8 E9 N* U3 X: c+ C. x4 g& ~DATA A ;
' u; `& I2 L; R' Q" [$ j  JINPUT aa ry;
cards;
11.8        58.72210953
' I4 ?; S( e0 R7 O* J20.1        88.0967142
$ ]& R" `5 a1 R& P20.1        81.24529722
20.1        84.29360695
4 P" l, ?$ T! Y20.1        81.10726644
1 S! T$ a. M* V5 U8 t9 D. [5 [7 f15.6        98.75298664
5 @* K' A* E: n31.6        68.86474304
4 W) U8 O7 X3 M3 F# E69.9        92.33568267
! I6 K9 z6 I, q$ `0 ?, O8.9        84.69947471
43.2        90.84231146
, ]7 s5 n! ]9 M, P5 c6 \+ f39.1        109.1894109
86.6        93.08064686
5 X- o, x# o9 o4 l! g! `9 Z4 g27.5        86.64394553
- T; U) u; g8 g8 N2 ^8 S65.1        76.04877217
% K% g% g' \% K5 o6 R7 Y144.8        91.16522565
60        78.80126767
% d6 F4 M" o. w& Z4 u: y' G3 _( A52        94.33156385
) J+ m' r  U( F8 T+ ?; s32.7        95.63820283
, n3 s1 t& h) f1 ]! P68.7        94.13253669
; t" O6 G$ c* D# S1 ]& ?101.4        89.24645622
84.25        101.196025
. ?7 {( Y# W+ k1 x( a" ?* h9 J133.1        85.42372881
75.6        85.5045225
139.5        96.59318637
5 A$ t% t/ G  |; o+ k56.4        107.7097506
7 @: `* N' M5 K" v- o7 \1 x) R112        94.01759531
$ r* M  {/ Q+ X45        97.63031155
76.6        85.12898331
14        92.92788563
264        84.72287086
1 _0 i$ e. d' }& r1 U11.8        27.16763006
8.9        66.1675435
8.4        69.86062718
7 D) V8 Z+ S/ J2 G9 l48.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
7 d% P  z; @% a3 P3 \* D5 |27.8        83.32801701
33.7        84.03827751
121.3        62.1226565
( u: Z/ z/ Z+ M5 P* u2 s199.8        79.47899302
270.5        72.63651182
206        87.08424454
" h. {, J. s0 S! R107.5        72.13753714
134.9        88.94919802
179        107.8421955
3.4        77.79831933
12.5        69.97430957
12.5        41.72274563
) M6 R4 H" g/ v+ N) U1 h12.5        42.59395726
) k; d2 m6 h( {: {& P21.5        96.77574591
18        80.52768162
17.5        91.15776212

$ U4 h+ o( T$ e7 O! ]6 g
& E* v3 O# A# [! U
# L! s- b+ w2 d5 s1 L6 f3 W' c5 J



0 H; \" C+ H& Q4 H( d
& ?, {! w1 D6 c: k
;
TITLE1 'Yield Response Data';
& H9 g6 s) ~" k7 B# `) bTITLE2 'LRP Model';
) O4 N" N. f) n# pPROC NLIN DATA=A BEST=3;
( W  G: t. z$ O: ]' o8 NPARMS A=27 to 112 by 5 B=0 to 30 by 1 K 3.4 to 271 by 2;
' f4 _, d% C* H2 Y5 v6 g& RIF AA<K THEN DO;
MODEL RY = A + B*AA;
8 G$ @. [' t2 S6 _6 ]" BDER.A = 1;
DER.B = AA;
DER.K = 0;
END;
5 z! Y: D! d' B  R' r4 PELSE DO;
MODEL RY = A + B*K;
3 |8 l: j, o/ s9 x& T* W4 G0 lDER.A = 1;
4 B3 X; O! }3 g# e' ?( n. d6 P6 ?DER.B = K;
DER.K = B;
END;
4 y2 N4 r/ B8 ~% c* BFILE PRINT;
IF _OBS_ = 1 AND _MODEL_= 0 THEN DO;
/ v; N" V! \/ p& t+ b) l* v+ jPLATEAU = A + B*K;
PUT PLATEAU=;
END;
& W0 ^( E' V+ [6 D( x5 b* V- cPLATEAU = A + B*K;
4 v( B% P* a0 d6 _ID PLATEAU;
OUTPUT OUT=NEW P=PRY PARMS=A B K SSE=SSE;
+ C$ e1 O- r1 ARUN;
PROC PLOT;
PLOT RY*AA='+' PRY*AA='*'/OVERLAY;
0 I8 T- y" c9 W2 O. h  KRUN;
/*   For you information:
. n8 a9 Q8 i$ g# ?3 p     */
5 y0 h$ C9 D0 w2 ^5 Y) Q/*   This part of the program outputs Rsquares, significance,
1 P9 [: P5 ^# ~, S: N# F: Fplateau,  */
/*   and soil-test value at the joining point into a nice dataset.
$ D$ b+ u, c( @" [# r/ V    */
2 d. |: T7 X+ v/ o8 [PROC MEANS NOPRINT;
VAR RY SSE A B K PLATEAU;
5 ^4 ~1 {9 k8 e8 X" IOUTPUT OUT=NEW2 N=TDF
% `5 J9 B4 g3 T                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;
2 }+ I: \. C" S. ^SSR=CSST-SSE;
MSR=SSR/2;
MSE=SSE/EDF;
F=MSR/MSE;
- Z5 h0 y0 @+ p$ @7 UPROBF=1-(PROBF(F,2,EDF));
4 H0 V1 v" \$ e' p. W  O; r7 N7 {KEEP INTERCPT SLOPE JOINT PLATEAU RSQ F PROBF;
PROC PRINT;
' j. C. a) p" K' q6 {RUN;

神秘了一场雨

自生为役~

我要体力啊啊

lynli

老牦

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

071029

额，不清楚