UUPharmacometrics / xpose4

Model building aid for population analysis using NONMEM
http://uupharmacometrics.github.io/xpose4/
GNU Lesser General Public License v3.0
33 stars 11 forks source link

Upper limit of quantification plots with Xpose.vpc.both #21

Open SimonKoele opened 2 years ago

SimonKoele commented 2 years ago

Dear all,

Recently I was making a model of data that had LLOQ and ULOQ data present. Using the Xpose.vpc.both function I generated VPC plots, however the categorical subplot shows no data when plotting ULOQ data (see example below). image

The function works fine when using LLOQ data. I investigated the code of the function and found that the following part of the code still selects the LLOQ output table when ULOQ is present/selected:

for (LEVS in 1:n.levs) { ret.new[(1 + (LEVS - 1) n.idv.levs):(n.idv.levs LEVS), 1:5] <- tmp.table.2[, c(idv.index, real.index[[LEVS]], lower.index[[LEVS]], median.index[[LEVS]], upper.index[[LEVS]])] ret.new[(1 + (LEVS - 1) n.idv.levs):(n.idv.levs LEVS), "by.var"] <- level.names[LEVS] }

Therefore it presents the user with an empty table. For now I have created a new function that selects the right table columns to subset from to generate the plots, which does work. I'm using xpose4_4.7.0 in R version 3.6.2. Does anyone else have this problem?

Kind regards, Simon Koele Radboudumc Pharmacometrics

andrewhooker commented 2 years ago

Hi Simon,

Thanks for the report. I’ll make the change and push it out later in a few days.

Best regards, Andy

Andrew Hooker, Ph.D. Professor of Pharmacometrics Dept. of Pharmacy Uppsala University Box 580, 751 23, Uppsala, Sweden Phone: +46 18 471 4355 https://katalog.uu.se/profile/?id=N4-631  https://katalog.uu.se/profile/?id=N4-631 http://www.farmbio.uu.se/research/researchgroups/pharmacometrics/

On 8 Dec 2021, at 11:16, SimonKoele @.***> wrote:

Dear all,

Recently I was making a model of data that had LLOQ and ULOQ data present. Using the Xpose.vpc.both function I generated VPC plots, however the categorical subplot shows no data when plotting ULOQ data (see example below). https://user-images.githubusercontent.com/93375358/145190802-baa3bb74-573c-4d86-add4-3ee2c6b78cd2.png The function works fine when using LLOQ data. I investigated the code of the function and found that the following part of the code still selects the LLOQ output table when ULOQ is present/selected:

for (LEVS in 1:n.levs) { ret.new[(1 + (LEVS - 1) n.idv.levs):(n.idv.levs LEVS), 1:5] <- tmp.table.2[, c(idv.index, real.index[[LEVS]], lower.index[[LEVS]], median.index[[LEVS]], upper.index[[LEVS]])] ret.new[(1 + (LEVS - 1) n.idv.levs):(n.idv.levs LEVS), "by.var"] <- level.names[LEVS] }

Therefore it presents the user with an empty table. For now I have created a new function that selects the right table columns to subset from to generate the plots, which does work. I'm using xpose4_4.7.0 in R version 3.6.2. Does anyone else have this problem?

Kind regards, Simon Koele Radboudumc Pharmacometrics

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/UUPharmacometrics/xpose4/issues/21, or unsubscribe https://github.com/notifications/unsubscribe-auth/ABJOMZCBTXXI4S2RRINWF3TUP4WBPANCNFSM5JTNJO6A. Triage notifications on the go with GitHub Mobile for iOS https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Android https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.

mariewijk commented 2 years ago

Dear Andrew,

I am having the same issue as Simon with trying to generate a VPC using xpose.vpc.both with ULOQ data. It gives the same empty ULOQ panel shown above.

Did you ever tend to this matter?

Best regards, Marie Wijk