Hello - I am estimating a zero truncated Poisson model in PROC NLMIXED I would like to get the LS Means of the interaction of group*time, but am unsure how to do this or if it is possible. If the model was not zero truncated, I could use PROC GLIMMIX. For example: proc glimmix data=mydata; class id time group var1; model y=time group var1 group*time group*var1 group*time*var1/link=log dist=poisson; random _residual_/subject=id; lsmeans group*time/ilink pdiff; lsmeans group*time*var1/ilink pdiff; run; Is it possible to get similar output in NLMIXED using the Estimate statements? I've also seen the NLMeans macro, but I had a hard time deciphering how to output the data needed for the macro. proc nlmixed data=mydata; log_lambda = intercept +b_var1*var1 + b_time1*time1+ b_time2*time2+b_group*group+ b_grouptime1*grouptime1+ b_grouptime2*grouptime2+ b_grouptime1var1*grouptime1var1+ b_grouptime2var1*grouptime2var1+u; lambda = exp(log_lambda); ll = y*log_lambda - lambda - log(1-exp(-lambda)) - lgamma(y+1); model y ~ general(ll); random u ~ normal(0,s2u) subject=id; run; Thanks.
... View more
