NONMEM Users Network Archive

Hosted by Cognigen

[NMusers] Incomplete circles drawn in VPC plotted using R scripts generated by PsN

From: HUI, Ka Ho <matthew.hui_at_link.cuhk.edu.hk>
Date: Thu, 19 Nov 2015 14:56:27 +0000

Dear all,

I was trying to generate VPC using the R scripts generated by PsN with the =
following command line:

vpc
-samples=200
-stratify_on=...
-rplots=1
-idv=TIME
-auto_bin=unique
-predcorr
-varcorr
-dir=vpc_FINAL_MODEL
FINAL_MODEL_vpc.mod

The R script generated looks like this:

#START OF AUTO-GENERATED PREAMBLE, WILL BE OVERWRITTEN WHEN THIS FILE IS US=
ED AS A TEMPLATE
#Created 2015-11-19 at 22:05

rplots.level <- 1
xpose.runno <- '1_vpc'
toolname <- 'vpc'
pdf.filename <- paste0('PsN_',toolname,'_plots.pdf')
pdf.title <- 'vpc diagnostic plots run 1_vpc'
working.directory<-'C:/...'
raw.results.file <- 'raw_results_FINAL_MODEL1_vpc.csv'
model.directory<-'C:/...'
model.filename<-'FINAL_MODEL1_vpc.mod'
subset.variable<-NULL
mod.suffix <- '.mod'
mod.prefix <- 'FINAL_MODEL'
tab.suffix <- ''
tool.results.file <- 'vpc_results.csv'
theta.labels <- c('T1','T2','T3','T4','T5','T6')
theta.fixed <- c(FALSE,FALSE,TRUE,FALSE,FALSE,FALSE)
omega.labels <- c('O1','O2')
omega.fixed <- c(FALSE,FALSE)
sigma.labels <- c('S1')
sigma.fixed <- c(FALSE)
n.eta <- 2
n.eps <- 1

vpctab <- 'vpctab'
have.loq.data <- FALSE
have.censored <- FALSE
is.categorical <- FALSE
is.tte <- FALSE
dv <- 'DV'
idv <- 'TIME2'

setwd(working.directory)

###########################################################################=
#
#END OF AUTO-GENERATED PREAMBLE
#WHEN THIS FILE IS USED AS A TEMPLATE THIS LINE MUST LOOK EXACTLY LIKE THIS


library(xpose4)

pdf(file=pdf.filename,width=10,height=7,title=pdf.title)

done <- FALSE
if (is.tte){
     #data is in the model directory, go there to read input
        setwd(model.directory)
        xpdb <- xpose.data(xpose.runno)
        plots <- kaplan.plot(object=xpdb,VPC=T)
        #go back to vpc directory
        setwd(working.directory)
        done <- TRUE
}

if (is.categorical & (!done)){
    plots<-xpose.VPC.categorical(vpc.info=tool.results.file,vpctab=vpct=
ab)
        done <- TRUE
}

if ((have.loq.data | have.censored) & (!done) ){
    plots<-xpose.VPC.both(vpc.info=tool.results.file,vpctab=vpctab)
        done <- TRUE
}

if (!done){
        plots<-xpose.VPC(vpc.info=tool.results.file,vpctab=vpctab)
        done <- TRUE
}

print(plots)

dev.off()

However the "circles" (indicating raw Cp data) are incomplete and take the =
shape of a "crescent moon" instead. When run with the -lloq option, the 1st=
 stratification has complete circles but the remaining 2 also has the cresc=
ent-moon circles instead. Other parts of the plots (Shades and plots of CI)=
 look normal.
Can anyone suggest any possible cause and solution?

Thanks and regards,
Matthew Hui

Received on Thu Nov 19 2015 - 09:56:27 EST

The NONMEM Users Network is maintained by ICON plc. Requests to subscribe to the network should be sent to: nmusers-request_at_iconplc.com. Once subscribed, you may contribute to the discussion by emailing: nmusers@globomaxnm.com.