Hello all, I'm trying to code plant measurement data (i.e. % cover, health rating 0-5, etc). into glimmix, however, I'm unsure how to go about it. I've had no issues doing soils data. Here is my code for my soils data: data soildata; set soildata;
if horizon="US" then delete;
if horizon="LS" then delete;
run;
proc print data=soildata;
RUN;
proc glimmix data=soildata;
title "TS EC";
Class rxn; /* rxn is my treatment prescription */
Model ec=rxn/ddfm=satterthwaite;
lsmeans rxn/ adjust=tukey plot=mean lines;
run; What do I need to add to deal with non continuous values? I know I need to include Dist=multinomial, but is there anything else? I suspect I also need a link function, but don't understand what those are, or what they do. Any help here would be great.
... View more