Update to a program for saving a model fit as a dataset



Stata Technical Bulletin

17


Koopman’s method

Let X and Y be two independent binomial variates based on sample sizes m and n and parameters p1 and p2, respectively.
Let
θ = pι∕p2∙ Koopman (1984) proposed a method for constructing confidence intervals for θ based on a chi-squared test.
oopman’s method has been widely used in medical research for evaluating drug efficacy and treatment effects.

Assume we test for H0 : θ = Θq against Ha : O ≠ 0o∙ For this problem there is no uniformly most powerful test as extreme
values may occur in the sample, but a chi-squared test seems a reasonable choice. The test statistic
Ug0 is then given by

ττ t ч {x - τnp1}2 (y - np2}2

θo x,y mp1 (I-Pi) np2(1 - p2)

where pɪ and p2 are the maximum likelihood estimates under the restriction θ = 0o, It can be proved that

ʌ _ θ0(m + y) + x + n- [{0o(m + ?/) + x + n}2 - ⅛θ0(m + n)(x + y)]1/2

2(m + n)

andp2 = pi∕0o.

For 0 = 1, the statistic Ug0(x,y) is the traditional Pearson chi-square. Rearranging Ug0(x,y) results in

rr ( ʌ (ж-Wi)2 ʃ m(0θ-Pι)l

This shows that under H0, Ug0{x, y) has asymptotically for m —> ∞ and n —> ∞, a chi-squared distribution with 1 degree
of freedom independent of
Θq (Bishop et al. 1977). Hence, an approximate 1 — α two-sided confidence region for θ is given by

{Uθo(x,y) < X1,1-α}

where χ2 1-α is the 1 — α fractile of the chi-squared distribution with 1 degree of freedom. Since U is a convex function of θ,
this is an asymmetric interval (0∕,0∙u), where

l(x,y) = UθJx,y} = χ2 ι-o.

and

θι < θu

As Ugl (ж, ye reduces to the usual chi-squared when 0 = 1, this interval will always agree with the chi-squared test.

Because there is no explicit expression for the inverse function of U, the values of 0/ and θu have to be solved by numerical
procedures. The main concern of the command koopman is to obtain 0/ and
θu by using repeated bisection as suggested by
Koopman (1984).

Syntax

koopman vαrxevent var.group [weihht [if ixp [in range [, level (#) ]

koopmani #x #m #y #n [, level (#) ]

koopman allows fweights.

Description

koopman computes confidence intervals for the ratio of two binomial proportions based on two independent binomially
distributed random variables using Koopman’s method. Point estimates and confidence intervals for the odds ratio are calculated.
ivent-var contains a one if the observation represents an event and 0 otherwise. group~var indicates the group to which each
observation belongs. The variable must have only two values. Observations with missing values are not used.

koopmani is the immediate form of koopman.



More intriguing information

1. The name is absent
2. The name is absent
3. Estimating the Impact of Medication on Diabetics' Diet and Lifestyle Choices
4. The name is absent
5. The name is absent
6. The name is absent
7. The name is absent
8. Higher education funding reforms in England: the distributional effects and the shifting balance of costs
9. The name is absent
10. Short Term Memory May Be the Depletion of the Readily Releasable Pool of Presynaptic Neurotransmitter Vesicles
11. Opciones de política económica en el Perú 2011-2015
12. On the Desirability of Taxing Charitable Contributions
13. he Effect of Phosphorylation on the Electron Capture Dissociation of Peptide Ions
14. Experimental Evidence of Risk Aversion in Consumer Markets: The Case of Beef Tenderness
15. BUSINESS SUCCESS: WHAT FACTORS REALLY MATTER?
16. The Tangible Contribution of R&D Spending Foreign-Owned Plants to a Host Region: a Plant Level Study of the Irish Manufacturing Sector (1980-1996)
17. Three Policies to Improve Productivity Growth in Canada
18. The name is absent
19. Globalization, Redistribution, and the Composition of Public Education Expenditures
20. Non Linear Contracting and Endogenous Buyer Power between Manufacturers and Retailers: Empirical Evidence on Food Retailing in France