Normal distribution random variable
1 view (last 30 days)
Show older comments
If normal random variable have parameters N(0,1), how can I find results from possible 1000 which I expect to show in interval 1-2
2 Comments
John D'Errico
on 18 May 2018
Edited: John D'Errico
on 30 May 2018
When you delete the question and even the subject, you insult those who bothered to spend the time to answer your question. You hurt the answers site.
If you cannot leave your question as it is, then you should not post questions.
Answers (2)
John D'Errico
on 17 May 2018
This is an extremely confusing question. Usually, that implies the person asking it is also confused.
It sounds like you have a population of 1000 samples from a normal random variable. For example,
z = randn(1,1000);
So the vector z will have 1000 samples. Now, do you want to do something with only those samples that lie in the interval [1,2]? If so, then you could just use find, testing if z>=1 and z<=2. Even simpler is just this:
z12 = z(z >= 1 & z <= 2);
If that is not what you are looking to see, then you need to VERY clearly explain your problem.
1 Comment
John D'Errico
on 17 May 2018
Edited: John D'Errico
on 17 May 2018
Actually, the theoretically expected frequency is
format long g
diff(normcdf([1 2]))
ans =
0.135905121983278
NOT 0.1365. So I'm not sure what you used to get that estimate.
Anyway, what is the problem? Do you seriously expect to see EXACTLY 136 samples in those bounds? If you do, when why not expect EXACTLY 135.905121983278 samples in that interval? Surely you see the problem in assuming the frequency to be exact.
x = randn(1,1000);
sum(x >= 1& x <= 2)
ans =
111
Which is reasonably close. Remember, these are samples taken from a normal distribution. The larger the sample size, the closer you will expect the fraction to be close to 13.59...%. That is just basic statistics.
x = randn(1,1000000);
sum(x >= 1& x <= 2)/numel(x)*100
ans =
13.607
x = randn(1,10000000);
sum(x >= 1& x <= 2)/numel(x)*100
ans =
13.5921
At the other extreme, if I sample from a Normal distribution, 50% of the time, the number "should" be greater than 0. Surely you agree with that? But what if I take only ONE sample? Sometimes it will be positive, sometimes negative. But a limited sample population will never show the EXACT distribution frequency, based on theory. So even though I expect 50% of my samples to exceed 0, oops! Look what happened:
randn(1)
ans =
-0.633677491587928
100% of my samples from this latest set were less than 0.
Again, this is all just basic statistics, and understanding random variables.
See Also
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!