cinid.plot {CINID}

Print the observed distribution and the estimated one

Description

This function computes a histogram of the observed headcapsule widths and the Gaussian curves of the probability densities of the four larval instars, defined formerly by the cinid.table function (you must use cinid.table before to use cinid.plot).

Usage

cinid.plot(cinid.out, breaks = 50, xlab = "Headcapsules width", ylab1 = "Density", 
ylab2 = "Number of Larvae", main = "", ...)

Arguments

cinid.out

the returned two-sized list by cinid.table function.

breaks

the breakpoints or the number of cells for the histogramm. See 'breaks' in 'hist'. See also 'Details'.

xlab

the X axis label. The default value is 'Headcapsule width'.

ylab1

the left Y axis label. The default value is 'Density'.

ylab2

the right Y axis label. The default value is 'Number of Larvae'.

main

the main title. The field is empty by default.

...

further arguments and graphical parameters passed to 'hist'

Details

The default for 'breaks' is '50' because it is usually the best way to cluster HCW data.

Value

An object of class 'histogram'. See 'hist'.

Author(s)

Aurelie Siberchicot [e-mail], Adrien Merville, Marie-Claude Bel-Venner and Samuel Venner

References

Merville A., Vallier A., Venner S., Siberchicot A., Fouchet D., Heddi A. and Bel-Venner M.-C. (2014). Determining the instar of a weevil larva (Coleoptera:Curculionidae) using a parsimonious method. European Journal of Entomology. 111(4): 567-573. doi: 10.14411/eje.2014.056

See Also

'hist'

Examples

data(s.oryzae)
soryzae <- s.oryzae$data
mu4.soryzae <- s.oryzae$mu4
sd4.soryzae <- s.oryzae$sd4

HCW.soryzae <- cinid.table(soryzae, mu4 = mu4.soryzae, sd4 = sd4.soryzae)
HCW.soryzae.plot <- cinid.plot(HCW.soryzae, breaks = 50)