Model-agnostic protocol
Use a per-subject likelihood or a batched program with automatic differentiation.
Bring a conditional model; estimate random effects and population parameters.
NLMEEstimation.jl defines estimation methods around a conditional likelihood supplied by a model. The package supplies the Gaussian random-effect prior and its covariance parameterization. Start with the model protocol, then compare the available estimators in the PK example.