求sas小程序！急！

祸劫

本人现在急求sas的几个小程序，请各位高手帮忙！多谢啦！！
主要想要下面几种：
1 过原点的线性加平台，求拐点
2 过原点的二次加平台，求拐点
" M; j. d  X' T+ ~: Q3 对数加平台，求拐点
4 ~& t5 H, Y+ \* Z# O4 u% n3 B
: ?2 O. ^" ]* f% p3 N下面是一个线性加平台的程序，但是不过原点，贴上来，希望能有个参考。我想要拟合的数据也写在里面了。
1 w7 N. P$ L6 Q5 o- J: E8 \
DATA A ;
! [" v' h. z: w! v5 LINPUT aa ry;
cards;
+ o$ `0 O$ @$ m* Q3 H+ h3 p11.8        58.72210953
20.1        88.0967142
20.1        81.24529722
20.1        84.29360695
20.1        81.10726644
, ~$ l4 O5 n2 @7 Q- v6 Q15.6        98.75298664
31.6        68.86474304
5 Y* }. `$ t" e69.9        92.33568267
8.9        84.69947471
$ u) O7 Q' N4 e43.2        90.84231146
) \$ w8 o) j( \- k39.1        109.1894109
86.6        93.08064686
27.5        86.64394553
65.1        76.04877217
144.8        91.16522565
60        78.80126767
, `2 M7 C2 ~, S) u& q8 F- [52        94.33156385
/ R# p$ I& C3 q# ^32.7        95.63820283
5 F2 ?3 m- [: C0 E. {) X68.7        94.13253669
& p$ B( T0 j) W" ~3 L101.4        89.24645622
84.25        101.196025
) h% ?* A+ A% M7 K133.1        85.42372881
/ w1 F9 g0 s% j' g/ H. ]75.6        85.5045225
139.5        96.59318637
6 {0 o2 S( D2 U56.4        107.7097506
9 B1 Y3 J$ O  Y: c3 d112        94.01759531
3 A9 l' ~/ X3 z! z) x45        97.63031155
6 C6 D* E4 F# Z76.6        85.12898331
% O* F' B8 S& P5 v  R. I5 n14        92.92788563
' T, S3 X) l5 z6 j" `& k264        84.72287086
11.8        27.16763006
8.9        66.1675435
8.4        69.86062718
! H% p5 [, R. A4 Z# v% |  Q' v% i48.2        95.39295393
+ h$ K" n- k; C! I# e85        90.46265185
9 i4 ]2 D# `8 @0 `% ], u+ T56.66666667        96.16099071
4 }- E6 W  o! R7 u  Q5 h21.5        89.88411761
, k9 l. M/ B% Q' c% u+ Q0 T: X167        87.70897566
+ b5 u3 t6 q8 C! x' B' h$ h130.5        97.58812616
167.5        99.29691725
27.5        94.3345668
) j/ _# L& g4 A& o27.8        83.32801701
33.7        84.03827751
121.3        62.1226565
199.8        79.47899302
270.5        72.63651182
  T  i7 Q/ R  B' ^6 v3 D; K206        87.08424454
107.5        72.13753714
. T( i$ p- M6 S  F134.9        88.94919802
179        107.8421955
1 M: e2 n& w( k3 b! n3.4        77.79831933
12.5        69.97430957
6 n# v( D$ M- J- K0 o+ z  q12.5        41.72274563
12.5        42.59395726
4 W' t& \3 i6 C  e7 E21.5        96.77574591
1 b: B9 i6 f% L' j* A! Y18        80.52768162
, f' g6 U* G: G3 F# q17.5        91.15776212


, ^' R( [- v6 I# L8 G5 h; x4 d% b
  S& @% U/ y( D8 x( S6 a; U

( g; o$ D0 s3 n* W& d1 X* V0 ~: o+ |



;
7 F( c/ j7 d$ V' |8 f5 h0 a7 QTITLE1 'Yield Response Data';
TITLE2 'LRP Model';
* n8 `4 O7 h$ s/ S# GPROC 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;
0 a% o+ e: K3 x, w* y5 V  E0 a' V9 uDER.B = AA;
, F: O% Q2 D# |4 F3 XDER.K = 0;
9 a. K7 ?& D1 Q; s  ?% G& I% c0 e; lEND;
# H# U  B6 c& o$ N, M, m5 _ELSE DO;
MODEL RY = A + B*K;
DER.A = 1;
DER.B = K;
' N, P, z! k, `9 F( oDER.K = B;
END;
FILE PRINT;
IF _OBS_ = 1 AND _MODEL_= 0 THEN DO;
0 q; K2 a) @6 V' P: z. P. rPLATEAU = A + B*K;
PUT PLATEAU=;
. U( d5 z8 t* qEND;
PLATEAU = A + B*K;
ID PLATEAU;
  P1 W! e0 r1 p! J. FOUTPUT OUT=NEW P=PRY PARMS=A B K SSE=SSE;
RUN;
! i8 X0 \" x  rPROC PLOT;
* l8 }: H  U/ ]+ z( D# {PLOT RY*AA='+' PRY*AA='*'/OVERLAY;
RUN;
/*   For you information:
     */
/*   This part of the program outputs Rsquares, significance,
plateau,  */
+ K2 i+ k! h9 p/*   and soil-test value at the joining point into a nice dataset.
    */
9 h1 K, I" [/ B- o: E: k4 F2 _% bPROC MEANS NOPRINT;
VAR RY SSE A B K PLATEAU;
OUTPUT OUT=NEW2 N=TDF
                MEAN=RY SSE A B K PLATEAU
6 K# b" r! L/ @6 I- e. W                CSS=CSST;
DATA NEW3; SET NEW2;
* A# h) d+ M; U# A& Z" Y6 }! [INTERCPT=A; SLOPE=B; JOINT=K;
RSQ=(CSST-SSE)/CSST;
EDF=TDF-3;
SSR=CSST-SSE;
0 E% G2 K6 a( H3 cMSR=SSR/2;
+ X. W6 L5 d, P0 E) |MSE=SSE/EDF;
0 ^; A* C- ~1 |. e* n* k, g) ZF=MSR/MSE;
PROBF=1-(PROBF(F,2,EDF));
KEEP INTERCPT SLOPE JOINT PLATEAU RSQ F PROBF;
PROC PRINT;
1 \& m( j  z& X9 u3 b. X+ n0 eRUN;

神秘了一场雨

自生为役~

我要体力啊啊

lynli

老牦

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

071029

额，不清楚