Bioequivalence and Bioavailability Forum • ‘Deep’ bug

‘Deep’ bug [Software]

posted by Helmut – Vienna, Austria, 2018-01-05 01:37 (2301 d 20:59 ago) – Posting: # 18147
Views: 10,958

Hi Mittyri,

❝ As I see the problem is in quantile function, looks like approximation is not correct for low DF.

Wow!

R:

for (df in 1:7) cat(sprintf("%9.6f", qt(0.025, df, lower.tail=FALSE)), "\n")

12.706205

 4.302653

 3.182446

 2.776445

 2.570582

 2.446912

 2.364624

Spreadsheets:

M$-Excel 2000   OO Calc 4.1.1   Gnumeric 1.12.15

  12.706150       12.706205         12.706205

   4.302656        4.302653          4.302653

   3.182449        3.182446          3.182446

   2.776451        2.776445          2.776445

   2.570578        2.570582          2.570582

   2.446914        2.446912          2.446912

   2.364623        2.364624          2.364624

tinv(x, y) was introduced as a custom function in PHX/WNL v6.4. Hence, no checking possible in earlier versions. However:

   v6.4       v7.0       v8.0  

12.706205  12.706205  12.706205

 4.302653   4.302653   4.302653

 3.182435   3.182435   3.182435

 2.776445   2.776445   2.776445

 2.570582   2.570582   2.570582

 2.446912   2.446912   2.446912

 2.364624   2.364624   2.364624

Since in my example both the estimate and its standard error are correct, something must be wrong ‘behind the curtains’ in PHX/WNL when calculating the CI.

—
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:

assessment of steady state in WinNonLin pkpdpkpd 2007-04-04 22:52
- 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 (example) apapanas 2018-01-03 14:44
    - steady state in WinNonlin (example) Helmut 2018-01-04 17:11
      - qt in WinNonlin mittyri 2018-01-04 23:13
        
        ‘Deep’ bugHelmut 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