I am trying to run a repeated measures model in proc mixed with relative seed production as my dependent variable and year (the repeated measure), genotype, and size as independent variables for each plant in my study. The model runs well if I include a single interaction term such as year*size, but when I try to include year*genotype*size I get the "stopped because of infinite likelihood error." My code is below. I've tried using type =un, type= AR(1) and type = cs, but nothing works. Any suggestions?
class genotype plant year;
model rel_seed_prod = year genotype size year*genotype year*size genotype*size year*genotype*size;
repeated year / subject = plant type = un;
You might try initializing the covariance structure parameters according to the values which were estimated when you have only the year*size interaction in the model. Having reasonable initial estimates of the parameters of the covariance structure can help one avoid some estimation problems.