bow TIE for NTIDs [Power / Sample Size]
Hi Astea,
Another function under revision… Starting with 1.4.6.9000 subject simulations are performed which take substantially longer than with the previous code. The approach used in ≤1.4.6 (based on the 4-dimensional t-distribution) was flawed. 
In the new version with the argument
Welcome to the club!
Where did you find that? 90.00–111.11% for AUC0–t of sirolimus, AUC0–72 of tacrolimus and everolimus.
I’m not willing to accept 111.11%. Heck, with an acceptable ∆ of 10% we get \(100(1-\Delta)^{-1}=111.1\dot{1}\) – nothing else! I will not eradicate my second-grade math only cause this number was proclaimed ex cathedra as truth by the omniscient almighty oracle. Had to swallow already rounding of the CI. Double rounding? Gimme a break!
❝ I tried this:
❝
❝ type1error.2TOST(CV=c(0.3,0.25), n=122, theta1=c(0.8, 0.9), theta2=c(1.25, 1.11), rho=0.9, details = FALSE)
❝ [1] 0.05000015
❝ Warning:
❝ Result of maximization over nullset may not be reliable.
![[image]](img/uploaded/image4.gif)

In the new version with the argument
details=TRUE
showing the TIE for all eight intersection null sets:type1error.2TOST(CV=c(0.3, 0.25), n=122, theta1=c(0.8, 0.9),
theta2=c(1.25, 1/0.9), rho=0.9, details=TRUE)
1e+06 simulations. Time consumed (secs)
user system elapsed
165.97 25.16 191.57
Intersection null P(Type I Error) theta0 #1 theta0 #2
1 H_A01 n H_Ca 0.049639 0.8000000 0.9494219
2 H_A02 n H_Ca 0.049800 1.2500000 1.0630689
3 H_Aa n H_C01 0.050149 0.9948235 0.9000000
4 H_Aa n H_C02 0.050005 0.9775390 1.1111111
5 H_A01 n H_C01 0.031504 0.8000000 0.9000000
6 H_A01 n H_C02 0.000000 0.8000000 1.1111111
7 H_A02 n H_C01 0.000000 1.2500000 0.9000000
8 H_A02 n H_C02 0.031817 1.2500000 1.1111111
❝ The best statisticians I've ever met told me they are not actually statisticians. Neither am I...
Welcome to the club!
❝ P.S. 111.00 instead of logical 111.11 because it was written in the mentioned guideline.
Where did you find that? 90.00–111.11% for AUC0–t of sirolimus, AUC0–72 of tacrolimus and everolimus.
❝ Increadible rounding methods!
I’m not willing to accept 111.11%. Heck, with an acceptable ∆ of 10% we get \(100(1-\Delta)^{-1}=111.1\dot{1}\) – nothing else! I will not eradicate my second-grade math only cause this number was proclaimed ex cathedra as truth by the omniscient almighty oracle. Had to swallow already rounding of the CI. Double rounding? Gimme a break!
—
Dif-tor heh smusma 🖖🏼 Довге життя Україна!![[image]](https://static.bebac.at/pics/Blue_and_yellow_ribbon_UA.png)
Helmut Schütz
![[image]](https://static.bebac.at/img/CC by.png)
The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes
Dif-tor heh smusma 🖖🏼 Довге життя Україна!
![[image]](https://static.bebac.at/pics/Blue_and_yellow_ribbon_UA.png)
Helmut Schütz
![[image]](https://static.bebac.at/img/CC by.png)
The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes
Complete thread:
- Sample size for 4-period 4-sequence crossover BE study Bryony Simmons 2018-02-01 12:16 [Power / Sample Size]
- function sampleN.TOST of package PowerTOST Helmut 2018-02-01 13:01
- Alpha adjustment in higher order crossover d_labes 2018-02-01 14:01
- Deficiencies Helmut 2018-02-01 15:48
- Deficiencies nobody 2018-02-01 17:10
- Deficiencies d_labes 2018-02-01 18:57
- Deficiencies Relaxation 2018-02-02 11:12
- Deficiencies nobody 2018-02-02 12:49
- Deficiencies Helmut 2018-02-02 16:14
- Deficiencies Relaxation 2018-02-02 19:41
- alpha... where is omega? Astea 2018-02-02 21:45
- α and no ω Helmut 2018-02-02 23:39
- TIE for NTIDs d_labes 2018-02-04 12:40
- TIE for NTIDs Astea 2018-02-04 20:04
- TIE for NTIDs Helmut 2018-02-05 01:01
- TIE for NTIDs d_labes 2018-02-05 16:40
- TIE for NTIDs Helmut 2018-02-05 17:49
- TIE for NTIDs d_labes 2018-02-05 22:17
- TIE for NTIDs Helmut 2018-02-06 12:34
- TIE for NTIDs d_labes 2018-02-05 22:17
- TIE for NTIDs Helmut 2018-02-05 17:49
- TIE for NTIDs d_labes 2018-02-05 16:40
- TIE for NTIDs d_labes 2018-02-05 16:35
- bow TIE for NTIDs Astea 2018-02-05 17:52
- bow TIE for NTIDsHelmut 2018-02-05 18:10
- 111.11 for NTIDs Astea 2018-02-05 19:27
- 111.11 for NTIDs Helmut 2018-02-06 00:12
- 111.11 for NTIDs Astea 2018-02-05 19:27
- bow TIE for NTIDs d_labes 2018-02-05 22:33
- bow TIE for NTIDsHelmut 2018-02-05 18:10
- bow TIE for NTIDs Astea 2018-02-05 17:52
- TIE for NTIDs Helmut 2018-02-05 01:01
- TIE for NTIDs Astea 2018-02-04 20:04
- TIE for NTIDs d_labes 2018-02-04 12:40
- α and no ω Helmut 2018-02-02 23:39
- alpha... where is omega? Astea 2018-02-02 21:45
- Deficiencies Relaxation 2018-02-02 19:41
- Deficiencies Relaxation 2018-02-02 11:12
- Deficiencies Helmut 2018-02-01 15:48
- Alpha adjustment in higher order crossover d_labes 2018-02-01 14:01
- function sampleN.TOST of package PowerTOST Helmut 2018-02-01 13:01