poped - evaluate_power error in example 11

[TheoPapath]

evaluate_power gives an error

This happens even though parameter number 5 is not fixed - likely a poped error and nothing can be done

Also in this example, the pediatric FIM does not match 1:1 with the poped examples. Not sure why that is.

[mattfidler]

Hi @TheoPapath

The reason the pediatric FIM does not match is because you mis-translated the model using:

    KA<-tKa*exp(eta.ka)* (pedCL^(isPediatric)) # add covariate for pediatrics

Once you correct to:

    CL<-tCl*exp(eta.cl)* (pedCL^(isPediatric)) # add covariate for pediatrics

The FIM matches. I have updated the example in your branch (you will need to pull the changes).

Also in the babelmixr2 translation the model does not have a parameter 5, so this is the reason this does not exist. If you update this to index 4 it works just fine.

However, this is really not that user-friendly, so I have added the function:

babelBpopIdx(babel.db.all, "pedCL")

Which will query the babel.db.all for the "pedCL" parameter and give the correct index. This is in #136

[mattfidler]

Thank you for working this example, I was unaware of the power of estimating covariates in PopED, this is something I am interested in 😄