求sas小程序！急！

祸劫

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

下面是一个线性加平台的程序，但是不过原点，贴上来，希望能有个参考。我想要拟合的数据也写在里面了。
" }8 o: Q: H7 t( b: F+ _* l
DATA A ;
INPUT aa ry;
cards;
11.8        58.72210953
. D$ a6 k8 D5 b% \9 J6 i20.1        88.0967142
20.1        81.24529722
5 Y- K! C+ ?! n5 ?; g20.1        84.29360695
# [- [, G/ C; c$ \, T6 o20.1        81.10726644
9 \4 ^3 X2 ?2 T+ \15.6        98.75298664
31.6        68.86474304
! r' v8 P+ ~" v  c1 a69.9        92.33568267
8.9        84.69947471
5 u6 L7 `% x, O9 _( C" t( x# J43.2        90.84231146
- b) Y' m5 [  n9 d3 ?7 E39.1        109.1894109
$ X) q& ?: |2 O+ A86.6        93.08064686
, w( }% ?8 u' n9 ]1 m% O1 ?27.5        86.64394553
. h" o7 ~1 X8 A. o; j+ N' l65.1        76.04877217
144.8        91.16522565
* S3 L$ k7 C! @9 |5 x) T. Y$ Q9 S6 ?4 m60        78.80126767
8 _* q0 V; _  I( p: p. ?$ R52        94.33156385
  X; n8 C) ~" u" S8 _7 u' U# @! g% P32.7        95.63820283
68.7        94.13253669
101.4        89.24645622
; F" C1 l, ^. O6 K5 ?84.25        101.196025
5 @  ^9 {' p5 k, p5 M- L133.1        85.42372881
' [9 _- h/ o* v4 i75.6        85.5045225
139.5        96.59318637
56.4        107.7097506
112        94.01759531
  ~, I* ]5 W2 [( g7 ]45        97.63031155
76.6        85.12898331
14        92.92788563
264        84.72287086
11.8        27.16763006
7 F  n3 x- b( d; M3 H! l: T8.9        66.1675435
4 b; m! w0 q" g% H& K8.4        69.86062718
48.2        95.39295393
* Q- T9 [9 d+ l2 b1 ~- z0 N8 t85        90.46265185
, J# L0 o9 w# F0 B  [( O& M56.66666667        96.16099071
2 ?* a8 n& Z% L6 R21.5        89.88411761
167        87.70897566
130.5        97.58812616
. A  y. |# N0 ^8 m% b' v167.5        99.29691725
) \4 k( c* G: [5 U' F27.5        94.3345668
27.8        83.32801701
33.7        84.03827751
121.3        62.1226565
) l: ~; V1 @9 O' H# `! F% l199.8        79.47899302
& X' }2 h. {+ G  [- t270.5        72.63651182
206        87.08424454
% i: |: i7 s2 j. ?( F107.5        72.13753714
134.9        88.94919802
5 g1 m- ?4 U/ ^' m& n% g6 _179        107.8421955
  S  I6 e1 O; T. X6 {) j  s6 {3.4        77.79831933
12.5        69.97430957
12.5        41.72274563
12.5        42.59395726
21.5        96.77574591
( l4 Z' H" ?. @0 H' _) |4 e18        80.52768162
3 y( t  L/ D# b* `* |17.5        91.15776212

: W: F4 X8 U0 Z$ C6 x5 Y0 O
8 N/ d! U$ \% W" l
' S8 \2 ]+ Y% x! N0 k' ^9 w* ?
, K- w% _% `0 W  x3 y& d


5 B: ~: B' L: b! L! F1 ?0 O6 n

;
TITLE1 'Yield Response Data';
TITLE2 'LRP Model';
  [# W9 T& d  @PROC NLIN DATA=A BEST=3;
5 o- B! ]( b5 D' Q! W2 YPARMS A=27 to 112 by 5 B=0 to 30 by 1 K 3.4 to 271 by 2;
! k9 }+ ?9 k; O- w8 CIF AA<K THEN DO;
MODEL RY = A + B*AA;
DER.A = 1;
DER.B = AA;
DER.K = 0;
END;
- i& u  s/ ]' a: W' mELSE DO;
7 \( [, y& H* ^6 U- q, C+ b" E6 PMODEL RY = A + B*K;
1 L6 r, ^3 M0 QDER.A = 1;
3 g' [9 k4 M& a; q: jDER.B = K;
' k4 T0 y/ }7 r9 y( n& Q+ E9 m- \DER.K = B;
END;
- n& T/ ^# {- ?  D( f* vFILE PRINT;
IF _OBS_ = 1 AND _MODEL_= 0 THEN DO;
PLATEAU = A + B*K;
! z1 O. {0 l8 d& iPUT PLATEAU=;
END;
PLATEAU = A + B*K;
6 j8 V$ Q5 w; [  Z" [/ W8 `ID PLATEAU;
OUTPUT OUT=NEW P=PRY PARMS=A B K SSE=SSE;
( K# k- C6 w- t* `) HRUN;
, [# N9 i1 G) |  g. wPROC PLOT;
PLOT RY*AA='+' PRY*AA='*'/OVERLAY;
; w( W: L3 f5 z$ Q6 hRUN;
4 m  B9 \' H) L7 q5 c/*   For you information:
# v/ j, Q+ s9 p     */
# B* p9 E2 w5 R% `: j6 b8 q0 b/*   This part of the program outputs Rsquares, significance,
( {+ j; l" q+ C1 l- Dplateau,  */
1 R7 B" H0 r, P. c/*   and soil-test value at the joining point into a nice dataset.
  i7 D" }% |$ Z: g    */
PROC MEANS NOPRINT;
0 F( J5 g( X; G+ ]. j( I- XVAR RY SSE A B K PLATEAU;
OUTPUT OUT=NEW2 N=TDF
3 _' g4 g1 N9 [% s  s* z( M                MEAN=RY SSE A B K PLATEAU
                CSS=CSST;
DATA NEW3; SET NEW2;
INTERCPT=A; SLOPE=B; JOINT=K;
RSQ=(CSST-SSE)/CSST;
& c; F' }3 T1 S# p8 O2 V1 c( XEDF=TDF-3;
SSR=CSST-SSE;
9 U. I: X( ?' nMSR=SSR/2;
3 Q; W! U) O8 K3 R7 ?! f$ |" {& m* {" {0 VMSE=SSE/EDF;
/ @3 `4 {3 P4 @- P+ z# i' NF=MSR/MSE;
1 Y% f, F# v7 ]9 jPROBF=1-(PROBF(F,2,EDF));
" D! u6 i$ S8 t6 T: a' r7 ~* c8 v# IKEEP INTERCPT SLOPE JOINT PLATEAU RSQ F PROBF;
" z2 B5 J8 x7 }- f# M8 ]$ _; _PROC PRINT;
RUN;
; ~/ U* f, g( M2 n' O0 c' F( J/ a

神秘了一场雨

自生为役~

我要体力啊啊

lynli

老牦

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

071029

额，不清楚