VPC question

________________________________ From: [email protected] [mailto:[email protected]] On Behalf Of Liu, Rong (BARDS) Sent: Wednesday, September 05, 2012 3:13 PM To: [email protected] Subject: [NMusers] VPC question Dear NM users, I have a basic VPC question. I run the following code to generate VPC plot using cluster: vpc run63.mod -lst=run63.lst -samples=1000 -idv=LAVES PlotVPC 7 I got the following error message. Reading /pkpd_ms/M4444/PPEN/npc_dir7/vpc_results.csv Error in xpose.VPC(vpc.file, vpctab = vpctab, PI = "both", PI.up.lty = 0, : trying to get slot "Data" from an object of a basic class ("NULL") with no slots Execution halted Could anybody kindly tell me what went wrong? Thank you very much, Rong Liu Associate Principle Scientist at Merck & Co., Inc Notice: This e-mail message, together with any attachments, contains information of Merck & Co., Inc. (One Merck Drive, Whitehouse Station, New Jersey, USA 08889), and/or its affiliates Direct contact information for affiliates is available at http://www.merck.com/contact/contacts.html) that may be confidential, proprietary copyrighted and/or legally privileged. It is intended solely for the use of the individual or entity named on this message. If you are not the intended recipient, and have received this message in error, please notify us immediately by reply e-mail and then delete it from your system. Notice: This e-mail message, together with any attachments, contains information of Merck & Co., Inc. (One Merck Drive, Whitehouse Station, New Jersey, USA 08889), and/or its affiliates Direct contact information for affiliates is available at http://www.merck.com/contact/contacts.html) that may be confidential, proprietary copyrighted and/or legally privileged. It is intended solely for the use of the individual or entity named on this message. If you are not the intended recipient, and have received this message in error, please notify us immediately by reply e-mail and then delete it from your system.

From: Norman Z Sent: Tuesday, October 09, 2012 7:30 AM To: [email protected] Subject: [NMusers] VPC question Dear NMusers, I am using PsN to conduct VPC. A “vpc_results.csv” was generated by PsN, where the percentiles of real and simulated data were reported. Generally, these percentiles plus the observed data were used for the VPC plot. My question is how to calculate number of observation outside the 90% prediction interval (PI). I am thinking about two way to do this. First, I may use the NPC results provided by the “vpc_results.csv” file, e.g. NPC results bin 3: 15 observations points below PI (count) points below PI (%) ... 0% PI 7 46.66667 ... 40% PI 5 33.33333 ... 80% PI 1 6.66667 ... 90% PI 0 0 ... 95% PI 0 0 ... Should the numbers in this table be used to obtain the observation outside the 90% PI? If so, how can they be used? If not, shall I calculate the observation outside the 90% PIs of each bin using R? Second, I can calculate the 90% PIs for each time point (e.g. 1.1h, 1.3h, etc) from the raw simulated data and compare the 90% PIs with the observation at one particular time point. Would the second strategy be more accurate than calculate the observation outside the 90% PIs of each binned time interval (e.g. 1-2 h, 2-3 h)? If this is the case, where can I find the raw simulated data? (I suspect the “DV_matrix.csv” and “vpc_simulation.1.npctab.dta” in “m1” folder are the raw simulated data, but want to refer to relevant document to understand how these files are organized.) Thanks for your help. Norman

From: [email protected] [mailto:[email protected]] On Behalf Of Norman Z Sent: 09 October 2012 16:31 To: [email protected] Subject: [NMusers] VPC question Dear NMusers, I am using PsN to conduct VPC. A "vpc_results.csv" was generated by PsN, where the percentiles of real and simulated data were reported. Generally, these percentiles plus the observed data were used for the VPC plot. My question is how to calculate number of observation outside the 90% prediction interval (PI). I am thinking about two way to do this. First, I may use the NPC results provided by the "vpc_results.csv" file, e.g. NPC results bin 3: 15 observations points below PI (count) points below PI (%) ... 0% PI 7 46.66667 ... 40% PI 5 33.33333 ... 80% PI 1 6.66667 ... 90% PI 0 0 ... 95% PI 0 0 ... Should the numbers in this table be used to obtain the observation outside the 90% PI? If so, how can they be used? If not, shall I calculate the observation outside the 90% PIs of each bin using R? Second, I can calculate the 90% PIs for each time point (e.g. 1.1h, 1.3h, etc) from the raw simulated data and compare the 90% PIs with the observation at one particular time point. Would the second strategy be more accurate than calculate the observation outside the 90% PIs of each binned time interval (e.g. 1-2 h, 2-3 h)? If this is the case, where can I find the raw simulated data? (I suspect the "DV_matrix.csv" and "vpc_simulation.1.npctab.dta" in "m1" folder are the raw simulated data, but want to refer to relevant document to understand how these files are organized.) Thanks for your help. Norman