Computes statistics for inferences regarding the population proportion and total given proportion data from a stratified random sample.
NINTS — Vector of length NSTRAT containing the observed number of units in each stratum from the class of interest. (Input)
NSAMPS — Vector of length NSTRAT containing the sample size in each stratum. (Input)
NPOPS — Vector of
length NSTRAT
containing the population in the strata. (Input)
If the
population strata sizes are not known, estimates must be entered in their
place.
PROPOR — Vector of length NSTRAT containing the within-strata proportion estimates. (Output)
STAT — Vector of
length 10 containing the resulting statistics. (Output)
These
are:
I |
STAT(I) |
1 |
Estimate of the proportion. |
2 |
Estimate of the total. |
3 |
Variance estimate of the proportion estimate. |
4 |
Variance estimate of the total estimate. |
5 |
Lower confidence limit for the proportion. |
6 |
Upper confidence limit for the proportion. |
7 |
Lower confidence limit for the total. |
8 |
Upper confidence limit for the total. |
9 |
Estimate (expressed as a percentage) of the coefficient of variation of the total estimate. |
10 |
Variance estimate of the proportion estimate assuming that sampling was simple random instead of stratified random. |
NSTRAT — Number
of strata into which the sample is divided. (Input)
In the
vectors of length NSTRAT, the elements
are all ordered in the same way.
Default: NSTRAT = size (NINTS,1).
CONPER —
Confidence level for two-sided interval estimate, in percent.
(Input)
A CONPER percent
confidence interval is computed; hence, CONPER must be greater
than or equal to 0.0 and less than 100.0. CONPER is often 90.0,
95.0, or 99.0. For a
one-sided confidence interval with confidence level
ONECL, set
CONPER = 100.0 − 2.0 * (100.0 − ONECL).
Default:
CONPER =
95.0.
Generic: CALL SMPPS (NINTS, NSAMPS, NPOPS, PROPOR, STAT [,…])
Specific: The specific interface names are S_SMPPS and D_SMPPS.
Single: CALL SMPPS (NSTRAT, NINTS, NSAMPS, NPOPS, CONPER, PROPOR, STAT)
Double: The double precision name is DSMPPS.
Routine SMPPS computes point and interval estimates for the population proportion and total from a stratified random sample. If the strata are formed so that the proportions differ greatly from one stratum to the next, considerable gain in statistical efficiency can be realized by use of stratified sampling (see Cochran 1977, page 107).
Let Nh be the number in the population in the h-th stratum, let nh be the number in the sample from the h-th stratum, let ah be the number of the class of interest in the sample from the h-th stratum, let N be the population size (Σ Nh), let ph be the proportion in the h-th stratum, ah/nh, and let L be the number of strata. Then, the estimate of the proportion is
and the estimate of the variance is
The confidence intervals are computed using a normal approximation.
This example is an artificial modification of an example used in routine SMPPR, which is from Cochran (1977, page 52). A list of 3042 names and addresses was built by an experienced secretary and a part-time student worker. The secretary entered 1838 names and addresses, and the student entered the remainder. Samples of size 100 were taken from the names entered by each. Verification of the addresses in the sample from the secretary’s work showed 12 to be wrong, and verification of the student’s sample showed 26 to be wrong. The objective is to estimate the total number of incorrect addresses.
USE UMACH_INT
USE SMPPS_INT
IMPLICIT NONE
INTEGER NSTRAT
PARAMETER (NSTRAT=2)
!
INTEGER NINTS(NSTRAT), NOUT, NPOPS(NSTRAT), NSAMPS(NSTRAT)
REAL CONPER, PROPOR(NSTRAT), SQRT, STAT(10), STDP, STDSRS, &
STDT
INTRINSIC SQRT
!
CALL UMACH (2, NOUT)
NINTS(1) = 12
NINTS(2) = 26
NSAMPS(1) = 100
NSAMPS(2) = 100
NPOPS(1) = 1838
NPOPS(2) = 1204
CONPER = 0.0
!
CALL SMPPS (NINTS, NSAMPS, NPOPS, PROPOR, STAT, CONPER=CONPER)
!
STDP = SQRT(STAT(3))
STDT = SQRT(STAT(4))
STDSRS = SQRT(STAT(10))
!
WRITE (NOUT,99999) STAT(1), STAT(2), STDP, STDT, STAT(9), STDSRS
99999 FORMAT (' Estimate of proportion bad: ', F7.3, /, &
' Estimate of total bad: ', F4.0, /, &
' Standard deviation estimate, proportion: ', F7.3, /, &
' Standard deviation estimate, total: ', F5.1, /, &
' Coefficient of variation: ', F5.1, &
'%', /, ' Std. dev. under simple random sampling: ', &
F7.3)
END
Estimate of proportion
bad:
0.175
Estimate of total
bad:
534.
Standard deviation estimate, proportion: 0.025
Standard
deviation estimate, total:
77.4
Coefficient of
variation:
14.5%
Std. dev. under simple random sampling: 0.027
PHONE: 713.784.3131 FAX:713.781.9260 |