I am using SAS 9.4. I have a factorial design for my experiment to try and establish common milkweed into a landscape using two main effects of: mowed plots vs unmowed plots, each of the main plots is divided into three sub-plots using three different rates of herbicide to suppress the grass in the landscape: a control rate, a low rate, and a high rate of glyphosate, and within the subplots I have four sub-subplots each with different milkweed seed densities: a control, a low density, a low density + mid-june mowing, and a high density. I am wanting to analyze the presence of common milkweed seedlings based seedlings being present or absent only. It runs fine up through the proc print data code; however once I try to run the Proc glimmix portion of the code it does not converge and I am needing help with trying to fix this so it converges. I attached my data excel file if that is useful in trying to help me.
Here is the code I have been trying to run:
* Import the data;
proc import out = seedling_counts
datafile = '/folders/myfolders/AES Consulting/Lizotte-Hall, Sydney/Data/Aug2016counts.xlsx'
dbms = xlsx
replace;
getnames = yes;
datarow = 2;
run;
* Create a absense/presense variable for seedlings;
data seedling_counts;
set seedling_counts;
if seedling_counts > 0 then seedling_present = 1;
else seedling_present = 0;
run;
* Print the data;
proc print data = seedling_counts;
run;
proc glimmix data = seedling_counts plots = residualpanel;
class mowing herb seed rep;
model seedling_present = mowing herb seed rep / dist = binary link = logit;
random rep*mowing*herb;
lsmeans herb / diff;
run;
