Arcsine distribution ($\mathrm{Beta}(\frac{1}{2},\frac{1}{2})$)

The standard arcsine distribution is a continuous probability distribution on $[0,1]$ with probability density functionIt is the beta distribution with both parameters equal to $1/2$. The density grows without bound toward both endpoints, yet neither endpoint has any probability mass. Its name comes from the inverse sine in its cumulative distribution function.^[ref]

A useful way to see the shape is to choose a point uniformly around a circle and project it onto a diameter. Equal angular intervals do not give equal horizontal distances: the projection moves slowly near either end. An ideal harmonic oscillator sampled at a uniformly chosen time in one complete period has the same endpoint concentration after rescaling its position. The distribution translates a simple geometric slowdown into a precise probability calculation.

For an original comparison, take equal-length regions: the middle half $[1/4,3/4]$, and the two outside quarters together. ThenA uniform distribution on $[0,1]$ instead assigns each region probability $1/2$. Both distributions have mean $1/2$; the mean alone hides their different concentration.

Choose $\Theta$ with uniform distribution on $(0,\pi /2)$ and set $X={\sin }^2\Theta$. Since the sine is increasing on this interval,Differentiation gives $f(x)=1/(\pi \sqrt{x(1-x)})$, and $F(0)=0$, $F(1)=1$ prove normalization. This also evaluates the density's improper integral without treating its divergent endpoint values as infinite probability.

For the circle construction, instead choose $\Phi$ with uniform distribution on $[0,2\pi )$ and put $Y=\cos \Phi$. Each $-1<y<1$ has two preimages in one revolution. The change of variables formula adds both contributions:The affine change of variables $X=(1+Y)/2$ gives the standard arcsine distribution. A single monotone-branch calculation without the second preimage would give only half the total probability.

The same geometry explains the ideal harmonic oscillator $y(t)=A\cos (\omega t+\varphi )$, with amplitude $A>0$ and angular frequency $\omega >0$. Choose $T$ uniformly in any complete period of duration $2\pi /\omega$. Its oscillation phase modulo $2\pi$ is uniform, and $Y=y(T)$ has densityNear a turning point, $|dy/dt|=\omega \sqrt{A^2-y^2}$ is small, so a fixed position interval takes more time to cross. For the chosen amplitude $A=2\mathrm{cm}$, exactly one third of the period is spent in $|y|\le 1\mathrm{cm}$, and two thirds outside that central interval. These are calculations for an ideal sinusoid, rather than measured oscillator data.

The sampling assumption matters. Taking one observation every full period returns the same position each time. It does not produce this probability distribution. Uniformly random observation times and a fixed regularly spaced sampling scheme are different models; aliasing can conceal the motion altogether.

Paul Lévy used a related geometric construction in his work on stochastic processes: project a uniformly chosen point on a circle onto its diameter to divide a segment. The projection description appears explicitly in §5 of his paper, following its arcsine density calculation.^[ref]

An unbounded probability density function does not imply an atom. HereThe last equivalence follows from $\arcsin u/u\to 1$. The endpoint interval has much more probability than the same interval under a uniform distribution, whose value is $\epsilon$, but its probability still tends to zero. This is an integrable singularity.

For example, the first 1% of the unit interval contains about $6.377\%$ of the probability; the last 1% contains the same amount because the arcsine distribution is a symmetric probability distribution. Each exact endpoint still has zero probability. Density and interval probability answer different questions.

A numerical integration routine must respect these endpoint integrable singularities. The substitution $x={\sin }^2\theta$ givesIt turns a difficult-looking density integral into a constant measure in angle. Integrating finite bins is also a faithful way to draw the distribution near its endpoints; sampling $f(0)$ to make a bar height is not.

Inverting $F(x)=p$ gives the quantile functionConsequently, inverse transform sampling uses $X={\sin }^2(\pi U/2)$ with $U$ uniform on $(0,1)$. The quartiles are $(2-\sqrt{2})/4$, $1/2$ and $(2+\sqrt{2})/4$.^[ref]

The first and third quartiles are about $0.1464$ and $0.8536$, so the central half of the probability spans about $0.7071$ of the unit interval. Under a uniform distribution, its width would be $0.5$.

The symmetric probability distribution gives $\mathbb{E}[X]=1/2$. The angle representation gives a short calculation of the second moment:Use ${\sin }^4\theta =(3-4\cos 2\theta +\cos 4\theta )/8$ to evaluate the integral. A uniform distribution on $[0,1]$ has variance $1/12$; the larger variance $1/8$ quantifies the greater spread toward the ends.

More generally, for a nonnegative integer $n$, the beta function yieldsFor $n=0$ this is one; for $n=1,2$ it reproduces the values above. The central binomial coefficient thus appears in an exact moment calculation for a continuous probability distribution.

Specify whether two parameters mean endpoints or location and width. With endpoints $a<b$, define $Y=a+(b-a)X$ for a standard arcsine distribution variable $X$. ThenThis is a location-scale family with location $a$ and width $b-a$. Its mean is $(a+b)/2$ and its variance is $(b-a{)}^2/8$. A convention using parameters $a,w$ for location and width instead has right endpoint $a+w$, so it must not be read as having endpoints $a,w$.^[ref]

For example, a chosen interval $[-2,6]$ gives mean $2$ and variance $8$, while its density still concentrates toward both endpoints. It is a bounded probability distribution; its endpoint singularities do not make its support infinite.