R Graphical Manual

Browse All

Last data update: 2014.03.03

R: Contour plot of spherical data using a von Mises-Fisher...

vmf.kerncontour

R Documentation

Contour plot of spherical data using a von Mises-Fisher kernel density estimate

Description

Contour plot of spherical data using a von Mises-Fisher kernel density estimate.

Usage

vmf.kerncontour(u, thumb = "none")

Arguments

`u`	A two column matrix. The first coolumn is the latitude and the second is the longitude.
`thumb`	This is either 'none' (defualt), or 'rot' for the rule of thumb suggested by Garcia-Portugues (2013). If it is "none" it is estimated via cross validation, with the fast function "vmfkde.tune_2".

Details

It calcculates the contour plot using a von Mises-Fisher kernel for spherical data only.

Value

A plot with the data in latitude and longitude along with their kernel density estimate contour lines.

Author(s)

Michail Tsagris

R implementation and documentation: Michail Tsagris <mtsagris@yahoo.gr> and Giorgos Athineou <athineou@csd.uoc.gr>

References

Garcia Portugues, E. (2013). Exact risk improvement of bandwidth selectors for kernel density estimation with directional data. Electronic Journal of Statistics, 7, 1655–1685.

Examples

x <- rvmf(100, rnorm(3), 15)
x <- euclid.inv(x)
par( mfrow = c(1, 2) )
vmf.kerncontour(x, "rot")
vmf.kerncontour(x, "none")

Results


R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(Directional)
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/Directional/vmf.kerncontour.Rd_%03d_medium.png", width=480, height=480)
> ### Name: vmf.kerncontour
> ### Title: Contour plot of spherical data using a von Mises-Fisher kernel
> ###   density estimate
> ### Aliases: vmf.kerncontour
> ### Keywords: Contour plot von Mises-Fisher kernel
> 
> ### ** Examples
> 
> x <- rvmf(100, rnorm(3), 15)
> x <- euclid.inv(x)
> par( mfrow = c(1, 2) )
> vmf.kerncontour(x, "rot")
> vmf.kerncontour(x, "none")
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>