Dear Victor,
breedR's function remlf90() actually wraps two different REML algorithms (AI and EM), depending on the 'method' argument.
EM is guaranteed to converge, by construction. AI, instead, is typically faster, but more prone to divergences.
One practical way of checking convergence is thus, to check whether both methods yield similar results. Although this is not always possible, as some breedR models only work with EM.
You can check the evolution of the likelihood by inspecting the contents of the element 'reml.out' in the object resulting from the remlf90 call. This contains the output from the algorithm in text form and includes many details, including the log-likelihood at each step. You can parse this output to extract the values for inspection.
Best,
ƒacu.-
--
Report issues at https://github.com/famuvie/breedR/issues
---
You received this message because you are subscribed to the Google Groups "breedR" group.
To unsubscribe from this group and stop receiving emails from it, send an email to breedr+un...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/breedr/dcf9a59e-72e5-43eb-ae41-94afc70528e0n%40googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/breedr/2836b273-a1ae-e0a6-2a3b-d9576238ecbb%40cirad.fr.