R script published pubmed site not work [🇷 for BE/BA]

posted by Weidson – Brazil, 2021-12-22 19:47 (827 d 03:21 ago) – Posting: # 22709
Views: 2,150

Dear Friends,

I would like to understand why the code bellow is not working. This code I was take out on the page: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3817376/. Would it be possible for any BEBAC member to adapt it to work? I need of the a graphics tool (for example: ggplot2) that to make graphics dinamic way as we can make at Phoenix WinNonlin.

(ggplot2)
library(nlme)
head(Theoph)

ggplot(data=Theoph, aes(x=Time, y=conc, group=Subject)) + geom_line() + labs(x=“Time (hr)”, y=“Concentration (mg/L)”)

p <- ggplot(data=Theoph, aes(x=Time, y=conc, group=Subject)) + geom_line() + labs(x=“Time (hr)”, y=“Concentration (mg/L)”) + stat_summary(fun.y=median, geom=“line”,aes(x=ntpd, y=conc, group=1), color=“red”, size=1)

print(p) # “p” is a ggplot object

# create a flag for body weight

Theoph$WT <- ifelse(Theoph$Wt<70, “WT < 70kg”, “WT >= 70kg”)
p + facet_grid(.~WT)""t>


In true, I need of a dinamic code (for example: ggplot2) for to plot a graphic of Individual plasma concentration profiles sorted by treatment and period with whih the mean. Example:

[image]

Complete thread:

UA Flag
Activity
 Admin contact
22,957 posts in 4,819 threads, 1,636 registered users;
72 visitors (0 registered, 72 guests [including 11 identified bots]).
Forum time: 23:08 CET (Europe/Vienna)

Nothing shows a lack of mathematical education more
than an overly precise calculation.    Carl Friedrich Gauß

The Bioequivalence and Bioavailability Forum is hosted by
BEBAC Ing. Helmut Schütz
HTML5