qt in WinNonlin [Software]
Hi Antigoni and Helmut,
That's really weird.
As I see the problem is in quantile function, looks like approximation is not correct for low DF.
R:
WNL 502:
For compliance the same solution could be implemented as R did:
For central qt, a C translation of
Hill, G. W. (1970) Algorithm 396: Student's t-quantiles. Communications of the ACM, 13(10), 619–620.
altered to take account of
Hill, G. W. (1981) Remark on Algorithm 396, ACM Transactions on Mathematical Software, 7, 250–1.
That's really weird.
As I see the problem is in quantile function, looks like approximation is not correct for low DF.
R:
qt(0.025, 1:7)
[1] -12.706205 -4.302653 -3.182446 -2.776445 -2.570582
[6] -2.446912 -2.364624
WNL 502:
-16.50068
-4.327589
-3.18348
-2.7764075
-2.5705446
-2.4469136
-2.3646432
For compliance the same solution could be implemented as R did:
For central qt, a C translation of
Hill, G. W. (1970) Algorithm 396: Student's t-quantiles. Communications of the ACM, 13(10), 619–620.
altered to take account of
Hill, G. W. (1981) Remark on Algorithm 396, ACM Transactions on Mathematical Software, 7, 250–1.
—
Kind regards,
Mittyri
Kind regards,
Mittyri
Complete thread:
- assessment of steady state in WinNonLin pkpdpkpd 2007-04-04 22:52 [Software]
- steady state in WinNonlin (example) Helmut 2007-04-05 02:30
- steady state in WinNonlin (example) pkpdpkpd 2007-04-05 16:33
- steady state in WinNonlin (example) Helmut 2007-04-05 17:16
- steady state in WinNonlin (example) pkpdpkpd 2007-04-05 17:39
- steady state in WinNonlin… Helmut 2007-04-05 18:40
- steady state in WinNonlin… pkpdpkpd 2007-04-05 21:37
- Paired Design in WinNonlin Helmut 2007-04-06 00:23
- steady state in WinNonlin… Helmut 2007-04-05 18:40
- steady state in WinNonlin (example) pkpdpkpd 2007-04-05 17:39
- steady state in WinNonlin (example) Helmut 2007-04-05 17:16
- steady state in WinNonlin (example) apapanas 2018-01-03 14:44
- steady state in WinNonlin (example) Helmut 2018-01-04 17:11
- qt in WinNonlinmittyri 2018-01-04 23:13
- ‘Deep’ bug Helmut 2018-01-05 00:37
- Workarounds for v6.3+ Helmut 2018-01-08 00:40
- Workarounds for v6.3+ apapanas 2018-01-10 12:47
- qt in WinNonlinmittyri 2018-01-04 23:13
- steady state in WinNonlin (example) Helmut 2018-01-04 17:11
- steady state in WinNonlin (example) pkpdpkpd 2007-04-05 16:33
- steady state in WinNonlin (example) Helmut 2007-04-05 02:30