[R] qt with ncp>37.62

data at tulipany.cz data at tulipany.cz
Sat Jun 14 22:44:02 CEST 2008


help(qt) states that: 
"ncp 	non-centrality parameter delta; currently except for rt(), only for
abs(ncp) <= 37.62"

so I would expect that calling qt with non-centrality parameter exceeding
37.62 should fail, instead e.g. calling

> mapply(function(x) qt(p = 0.9, df = 55, ncp = x),35:45)

gives:

 [1] 40.21448 41.35293 42.49164 43.68862 44.82945 45.97048 47.11170 48.25310
 [9] 49.39467 50.53639 51.67826
Warning messages:
1: In qt(p = 0.9, df = 55, ncp = x) :
  full precision was not achieved in 'pnt'
2: In qt(p = 0.9, df = 55, ncp = x) :
  full precision was not achieved in 'pnt'
3: In qt(p = 0.9, df = 55, ncp = x) :
  full precision was not achieved in 'pnt'

so it seems calculation for (according to what is written in documentation)
allowed values of ncp, i.e. in this case 35,36 and 37 is done and precision is
checked, whereas calculation for the rest may be completely incorrect?
Or was there any update of code (in pnt.c ?), allowing calculation of pt with
higher ncp, not followed by documentation update?

Thank you,
Nikola Kaspříková



More information about the R-help mailing list