On Jul 21, 2010, at 1:18 PM, Knut Krueger wrote:

David Winsemius schrieb:


And exactly why did you think I offered

?Chisquare
I was completely on the wrong way, and tried to find a solution with the formula instead to substitute the formula.

So I tried to implement pchisq into the formula - and of course I got wrong values ...

> p <- function(x) ((0.5^(5/2))/gamma(5/2))*(x^((5/2)-1))* (2.718281828459^(-x/2))
> integrate(p, 25.508, Inf)
0.0001111 with absolute error < 2.7e-05

Check result>>>>

> pchisq(25.508,5)
[1] 0.9999
> 1-pchisq(25.508,5)
[1] 0.0001111


Thank's Knut

David Winsemius, MD
West Hartford, CT

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to