[Rd] R stat functions do not work as stated on the mannual ( PR#7419)

Liaw, Andy andy_liaw at merck.com
Tue Dec 14 12:41:26 CET 2004


> From: casadoj at ecc.es
> 
> 
> Dear R Developers:
> 
> I have been playing with R, release 2.0.1 for a week now and 
> have detected =
> that all stat functions related to distribution probabilities 
> have the same=
>  problem:
> 
> 1.- According to the manual the log.p parameter is always the 
> last one.
> 2.- When you use the software, the last parameter seems to be 
> lower.tail
> 
> Example:
> 
> > pt (1.1, 5)
> [1] 0.8392746
> > pt (1.1, 5, F)
> [1] 0.8392746
> > pt (1.1, 5, F, T)
> [1] 0.8392746
> >=0D
> 
> On this example, I have used the Student T distribution. The 
> result of this=
>  example has been tested with the stat calculator at 
http://calculators.sta=
t.ucla.edu/.
1.- This first line shows that when this function has two arguments, the lo=
g.p default value is false, and the lower.tail is true.
2.- The second line shows that when this function has 3 arguments, the last=
 one is the log.p argument and the lower.tail is taken by default to true.
3.- The third line shows that when this function has 4 arguments, the third=
 one keeps on being the log.p argument and the lower.tail is the last one.
4.- Acording to the mannual, the lower.tail should be the third argument an=
d not the last one.

Which manual did you read?  ?pt has:

Usage
[...]
pt(q, df, ncp=0, lower.tail = TRUE, log.p = FALSE)

So none of the calls you showed actually passed any thing for log.p to pt():
It's the 5th argument.  While it seems to be a convention that log.p is the
last argument, the actual position can change from distribution to
distribution as the number of parameters can differ.

Andy

Best regards,
=0D
Jos=E9 Luis Casado Mart=EDnez
------------------------------------------------------------------
European Computing Consultants
C/ Hermanos Garc=EDa Noblejas, N=BA 39, 5=AA, N 1
28037 Madrid
Telf.: 34-91-406 19 15. Fax: 34-91-406 19 16
Movil: 34-607-750 316
------------------------------------------------------------------

The last line shows that the sum of probabilities of 1.1 using 5 degrees of=
 freedom adding the lower tail (pt (1.1, 5, F, T)) and not adding the lower=
 tail ( pt (1.1, 5, F, F)) is 1. This is not the case with log (p).



_____________________________________________________________________
Mensaje analizado y protegido, tecnologia antivirus www.trendmicro.es
	[[alternative HTML version deleted]]

______________________________________________
R-devel at stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel



More information about the R-devel mailing list