Proposal: Simulation from PK/PD and systems pharmacology models in R with mrgsolve.

[kylebaron]

ISoP_mrgsolve_intro_13May2016.zip

• General Topic: Introduction to the mrgsolve R package
• YOUTUBE Live Feed Link: http://www.youtube.com/watch?v=ts3Lj3hVyiw
• Potential Presenter(s): Kyle Baron @kylebmetrum
• Potential Date(s): 13 May 2016
• Time: 12:00 PM EDT
• Background Resources for the Audience: https://github.com/metrumresearchgroup/mrgsolve https://mrgsolve.github.io/user_guide
• Get the modmrg package here https://github.com/metrumresearchgroup/modmrg
• YouTube Video: https://youtu.be/ts3Lj3hVyiw

• R code and .html files attached as .zip ISoP_mrgsolve_intro_code_13May2016.zip ISoP_mrgsolve_intro_documents_13May2016.zip

  archive

[jacksonan]

I have everything loaded correctly except the compilers.

Where on Github can I post this question to get clarified instructions for downloading compilers?

[kylebaron]

You can post questions here: https://github.com/metrumresearchgroup/mrgsolve/issues

The installation guide is here: https://github.com/metrumresearchgroup/mrgsolve/wiki/mrgsolve-Installation

[Anathawa]

Adding mod <- mread("pk1cmt", modlib()) to the description? Loving the video presentation. :)

This is not an important issue, but I get a result that differs from the video, despite trying to do the exact same steps. :) Not sure this q really fits under issues undermrgsolve?

I realized that the problem is that the video actually shows the output from the code in the supplied code, and not that shown in the video:

In the documentation: data <- expand.ev(amt=1000, ii=24, addl=100, rate = 100, cmt=2, CL = seq(0.5,2.5,0.25))

In the video:

data <- expand.ev(ID=1,amt=c(1000), cmt=2, ii=24, addl=100, CL = seq(0.5,2.5,0.25)) data %>% mutate(rate=amt/10)

The rate gets defined by both rate=100, and mutate(rate=amt/10), it gets clicked but is edited out, so the plot doesn't correspond to the visible code. Just leaving it as a comment if anyone else find it worrying. :)

[kylebaron]

Thanks @Anathawa for the report. The modmrg package was deprecated and we just rolled all of those models into the mrgsolve package itself and they can be used with the same mread call that you'd use for any other model; just make project argument point to the value returned by modlib(). You can see a listing of the models with modlib(list=TRUE) or ?modlib help topic.

Apologies for the discrepancy in the coding. I was torn between working from a pre-coded script (less stressful) and coding on the fly (more stressful). I might have mixed it up a bit there.