PRsurface {adaptIVPT}  R Documentation 
This function plots the power (passingrate) curve and power (passingrate) surface of the mixed scaling (MS) approach. A power curve shows the statistical power across different effect sizes. In IVPT studies, the effect size is captured by the difference between the means of logmeasurements of the test and reference products (i.e., logGMR). For the passingrate surface, the corresponding function considers different values of the standard deviation.
PRsurface(
n,
r,
observed_GMR = 0.95,
observed_sigmaWR = 0.294,
GMR_grid = seq(0.75, 1.3, length.out = 100),
sigmaWR_grid = seq(0.2, 1, length.out = 100),
params = list(),
nsim = 1000,
ncores = NULL,
verbose = FALSE,
plot = TRUE
)
n 
The number of donors in each simulation. 
r 
The number of replicates from each donor for each simulated dataset. 
observed_GMR 
The observed (estimated) GMR of the user's data. Along with the observed sigmaWR, the corresponding passing rate will be displayed in the 3D plot as a vertical line parallel to the zaxis. 
observed_sigmaWR 
The observed (estimated) sigmaWR of the user's data. Along with the observed GMR, the corresponding passing rate will be displayed in the 3D plot as a vertical line parallel to the zaxis. 
GMR_grid 
The grid of GMR values to be used for plotting the 3D surface of passing rates. 
sigmaWR_grid 
The grid of sigmaWR values to be used for plotting the 3D surface of passing rates. 
params 
(Optional) The list of true parameters to be assumed in data generation.

nsim 
(Optional) The number of total simulations to be conducted. Defaults to 1,000. 
ncores 
(Optional) The number of CPU cores to use for parallel processing (OpenMP). If R hasn't been installed with OpenMP configured, this will not take effect. When OpenMP is available, it should not exceed the number of existing cores. If unspecified, it will default to 2 cores or the number of existing cores, whichever is smaller. 
verbose 
(Optional) A logical value ( 
plot 
(Optional) A logical value ( 
A list
GMR
 A list of true parameter settings.
passing_rate
 The estimated passing rate.
runtime
 The total elapsed time charged for the execution of the program.
Daeyoung Lim, daeyoung.lim@uconn.edu
Davit, B. M., Chen, M. L., Conner, D. P., Haidar, S. H., Kim, S., Lee, C. H., Lionberger, R. A., Makhlouf, F. T., Nwakama, P. E., Patel, D. T., Schuirmann, D. J., & Yu, L. X. (2012). Implementation of a referencescaled average bioequivalence approach for highly variable generic drug products by the US Food and Drug Administration. The AAPS journal, 14(4), 915924.
out < PRsurface(6, 3, GMR_grid = c(0.90, 1), sigmaWR_grid = c(0.2, 0.5), nsim = 2, plot = FALSE)