‘Deep’ bug [Software]

posted by Helmut Homepage – Vienna, Austria, 2018-01-05 01:37 (2275 d 07:47 ago) – Posting: # 18147
Views: 10,853

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.1824
35   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 🖖🏼 Довге життя Україна! [image]
Helmut Schütz
[image]

The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes

Complete thread:

UA Flag
Activity
 Admin contact
22,957 posts in 4,819 threads, 1,636 registered users;
97 visitors (0 registered, 97 guests [including 8 identified bots]).
Forum time: 09:24 CET (Europe/Vienna)

With four parameters I can fit an elephant,
and with five I can make him wiggle his trunk.    John von Neumann

The Bioequivalence and Bioavailability Forum is hosted by
BEBAC Ing. Helmut Schütz
HTML5