Displaying 1 result from an estimated 1 matches for "expqq".
Did you mean:
exp
2004 Mar 08
0
Graphical Test
...advance for your help.
:-)
Jei
> ### Graphical test of suitability of the
exponential and gamma model for the aircondit data
>
>
> library(boot)
>
> data(aircondit)
>
>
>
> # Testing if the aircondit data comes from an
exponential distribution
>
> expqq.fun <- function(data, q) sort(data)/mean(data)
> exp.gen <- function(data, mle) rexp(length(data),
mle)
> n <- nrow(aircondit)
> qq <- qexp((1:n)/(n+1))
> exp.boot <- boot(aircondit$hours, expqq.fun,R=999,
sim="parametric",
+ ran.gen=exp.gen, mle=1/mean(airc...