Description Usage Arguments Value Author(s) References Examples
The projected normal distribution provides a flexible distribution for circular data, e.g., asymmetry and possible bimodality.
1 2 |
x |
a vector. The |
n |
number of observations. |
mu |
the mean vector of the bivariate normal. |
sigma |
the 2x2 variance and covariance matrix of the bivariate normal. |
log |
logical. If |
control.circular |
the attribute of the resulting object. |
dpnorm
gives the density, rpnorm
generates random deviates.
Claudio Agostinelli
S.R. Jammalamadaka and A. SenGupta (2001). Topics in Circular Statistics, Section 2.2.4, World Scientific Press, Singapore. K.V. Mardia (1972). Statistics of Directional Data. Academic Press. London and New York. F. Wang and A.E. Gelfand (2013). Directional data analysis under the general projected normal distribution. Stat Methodol. 10(1):113-127. doi:10.1016/j.stamet.2012.07.005.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 | data1 <- rpnorm(100, mu=c(0,0), sigma=diag(2), control.circular=list(units="degrees")) # Uniform on the circle
plot(data1)
ff <- function(x) dpnorm(x, mu=c(0,0), sigma=diag(2)) # Uniform on the circle
curve.circular(ff, join=TRUE,
main="Density of a Projected Normal Distribution \n mu=(0,0), sigma=diag(2)")
ff <- function(x) dpnorm(x, mu=c(1,1), sigma=diag(2)) # Unimodal
curve.circular(ff, join=TRUE, xlim=c(-1, 2.3),
main="Density of a Projected Normal Distribution \n mu=(1,1), sigma=diag(2)")
sigma <- matrix(c(1,0.9,0.9,1), nrow=2)
ff <- function(x) dpnorm(x, mu=c(0.5,0.5), sigma=sigma) # Bimodal
curve.circular(ff, join=TRUE, xlim=c(-1, 2.3),
main="Density of a Projected Normal Distribution \n mu=(0.5,0.5), rho=0.9")
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.