comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » Re: error on stddev estimate
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Switch to threaded view of this topic Create a new topic Submit Reply
Re: error on stddev estimate [message #61742] Tue, 29 July 2008 10:30
R.G. Stockwell is currently offline  R.G. Stockwell
Messages: 363
Registered: July 1999
Senior Member
"Vince Hradil" <hradilv@yahoo.com> wrote in message
news:570eec93-2c58-4a2a-b615-b5fb6238f8f4@l42g2000hsc.google groups.com...
On Jul 29, 11:44 am, pgri...@gmail.com wrote:
> R.G. Stockwell wrote:
>> I am analyzing some data, and I am making a point about
>> how the stddev of a time series is changing over time.
>
>> After a quick look through the help, there doesn't seem to be
>> an estimate on the error of a standard deviation estimate (in the code).
>
>> Anyone have that code handy?
>
> Numerical recipes (3d ed.. page 722) says that for gaussians,
> the standard error of the measured standard deviation sigma
> is given by sigma/sqrt(2*n) for a n-element sample.
>
> Ciao,
> Paolo
>
>
>
>> As an example, my data is like 2006 had a mean of 10 and a std dev of 2.
>> 2007 had a mean of 10 and a std dev of 4. I'd like to get estimates
>> on the error in my results of std devs (of the 2 and the 4).
>
>> Cheers,
>> bob
>
>

> You could bootstrap to find the standard error:
> 1-resample, with replacement
> 2-determine standard deviation
> 3-repeat lots of times
> 4-standard error of sd is the standard deviation of the determined
> values.

Thanks Vince and Paolo,
yeah I was thinking bootstrapping would be the way to go.
Paolo, thanks for the tip, i'll check out NR.

Cheers,
bob
Re: error on stddev estimate [message #61743 is a reply to message #61742] Tue, 29 July 2008 10:01 Go to previous message
Vince Hradil is currently offline  Vince Hradil
Messages: 574
Registered: December 1999
Senior Member
On Jul 29, 11:44 am, pgri...@gmail.com wrote:
> R.G. Stockwell wrote:
>> I am analyzing some data, and I am making a point about
>> how the stddev of a time series is changing over time.
>
>> After a quick look through the help, there doesn't seem to be
>> an estimate on the error of a standard deviation estimate (in the code).
>
>> Anyone have that code handy?
>
> Numerical recipes (3d ed.. page 722) says that for gaussians,
> the standard error of the measured standard deviation sigma
> is given by sigma/sqrt(2*n) for a n-element sample.
>
> Ciao,
> Paolo
>
>
>
>> As an example, my data is like 2006 had a mean of 10 and a std dev of 2.
>> 2007 had a mean of 10 and a std dev of 4.  I'd like to get estimates
>> on the error in my results of std devs (of the 2 and the 4).
>
>> Cheers,
>> bob
>
>

You could bootstrap to find the standard error:
1-resample, with replacement
2-determine standard deviation
3-repeat lots of times
4-standard error of sd is the standard deviation of the determined
values.
Re: error on stddev estimate [message #61744 is a reply to message #61743] Tue, 29 July 2008 09:44 Go to previous message
pgrigis is currently offline  pgrigis
Messages: 436
Registered: September 2007
Senior Member
R.G. Stockwell wrote:
> I am analyzing some data, and I am making a point about
> how the stddev of a time series is changing over time.
>
> After a quick look through the help, there doesn't seem to be
> an estimate on the error of a standard deviation estimate (in the code).
>
> Anyone have that code handy?

Numerical recipes (3d ed.. page 722) says that for gaussians,
the standard error of the measured standard deviation sigma
is given by sigma/sqrt(2*n) for a n-element sample.

Ciao,
Paolo


>
> As an example, my data is like 2006 had a mean of 10 and a std dev of 2.
> 2007 had a mean of 10 and a std dev of 4. I'd like to get estimates
> on the error in my results of std devs (of the 2 and the 4).
>
> Cheers,
> bob
  Switch to threaded view of this topic Create a new topic Submit Reply
Previous Topic: error on stddev estimate
Next Topic: how to save the console output as a ascii or text file in IDL

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Wed Oct 08 13:53:02 PDT 2025

Total time taken to generate the page: 0.00774 seconds