求sas小程序！急！

祸劫

本人现在急求sas的几个小程序，请各位高手帮忙！多谢啦！！
- M8 U6 w: T3 E2 @- g主要想要下面几种：
1 过原点的线性加平台，求拐点
/ O" @9 m0 t$ d2 g6 R3 p5 j2 过原点的二次加平台，求拐点
  E( V' g1 e5 z5 N0 }3 对数加平台，求拐点

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

" n# a, V4 D7 m& d% oDATA A ;
INPUT aa ry;
cards;
11.8        58.72210953
$ m6 b* z8 G1 b5 M20.1        88.0967142
20.1        81.24529722
! U2 o+ O# _) M% T2 B+ {8 f- ?6 h- P20.1        84.29360695
20.1        81.10726644
15.6        98.75298664
31.6        68.86474304
4 ?+ c8 b- {0 O) Y, K. h+ v5 }69.9        92.33568267
8.9        84.69947471
43.2        90.84231146
# G- U. ~; t  q) _7 x6 V7 x39.1        109.1894109
86.6        93.08064686
2 P4 U, |# p5 V0 h' p6 ?27.5        86.64394553
65.1        76.04877217
8 ~- w. d: c$ ~# u, |144.8        91.16522565
: V6 q2 G5 {5 [; `60        78.80126767
; n. i8 J; `1 {. p. v  ]52        94.33156385
/ S, f+ A# \+ m7 T8 b3 |32.7        95.63820283
68.7        94.13253669
101.4        89.24645622
! s0 b  g  ~6 f84.25        101.196025
133.1        85.42372881
" L, k( J: c% X0 n- e" M* p3 ^75.6        85.5045225
4 B1 E( Q) W  L% ]4 r' `" H139.5        96.59318637
) c; K& q: Z% N56.4        107.7097506
3 j; V* a' e9 y, P6 E1 @% y112        94.01759531
45        97.63031155
76.6        85.12898331
1 ~' K2 {2 k9 d. U9 u' X! r5 K14        92.92788563
; n9 ]+ l6 l' i; J8 {264        84.72287086
/ U; {/ p. [/ z! K11.8        27.16763006
8.9        66.1675435
/ h2 P" O" L$ Y$ u3 ?# F9 q& f! i8.4        69.86062718
& \* G/ K( C7 ^& C: c4 M48.2        95.39295393
85        90.46265185
56.66666667        96.16099071
' [0 ?# N( V( S21.5        89.88411761
! h6 t" P  S. v5 v! ~/ y! p167        87.70897566
130.5        97.58812616
167.5        99.29691725
+ L! w8 c  b$ H4 h% z6 f0 B27.5        94.3345668
; C! t# X( F$ C3 L+ I27.8        83.32801701
33.7        84.03827751
121.3        62.1226565
199.8        79.47899302
5 I# R+ a+ @8 K% r! t270.5        72.63651182
206        87.08424454
107.5        72.13753714
# V; ]9 B+ N8 s9 @5 c" v7 p1 ]134.9        88.94919802
179        107.8421955
/ `, q8 N# x( u  C3.4        77.79831933
12.5        69.97430957
12.5        41.72274563
  n( T# U1 Z8 c  }# O12.5        42.59395726
21.5        96.77574591
5 d; ]) `  y! V# U0 I2 V2 P5 E+ h18        80.52768162
; i2 Z9 i: K; Q5 q7 N17.5        91.15776212



' G7 X/ N6 [0 l1 ?) k
9 G' q2 P6 C; V  x# m7 p

% x6 m4 f# {, \8 C1 Q
0 z- r) A# n' h' l; I4 I2 a

  q' d! P$ r* ]8 g: V9 b/ S;
TITLE1 'Yield Response Data';
TITLE2 'LRP Model';
* O% o3 \( x6 b: D4 bPROC 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;
9 n* E. @, j( G3 v9 k5 E4 O1 ~MODEL RY = A + B*K;
! K  [! _- u2 cDER.A = 1;
1 P( [, j  P' k. v% N5 ~DER.B = K;
: J3 e' Q5 b% E: _& [( wDER.K = B;
. H% ~! _. Q) dEND;
' \/ L) x& ]+ h* G, VFILE PRINT;
IF _OBS_ = 1 AND _MODEL_= 0 THEN DO;
PLATEAU = A + B*K;
PUT PLATEAU=;
  a% u% B% {! W7 T1 DEND;
PLATEAU = A + B*K;
ID PLATEAU;
# `% a* c0 ~: u3 G" ^' tOUTPUT OUT=NEW P=PRY PARMS=A B K SSE=SSE;
5 x* k# D' a. e  T% X  `! J* [) xRUN;
0 Z2 v. P" ~. m0 D+ j% OPROC PLOT;
  F" R8 _) t# hPLOT RY*AA='+' PRY*AA='*'/OVERLAY;
RUN;
3 b0 b" ]. ~8 |$ X1 W9 G/*   For you information:
# E) \; |, Q# T( R8 w; I! R# }* q     */
/*   This part of the program outputs Rsquares, significance,
plateau,  */
/*   and soil-test value at the joining point into a nice dataset.
9 v, y, s' }9 w1 k3 A9 J# H: p& p    */
PROC MEANS NOPRINT;
VAR RY SSE A B K PLATEAU;
OUTPUT OUT=NEW2 N=TDF
% q$ T% L/ u. _0 q- e4 f& J                MEAN=RY SSE A B K PLATEAU
                CSS=CSST;
1 g( S, w* }, q! M0 VDATA NEW3; SET NEW2;
8 d( V3 X9 x* @0 c  m" l, SINTERCPT=A; SLOPE=B; JOINT=K;
RSQ=(CSST-SSE)/CSST;
3 j# `" R6 R8 s1 G# JEDF=TDF-3;
6 K$ P- O4 C6 o1 W9 v! gSSR=CSST-SSE;
MSR=SSR/2;
MSE=SSE/EDF;
* b" n) v. X/ y  I2 S: x- u! \F=MSR/MSE;
7 P$ U8 @' g) N6 X. a) E3 bPROBF=1-(PROBF(F,2,EDF));
KEEP INTERCPT SLOPE JOINT PLATEAU RSQ F PROBF;
% @! ?& N  o7 D0 H. zPROC PRINT;
" m7 f) j9 e/ Z. [, Y5 |  URUN;

神秘了一场雨

自生为役~

我要体力啊啊

lynli

老牦

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

071029

额，不清楚