I have a process 𝑃 generating random variables X_1, ... X_n
. From each of these I've sampled a set of samples S_i = (s_i1, ... s_im)
, respectively.
I want to test whether the process/real-world dynamic generating the variables X_i
is consistently (or pretty consistently) generating normally distributed variables (not identically distributed, though! They can have different μ and σ).
One suggestion I got was to think of all samples (S_1, ..., S_n)
as coming from a gaussian mixture. However, it didn't come with a suggestion for a statistical test, and I couldn't find one.
My original thought was to test each sample separately for normality (with something like Shapiro-Wilks), control the False Discovery Rate in the set of tests (using something like Benjamini-Yekutieli; I can't assume independence), and then see whether the rate of rejections (in the case of Shapiro-Wilks, this means non-normality) is under/over some predetermined threshold (like 10% or 20%), in which case I deem the mysterious process as one not generating gaussians consistently enough.
However, I'm no statistician, so I've come to ask for help. Thank you! :)
Clarification:
One shouldn't think of X_i
as if they were collected from the same process. 𝑃 is a process generating processes, if you want. Or distributions, or populations.
Think of 𝑃 as "my marketing team" and different X as different marketing campaings they output, where each such campaign generates different revenue streams, or engagement streams. One might be a tv campaign in Belfast, another an online campaign in Kenya. We thus do not expect the revenue/engaement/whatever distribution for all of them to have the same μ and σ. Therefore, the question is whether revenue streams from campaigns generally distribute normally or not.