Enkytræer og tørke

data a;set kmalc50;
treat=-treat;
mort=n0-nt;mortrate=mort/n0;
run;

proc nlmixed data=a df=264;
parms
    a1= -3.4472  b1= 0.1396 c1= 0.009502
    a2= -2.9870  b2= 0.2174 c2= 0.1012
    a3= -1.9326  b3= 0.1386
    a4= -2.1757  b4= 0.1765 c4= 0.06954
    a5=  1.4405             c5= 0.01358
    a6= -4.6539  b6= 2.3555 c6= 0.04574 ;

pr_response=
(tid=  4 and spp='ALI ')*(c1+(1-c1)* (probnorm(a1+b1*treat)))+

(tid=  4 and spp='CRYP')*(c2+(1-c2)* (probnorm(a2+b2*treat)))+

(tid=  7 and spp='ALI')*(.25+(1-.25)*(probnorm(a3+b3*treat)))+

(tid=  7 and spp='CRYP')*(c4+(1-c4)* (probnorm(a4+b4*treat)))+

(tid= 28 and spp='CRYP')*(c5*(exp(a5*treat)))+

(tid= 56 and spp='ALI ')*(c6+(1-c6)*  (probnorm(a6+b6*treat)))
;

model mort ~ binomial(n0,pr_response);

LC_A4  =(probit(0.5)-a1)/b1;
LC_C4  =(probit(0.5)-a2)/b2;
LC_A7  =(probit(0.5)-a3)/b3;
LC_C7  =(probit(0.5)-a4)/b4;
LC_C28 =a5**-1*log((c5+(1-c5)*.5)/c5);
LC_A56 =(probit(0.5)-a6)/b6;

estimate 'LC50  4 ALI ' LC_A4  ;
estimate 'LC50  4 CRYP' LC_C4  ;
estimate 'LC50  7 ALI ' LC_A7  ;
estimate 'LC50  7 CRYP' LC_C7  ;
estimate 'LC50 28 CRYP' LC_C28 ;
estimate 'LC50 56 ALI ' LC_A56 ;

contrast 'LC_A4  vs. LC_C4 ' LC_A4 -LC_C4 ;
contrast 'LC_A4  vs. LC_A7 ' LC_A4 -LC_A7 ;
contrast 'LC_A4  vs. LC_C7 ' LC_A4 -LC_C7 ;
contrast 'LC_A4  vs. LC_C28' LC_A4 -LC_C28;
contrast 'LC_A4  vs. LC_A56' LC_A4 -LC_A56;
contrast 'LC_C4  vs. LC_A7 ' LC_C4 -LC_A7 ;
contrast 'LC_C4  vs. LC_C7 ' LC_C4 -LC_C7 ;
contrast 'LC_C4  vs. LC_C28' LC_C4 -LC_C28;
contrast 'LC_C4  vs. LC_A56' LC_C4 -LC_A56;
contrast 'LC_A7  vs. LC_C7 ' LC_A7 -LC_C7 ;
contrast 'LC_A7  vs. LC_C28' LC_A7 -LC_C28;
contrast 'LC_A7  vs. LC_A56' LC_A7 -LC_A56;
contrast 'LC_C7  vs. LC_C28' LC_C7 -LC_C28;
contrast 'LC_C7  vs. LC_A56' LC_C7 -LC_A56;
contrast 'LC_C28 vs. LC_A56' LC_C28-LC_A56;

run;