fitdistribution: Fits a leaf angle distribution
In LeafAngle: Analysis and Visualization of Plant Leaf Angle Distributions

Description Usage Arguments Details Value Author(s) References See Also Examples

Fits one of nine leaf angle distributions to data, using either log-likelihood (the preferred method), or a method based on minimizing the chi-squared statistic. This latter method is included to be consistent with Wang et al. (2007), as well as others.

A leaf angle distribution is a distribution valid for 0 <= x <= 90.

1 2	fitdistribution(angles, distribution, fitmethod = c("loglik", "chisq"), ellipsmethod = 2, degrees = TRUE, ...)

`angles`	Vector of leaf angles (0-90 degrees, or 0 - pi/2 radians)
`distribution`	Name of the distribution. Currently included are 'ellipsoid', 'rotatedell', 'twoparbeta', 'spherical', 'planophile', 'extremophile', 'erectophile', 'uniform', and 'plagiophile'.
`fitmethod`	Method to fit the distribution, either "loglik" (log-likelihood) or "chisq" (chi-squared statistic).
`ellipsmethod`	If distribution='ellipsoid', and ellipsmethod=1, the method reported by Wang et al. 2007 is used. Not recommended.
`degrees`	If TRUE, the default, the sample of angles is in degrees, otherwise radians.
`...`	Other parameters passed to `ftheta`.

See drawsample to simulate from a fitted leaf angle distribution, and fitalldistributions for a convenience function to fit all built-in distributions to one dataset.

Returns an object of class 'angledist'. Methods exist for print, summary, and plot. See Examples.

Remko Duursma

Wang W.M., Li Z.L. & Su H.B. (2007). Comparison of leaf angle distribution functions: Effects on extinction coefficient and fraction of sunlit foliage. Agricultural and Forest Meteorology, 143, 106-122.

fitdistribution, fitalldistributions

# Some leaf angle data:
data(eteret)

# Fit the ellipsoidal distribution:
f <- fitdistribution(eteret, "ellipsoid")
f

# Standard plot, histogram with fitted density:
plot(f)