I'm comparing wind direction distributions using Python's scipy.stats.circmean
and scipy.stats.circstd
, and I encountered some unexpected results.
As illustrated for example in this answer, for 'normal' (non-circular) data, the mean of the difference between two samples is equivalent to the difference between their means. I blindly assumed this would be the same for circular statistics, but I'm getting different results and I'm suspecting that for circular statistics, things might be different. I just can't wrap my head around why this would be the case.