Computing optimal sampling designs for two-stage studies



26


Stata Technical Bulletin


STB-58


Description

The meanscor command performs a weighted logistic regression using the mean score method. This function requires the
complete covariate(s) to be categorical, and the default output contains the regression coefficient estimates and their standard
errors in odds-ratio form.

An important area of application of this function is in the analysis of data from a two-stage study. In this type of study,
some variables are incomplete due to only a subset of the study subjects being sampled at the second stage (Reilly 1996).

Options

f irst Vaarlist) specifies the complete covariates.

second(aαrlist) specifies the incomplete covariates.

odd(#) specifies whether the odds-ratio (odd = 1) or regression coefficients (odd = 0) format is reported. Default value is 1.

Methods and Formulas

The mean score estimates will maximize the weighted likelihood

Z nI(zi,yi)

Σ (1 + √⅛j) ⅛spβ(yih)

where          is the number of incomplete observations in each stratum defined by the different levels of response and

complete covariates Zi, and nc^i,yi^ is the number of complete observations in each stratum.

As the above equation indicates, the mean score method weights each complete observation according to the total number
of observations in the same stratum.

The asymptotic variance of the mean score estimate is given by

Var(3) = ɪ 11v)

where n is total number of observations, and I is the usual information matrix. V is estimated by the matrix

n(y, Z)n1^y'z'>

Σ  nc(y,z)  Var(⅜*)

(У,*)

where Var{Sβ{y, z)) is the variance-covariance matrix of the score in each {y,z) stratum.

We can regard the second term of the variance expression as a penalty for the incompletely observed observations. Hence,
the mean score estimates will have larger variance than the estimates obtained if all observations were complete but smaller
variance than the estimates from an analysis of complete cases only.

Examples

We begin with a simulated dataset. We generated 1,000 observations of a predictor variable x from the standard normal
distribution. The response variable
y was then generated as a Bernoulli random variable with p = ехр(ж)/{1 + ехр(ж)}. A
dichotomous surrogate variable for
x. called z, was generated as one for positive x and zero otherwise.

A random subsample of 500 observations had their x value deleted (set to missing). The dataset, called sim_miss.dta is
provided with this insert as an illustration and can be analyzed using the mean score method by

. use sim-miss

. meanscor у x,first(z) second(x)
meanscore estimates


I

odd-ratio

Std. Err.

z

P>z

[957. Conf.

Interval]

cons I

1.050643

.0751759

0.690

0.490

.9131663

1.208817

X I

2.770173

.282211

10.002

0.000

2.268772

3.382384



More intriguing information

1. The Role of Immigration in Sustaining the Social Security System: A Political Economy Approach
2. The name is absent
3. Confusion and Reinforcement Learning in Experimental Public Goods Games
4. CGE modelling of the resources boom in Indonesia and Australia using TERM
5. Existentialism: a Philosophy of Hope or Despair?
6. Credit Market Competition and Capital Regulation
7. Pass-through of external shocks along the pricing chain: A panel estimation approach for the euro area
8. Intertemporal Risk Management Decisions of Farmers under Preference, Market, and Policy Dynamics
9. Segmentación en la era de la globalización: ¿Cómo encontrar un segmento nuevo de mercado?
10. The name is absent