inst/examples/01_simple_linuplogdn/05_Ind_Log_Bolus/run.R

In Quantmetrics770/IQnca: Noncompartmental Analysis in R

library(IQnca)
rm(list=ls())

# -------------------------------------------------------------------------#
# Import data ----
# -------------------------------------------------------------------------#
data <- as.data.frame(Indometh)
data$Subject <- as.numeric(as.character(data$Subject))
names(data) <- c("USUBJID","NTIME","ACONC")
data$DOSE <- 25 # Same dose as in reference used
data$STUDYID <- "Indometh R dataset"
data$COMPOUND <- "Indometh"
data$ANALYTE <- "Indometh"
data$MATRIX <- "Plasma"
data$PROFILE <- "Single dose profile"
data$PROFTYPE <- "SD"
data$GROUP <- "25 mg SD"
data$GROUPN <- "25"
data$GROUPU <- "mg"
data$DAY <- 1
data$ATIME <- NA
data$TIMEUNIT <- "Hours"
data$CONCUNIT <- "ug/mL"
data$LLOQ <- 0
data$ADM <- "Bolus"
data$DOSEUNIT <- "mg"


# -------------------------------------------------------------------------#
# Create IQdataNCA dataset ----
# -------------------------------------------------------------------------#

cat("Indometh Log Bolus\n")

dataNCA <- IQdataNCA(data = data,FLAGTIME = "nominal",AUCMETHD = "LinearUp LogDown")

resNCA <- nca_IQdataNCA(data = dataNCA)

resNCA