Notice: On March 31, it was announced that Statalist is moving from an email list to a forum. The old list will shut down on April 23, and its replacement, statalist.org is already up and running.

# Re: st: RE: RE: calculating standard deviation of age at death from period life tables

 From Steve Samuels To statalist@hsphsun2.harvard.edu Subject Re: st: RE: RE: calculating standard deviation of age at death from period life tables Date Mon, 8 Mar 2010 12:42:23 -0500

```Tony is on the right track.  This requires a life-table solution.  The
SD for age at death in period (or "current") life tables  is that of
the ages at death for a synthetic cohort whose age-specific death
rates are those observed for the period, in this case 2004.  This
calculation requires information about the numbers at risk in each age
during the year, not just about the numbers who died at each age.

I've not constructed a period life table in Stata, but see:

http://data.princeton.edu/eco572/periodlt.html
http://www.stata.com/statalist/archive/2004-02/msg00127.html

Getting from these results to the SD will take some work: from the
table probabilities estimate the mean at death m; then estimate the
mean of (age -m)^2

-Steve

On Mon, Mar 8, 2010 at 11:57 AM, Lachenbruch, Peter
<Peter.Lachenbruch@oregonstate.edu> wrote:
> Am I missing something?  Doesn't the Greenwood formula give this to you, and isn't it in the ltable command?
>

> -----Original Message-----
> From: owner-statalist@hsphsun2.harvard.edu [mailto:owner-statalist@hsphsun2.harvard.edu] On Behalf Of Martin Weiss

> Either use frequency weights (first example) or -expand- to the required
> number:
>
>
> I want to calculate the standard deviation of age at death from period
> life tables:
> So I know how many people died at age x, where x runs from 0,1,2...110:
> For example, in 2004 in the US people 2000 died at age 0, 1900 at age 1,
> 600 at age 2, etc.
>
> What I think is the right way of calculating is:
>
> To generate a new variable that contains  2000 times 0, 1900 times 1,
> 600 times 2 etc...
> But I don't know how I can generate this new variable in Stata most
> efficiently...
> Thanks a lot for any hint!
> Regards,
> Johannes
>
> --
>
>
Steve Samuels
sjsamuels@gmail.com
18 Cantine's Island
Saugerties NY 12477
USA
845-246-0774

*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/
```