summary.quax | R Documentation |
This function is the summary
method for class quax
objects as returned
by quax
.
## S3 method for class 'quax'
summary(object, ...)
object |
The function returned by |
... |
not in use here |
The value
component of the result can be used to compare the
quality of the fit of different dispersal kernels for the same quantile
to the same data.
A list with the following components:
coefficients
The parameters of the estimated dispersal kernel.
value
The attained value of the objective function that is minimised in the quantile regression.
## Prepare artificial data:
set.seed(0)
r <- rgamma(200, shape=2, scale=150)
simulated.data <- data.frame(distance = r, density =
rpois(length(r), k_lognormal(r, par=c(6,0), N=1000000, d=2)))
plot(density ~ distance, simulated.data)
## Fit a log-normal and a power-law dispersal kernel to the data:
f1 <- quax(density ~ distance, simulated.data,
tau = 0.9, fun = k_lognormal)
f2 <- quax(density ~ distance, simulated.data,
tau = 0.9, fun = k_power)
## Compare both fits:
summary(f1)
summary(f2)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.