Hi, My question deals with results in Glimmix that do not show all LSMEANS with my data. Here's my code Data Cards; input nozzle$ program$ sprayer$ @@; input yield@; output; Cards; ...... PROC Glimmix; class nozzle program sprayer; model yield = nozzle|program|sprayer / DDFM=KR; lsmeans nozzle|program|sprayer / pdiff cl adjust=sidak; lsmeans nozzle|program|sprayer / diff lines ADJUST=sidak; *ods output diffs=ppp lsmeans=mmm; *ods exclude diffs lsmeans; run; I have 5 nozzles, 3 programs, 2 sprayers and four plots with yield data for each combination. This means with four replicates, I have 120 data points, and 8 data points for the untreated plots. An example would look like XR B T 120 - Where XR is my nozzle, B is my application program, T is my sprayer type, and 120 is my yield. My issue is that my untreated plots (four for each sprayer) do not have a nozzle type or program, so the LSMEANS is not showing for comparisons of nozzle type or program. What can I do to make my code allow for keeping untreated plot data but be able to see results for LSMEANS of program and nozzle type? I appreciate your help and look forward to hearing back soon! Thanks!!!
... View more