Re: WHY IS SOUND to get C.V.´s from Mont Carlo
- From: Jack Tomsky <jtomsky@xxxxxxxxxxxxx>
- Date: Thu, 29 Mar 2007 19:32:50 EDT
Let be the Hypotheses Test:
_____H0: m=0 , Ha: m =/0
for means m of normal samples WITH SD=1 with size n
----N(m,1):n
It should be noted that the source population as a
KNOWN MEAN that is equal to 1.
Procedure:
Because throughout the test I must admit that H0 is
true I construct an huge number N =1 million (or 4
millions, . . . ), of N(0,1):n samples and I evaluate
the mean of each one.
The set of N means is very similar to the analytical
Distribution Function of the statistics (the
Givenko-Cantelli theorem assures it).
Then I evaluate the quantiles
_____u = prob(mean)< alpha/2
_____v = prob (mean) > 1- alpha/2
that bounds the rejection region FOR ALL MEANS OF THE
SAMPLES OF SIZE n DRAWN FROM N(m,1).
I evaluate (other thread) this bounds for the
confidence levels 95% and 99%. They were:
__Size____
__15_____0.354, 0.646_______0.310, 0.689
__16_____0.359, 0.641_______0.316, 0.684
__17_____0.363, 0.637_______0.322, 0.678
__18_____0.367, 0.633_______0.327, 0.674
__19_____0.370, 0.630_______0.331, 0.669
__20_____0.374, 0.626_______ 0.335, 0.665
This values are ABSOLUTELLY EQUAL TO THOSE OBTAINED
BY THEORY.
____licas (Luis A. Afonso)
These "confidence intervals" are not equal to those obtained by theory. The exact 1-alpha level confidence intervals on the population mean for samples of size n with sigma known and equal to one are
xbar +/- z(1-alpha/2)/sqrt(n).
Jack
.
- References:
- Re: WHY IS SOUND to get C.V.´s from Mont Carlo
- From: Luis A. Afonso
- Re: WHY IS SOUND to get C.V.´s from Mont Carlo
- Prev by Date: Re: WHY IS SOUND to get C.V.´s from Mont Carlo
- Next by Date: Re: legally liable
- Previous by thread: Re: WHY IS SOUND to get C.V.´s from Mont Carlo
- Next by thread: Re: WHY IS SOUND to get C.V.´s from Mont Carlo
- Index(es):