ods listing;
ods html path='\\NETAPP8\fid_forageAgronomy\home\jjpittman\PHD\Class Spring 2013\SOIL5112\LPHW';

DATA one;
input x y;
cards;
0	2200
50	3160
75	3760
100	4680
150	4560
200	5440
;

proc nlin data = one best = 3;
  parms b0=1000 to 15000 by 500 b1=10 to 100 by 10 njoint=50 to 80 by 5;
  if x<njoint then do;
  model y = b0 + b1*x;
der.b0=1;
der.b1=x;
der.njoint=0;
end;
else do;
model y=b0+b1*njoint;
der.b0=1;
der.b1=njoint;
der.njoint=b1;
end;
file print;
if _obs_ =1 and _model_ =0 then do;
  plateau = b0 + b1*njoint;
  put plateau=;
end;
  plateau=b0+b1*njoint;
id plateau;
  output out = new p = pry parms=b0 b1 njoint sse=sse;
run;
proc plot;
plot y*x='+' pry*x='*'/overlay;
run;
proc means noprint;
var y sse b0 b1 njoint plateau;
output out = new2 n = tdf
mean = y sse b0 b1 njoint plateau
css=csst;
data new3; set new2;
intercpt=b0; slope=b1; joint=njoint;
rsq=(csst-sse)/csst;
edf=tdf-3;
ssr=csst-sse;
msr=ssr/2;
mse=sse/edf;
f=msr/mse;
probf=1-(probf(f,2,edf));
keep intercpt slope joint plateau rsq f probf;
proc print;
run;