note file PRE_0

switch b10
  case 1,2: note poisson or negbin so start with loge(Y)
            yvar b5
            chan 0 0.001 cb5 0.15 c1189
            switch b16
              case 0: note not mixed model
                 loge c1189 c1189
                 leave
              case 1: note mixed model
                 note mask is 1 if discrete, 0 if continuous
                 calc c1189=c1189*'mask'+1
                 calc c1189=loge(c1189)+cb5*(1-'mask')
                 leave
         
            ends
			note AL subtract the offset if it exists
            switch b20
              case 1: note offset exists
           note      calc c1189=c1189-'offs'
                 leave
            ends
   resp c1189
ends

note remove any bottom level distributional constraints
obey nobvar

nlev b4
loop b1 b15 b4
   offs b1
endloop