Title  h=1.5 "Likelihood Function, Germination Example"; 
goptions reset=all; 

Data Likelihood;  array temp(15); 
do i=1 to 15; input temp(i) @@; end; 
  do A = -3 to 6 by .2; 
    Do B = -.3 to 2.5 by .05; 
	  L=1; 
*yes;	Do i= 1 to 9; line = A+B*(temp(i)-70); L=L*exp(line)/(1+exp(line)); end; 
*no ;   Do i = 10 to 15; line =A+B*(temp(i)-70); L=L/(1+exp(line)); end; 
      Log_L = log(L); keep A B L log_L;
	output; end; end; 
datalines;  
	     64 67   70 71 73      77 78 82 85    
50 58 63    67 62         75
 ; 
 proc g3d data=Likelihood; plot A*B=L; plot A*B=log_L/rotate = 10; 
 run;