mmds (version 1.1)

plot.ds.mixture: Plot the results of a mixture model detection function fit.

Description

Plots the detection function (or pdf) of a fitted mixture model detection function, optionally overlayed on a histogram of the observed data.

Usage

"plot" (x, style = "", main = "", breaks = "Sturges", ylim = NULL, xlim = NULL, pdf = FALSE, plot.formula = NULL, hide.hist = FALSE, nomf = FALSE, x.axis = NULL, xlab = "Distance", ylab = NULL, ...)

Arguments

x

a ds.mixture object.

style

If set to "comp", composite plots of the detection function will be shown for the detection function (the averaged detection function in the covariate case).

main

(A vector of) title(s) for the plot(s). By default these are set by the function (and are fairly ugly but descriptive).

breaks

Breaks to be used for the histogram. This can be a vector of numbers or any of the permissable options used in hist. Defaults to "Sturges".

ylim

Used to manually set the y limit of the plot. Defaults to NULL.

xlim

Used to manually set the x limit of the plot. Defaults to NULL.

pdf

Should the pdf be plotted rather than detection function be plotted? Only really useful with point transect data. Defaults to FALSE.

plot.formula

Formula of covariates to be plotted. Defaults to NULL, which plots all covariates. No effect with non-covariate models.

hide.hist

Should the histogram be hidden, leaving only the detection function (or pdf) to be plotted? Defaults to FALSE.

nomf

Should the mfrow value be altered? Useful when creating custom plots for publication. Defaults to FALSE (yes, change the mfrow value).

x.axis

Set the x axis labels. Again, useful for publication plots. Defaults to NULL, which uses the default R values.

xlab

Label for the x axis.

ylab

Label for the y axis.

...

not used at the moment

Value

a plot!

Details

For covariate models, all the levels of factor variables are plotted or the 25, 50 and 75th percentiles of continuous variables are plotted averaged over the values of the other covariates.

Examples

library(mmds)
set.seed(0)
## simulate some line transect data from a 2 point mixture
sim.dat<-sim.mix(c(-0.223,-1.897,inv.reparam.pi(0.3)),2,100,1)
## fit the model
fit.sim.dat<-fitmix(sim.dat,1,2)
## plot
plot(fit.sim.dat)