cancel
Showing results for 
Show  only  | Search instead for 
Did you mean: 
Check out the JMP® Marketplace featured Capability Explorer add-in
Choose Language Hide Translation Bar
KernDawg
Level I

Mixed Model Trouble

What is the issue with this data? I can't get any output for the fixed effects tests.

I am performing a mixed model of a randomized complete block design. Rep is a random factor while S Rate and N Rate are continuous variables. I would like to see if S Rate or N Rate impact yield. Yield is also a continuous variable.  

 

I can get output when I run it as a standard least squares with the reml method, but when I change it to mixed model I don't get any numbers for the fixed effects test. 

3 REPLIES 3
Phil_Kay
Staff

Re: Mixed Model Trouble

Hi,

I am not sure what is happening here. In JMP Pro 16.2.0 I get the same problem as you have. By chance I ran the report in an early adopter version of JMP Pro 17 (EA6) and in that case it gave estimates for the fixed effects. I suggest that you contact technical support (email: support@jmp.com) so they can investigate.

All of that aside, it appears that your factors have little or no effect on the response. I suspect that is part of the reason for the problem that you are seeing.

I hope that helps,

Phil

Re: Mixed Model Trouble

Hello @KernDawg,  The convergence criteria and optimization algorithm for Fit Mixed and Fit Least Squares (with a random effect specified) are going to be different; which I suspect is also part of the problem here.  Another issue with this dataset is the covariance structure (high confounding risk), which can cause convergence problems in my experience: 

 PatrickGiuliano_1-1646714315406.png

Note that as a rule of thumb, we might consider a Pearson correlation of ~ 0.3 to indicate moderate risk, with ~0.5 or higher to indicate higher risk (but this will depend on subject-matter-knowledge and your own personal risk-definition).  

 

The Fit Mixed platform may be more sensitive to this confounding risk and hence, what appears to be a convergence problem (where effect estimates are not reported).  Agree with @Phil_Kay that low signal to noise ratio in this model may contribute to the problem we are seeing.

 

The Fit Least Squares personality with a random effect specified looks like a perfectly acceptable model specification for your modeling scenario here; in which case, none of your effects are statistically significant on your Yield Actual response (except for the Intercept term).   

 

PatrickGiuliano_2-1646714518880.png

 

The fact that only the intercept is "significant" is likely another clue.  You might want to look at a better-fitting model here, or consider why your expected model is not fitting your data particularly well (Adj-R-squared = 6% and is less than 1/2 of R-squared = 13%): 

 

PatrickGiuliano_3-1646714932424.png

 

-@PatrickGiuliano 

 

Re: Mixed Model Trouble

@KernDawg: I'd like to add one more thing here as food for thought. If we pay attention to the data structure in the "S Rate" and the "N Rate" variables, we can see that they are both rather discrete (not strictly continuous), even if they are modeled in your data table this way.  You can see this in the pattern of the data as shown in the normal quantile plots.  This type of pattern is often indicative of limited measurement resolution (in the context of a measurement system used to collect data).

 

PatrickGiuliano_0-1646800397855.png

PatrickGiuliano_4-1646800749555.png

 

If we fit an exponential distribution to these data (using the Distribution Platform and selecting this distribution fit for the purposes of demonstration and because Johnson Sb risks overfitting to the data), we can ask JMP to save a simulated continuous column for these fits for each of these two variables: 

 

PatrickGiuliano_5-1646801443362.png

 

PatrickGiuliano_2-1646800567368.png

 

When we do this, and we try fitting the same model in Fit Mixed (using these variables instead of the original ones), the model converges and we successfully obtain Fixed Effects Estimates: 

PatrickGiuliano_3-1646800653348.png