Hi everyone,
recently questions came up with regards to statistics on the unit sphere in scipy.
Issue 12041 raises the question if spherical mean/variance should be added as descriptive statistics similar to
circmean/circvar . The canonical reference for all things related to statistics on the unit sphere is "Directional Statistics" by Mardia & Jupp.
Another important basic functionality is the generation of random samples on the unit sphere: discussion was started in
this issue .
Any strong opinion against these? I think at least the random number generation is reimplemented every day around the world (done so myself at least 2 times) and would strongly benefit from being part of scipy.
Best