Statalist The Stata Listserver

[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

Re: st: jacknife 2 replication method

From (Jeff Pitblado, StataCorp LP)
Subject   Re: st: jacknife 2 replication method
Date   Tue, 16 May 2006 14:41:12 -0500

Leslie R Hinkson (lhinkson@Princeton.EDU) needs to use the jackknife 2
variance replication method:

> I have a dataset that already has 62 replicate weights as well as a main
> weight that I thought was an iweight but am beginning to suspect may be a
> pweight. I was told that in order to do any analysis (including simply
> generating means and standard deviations) that I'd have to perform a
> jacknife 2 replication method. I tried to use the jknife command in Stata
> but couldn't get anywhere. Is this because I'm a qualitative researcher
> struggling with Stata or is it because I'm using the wrong command? 

Leslie didn't mention Stata 8 or Stata 9.  For Stata 8, Leslie will have to
use Nick Winter's -svr- suite of survey commands.

	. ssc install svr
	. help svr

Jackknife variance estimation using replication weight variables is a new
feature in Stata 9.  Assuming Leslie's sampling weight variable is named
-sampwgt-, and the replicate weight variables are named -jkrw1-, -jkrw2-, ...,
-jkrw62-; the -svyset- command is as follows:

	. svyset [pw=sampwgt], jkrweight(jkrw*, multiplier(1)) vce(jackknife)

(note: the shortcut 'jkrw*' will expand to 'jkrw1 jkrw2 ... jkrw62')

Then Leslie can compute means using -svy: mean-.  Note that I specified the
suboption -multiplier(1)- in the -jkrweight()- option, it specifies the value
of the m_h term in the jackknife variance formula.

	By default, -svy jackknife- assumes the 'jackknife n' replication
	method, using the formula

		m_h = (n_h-1)/n_h

	for his multiplier (n_h is the stratum sample size), but is not
	correct for the 'jackknife 2' replication method.

Leslie might also want to specify the -mse- option on -svyset- so that -svy
jackknife- uses the mean squared error formula for the variance calculation.
See [SVY] vce and [SVY] svyset for more details about jackknife variance
estimation and survey design settings.

*   For searches and help try:

© Copyright 1996–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index