Contour plot of the skew skew-normal distribution in S^2
Contour plot of the skew skew-normal distribution in S^2.
skewnorm.contour(x, type = "alr", n = 100, appear = TRUE)
x |
A matrix with the compositional data. It has to be a 3 column matrix. |
type |
This is either "alr" or "ilr", corresponding to the additive and the isometric log-ratio transformation respectively. |
n |
The number of grid points to consider over which the density is calculated. |
appear |
Should the available data appear on the ternary plot (TRUE) or not (FALSE)? |
The alr or the ilr transformation is applied to the compositional data at first. Then for a grid of points within the 2-dimensional simplex the bivariate skew skew-normal density is calculated and the contours are plotted along with the points.
A ternary diagram with the points (if appear = TRUE) and the bivariate skew skew-normal contour lines.
Michail Tsagris
R implementation and documentation: Michail Tsagris mtsagris@uoc.gr and Giorgos Athineou <gioathineou@gmail.com>
Azzalini A. and Valle A. D. (1996). The multivariate skew-skewnormal distribution. Biometrika 83(4):715-726.
Aitchison J. (1986). The statistical analysis of compositional data. Chapman \& Hall.
x <- as.matrix(iris[51:100, 1:3]) x <- x / rowSums(x) skewnorm.contour(x)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.