Sample uniformly distributed circular and spherical data
Sample uniformly distributed circular and spherical data
Simulation of the uniform distribution on [0,2π) and c("\n", "Sp−1:=xinRp:∣∣x∣∣=1"), p≥2.
r_unif_cir(n, M =1L, sorted =FALSE)r_unif_sph(n, p, M =1L)
Arguments
n: sample size.
M: number of samples of size n. Defaults to 1.
sorted: return each circular sample sorted? Defaults to FALSE.
p: integer giving the dimension of the ambient space Rp that contains Sp−1.
Returns
r_unif_cir: a matrix of size c(n, M) with M random samples of size n of uniformly-generated circular data on [0,2π).
r_unif_sph: an array of size c(n, p, M) with M random samples of size n of uniformly-generated directions on Sp−1.
Examples
# A sample on [0, 2*pi)n <-5r_unif_cir(n = n)# A sample on S^1p <-2samp <- r_unif_sph(n = n, p = p)samp
rowSums(samp^2)# A sample on S^2p <-3samp <- r_unif_sph(n = n, p = p)samp
rowSums(samp^2)