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.
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
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 🖖🏼 Довге життя Україна!
Helmut Schütz
The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes
Dif-tor heh smusma 🖖🏼 Довге життя Україна!
Helmut Schütz
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
- 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