Bioequivalence and Bioavailability Forum

Dear Helmut,

❝ Sometimes people don’t know which version they (or their IT guys) installed. Maybe it’s safer to replace

❝ ...

❝ by

❝     if(packageVersion("PowerTOST") < "1.1.13" & des == "paired") {


❝       pwrN <- power.TOST(CV=CV, n=Ns[j], theta0=GMR, design=des)
❝     } else {
❝       pwrN <- power2.TOST(CV=CV, n=n, theta0=GMR, design=des)
❝     }

This is only necessary for a function definition outside package. Within the next release of PowerTOST, which will introduce the sensitivity analysis, the correct version of power2.TOST() is of course available.

❝ Most input errors are already trapped by PowerTOST. At the beginning of Sens() add:

❝   if(pwr.min < 0.05) stop("Minimum acceptable power must not be <0.05!")

❝ … otherwise uniroot() throws an error.

Done. Do you know the real reason behind that error?

❝ Consider adding ... to transfer optional parameters (theta1,…).

Very good idea . But ran into trouble doing so with some inexplicable (for me) errors. Seems have to study how to work with the argument '...'.

❝

---

Sneak Preview

Funny .
Especially the deviation for the CV!
Seems we don't have to care for any deviation here. In the first look surprising.

---

Will send you the code of the package version of your function Sens() soon.
Some adaption of the argument names seems necessary to me to be in accordance with the naming within PowerTOST up to know (f.i. pwr.target=targetpower, pwr.min=minpower, GMR=theta0 ...). Although I must confess that I find your naming more appealing.

In the output you name the achieved power 'expected' power. But this term is already exhausted for the output of exppower.TOST().

I would suggest the function name to be something like psa.ABE() (power sensitivity analysis for ABE) or psa.TOST() to account for the upcoming things seen in your Sneak Preview.