Heller_pdproject_dec-2

Heller title: “PD Project Dec 2, 2024_1” output: html_document — #libraries

library(ggplot2)
library(dplyr)
library(knitr)

#theme

my_theme<-function(x){theme_bw()+
    theme(text = element_text(size=20))+
    theme(axis.line.y = element_line(size = 2.0))+
    theme(axis.line.x = element_line(size = 2.0))+
    theme(axis.ticks = element_line(size = 1.5,colour="black"))+
    theme(axis.ticks.length=  unit(0.45, "cm"))+
    theme(axis.title.y =element_text(vjust=1.2))+
    theme(axis.title.x =element_text(vjust=-0.2))+
    theme(axis.text=element_text(colour="black"))+
    theme(panel.background = element_rect(fill ="white"))}

#40 mg dose

#sim40<-read.table("/work/ac0837/PHAR7380/week11/project/q2w40mg.res",skip=1,header=T)
sim40<-read.table("C:\\Heller\\PHAR7380\\pdproject2\\q2w40mg.res",skip=1,header=T)
#PK Profile
ggplot(data=sim40%>%filter(CMT==2),aes(TIME/7,PRED))+
geom_line(size=0.5)+
theme_bw()+
my_theme()+
labs(x="Time after dose (week)",y="UNT763 Concentration (mg/L)")+
    scale_x_continuous(breaks = c(0,2,4,6,8,10,12,14,16),limits = c(0,16))

## Warning: Using `size` aesthetic for lines was deprecated in ggplot2 3.4.0.
## i Please use `linewidth` instead.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.

## Warning: The `size` argument of `element_line()` is deprecated as of ggplot2 3.4.0.
## i Please use the `linewidth` argument instead.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.

#PD Profile
ggplot(data=sim40%>%filter(CMT==3),aes(TIME/7,PRED))+
geom_line(size=0.5)+
theme_bw()+
my_theme()+
labs(x="Time after dose (week)",y="CRP (mg/L)")+
geom_hline(yintercept=8, linetype="dashed", color = "red")+
scale_x_continuous(breaks = c(0,2,4,6,8,10,12,14,16),limits = c(0,16))

#80 mg dose

#sim80<-read.table("/work/ac0837/PHAR7380/week11/project/q2w80mg.res",skip=1,header=T)
sim80<-read.table("C:\\Heller\\PHAR7380\\pdproject2\\ld80mg.res",skip=1,header=T)
#PK Profile
ggplot(data=sim80%>%filter(CMT==2),aes(TIME/7,PRED))+
geom_line(size=0.5)+
theme_bw()+
my_theme()+
labs(x="Time after dose (week)",y="UNT763 Concentration (mg/L)")+
    scale_x_continuous(breaks = c(0,2,4,6,8,10,12,14,16),limits = c(0,16))

#PD Profile
ggplot(data=sim80%>%filter(CMT==3),aes(TIME/7,PRED))+
geom_line(size=0.5)+
theme_bw()+
my_theme()+
labs(x="Time after dose (week)",y="CRP (mg/L)")+
geom_hline(yintercept=8, linetype="dashed", color = "red")+
scale_x_continuous(breaks = c(0,2,4,6,8,10,12,14,16),limits = c(0,16))

#80 mg loading dose, then q2w 40 mg dose md

#sim8040<-read.table("/work/ac0837/PHAR7380/week11/project/ld80mgq2w40mgmd.res",skip=1,header=T)
sim8040<-read.table("C:\\Heller\\PHAR7380\\pdproject2\\ld80mgq2w40mgmd.res",skip=1,header=T)
#PK Profile
ggplot(data=sim8040%>%filter(CMT==2),aes(TIME/7,PRED))+
geom_line(size=0.5)+
theme_bw()+
my_theme()+
labs(x="Time after dose (week)",y="UNT763 Concentration (mg/L)")+
    scale_x_continuous(breaks = c(0,2,4,6,8,10,12,14,16),limits = c(0,16))

#PD Profile
ggplot(data=sim8040%>%filter(CMT==3),aes(TIME/7,PRED))+
geom_line(size=0.5)+
theme_bw()+
my_theme()+
labs(x="Time after dose (week)",y="CRP (mg/L)")+
geom_hline(yintercept=8, linetype="dashed", color = "red")+
scale_x_continuous(breaks = c(0,2,4,6,8,10,12,14,16),limits = c(0,16))

#160 mg loadig dose, then q2w 40 mg dose md

#sim16040<-read.table("/work/ac0837/PHAR7380/week11/project/ld160mgq2w40mgmd.res",skip=1,header=T)
sim16040<-read.table("C:\\Heller\\PHAR7380\\pdproject2\\ld160mgq2w40mgmd.res",skip=1,header=T)
#PK Profile
ggplot(data=sim16040%>%filter(CMT==2),aes(TIME/7,PRED))+
geom_line(size=0.5)+
theme_bw()+
my_theme()+
labs(x="Time after dose (week)",y="UNT763 Concentration (mg/L)")+
    scale_x_continuous(breaks = c(0,2,4,6,8,10,12,14,16),limits = c(0,16))

#PD Profile
ggplot(data=sim16040%>%filter(CMT==3),aes(TIME/7,PRED))+
geom_line(size=0.5)+
theme_bw()+
my_theme()+
labs(x="Time after dose (week)",y="CRP (mg/L)")+
geom_hline(yintercept=8, linetype="dashed", color = "red")+
scale_x_continuous(breaks = c(0,2,4,6,8,10,12,14,16),limits = c(0,16))