5
$\begingroup$

I have some data and I am trying to identify its distribution. The nearest I can get is a skewed-Gaussian distribution, but I don't think it is. The data itself consists of 130000 points and is binned with the Freedman–Diaconis rule.

I also tried Poisson, Normal-Log, gamma and chi-squared distributions which have the right shape but the parameters never match the data. Here is a plot of the data: enter image description here

The black curve is the best approximation I can get -- the skewed Gaussian, however if I try and generate a set of artificial data with the fit results (using Mathematica's SkewNormalDistribution[...]) it doesn't match the original set at all.

I'm adding some further information here. The data shown here is the background noise of multiple spectra acquired from experiment. I want to understand the characteristics of this noise, so that I may reproduce it in simulation.

In order to do this I have attempted to fit the distributions to the histogram I have shown here to try and determine the distribution this spectral noise takes. If I have a successful fit I can use the extracted fit parameters to generate the simulated data. For example with a Skew-Normal distribution I can extract $\mu$, $\alpha$, and $\sigma$ and then use these to generate my simulated distribution.

Here is a crude probability plot I have made quickly (I have little contact with these kind of plots so am unable to do anything more clever):

enter image description here

I am adding some further information as the source of the data seems to be relevant. The data shown in the histogram comes from the amplitude of the noise floor of some FFT spectra. The unit the data is recorded in is $\rm{dBV_{pk}}$, which is $20\log_{10}(V_{pk})$ I have extracted $V_{pk}$ and multiplied by $10^6$ for the purposes of fitting (it's usually easier to fit scaled data). The voltages I am dealing with are consequently at the $\rm{\mu V}$ scale. Hence everything exists between $0$ and $1$ in the histogram.

I would have EXPECTED this distribution to be Gaussian white. Which is what Johnson-Nyquist noise is (at least in the regime I am measuring in). But it clearly isn't. There is something going on here, and THIS is why I want to know the distribution. Is it Gaussian convoluted with a Boltzmann? Possibly but unlikely the temperature gradient would be FAR too steep. Is this a consequence of some filter at the input of my FFT spectrum analyser? I deal with distributions a lot, but I've never come across anything like this -- so hence I am asking the stats experts!

So I think it is probably necessary to show what I am dealing with from the beginning. In the following graph we see an FFT spectrum:

enter image description here

This is an FFT of some transient data (I don't have access to the transient itself). The red points are the peak feature -- I DROP these from the data set for the purposes of these discussions and only take the blue circles. I am interested in the distribution of these circles, for many many spectra. So there is an underlying signal which is noised, but the linewidth of the feature is so small ($\rm{mHz}$ level) I expect the remaining data to not be biased by the peak feature. The amplitude data has been linearised originally in $20\log_{10}(V_{pk})$ units. Looking at the spectra one can already see that it isn't exactly white but that's why I want to learn about the distribution.

In case it helps here is the same data but left in its original $\rm{dB_{V_{pk}}} = 20\log_{10}(V_{pk})$ form. This DOES look like a skewed normal distribution, so I suppose I could always work from this angle and convert the results back to my linear units for use in simulation.

enter image description here

|cite|improve this question
$\endgroup$
  • 4
    $\begingroup$ Histograms do not provide sufficient detail to do anything more than (a) rule out obviously poor candidates and (b) permit intelligent guesses from among a huge array of possibilities. What you need to provide to make this question answerable are (a) information about why you are fitting a distribution and (b) better descriptions of the data, such as probability plots, N-letter summaries, etc. $\endgroup$ – whuber♦ Jun 3 at 16:41
  • $\begingroup$ @whuber I'be added a little more detail. I hope this is enough. I can provide a dump of the data if this is helpful. $\endgroup$ – QuantumPenguin Jun 3 at 17:03
  • 1
    $\begingroup$ I am reading that the data are positive and continuous. That rules out Gaussian and Poisson. $\endgroup$ – Nick Cox Jun 3 at 17:59
  • 2
    $\begingroup$ How about just sampling randomly from your 130k noise observations when you need to add noise to what you simulate. $\endgroup$ – BruceET Jun 3 at 19:37
  • 2
    $\begingroup$ I'm still trying to understand the sense in which non-negative values represent "noise." Almost by definition, "noise" is a random component in a system that, on the average, does not change the values. Thus, if it's additive noise it must have a mean of zero and if it's multiplicative noise it must be centered at $1.$ This histogram is neither. I begin to suspect it might be a histogram of some measure of the variance or dispersion of a "noise" component in a signal, but the lack of details prevents further inference of that sort. What does this histogram really represent?? $\endgroup$ – whuber♦ Jun 3 at 19:39
2
$\begingroup$

I collected $2^{20}$ values from a unit normal process, did the FFT and binned the magnitudes. Then overplotted with a Rayleigh distribution:

Hist & Rayleigh distribution

I did no scaling on anything, because I was working fast, but I will go back and do it.

|cite|improve this answer
$\endgroup$
  • $\begingroup$ That looks really super promising! I've not come across a Rayleigh distribution before! $\endgroup$ – QuantumPenguin Jun 5 at 8:29
  • $\begingroup$ Exactly perfect. And reading the literature, the distribution matches perfectly to what we see -- a distribution for positive values only. $\endgroup$ – QuantumPenguin Jun 5 at 9:22
  • $\begingroup$ Very happy this turned out well, but a bit annoyed at myself: after the hugh clue that the FFT spectrum analyzer output was absolute values, I should have known immediately that they must be Rayleigh distributed. $\endgroup$ – Ed V Jun 5 at 21:06
  • $\begingroup$ Somewhat the same! I was previously unaware of a Rayleigh distribution but it didn't occur to me quick enough that that the sign and phase information is quite obviously destroyed...anyway thanks again old timer! $\endgroup$ – QuantumPenguin Jun 5 at 22:07
1
$\begingroup$

So I have understood my problem. This is basically a consequence of taking the FFT of time transient data and taking the absolute value of it. The FFT spectrum analyser device actually spits out the absolute value -- so the phase and sign information of the original transient is LOST.

You can prove this simply by generating a random list of numbers, normal distributed, and FFT it. Then take the absolute value and plot as a histogram. You get exactly the same distribution as I have shown in my question.

It would still be nice to know what the actual distribution of this data is -- as in the shape of it. But I can basically reconstruct my original noise distribution and verify that it is indeed Gaussian distributed.

|cite|improve this answer
$\endgroup$
  • $\begingroup$ I just did what you suggested in the second paragraph and confirmed it. As you say, it would still be nice to know the actual distribution. $\endgroup$ – Ed V Jun 4 at 15:52
  • 1
    $\begingroup$ Please look at this link near the bottom : dsp.stackexchange.com/a/44524/41790 $\endgroup$ – Ed V Jun 4 at 16:03
  • $\begingroup$ @EdV Thanks for your input!! It's always good to have the input of the experienced! I'm working in something a little different to you. Trapped ion spectroscopy with RF integration. $\endgroup$ – QuantumPenguin Jun 4 at 17:47
  • 1
    $\begingroup$ Nice research topic! Also, although I have been of almost no help thus far, I am following a lead that suggests the noise is Rician distributed. I can easily test this tonight and will update ASAP! $\endgroup$ – Ed V Jun 4 at 19:43

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.