I am looking into the pros and cons of each normalisation technique for work and it got me thinking. What if I used trimmed means and the sqrt of Winsorized variances to compute the standardised data? Instead of
x<- rnorm(100, 1, 2)
y<- x-mean(x)/sd(x)
it becomes
y<- x- mean(x, trim= 0.2)/sqrt(winvar(x))
My thinking is that this won't make much difference to normally distributed data but in the case of non-normal data it might place the mean within a more accurate place in the distribution, such that the true number of points below or above the mean will be known.
This might all be rubbish but let me know what you think.