[R] Standard error of standard deviation: bootstrap or theoretical results?

Thomas W Blackwell tblackw at umich.edu
Wed Aug 6 15:04:13 CEST 2003


Huan  -

The difference between the empirical ("bootstrap') result and the
theoretical results shows evidence for autocorrelation in the time
series data.

-  tom blackwell  -  u michigan medical school  -  ann arbor  -

On Wed, 6 Aug 2003 huan.huang at bnpparibas.com wrote:

> This is more a statistical question rather than an R question. I'd
> appreciate it if you can give me some suggestions.
>
> I have a sample of a time series (sample size 500, fat tail in density). I
> am trying to calculate the Standard error of standard deviation of a
> sub-block-sample (sample size 250). I take 100 this kind of
> sub-block-sample, randomly. For these 100 subsamples, I use the following 3
> methods to calculate the standard error of standard deviation:
>
> 1. From book "Handbook of applicable mathematics", Walter Ledermann (chief
> editor) Volumn VI: Statistics, Part A, Lloyd, John Wiley & Sons.
>
> Page 30-32:
>
> var(S) = (mu4 - mu2^2)/(4 * mu2 * n)
> mu4 = E(X - mu)^4, mu2 = E(X - mu)^2, S^2 = sum(X - mu)^2/n
>
> The results are about: 0.00090
>
> 2. From   http://davidmlane.com/hyperstat/A19196.html
> The results are about 0.00066
>
> 3. From   http://mathworld.wolfram.com/StandardDeviationDistribution.html
> The results are about 0.00065
>
> Finally I calculate the standard deviation for each of the 100 subsamples
> and the standard error of those 100 standard deviations ( I reckon this is
> the bootstrap result for the standard error of the standard deviation I
> want).
> I get 0.00024
>
> I tried all above a couple of times and got similar results for each
> methods I used. The results from the first 3 methods are apparently higher
> than the bootstrap one. I am a bit confused. Do I miss anything? Which one
> do you believe?
>
> Huan




More information about the R-help mailing list