\name{Contour plot (on the sphere) of a mixture of elliptically symmetric distributions}
\alias{mixsespc.contour}
\alias{mixesag.contour}
\title{
Contour plot (on the sphere) of a mixture of elliptically symmetric  distributions
}

\description{
The contour plot (on the sphere) of a mixture of elliptically symmetric  distributions.
}

\usage{
mixsespc.contour(probs, mu, theta, bgcol = "snow", dat = NULL, col = NULL,
lat = 50, long = 50)
mixesag.contour(probs, mu, gam, bgcol = "snow", dat = NULL, col = NULL,
lat = 50, long = 50)
}

\arguments{
\item{probs}{
This is avector with the mixing probability of each group.
}
\item{mu}{
A matrix with the mean direction of each group.
}
\item{theta}{
A matrix with the \eqn{\theta} parameters of the SESPC distribution of each group.
}
\item{gam}{
A matrix with the \eqn{\gamma} parameters of the ESAG distribution of each group.
}
\item{bgcol}{
The color of the surface of the sphere.
}
\item{dat}{
If you have you want to plot supply them here. This has to be a numerical matrix with three columns,
i.e. unit vectors.
}
\item{col}{
If you supplied data then choose the color of the points. If you did not choose a color,
the points will appear in red.
}
\item{lat}{
A positive number determing the range of degrees to move left and right from the latitude center.
See the example to better understand this argument.
}
\item{long}{
A positive number determing the range of degrees to move up and down from the longitude center.
See the example to better understand this argument.
}
}

\details{
The goal of this function is for the user to see how the mixtures of SESPC or ESAG distributions look like.
}

\value{
A plot containing the contours of the mixture distribution.
}

\references{
Perdikis T., Alharbi N. and Tsagris M. (2026). Model--based clustering for spherical
and hyper--spherical data using elliptically symmetric distributions.

https://arxiv.org/abs/2605.27496
}

\author{
Michail Tsagris.

R implementation and documentation: Michail Tsagris \email{mtsagris@uoc.gr}.
}

%\note{
%%  ~~further notes~~
%}

\seealso{
\code{\link{dmixsespc}, \link{mixsespc.mle}
}
}

\examples{
\donttest{
probs <- c(0.5, 0.5)
mu <- matrix(rnorm(6, 10, 5), ncol = 3)
theta <- matrix(c(-2, 1, -2, 1), ncol = 2)
x <- rmixsespc(200, probs, mu, theta)$x
b <- dmixsespc(x, probs, mu, theta)
## the lat and long are decreased to 10. Increase them back to 50 to
## see the difference
mixsespc.contour(probs, mu, theta, lat = 10, long = 10)
}
}

