MixedModels residuals and post-hocs

Carlos Familia

unread,

Jun 17, 2016, 5:26:35 AM6/17/16

to julia-stats

Hello,

I was wondering if it is possible to take the residuals from a lmm model and if there is any lsmeans package for post-hocs in julia?

Best regards,

Carlos Família

Douglas Bates

unread,

Jun 17, 2016, 4:52:05 PM6/17/16

to julia-stats

On Friday, June 17, 2016 at 4:26:35 AM UTC-5, Carlos Familia wrote:

Hello,

I was wondering if it is possible to take the residuals from a lmm model and if there is any lsmeans package for post-hocs in julia?

There is no lsmeans package that I know of and it is highly unlikely that I will ever write one, but someone else may.

As for residuals, like many other properties of mixed-effects models, you need to be careful how you define them. You need to decide if the random effects are included in the fitted values or not. Perhaps it is better to define a fitted method for LinearMixedModel, which can then be the basis of a residuals method.

Could you open an issue in the MixedModels repository? It is probably better to have this discussion there.

The reason that I say it is unlikely I will write any lsmeans code is because lsmeans are a bizarre SAS construction from the days when we stored data on magnetic tape and programs on punched card decks. It would take so long to get a single model fit that users wanted every possible inference or diagnostic value calculated from the model, whether or not they made sense.

Carlos Familia

unread,

Jun 29, 2016, 1:53:49 PM6/29/16

to julia-stats

Thank you for your answer, and sorry for the late reply.

I will open an issue in the MixedModels repository.

Carlos Familia

unread,

Jul 1, 2016, 7:56:43 AM7/1/16

to julia-stats

So what options are there for evaluating the significance of the main effects, and for diagnostics within the julia's MixedModels package?

On Friday, June 17, 2016 at 9:52:05 PM UTC+1, Douglas Bates wrote:

Reply all

Reply to author

Forward