Questions tagged [circular-statistics]
Directional statistics (also called circular or spherical statistics) is the discipline of statistics that deals with directions (unit vectors in $\mathbb{R}^n$), axes (lines through the origin in $\mathbb{R}^n$) or rotations in $\mathbb{R}^n$.
284
questions
0
votes
0
answers
16
views
GLMM with interaction terms between two circular predictor variables?
I am running a GLMM to see if several weather and nest box covariates influence occupancy (binary linear response). I would like to include two circular predictors (wind direction and box entrance ...
4
votes
1
answer
176
views
Use of circular predictor in GLMM
I am developing a mixed-effects binomial logistic regression (using glmmTMB, family = binomial) where the response is presence-absence. One of my potential predictors is hour of day, which takes ...
3
votes
1
answer
70
views
Reference for Directional Statistics of Plane Orientation
I've got a project I'm working on where I've got the orientation (normal) vectors of planes. These vectors are all within a unit hemisphere where the $z$-coordinate is strictly positive. The ...
0
votes
0
answers
15
views
is there a closed form of the CF (characteristic function) of a bivariate von Mises distribution?
is there a closed form of the CF (characteristic function) of a bivariate von Mises distribution?
And if I have two parameters that follow von Mises distribution, but my two parameters will be mixed ...
0
votes
0
answers
20
views
How to cluster time-of-day event data?
Suppose I have data about patient sneezes, e.g.,:
Name Time Loudness
Alex 07:59 10
Bob 08:03 12
Charlie 17:06 9
Alex 08:09 13
...
You can see ...
0
votes
0
answers
11
views
Characterize conditions in which Taylor moment approximation is good
I am working with the Projected Gaussian, or Angular Gaussian distribution, which is given by $z = \frac{x}{||x||}$, where $x \sim \mathcal{N}(\mu, \Sigma)$. This is a distribution on the sphere in $\...
0
votes
1
answer
37
views
Predicted circular regression curve from a bpnr regression
I'm trying to understant how to interpretresults from the bpnreg package to do glm on a circular response variable (flight directions).
I understood from the paper that the interpretation can be done ...
2
votes
1
answer
87
views
Modelling the effect of a continuous (linear) predictor on a non-von Mises circular data
I am working with circular data (n=800+) representing the time of animal activity in radians, scaled between 0 and 2π. My objective is to determine the influence of a linear predictor, specifically ...
1
vote
0
answers
153
views
Analytically estimate joint Von Mises distribution parameters from multiple underlying distributions with arbitrary weights
Given a set of n one dimensional (circular) Von Mises distributions, it is possible to randomly sample each distribution (with a different weight, ...
0
votes
0
answers
30
views
Maximum Entropy distribution of a ticking clock
Say I have a clock that emits "ticks". An ideal clock looks like a dirac comb. It has:
perfect periodicity of ticks (there is a precise fixed time interval between any two consecutive ticks)...
1
vote
0
answers
34
views
I am looking for a regression formula to calculate a polar curve that goes through a set of points analogous to linear regression, how do I do this? [closed]
Suppose you have an arbitrary set of points in a two dimensional space, what is the "curve of best fit" not line of best fit that can be used to model the data?
1
vote
1
answer
53
views
Conclusions on modality from Rayleigh's Uniformity Test
Rayleigh's Uniformity Test (RUT), when run on a circular distribution C, returns p ~= 0. Then, we can confidently assert that ...
10
votes
2
answers
402
views
Formulas, approximations, or bounds for $\mathbb{E}\left( \frac{X}{\lVert X \rVert} \right)$, $X\sim N(\mu, \Sigma)$?
In another question, I asked for $\mathbb{E}\left( \frac{X}{\lVert X \rVert} \right)$, in the case where $X \in \mathbb{R}^d \sim N(\mu, I_{d})$. Somebody posted an exact formula based on the symmetry ...
2
votes
1
answer
110
views
Distribution of intersection of $n$-sphere and hyperplane projected onto $c$
Suppose $x$ is uniformly sampled from an n-sphere of radius $1$ restricted to points with $\langle a, x\rangle=b$ for a given unit vector $a$ and a constant $b$
Given a unit vector $c$, what is the ...
2
votes
2
answers
117
views
Significance test of the amplitude of sinusoidal regression
I am working on a multi-level logistic regression model with a binary outcome and a circular predictor. Thanks to this post, I know that I have to add sine and cosine of the circular variable as a ...
5
votes
1
answer
116
views
Overlapping circular bearing distributions on a plane
I have some directional hydrophones capable of recognizing transient signals/sound and estimating the circular probability density function of the bearing, or direction, that the sound came from. I ...
2
votes
1
answer
38
views
Distribution of Data Within Hypersphere That Maximizes the Nuclear Norm?
Suppose I have $N$ points in $\mathbb{R}^D$ such that for each point $x_n$, its L2 norm is at most 1:
$$||x_n||_2^2 = x_n^T x_n \leq 1$$
Assuming $N > D >> 0$, if I construct a matrix $X \in \...
3
votes
1
answer
119
views
What concentration $\kappa \in [0, \infty)$ maximizes the entropy of the von Mises-Fisher distribution?
I'd like to prove what concentration parameter $\kappa \in [0, \infty)$ maximizes the (differential) entropy of a von-Mises Fisher Distribution.
The differential entropy of of a von Mises-Fisher ...
0
votes
0
answers
17
views
Cubic spline with circular predictor [duplicate]
I have a set of observations $y_i$ for a set of values of the independent variable $x_i$. $x_i$ takes values of angles, so it is a circular variable. Is there some method to perform cubic splines or ...
5
votes
2
answers
161
views
Calculating summary statistics on a distribution of "day of year" data [R]
Say I have a dataset that consists of 1,000 integers between 1 and 365 that represent days of the year a certain event happened. I am trying to figure out how I could calculate summary statistics of ...
4
votes
1
answer
126
views
Non-Uniform Spherical Distributions
Suppose $X_i\overset{\text{iid}}{\sim} N(0,1)$, and define the random vector $\mathbf{X}=(X_1,\ldots,X_n)$. Then the normalized vector $\mathbf{Z}:=\frac{\mathbf{X}}{\|\mathbf{X}\|_2}$ is uniformly ...
5
votes
1
answer
199
views
How to plot angular time series?
I am trying to inspect a circular time series (a long time series of angular measures in 0-360°). The main aim would be to identify abrupt changes in the time series, but as a start I would like to ...
5
votes
1
answer
70
views
What is the entropy of a mixture of von Mises-Fisher distributions?
What results exist for computing (or approximating) the entropy of a mixture of von Mises-Fisher distributions?
1
vote
0
answers
37
views
What is the value of $\sum_i x_i^2 v_i$ for isotropic $x$?
Suppose $x$ is an isotropic random variable in $\mathbb{R}^d$ with $E[\|x\|^2]=d$ and $v$ is some vector. It appears that $\sum_i x_i^2 v_i \approx \sum_i v_i$ when $d \approx \infty$.
What is an easy ...
8
votes
1
answer
289
views
Donut-like Distribution in Cartesian Coordinates
I have a set of points $P_i$ which are described by an angle $\theta_i$ and a magnitude $r_i$.
$\theta_i$ follows a Uniform distribution $(\theta_i \sim U(0, 2\pi))$ and $r_i$ follows a chi-k ...
8
votes
1
answer
241
views
Radial axis transformation in polar kernel density estimate
Consider a kernel density estimate of a continuous, non-negative random variable defined over the unit circle with no discontinuity between 360 and 0 degrees.
Unlike in the most common KDE ...
0
votes
0
answers
59
views
How to choose priors for bounds on circular truncated distributions?
I am considering choices of priors for truncated distributions on a circle. Let's take the truncated normal distribution on the unit circle as an example. It has parameters $\mu \in [-\pi, \pi]$ and $\...
4
votes
1
answer
80
views
Possibility priors in Bayesian analysis?
A couple of trains of thought have come together for a model I am designing. Let's start with the first part: Bayesian inference doesn't update strongly enough.
One of the parameters $\theta$ is an ...
1
vote
1
answer
105
views
watson's two sample test with ties - does it exist on R?
I would like to compare two samples with the watson's two sample test in R (circular data). As I have a lot of ties in my samples, I followed the example given in the book from Zar (Biostatistical ...
2
votes
0
answers
168
views
How to calculate confidence intervals for circular data?
I managed to convert my data i circular/polar data (angle, r). Now I would like to compare two treatments to see if they are statistically different. I thought to use the Watson's U2 test but ...
2
votes
0
answers
37
views
Does it make sense to transform a feature containing hours (24h) into two features with xy-coordinates of each hour in the space? [duplicate]
I have a clustering problem that I might solve with an algorithm based on Euclidean distance (e.g. K-Means).
One potential feature is the "hour" at which each user began an interaction.
As ...
1
vote
0
answers
347
views
How do I interpret the sine and cosine terms in a linear-circular regression?
When using a circular variable as a predictor (X) in a linear regression, the consensus seems to be to split the variable into sine and cosine terms.
e.g. Y ~ cos(X) + sin(X) + etc
This creates two ...
1
vote
1
answer
153
views
Python / R package for multivariate gaussian process regression for circular data?
I have two circular predictor features (two angles between 0 and 360 degrees) and a circular outcome (another angle, between 0 and 360 degrees). I'd like to be able to fit a model and get predictions ...
3
votes
1
answer
573
views
Regression of circular variable with scikit-learn
I am trying to use Support Vector Regression on a (neurophysiological) dataset where the position of points on a circular manifold in N dimensions is correlated with a circular variable (phase of an ...
3
votes
0
answers
95
views
Clustering and grouping of rhythmic data (acrophases)
I am looking for suggestions towards a clustering method for rhythmic/oscillation data.
We performed cosinor regression (sin(2*pi*time/period)+cos(2*pi*time/period))...
2
votes
0
answers
116
views
Calculating the parameters of von Mises distribution
I would like to calculate the values of the concentration ($\kappa$) and mean direction ($\mu$) for a von Mises mixture model from the theta values given by the movMF() function in R. At the bottom of ...
1
vote
0
answers
49
views
Within a circle, Find a probability distribution to model a random variable with distance to center and angle information
Within a circle with radius R, a random variable with two types of information as below:
Distance to the circle center. The distance should satisfy uniform distribution~U[a R], where 0<=a<R is ...
2
votes
3
answers
417
views
Measuring dispersion in circular data
stats noob here.
I have a (circular) dataset with values in [0, 2pi]. I need some kind of a measure of how disperse or diverse the dataset is. I have looked a bit into non circular (regular) data, and ...
2
votes
0
answers
74
views
Is circular correlation maximized by points laying on a square grid on the plane?
Background
Circular correlation is given by:
$$R_{\operatorname{circular}} \triangleq \frac{\sum_{i=1}^m \sin (x_i - \bar x) \sin (y_i - \bar y)}{\sqrt{\sum_{i=1}^m \sin (x_i - \bar x)^2} \sqrt{\sum_{...
1
vote
1
answer
374
views
Definition of directional variance
We are currently adding some basic functionality for directional statistics to SciPy. Directional statistics refer to data whose magnitude does not matter such as unit vectors.
While implementing the ...
5
votes
0
answers
169
views
Running a poisson GLM with cyclical explanatory variable?
Im running a few poisson GLMs looking at count data of bats sightings in relation to lunar cycle.
We specifically want to look at how species are affected leading up to and away from new moon/full ...
0
votes
0
answers
36
views
How best to model day time dependet activity?
I aquired the median activity data of animals per hour and wanted to model them dependent on time. Now I have tried different models and came up with the following. (This only includes a small portion ...
1
vote
1
answer
346
views
How to calculate the value (wind direction) of "maximal effect" in a linear-circular correlation?
I am trying to evaluate the effect of wind direction (circular variable) on a dependent linear variable. I have used circular-linear regression to find the correlation coefficient between the two ...
3
votes
1
answer
353
views
Which statistical test should I use to see if the month of onset of a disease is random or clusters in the winter?
I have data showing which month a disease started for 119 subjects. The months are represented by 1-12 (Jan - December). We hypothesize the disease more often onsets during winter months. What test ...
1
vote
0
answers
110
views
How can I quantify this non-linear relationship where y-values peak within a narrow range of x-values?
I have plotted this data, and to me there appears to be a clear relationship between the variables (sediment accumulates when wind direction is between 150 and 250 degrees). Is there a appropriate ...
11
votes
1
answer
2k
views
Why doesn't mean square error work in case of angular data?
Suppose, the following are the first few lines from a dataset for solving a regression problem:
...
2
votes
0
answers
128
views
Machine learning regressor for cyclic values
I have a question regarding training a Machine Learning regressor for cyclic values. For instance, consider phase values. We cannot use conventional regression methods for this as for example the ...
7
votes
3
answers
2k
views
Why does the von Mises-Fisher distribution need two parameters?
The von Mises-Fisher distribution has two parameters: the mean $\mu \in \mathbb{R}^p$ and concentration $\kappa \geq 0$, where $\mu$ is constrained to have unit norm.
Why not instead define the ...
2
votes
0
answers
303
views
Variance of von Mises-Fisher Distribution
As a follow up to this previous question on the expectation of the von Mises-Fisher distribution, what is the variance of a von-Mises Fisher distribution as a function of the mean direction $\mu$ and ...
13
votes
2
answers
1k
views
Fitting sine wave with lm in R for circadian activity- frequencies?
I'm trying to fit a sine wave over some activity data, just like this post. I've managed to get a reasonable looking graph for one condition:
However, when I plot the other condition the wave looks ...