Closed ghost closed 1 year ago
Hmm, I'm worried that it's gaussianizing each row separately (i.e., each row is considered a single sample of some 1-d data). Can you try transposing the data matrix and see if it works?
The result is the same, but now the coefficients in object out are a list of 4001 tuples of size 3. What I keep seeing is that the 3rd value in the tupple is always equal to zero. Maybe this is the problem?
As you can see in the picture below, my data (x) are somewhat close to gaussian, but zero change (in y) seems strange to me.
Maybe I am missing something so I am going to be more clear: I have a timeseries of 4001 points and I want to gaussianize it. I expect the values close to 6 and 7 to contract and the values around 0 and 1 to dialate so that my above histogram will transform into more gaussian-like. Is this the intent of your class?
Hi Nikos, you are correct in interpreting the intent of the class. I'm afraid I haven't used it in a long time and don't remember much about the details. Let me make a few observations though.
Sorry I can't be of more help, this was kind of a weekend project from a few years ago and I haven't really used it much myself since then. If you do figure anything out, please comment back here in case it's useful for others trying to do this type of Gaussianization.
Hi @Nikos-T , FWIW looking at your histogram chart it does not seem like your data even has heavy tails (its a truncated distribution on the left, and on the right it's barely heavy tailed). I wouldnt be surprised if 'delta' is estimated at 0. What's the kurtosis / skewness of your data in the histogram?
Hello, unfortunately, I do not remember neither my use-case, nor its data. However re-reading my problem, I suppose you are correct and I probably started with a non-skewed distribution, or I misinterpreted the scope of this "gaussanization". I am closing the issue, thanks for the help.
I tried to transform my data to gaussian distribution and the output of the transform function is the same as the original data. I created this script to run some tests. Thank you in advance. Test script:
Output:
The fit of x_uni has different coefficients compared to x_norm but y_uni and y_norm are equal to x_uni and x_norm respectively. In my understanding the transformed signal should scale to resemple a gaussian distribution. Am I missing something or not running the transforms correctly?