Monolix

Data Sets

Info

Many different data sets were downloaded from Monolix's website
Structure of downloaded data sets follow rules with one exception (see below)

Original Form

Bolus
- DATASET_2CPT\\01_BOLUS_2CPT_VCL\\bolus_2cpt_VCL_data.txt
Infuse
- DATASET_2CPT\\03_INF_2CPT_VCL\\Infusion_2cpt_VCL_SS_data.txt
Oral
- DATASET_2CPT\\08_ORAL1_2CPT_KAVVMKM\\oral1_2cpt_KAVVMKM_data.txt

Massage Data

Remove '#' in the very first line before importing
Replace dataset["AMT"]=="." with 0

Final Form

Download Bolus(source file)

> head(bolus)
  ID TIME     DV   LNDV MDV    AMT EVID   DOSE    V1   CL    Q    V2
1  1 0.00     NA     NA   1 120000    1 120000 101.5 3.57 6.99 59.19
2  1 0.25 1040.7 6.9476   0      0    0 120000 101.5 3.57 6.99 59.19
3  1 0.50 1629.0 7.3957   0      0    0 120000 101.5 3.57 6.99 59.19
4  1 0.75  877.8 6.7774   0      0    0 120000 101.5 3.57 6.99 59.19
5  1 1.00 1247.2 7.1286   0      0    0 120000 101.5 3.57 6.99 59.19
6  1 1.50 1225.1 7.1107   0      0    0 120000 101.5 3.57 6.99 59.19

Download Infusion(source file)

> head(infuse)
  ID   TIME     DV   LNDV MDV   AMT  RATE EVID  DOSE   V1 CL    Q    V2 SS TAU TMSS
1  1 216.00     NA     NA   1 60000 60000    1 60000 67.2  4 5.38 54.27  1  24 0.00
2  1 216.25  489.2 6.1927   0     0    NA    0 60000 67.2  4 5.38 54.27 NA  NA 0.25
3  1 216.50  900.8 6.8032   0     0    NA    0 60000 67.2  4 5.38 54.27 NA  NA 0.50
4  1 216.75  999.7 6.9074   0     0    NA    0 60000 67.2  4 5.38 54.27 NA  NA 0.75
5  1 217.00 1360.5 7.2156   0     0    NA    0 60000 67.2  4 5.38 54.27 NA  NA 1.00
6  1 217.50 1154.3 7.0512   0     0    NA    0 60000 67.2  4 5.38 54.27 NA  NA 1.50

Download Oral(source file)

> head(oral)
  ID TIME    DV   LNDV MDV   AMT EVID  DOSE    V1    VM    KM    Q    V2    KA
1  1 0.00    NA     NA   1 20000    1 20000 94.18 606.5 249.8 6.84 77.56 0.884
2  1 0.25  38.6 3.6525   0     0    0 20000 94.18 606.5 249.8 6.84 77.56 0.884
3  1 0.50  77.2 4.3466   0     0    0 20000 94.18 606.5 249.8 6.84 77.56 0.884
4  1 0.75  99.6 4.6009   0     0    0 20000 94.18 606.5 249.8 6.84 77.56 0.884
5  1 1.00 115.4 4.7486   0     0    0 20000 94.18 606.5 249.8 6.84 77.56 0.884
6  1 1.50 124.8 4.8269   0     0    0 20000 94.18 606.5 249.8 6.84 77.56 0.884

Code

Code to alter each file

make.numeric <- function(data){
	column.names <- colnames(data)
	print(head(data))
	for( i in 1:ncol(data)){
#		cat("\nColumn = ",column.names[i],"\n")
#		print(head(data[,i]))
		data[,i] <- as.numeric(data[,i])
	}
	return(data)
} 

cat("\nbolus_tcpt_VCL_data\n")

bolus <- read.table(file.choose(),header=TRUE,stringsAsFactors = FALSE)

bolus[bolus["AMT"]==".","AMT"] <- 0
bolus[bolus=="."] <- NA
bolus <- make.numeric(bolus)


cat("\ninfuse\n")

infuse <- read.table(file.choose(),header=TRUE,stringsAsFactors = FALSE)

infuse[infuse["AMT"]==".","AMT"] <- 0
infuse[infuse=="."] <- NA
infuse <- make.numeric(infuse)


cat("\noral\n")

oral <- read.table(file.choose(),header=TRUE,stringsAsFactors = FALSE)

oral[oral["AMT"]==".","AMT"] <- 0
oral[oral=="."] <- NA
oral <- make.numeric(oral)