Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

questions about mrgsim.sa #4

Open
kevin8897 opened this issue Jun 7, 2023 · 2 comments
Open

questions about mrgsim.sa #4

kevin8897 opened this issue Jun 7, 2023 · 2 comments

Comments

@kevin8897
Copy link

what dose "eps" mean in the code, parameter change value extent or a ratio?
how the sensitivity indices was calculated, and how the parameters changed in default?

Examples
mod <- mrgsolve::house(delta=0.1)

par <- "CL,VC,KA"

var <- "CP"

dose <- ev(amt = 100)

fun <- function(mod, ...) mrgsolve::mrgsim_e(mod, dose, output="df")

out <- lsa(mod, par, var, fun)

@kylebaron
Copy link
Owner

Yeah ... it's a fractional change in the parameter value.

The source code is here: https://github.com/kylebaron/mrgsim.sa/blob/main/R/lsa.R#L85C1-L86

I need to write a vignette about this or get it in the documentation so you don't have to read the code for this.

@kevin8897
Copy link
Author

Yeah ... it's a fractional change in the parameter value.

The source code is here: https://github.com/kylebaron/mrgsim.sa/blob/main/R/lsa.R#L85C1-L86

I need to write a vignette about this or get it in the documentation so you don't have to read the code for this.

Hi, kylebaron, thanks a lot for your reply.
so the default eps value is 1e-08?
and based on the output, dose lsa only run a single side analysis?
in the source code, dose dvalue ="sens" in the output file?
dvalue <- function(sim,ref,scale) {
(sim[["dv_value"]] - ref[["dv_value"]])/scale}, what is scale? is it baseline ref[["dv_value"]]?

sorry for so many questiones, looking forward to the vignette, thanks!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants