Re: Question on nlme

On 9/19/2025 2:36 PM, sasha a wrote:
> Hello,
> 
> I am using the nlme package in R:
> 
> model <- lme(
>      response ~ response_type + x_y + x_z - 1,
>      random = ~ response_type - 1 | group,
>      weights = varIdent(form = ~ 1 | response_type),
>      data = data_long,
>      control = lmeControl(opt = "optim")
> )
> 
> Is it possible to know what estimation process this is using? For example, is it using EM? Or is it use Newton-Raphson with GLS/RMLE?
> 
> Thanks!
> 
> 
> 	[[alternative HTML version deleted]]
> 
> ______________________________________________
> [email protected] mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide https://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
Hello,

I am not sure this answers the question but nlme::lme.formula 
optimization function defaults to ?optim, which defaults to method 
Nelder-Mead.
Users can change this via ?lmeControl, with two possible values for the 
method, ?nlminb and optim.

Hope this helps,

Rui Barradas