R Graphical Manual

Browse All

Last data update: 2014.03.03

R: Contour plot of the Kent distribution without any data.

kent.datacontour

R Documentation

Contour plot of the Kent distribution without any data.

Description

The contour plot of the Kent distribution on the sphere for some data is produced.

Usage

kent.datacontour(x)

Arguments

`x`	A two column matrix, where the first column is the latitude and the second comlumn is the longitude. If the matrix has two columns, it is assumed to have unit vectors and in this case it is turned into latitude and longitude.

Details

MLE of the parameters of the Kent distribution are calculated, then the contour plot is plotted using these estimates and finally the data are also plotted.

Value

A plot containing the contours of the distribution along with the data.

Author(s)

Michail Tsagris

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

References

Kent John (1982). The Fisher-Bingham distribution on the sphere. Journal of the Royal Statistical Society, Series B, 44(1): 71-80.

Examples

x <- rvmf(100, rnorm(3), 10)
kent.mle(x)
y <- euclid.inv(x)
kent.datacontour(y)
vmf.kerncontour(y, thumb = "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/kent.datacontour.Rd_%03d_medium.png", width=480, height=480)
> ### Name: kent.datacontour
> ### Title: Contour plot of the Kent distribution without any data.
> ### Aliases: kent.datacontour
> ### Keywords: Contour plot Kent distribution
> 
> ### ** Examples
> 
> x <- rvmf(100, rnorm(3), 10)
> kent.mle(x)
$G
           mean      major      minor
[1,] -0.8503029 -0.3818506 -0.3621811
[2,]  0.4902817 -0.8249060 -0.2813430
[3,] -0.1913344 -0.4167976  0.8886343

$para
     kappa       beta 
10.3943405  0.3287205 

$logcon
[1] 9.893739

$loglik
[1] -50.10158

$runtime
   user  system elapsed 
  0.012   0.004   0.015 

> y <- euclid.inv(x)
> kent.datacontour(y)
> vmf.kerncontour(y, thumb = "none")
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>