I am trying to replicate a SAS call
PROC GLM DATA=MAD.MDA_plot_controls outstat=MAD.plot_control_anova_stat;
BY Year Location Trait;
CLASS Row Column;
MODEL Value = Row Column/SS3;
RANDOM Row Column;
LSMEANS Row Column;
RUN;
my first try ...which I am sure isnt correct is the following
x4<-GLM(Value ~ factor(Row) + factor(Column), MDA_plot_controls, eps=1e-8)
the parts of the SAS command I am not sure I am replicating are
BY Year Location Trait
.
RANDOM Row Column;
LSMEANS Row Column;
any ideas ?