The Nationwide Inpatient Sample (NIS) is one of a family of databases and software tools developed as part of the Healthcare Cost and Utilization Project (HCUP), a FederalStateIndustry partnership sponsored by the Agency for Healthcare Research and Quality (AHRQ). The NIS is the largest nationwide allpayer hospital inpatient care database in the U.S. Each year the NIS contains data from approximately seven to eight million hospital stays – all discharge data from nearly 1,000 hospitals selected from HCUP State Inpatient Databases (SID) data.
The HCUP NIS team developed the NIS to provide analyses of hospital utilization, charges, and quality of care across the United States. This report describes the NIS sample and weights, summarizes the contents of the 2003 NIS, and discusses data analysis issues. Previous NIS releases covered 1988 through 2002. This document highlights cumulative information for all previous years to provide a longitudinal view of the database. Once again, we have enhanced the nationwide representation of the sample by incorporating data from additional HCUP State Partners. The 2003 NIS includes data from 37 states, two more than the 2002 NIS.
The NIS sampling frame included all community, nonrehabilitation hospitals in the SID that could be matched to the corresponding American Hospital Association (AHA) Annual Survey data. Based on data from 37 states, there were 3,763 hospitals in the 2003 sampling frame, a 5.4% increase from the 2002 NIS. The target universe includes all acute care discharges from nonrehabilitation, community hospitals in the United States. There were 4,836 hospitals in the target universe in 2003.
The NIS is a stratified probability sample of hospitals in the frame, with sampling probabilities calculated to select 20% of the universe contained in each stratum. The overall objective was to select a sample of hospitals representative of the target universe. With this objective in mind, we defined NIS sampling strata based on the following five hospital characteristics contained in the AHA hospital files:
After stratifying the universe of hospitals, we randomly selected up to 20% of the total number of U.S. hospitals within each stratum. If a stratum contained too few frame hospitals, then all were selected for the NIS, subject to sampling restrictions specified by states. The resulting sample for 2003 included 994 hospitals, representing 20.6% of the total hospital universe of 4,836 hospitals.
Given the increase in the number of contributing states, the NIS team evaluated and revised the sampling and weighting strategy for 1998 and subsequent data years in order to best represent the U.S. These changes included:
Also, beginning with the 1998 NIS sampling procedures, all frame hospitals within a stratum have an equal probability of selection for the sample, regardless of whether they had appeared in prior NIS samples. This deviates from the procedure used for earlier samples, which maximized the longitudinal component of the NIS series. A full description of the evaluation and revision of the NIS sampling strategy for 1998 and subsequent data years can be found in the special report, Changes in NIS Sampling and Weighting Strategy for 1998. This document is available on the 2003 NIS Documentation CDROM and on the HCUP User Support Website at http://www.hcupus.ahrq.gov/db/nation/nis/nisrelatedreports.jsp.
The 2003 NIS sampling frame included data provided by 37 HCUP State Partners. On average, 97% of the hospital universe is included in the sampling frame for all but four of these states. Two State Partners – Hawaii and South Carolina – limited the number of state hospitals included in the frame to between 70 and 85 percent. (Restrictions from other states did not have an appreciable effect on the percentage of hospitals in the sampling frame.) One State Partner, Texas, supplied data from only 73% of the state’s hospitals because some Texas hospitals, mostly small rural facilities, are exempt from statutory reporting requirements. Finally, we dropped 33 Michigan hospitals that did not report total charges from the sampling frame, leaving 70% of Michigan hospitals in the frame.
While 20% of the hospitals from each region are selected for the NIS, the comprehensiveness of the sampling frame varies by region. In the Northeast, 91.9% of hospitals are included in the sampling frame, compared with 90.4% in the Midwest (a 14% increase from 2002), 78.6% in the West (an 8% increase from 2002); and 63.1% in the South. Because the NIS sampling frame has a disproportionate representation of the more populous states and includes hospitals with more annual discharges, its comprehensiveness in terms of discharges is higher. The states in the NIS sampling frame contain 97.6% of the population in the Northeast, 99.0% in the Midwest, 92.0% in the West, and 81.3% in the South. Overall, the 2003 NIS sampling frame comprises 77.8% of all U.S. hospitals and covers 90.8% of the U.S. population.
The final 2003 sample included 7,977,728 discharges from 994 hospitals selected from all 37 frame states. Hospitals were sampled throughout each region of the United States. In the Northeast and Midwest, where a higher proportion of states are represented, relatively fewer hospitals are sampled from each state than in the South and West, where the proportion of states in the NIS is lower. Since the inception of the original 1988 NIS , its scope has expanded across several dimensions:
The additional states have enhanced the representation of the nationwide population. The 2003 NIS includes data from 37 states – two more states than the 2002 NIS and 29 more states than the original 1988 NIS. The states added to the 2003 NIS have increased the percentage of population represented in the Midwest and the West. With the addition of Indiana, the percentage of Midwest population represented in the NIS increased from 90% for 2002 to 99% for 2003. With the return of Arizona, the percentage of Western population represented in the NIS increased from 84% to 92 percent.
Ideally, relationships among outcomes and their correlates estimated from the NIS should accurately represent all U.S. hospitals. However, when creating nationwide estimates, it is advisable to check these estimates against other data sources, if available. For example, the National Hospital Discharge Survey (http://www.cdc.gov/nchs/about/major/hdasd/nhds.htm) can provide benchmarks against which to verify national estimates for hospitalizations with more than 5,000 cases.
The NIS Comparison Report assesses the accuracy of NIS estimates. The most recent report is available on the NIS Documentation CDROM and provides a comparison of a previous year’s NIS with other data sources. The updated report for the current NIS will be posted on the HCUP User Support Website (http://www.hcupus.ahrq.gov/db/nation/nis/nisrelatedreports.jsp) as soon as it is completed.
Two nonoverlapping 10% subsamples of discharges were drawn from the NIS file for several reasons pertaining to data analysis. One reason for creating the subsamples was to reduce processing costs for selected studies that will not require the entire NIS. Another reason is that the two subsamples may be used to validate models and obtain unbiased estimates of standard errors. The subsamples were selected by drawing every tenth discharge, starting with two different, randomlyselected starting points. Having a different starting point for each of the two subsamples guaranteed that the resulting subsamples would not overlap.
It is necessary to incorporate sample weights to obtain nationwide estimates. Therefore, sample weights were developed separately for hospital and dischargelevel analyses. Within a stratum, each NIS sample hospital's universe weight is equal to the number of universe hospitals it represents during the year. Since 20% of the AHA universe hospitals in each stratum are sampled when possible, the hospital weights (HOSPWT) are usually near five. The calculations for dischargelevel sampling weights (DISCWT) are similar to the calculations for hospitallevel sampling weights. In the 10% subsamples, each discharge has a 10% chance of being drawn. Therefore, the discharge weights (DISCWT10) are multiplied by 10 for each of the subsamples. Because the 10% subsamples are based on samples of discharges, each hospital is represented in the subsamples. Thus, no adjustment is required for the hospital weight when using the subsamples.
To produce nationwide estimates, the discharge weights should be used to extrapolate sampled discharges in the Core file to the discharges from all U.S. community, nonrehabilitation hospitals. For the 2000 NIS, DISCWT should be used to create nationwide estimates for all analyses except those that involve total charges, and DISCWTCHARGE should be used to create nationwide estimates of total charges. For all other years of the NIS, including the 2003 NIS, DISCWTCHARGE is not required, and DISCWT (DISCWT_U prior to the 1998 NIS) should be used to create all estimates. For a 10% subsample file, use the corresponding subsample discharge weight, DISCWT10 (D10CWT_U prior to the 1998 NIS) or DISCWTCHARGE10.
Missing Values
Missing data values can compromise the quality of estimates. If the outcome for discharges with missing values is different from the outcome for discharges with valid values, then sample estimates for that outcome will be biased and will not accurately represent the discharge population. Also, when estimating totals for nonnegative variables with missing values, sums would tend to be underestimated because the cases with missing values would be omitted from the calculations. Several techniques are available to help overcome this bias. One strategy is to impute acceptable values to replace missing values. Another strategy is to use sample weight adjustments to compensate for missing values. Such data preparation and adjustment is outside the scope of this report. However, if necessary, it should be done before analyzing data with statistical procedures.
Variance Calculations
It may be important for researchers to calculate a measure of precision for some estimates based on the NIS sample data. Variance estimates must take into account both the sampling design and the form of the statistic. Standard formulas for a stratified, singlestage cluster sample without replacement may be used to calculate statistics and their variances in most applications.
The NIS database includes a Hospital Weights file with variables required to calculate finite population statistics. In addition to the sample weights described earlier, hospital identifiers (Primary Sampling Units or PSUs), stratification variables, and stratumspecific totals for the numbers of discharges and hospitals are included so that finitepopulation corrections (FPCs) can be applied to variance estimates. Examples of the use of SAS, SUDAAN, and STATA to calculate variances in the NIS are presented in the special report: Calculating Nationwide Inpatient Sample Variances. This report is available on the NIS Documentation CDROM and on the HCUP User Support Website at www.hcupus.ahrq.gov.
Longitudinal Analyses
All frame hospitals within a stratum have an equal probability of being selected for the sample, regardless of whether they have appeared in prior NIS samples. This deviates from the procedure used for earlier samples, prior to data year 1998, which maximized the longitudinal component of the NIS series. Hospitals that continue in the NIS for multiple consecutive years are a subset of the NIS hospitals for any one of those years. Consequently, longitudinal analyses of hospitallevel outcomes may be biased if they are based on any subset of NIS hospitals limited to continuous NIS membership. The analyses may be more efficient (e.g., produce more precise estimates) if they account for the potential correlation between repeated measures on the same hospital over time.
Studying Trends
When studying trends over time using the NIS, be aware that the sampling frame for the NIS changes almost annually, i.e., more states have been added over time. Estimates from earlier years of the NIS may be subject to more sampling bias than later years of the NIS. In order to facilitate analysis of trends using multiple years of NIS data, an alternate set of NIS discharge and hospital weights for the 19881997 HCUP NIS were developed. These alternative weights were calculated in the same way as the weights for the 1998 and later years of the NIS. The special report, Using the HCUP Nationwide Inpatient Sample to Estimate Trends, includes details regarding the alternate weights and other recommendations for trends analysis. Both the NIS Trends Report and the alternative weights are available on the HCUPUS web site under Methods Series (http://www.hcupus.ahrq.gov/reports/methods/methods_topic.jsp). The NIS Trends Report is also available on the NIS Documentation CDROM.
The Nationwide Inpatient Sample (NIS) is one of a family of databases and software tools developed as part of the Healthcare Cost and Utilization Project (HCUP), a FederalStateIndustry partnership sponsored by the Agency for Healthcare Research and Quality (AHRQ). The NIS is the largest nationwide allpayer hospital inpatient care database in the U.S. Each year the NIS contains data from approximately seven to eight million hospital stays – all discharge data from nearly 1,000 hospitals selected from HCUP State Inpatient Databases (SID) data.
The HCUP NIS team developed the NIS to facilitate analyses of hospital utilization, charges, and quality of care across the United States. Potential research issues focus on both discharge and hospitallevel outcomes. Discharge outcomes of interest include trends in inpatient treatments with respect to:
Hospital outcomes of interest include:
These and other outcomes are of interest for the nation as a whole and for policyrelevant inpatient subgroups defined by geographic regions, patient demographics, hospital characteristics, physician characteristics, and pay sources.
This report describes the NIS sample and weights, summarizes the contents of the 2003 NIS, and discusses data analysis issues. The 2003 NIS includes data for calendar year 2003, while previous NIS releases covered 1988 through 2002. This document highlights cumulative information for all previous years to provide a longitudinal view of the database.
Table 1 displays the number of states, hospitals, and discharges in each year and reveals the increase in the number of participating states over time. The two additional states in the 2003 NIS have enhanced the nationwide representation of the sample, making this the most comprehensive NIS to date. Note that one state that appeared in previous years of the NIS, Maine, was unable to provide data for the 2003 NIS.
Calendar Year  States in the Frame 
Number of States  Sample Hospitals  Sample Discharges (Millions) 

1988–1992  Arizona, California, Colorado, Florida, Iowa, Illinois, Massachusetts, New Jersey, Pennsylvania, Washington, and Wisconsin  8–11  759–875  5.2–6.2 
1993  Add Connecticut, Kansas, Maryland, New York, Oregon, and South Carolina  17  913  6.5 
1994  No new additions  17  904  6.4 
1995  Add Missouri and Tennessee  19  938  6.7 
1996  No new additions  19  906  6.5 
1997  Add Georgia, Hawaii, and Utah  22  1012  7.1 
1998  No new additions  22  984  6.8 
1999  Add Maine and Virginia  24  984  7.2 
2000  Add Kentucky, North Carolina, Texas, and West Virginia  28  994  7.5 
2001  Add Michigan, Minnesota, Nebraska, Rhode Island, and Vermont  33  986  7.5 
2002  Add Nevada, Ohio, and South Dakota; Drop Arizona  35  995  7.9 
2003  Add Arizona, Indiana and New Hampshire; Drop Maine  37  994  8.0 
The hospital universe is defined as all hospitals located in the U.S. open during any part of the calendar year and designated as community hospitals in the American Hospital Association (AHA) Annual Survey. The AHA defines community hospitals as follows: “All nonfederal shortterm general and other specialty hospitals, excluding hospital units of institutions.” Consequently, Veterans Hospitals and other federal facilities (Department of Defense and Indian Health Service) are excluded. Beginning with the 1998 NIS, community rehabilitation hospitals were excluded from the universe because the type of care provided and the characteristics of the discharges from these facilities were markedly different from other shortterm hospitals. Figure 1 displays the number of universe hospitals for each year based on the AHA Annual Survey. Between the years 19882001, a steady decline in the number of hospitals is evident. However, beginning in 2002, the number of universe hospitals has stabilized.
Figure 1: Hospital Universe, by Year^{1} (text version)
All U.S. hospital entities designated as community hospitals in the AHA hospital file, except rehabilitation hospitals, were included in the hospital universe. Therefore, when two or more community hospitals merged to create a new community hospital, the original hospitals and the newlyformed hospital were all considered separate hospital entities in the universe during the year they merged. Likewise, if a community hospital split, the original hospital and all newlycreated community hospitals were treated as separate entities in the universe during the year this occurred. Finally, community hospitals that closed during a given year were included in the hospital universe, as long as they were in operation during some part of the calendar year.
Given the increase in the number of contributing states, the NIS team evaluated and revised the sampling and weighting strategy for 1998 and subsequent data years, in order to best represent the U.S. This included changes to the definitions of the strata variables, the exclusion of rehabilitation hospitals from the NIS hospital universe, and a change to the calculation of hospital universe discharges for the weights. A full description of this process can be found in the special report on Changes in NIS Sampling and Weighting Strategy for 1998. This report is available on the 2003 NIS Documentation CDROM and on the HCUP User Support Website at www.hcupus.ahrq.gov. A description of the sampling procedures and definitions of strata variables used from 1988 through 1997 can be found in the special report: Design of the HCUP Nationwide Inpatient Sample, 1997. This report is available on the 1997 NIS Documentation CDROM and on the HCUP User Support Website.
The NIS sampling strata were defined based on five hospital characteristics contained in the AHA hospital files. Beg inning with the 1998 NIS, the stratification variables were defined as follows:
Region  States 

1: Northeast  Connecticut , Maine, Massachusetts, New Hampshire, New Jersey, New York, Pennsylvania, Rhode Island, Vermont 
2: Midwest  Illinois, Indiana, Iowa, Kansas, Michigan, Minnesota, Missouri, Nebraska, North Dakota, Ohio, South Dakota, Wisconsin 
3: South  Alabama, Arkansas, Delaware, District of Columbia, Florida, Georgia, Kentucky, Louisiana, Maryland, Mississippi, North Carolina, Oklahoma, South Carolina, Tennessee, Texas, Virginia, West Virginia 
4: West  Alaska, Arizona, California, Colorado, Hawaii, Idaho, Montana, Nevada, New Mexico, Oregon, Utah, Washington, Wyoming 
We chose the bed size cutoff points so that approximately onethird of the hospitals in a given region, location, and teaching status combination would fall within each bed size category (small, medium or large). We used different cutoff points for rural, urban nonteaching, and urban teaching hospitals because hospitals in those categories tend to be small, medium, and large, respectively. For example, a mediumsized teaching hospital would be considered a rather large rural hospital. Further, the size distribution is different among regions for each of the urban/teaching categories. For example, teaching hospitals tend to be smaller in the West than they are in the South. Using differing cutoff points in this manner avoids strata containing small numbers of hospitals.
We did not split rural hospitals according to teaching status, because rural teaching hospitals were rare. For example, in 2003, rural teaching hospitals comprised only 1% of the total hospital universe. We defined the bed size categories within location and teaching status because they would otherwise have been redundant. Rural hospitals tend to be small; urban nonteaching hospitals tend to be mediumsized; and urban teaching hospitals tend to be large. Yet it was important to recognize gradations of size within these types of hospitals. For example, in serving rural discharges, the role of "large" rural hospitals (particularly rural referral centers) often differs from the role of "small" rural hospitals.
To further ensure accurate geographic representation, implicit stratification variables included state and threedigit ZIP Code (the first three digits of the hospital's fivedigit ZIP Code). Within each stratum, we sorted hospitals by threedigit ZIP Code prior to systematic random sampling.
Location and Teaching Status  Hospital Bed Size  

Small  Medium  Large  
NORTHEAST  
Rural  149  5099  100+ 
Urban, nonteaching  1124  125199  200+ 
Urban, teaching  1249  250424  425+ 
MIDWEST  
Rural  129  3049  50+ 
Urban, nonteaching  174  75174  175+ 
Urban, teaching  1249  250374  375+ 
SOUTH  
Rural  139  4074  75+ 
Urban, nonteaching  199  100199  200+ 
Urban, teaching  1249  250449  450+ 
WEST  
Rural  124  2544  45+ 
Urban, nonteaching  199  100174  175+ 
Urban, teaching  1199  200324  325+ 
The universe of hospitals was established as all community hospitals located in the U.S. with the exception, beginning in 1998, of rehabilitation hospitals. However, some hospitals do not supply data to HCUP. Therefore, we constructed the NISsampling frame from the subset of universe hospitals that released their discharge data for research use. When the 2003 sample was drawn, the Agency for Healthcare Research and Quality (AHRQ) had agreements with 37 HCUP State Partner organizations to include their data in the NIS. The number of State Partners contributing data to the NIS has expanded over the years, as shown in Table 1. As a result, the number of hospitals included in the NIS sampling frame has also increased over the years, as displayed in Figure 3.
The list of the entire frame of hospitals was composed of all AHA community hospitals in each of the frame states that could be matched to the discharge data provided to HCUP. If an AHA community hospital could not be matched to the discharge data provided by the data source, it was eliminated from the sampling frame (but not from the target universe).
Figure 3: NIS Hospital Sampling Frame, by Year (text version)
Figure 4 illustrates the number of hospitals in the universe, frame, and sample and the percentage of universe hospitals in the frame for each state in the sampling frame for 2003. In most cases, the difference between the universe and the frame represents the difference in the number of community, nonrehabilitation hospitals in the 2003 AHA Annual Survey of Hospitals and the hospitals for which data were supplied to HCUP. For example, for Connecticut, Massachusetts, Minnesota, North Carolina, and Texas, the data organization supplied fewer hospitals than report to the AHA.
The largest discrepancy between HCUP data and AHA data is in Texas. As is evident in Figure 4, only 303 out of 414 Texas community, nonrehabilitation hospitals supplied data to HCUP for 2003. Certain Texas statelicensed hospitals are exempt from statutory reporting requirements. Exempt hospitals include:
The Texas statute that exempts rural providers from the requirement to submit data defines a hospital as a rural provider if it:
These exemptions apply primarily to smaller rural public hospitals and, as a result, these facilities are less likely to be included in the sampling frame than other Texas hospitals. While the number of hospitals omitted appears sizable, those available for the NIS include 91.8% of inpatient discharges from Texas universe hospitals.
However, for Georgia, Hawaii, Indiana, Michigan, Nebraska, South Carolina, and South Dakota, we had to drop several HCUP hospitals from the frame, as described below.
The Georgia frame contains two fewer hospitals than the state universe. One hospital was excluded because of sampling restrictions stipulated by the State Partner, and one hospital identified in AHA data was not included in the data supplied to HCUP.
The Hawaii frame contains seven fewer hospitals than the state universe. Four hospitals were excluded because of sampling restrictions stipulated by the State Partner, and three hospitals identified in AHA data were not included in the data supplied to HCUP.
Similarly, the Indiana frame contains four fewer hospitals than the state universe. One hospital was excluded because of sampling restrictions stipulated by the State Partner, and three hospitals identified in AHA data were not included in the data supplied to HCUP.
The Michigan frame contains 43 fewer hospitals than the state universe. The NIS team decided to drop 33 hospitals from the frame that did not provide total charges. Our reasoning is that charges represent a critical outcome variable in the NIS. By dropping these hospitals, we avoid having to adjust the weights or create another weighting variable specifically for total charges. These hospitals are fairly evenly distributed by hospital type. There are no sampling strata in the state containing only hospitals without charges. The total charge data reported for Michigan is similar to total charge data reported by other Midwestern states. Thus, there does not seem to be an obvious bias in the type of cases for which charges are reported. The stratification and weighting scheme will adjust for the hospitals that are being dropped. In addition, 10 Michigan hospitals identified in AHA data were not included in the data supplied to HCUP.
The Nebraska frame contains five fewer hospitals than the state universe. One hospital was excluded because of sampling restrictions stipulated by the State Partner. We dropped three additional hospitals from the sampling frame because they had incomplete data and were missing a high percentage of Medicare Discharges. One hospital identified in AHA data was not included in the data supplied to HCUP.
The South Carolina frame contains nine fewer hospitals than the state universe. Seven hospitals were excluded because of sampling restrictions stipulated by South Carolina, and two hospitals identified in AHA data were not included in the data supplied to HCUP.
Likewise, the South Dakota frame contains five fewer hospitals than the South Dakota universe. Two hospitals were excluded because of sampling restrictions stipulated by South Dakota, while three hospitals identified in AHA data were not included in the data supplied to HCUP.
Figure 4: Number of Hospitals in the 2003 Universe, Frame, and Sample for Frame States (text version)
Part A: Arizona – North Carolina
Figure 4: Number of Hospitals in the 2003 Universe, Frame, and Sample for Frame States (text version)
Part B: Nebraska – West Virginia
The NIS is a stratified probability sample of hospitals in the frame, with sampling probabilities calculated to select 20% of the universe of U.S. community, nonrehabilitation hospitals contained in each stratum. This sample size was determined by AHRQ based on their experience with similar research databases. The overall design objective was to select a sample of hospitals that accurately represents the target universe, which includes hospitals outside the frame (i.e., having zero probability of selection). Moreover, this sample was to be geographically dispersed, yet drawn only from data supplied by HCUP Partners.
It should be possible, for example, to estimate DRGspecific average lengths of stay across all U.S. hospitals using weighted average lengths of stay, based on averages or regression coefficients calculated from the NIS. Ideally, relationships among outcomes and their correlates estimated from the NIS should hold across all U.S. hospitals. However, the 2003 NIS includes data from only 37 states. Therefore, it is advisable to verify your estimates against other data sources, if available. For example, the National Hospital Discharge Survey (http://www.cdc.gov/nchs/products/pubs/pubd/series/sr13/ser13.htm) can provide benchmarks against which to check your national estimates for hospitalizations with more than 5,000 cases. The NIS Comparison Report assesses the accuracy of NIS estimates. The most recent report is available on the NIS Documentation CDROM and provides a comparison of a previous year’s NIS with other data sources. The updated report for the current NIS will be posted on the HCUP User Support Website (http://www.hcupus.ahrq.gov/db/nation/nis/nisrelatedreports.jsp) as soon as it is completed.
The NIS team considered alternative stratified sampling allocation schemes. However, allocation proportional to the number of hospitals was preferred for several reasons:
After stratifying the universe of hospitals, we randomly selected up to 20% of the total number of U.S. hospitals within each stratum. If too few frame hospitals appeared in a cell, then we selected all frame hospitals for the NIS, subject to sampling restrictions specified by states. To simplify variance calculations, we drew at least two hospitals from each stratum. If fewer than two frame hospitals were available in a stratum, then we merged it with an "adjacent" cell containing hospitals with similar characteristics.
We drew a systematic random sample of hospitals from each stratum, after sorting hospitals by stratum, then by the threedigit ZIP Code (the first three digits of the hospital's fivedigit ZIP Code) within each stratum, and then by a random number within each threedigit ZIP Code. These sorts ensured further geographic generalizability of hospitals within the frame states, as well as random ordering of hospitals within threedigit ZIP Codes.
Generally, threedigit ZIP Codes that are proximal in value are geographically near one another within a state. Furthermore, the U.S. Postal Service locates regional mail distribution centers at the threedigit level. Thus, the boundaries tend to be a compromise between geographic size and population size.
We drew two nonoverlapping 10% subsamples of discharges from the NIS file for each year. The subsamples were selected by drawing every tenth discharge starting with two different starting points (randomly selected between 1 and 10). Having a different starting point for each of the two subsamples guaranteed that they would not overlap. Discharges were sampled so that 10% of each hospital's discharges in each quarter were selected for each of the subsamples. The two samples can be combined to form a single, generalizable 20% subsample of discharges.
Beginning with the 1998 NIS sampling procedures, all frame hospitals within a stratum have an equal probability of selection for the sample, regardless of whether they appeared in prior NIS samples. This deviates from the procedure used for earlier samples, which maximized the longitudinal component of the NIS series.
Further description of the sampling procedures for earlier releases of the NIS can be found in the special report: Design of the HCUP Nationwide Inpatient Sample, 1997. This report is available on the 1997 NIS Documentation CDROM and on the HCUP User Support Website at www.hcupus.ahrq.gov. For a description of the development of the new sample design for 1998 and subsequent data years, see the special report: Changes in NIS Sampling and Weighting Strategy for 1998. This report is available on the 2003 NIS Documentation CDROM and on the HCUP User Support Website at http://www.hcupus.ahrq.gov/db/nation/nis/nisrelatedreports.jsp.
Beginning with the 1993 NIS, the NIS samples no longer contain zeroweight hospitals. For a description of zeroweight hospitals in the 19881992 samples, see the special report: Design of the HCUP Nationwide Inpatient Sample, Release 1. This report is available on the 19881992 NIS Documentation CDROM.
Figure 5 depicts the numbers of hospitals sampled each year, and Figure 6 presents the numbers of discharges in each year of the NIS. For the 19881992 NIS, zeroweight hospitals were maintained to provide a longitudinal sample. Therefore, two figures exist for each of these years: one number for the regular NIS sample and another number for the total sample.
Figure 7 displays the weighted number of discharges sampled each year. Note that this number decreased from 35,408,207 in 1997 to 34,874,001 in 1998, a difference of 534,206 (1.5%). This slight decline is associated with two changes to the 1998 NIS design: the exclusion of community, rehabilitation hospitals from the hospital universe, and a change to the calculation of hospital universe discharges for the weights. Prior to 1998, we calculated discharges as the sum of total facility admissions (AHA data element ADMTOT), which includes longtermcare admissions, plus births (AHA data element BIRTHS) reported for each U.S. community hospital in the AHA Annual Survey.
Beginning in 1998, we calculate discharges as the sum of hospital admissions (AHA data element ADMH) plus births for each U.S. community, nonrehabilitation hospital. This number is more consistent with the number of discharges we receive from the state data sources. We also substitute total facility admissions if the number of hospital admissions is missing. Without these changes, the weighted number of discharges for 1998 would have been 35,622,743. The exclusion of community, rehabilitation hospitals reduced the number of universe hospitals by 177 and the number of weighted discharges by 214,490. The change in the calculation of discharges reduced the weighted number of discharges by 534,252.
Figure 5: Number of Hospitals Sampled, by Year (text version)
Figure 6: Number of NIS Discharges, Unweighted, by Year (text version)
Figure 7: Number of NIS Discharges, Weighted, by Year (text version)
Figure 8 presents a summary of the 2003 NIS hospital sample by geographic region and the number of:
For example, in 2003, the Northeast region contained 657 hospitals in the universe. It also included 604 hospitals in the frame, of which 134 were drawn for the sample. This was three more than the target sample size of 131 hospitals, resulting in a surplus of three hospitals beyond the target. The total sample exceeded the target by 27 hospitals, with a resulting sample of 20.6% of the total hospital universe. We sampled more than the target number of hospitals in each region because we rounded the target sample size for each stratum up to the next highest integer whenever it was not an integer.
Figure 9 summarizes the estimated U.S. population by geographic region on July 1, 2003^{3}. For each region, the figure reveals:
For example, the estimated population of the Northeast region on July 1, 2003 was 54,426,252. On that same date, the estimated population of states in the Northeast region that were included in the 2003 NIS was 53,117,047. This represents 97.6% of the total Northeast region’s population. The percentage of estimated U.S. population included in states in the 2003 NIS was lower in the West (92.0%) and in the South (81.3%). However, the states newly added to the 2003 NIS have substantially increased the percentage of the Midwest and West populations represented. With the addition of Indiana, the Midwest population that is represented grew from 58,308,004 in 2002 to 64,795,510 in 2003 – an increase of 11 percent. The West region experienced a similar increase; with the return of Arizona, the represented population rose from 54,817,359 to 61,128,276 – an increase of 11.5 percent. Although New Hampshire was added to the Northeast region, those gains were offset by the loss of Maine from this year’s NIS. Overall, the states in the 2003 NIS include an estimated 90.8% of the entire U.S population, up 4% from 2002.
Figure 10 depicts the number of discharges in the 2003 sample for each state. The number of sampled discharges from each state ranges from 7,141 (South Dakota) to 871,681 (California).
Figure 8: Number of Hospitals in 2003 Universe, Frame, Sample, Target, and Surplus, by Region (text version)
Figure 9: Percentage of U.S. Population in 2003 NIS States, by Region (text version)
Figure 10: Number of Discharges in the 2003 NIS, by State (text version)
To obtain nationwide estimates, we developed discharge weights using the AHA universe as the standard. These were developed separately for hospital and dischargelevel analyses. Hospitallevel weights were developed to extrapolate NIS sample hospitals to the hospital universe. Similarly, dischargelevel weights were developed to extrapolate NIS sample discharges to the discharge universe.
Hospital weights to the universe were calculated by poststratification. For each year, hospitals were stratified on the same variables that were used for sampling: geographic region, urban/rural location, teaching status, bed size, and control. The strata that were collapsed for sampling were also collapsed for sample weight calculations. Within each stratum s, each NIS sample hospital's universe weight was calculated as:
W_{s}(universe) = N_{s}(universe) ÷ N_{s}(sample)
where W_{s}(universe) was the hospital universe weight, and N_{s}(universe) and N_{s}(sample) were the number of community hospitals within stratum s in the universe and sample, respectively. Thus, each hospital's universe weight (HOSPWT) is equal to the number of universe hospitals it represents during that year. Because 20% of the hospitals in each stratum were sampled when possible, the hospital weights are usually near five.
The calculations for dischargelevel sampling weights were similar to the calculations for hospitallevel sampling weights. The discharge weights are usually constant for all discharges within a stratum. The only exceptions are for strata with sample hospitals that, according to the AHA files, were open for the entire year but contributed less than a full year of data to the NIS. For those hospitals, we adjusted the number of observed discharges by a factor of 4 ÷ Q, where Q was the number of calendar quarters for which the hospital contributed discharges to the NIS. For example, when a sample hospital contributed only two quarters of discharge data to the NIS, the adjusted number of discharges was double the observed number. This adjustment was done only for weighting purposes. The NIS data set includes only the actual (unadjusted) number of observed discharges.
With that minor adjustment, each discharge weight is essentially equal to the number of AHA universe discharges that each sampled discharge represents in its stratum. This calculation was possible because the number of total discharges was available for every hospital in the universe from the AHA files. Each universe hospital's AHA discharge total was calculated as the sum of newborns and hospital discharges.
Discharge weights to the universe were calculated by poststratification. Hospitals were stratified just as they were for universe hospital weight calculations. Within stratum s, for hospital i, each NIS sample discharge's universe weight was calculated as:
DW_{is}(universe) = [DN_{s}(universe) ÷ ADN_{s}(sample)] * (4 ÷ Q_{i})
where DW_{is}(universe) was the discharge weight; DN_{s}(universe) represented the number of discharges from community hospitals in the universe within stratum s; ADN_{s}(sample) was the number of adjusted discharges from sample hospitals selected for the NIS; and Q_{i} was the number of quarters of discharge data contributed by hospital i to the NIS (usually Q_{i} = 4). Thus, each discharge's weight (DISCWT) is equal to the number of universe discharges it represents in stratum s during that year. Because all discharges from 20% of the hospitals in each stratum were sampled when possible, the discharge weights are usually near five.
To produce nationwide estimates, use one of the following discharge weights to extrapolate discharges in the NIS Core file to the discharges from all U.S. community, nonrehabilitation hospitals. When using one of the 10% subsample files, use the subsample discharge weight (the discharge weight multiplied by 10). When using the hospital weights with the subsample files, there is no need to multiply the hospital weights because all hospitals will be represented in the subsample files. Thus, the same hospital weight (HOSPWT) can be used for the full NIS and for the subsample files.
NIS Year  Name of Discharge Weight on the Core File to Use for Creating Nationwide Estimates  Name of Discharge Weight on the 10% Subsample File to Use for Creating Nationwide Estimates 

20012003 


2000 


19981999 


19881997 


Missing data values can compromise the quality of estimates. If the outcome for discharges with missing values is different from the outcome for discharges with valid values, then sample estimates for that outcome will be biased and inaccurately represent the discharge population. There are several techniques available to help overcome this bias. One strategy is to use imputation to replace missing values with acceptable values. Another strategy is to use sample weight adjustments to compensate for missing values^{4}. Such data preparation and adjustment is beyond the scope of this report. However, if necessary, imputation or adjustments should be done before analyzing data using statistical procedures.
On the other hand, if the cases with and without missing values are assumed to be similar with respect to their outcomes, no adjustment may be necessary for estimates of means and rates. This is because the nonmissing cases would be representative of the missing cases. However, some adjustment may still be necessary for the estimates of totals. Sums of data elements containing missing values would be incomplete because cases with missing values would be omitted from the calculations.
It may be important for researchers to calculate a measure of precision for some estimates based on the NIS sample data. Variance estimates must take into account both the sampling design and the form of the statistic. The sampling design consisted of a stratified, singlestage cluster sample. A stratified random sample of hospitals (clusters) was drawn and then all discharges were included from each selected hospital.
If hospitals inside the frame are similar to hospitals outside the frame, the sample hospitals can be treated as if they were randomly selected from the entire universe of hospitals within each stratum. Standard formulas for a stratified, singlestage cluster sample without replacement could be used to calculate statistics and their variances in most applications.
A multitude of statistics can be estimated from the NIS data. Several computer programs are listed below that calculate statistics and their variances from sample survey data. Some of these programs use general methods of variance calculations (e.g., the jackknife and balanced halfsample replications) that take into account the sampling design. However, it may be desirable to calculate variances using formulas specifically developed for some statistics.
These variance calculations are based on finitesample theory, which is an appropriate method for obtaining crosssectional, nationwide estimates of outcomes. According to finitesample theory, the intent of the estimation process is to obtain estimates that are precise representations of the nationwide population at a specific point in time. In the context of the NIS, any estimates that attempt to accurately describe characteristics and interrelationships among characteristics of hospitals and discharges during a specific year should be governed by finitesample theory. Examples of this would be estimates of expenditure and utilization patterns or hospital market factors.
Alternatively, in the study of hypothetical population outcomes not limited to a specific point in time, the concept of a “superpopulation” may be useful. Analysts may be less interested in specific characteristics from the finite population (and time period) from which the sample was drawn, than they are in hypothetical characteristics of a conceptual "superpopulation" from which any particular finite population in a given year might have been drawn. According to this superpopulation model, the nationwide population in a given year is only a snapshot in time of the possible interrelationships among hospital, market, and discharge characteristics. In a given year, all possible interactions between such characteristics may not have been observed, but analysts may wish to predict or simulate interrelationships that may occur in the future.
Under the finitepopulation model, the variances of estimates approach zero as the sampling fraction approaches one. This is the case because the population is defined at that point in time, and because the estimate is for a characteristic as it existed at the time of sampling. This is in contrast to the superpopulation model, which adopts a stochastic viewpoint rather than a deterministic viewpoint. That is, the nationwide population in a particular year is viewed as a random sample of some underlying superpopulation over time. Different methods are used for calculating variances under the two sample theories. The choice of an appropriate method for calculating variances for nationwide estimates depends on the type of measure and the intent of the estimation process.
The hospital weights are useful for producing hospitallevel statistics for analyses that use the hospital as the unit of analysis, while the discharge weights are useful for producing dischargelevel statistics for analyses that use the discharge as the unit of analysis. The discharge weights may be used to estimate nationwide population statistics.
In most cases, computer programs are readily available to perform these calculations. Several statistical programming packages allow weighted analyses^{5}. For example, nearly all SAS (Statistical Analysis System) procedures incorporate weights. In addition, several statistical analysis programs have been developed to specifically calculate statistics and their standard errors from survey data. Version eight or later of SAS contains procedures (PROC SURVEYMEANS and PROC SURVEYREG) for calculating statistics based on specific sampling designs. STATA and SUDAAN are two other common statistical software packages that perform calculations for numerous statistics arising from the stratified, singlestage cluster sampling design. Examples of the use of SAS, SUDAAN, and STATA to calculate NIS variances are presented in the special report: Calculating Nationwide Inpatient Sample Variances. This report is available on the 2003 NIS Documentation CDROM and on the HCUP User Support Website at www.hcupus.ahrq.gov. For an excellent review of programs to calculate statistics from survey data, visit the following Website: http://www.hcp.med.harvard.edu/statistics/surveysoft/ .
The NIS database includes a Hospital Weights file with variables required by these programs to calculate finite population statistics. In addition to the sample weights described earlier, hospital identifiers (Primary Sampling Units or PSUs), stratification variables, and stratumspecific totals for the numbers of discharges and hospitals are included so that finitepopulation corrections (FPCs) can be applied to variance estimates.
In addition to these subroutines, standard errors can be estimated by validation and crossvalidation techniques. Given that a very large number of observations will be available for most analyses, it may be feasible to set aside a part of the data for validation purposes. Standard errors and confidence intervals can then be calculated from the validation data.
If the analytic file is too small to set aside a large validation sample, crossvalidation techniques may be used. For example, tenfold crossvalidation would split the data into ten equalsized subsets. The estimation would take place in ten iterations. In each iteration, the outcome of interest is predicted for onetenth of the observations by an estimate based on a model fit to the other ninetenths of the observations. Unbiased estimates of error variance are then obtained by comparing the actual values to the predicted values obtained in this manner.
Finally, it should be noted that a large array of hospitallevel variables are available for the entire universe of hospitals, including those outside the sampling frame. For instance, the variables from the AHA surveys and from the Medicare Cost Reports are available for nearly all hospitals. To the extent that hospitallevel outcomes correlate with these variables, they may be used to sharpen regional and nationwide estimates.
As a simple example, each hospital's number of cesarean sections would be correlated with their total number of deliveries. The figure for cesarean sections must be obtained from discharge data, but the number of deliveries is available from AHA data. Thus, if a regression model can be fit predicting cesarean sections from deliveries based on the NIS data, that regression model can then be used to obtain hospitalspecific estimates of the number of cesarean sections for all hospitals in the universe.
Hospitals that continue in the NIS for multiple consecutive years are a subset of the hospitals in the NIS for any one of those years. Consequently, longitudinal analyses of hospitallevel outcomes may be biased, if they are based on any subset of NIS hospitals limited to continuous NIS membership. In particular, such subsets would tend to contain fewer hospitals that opened, closed, split, merged, or changed strata. Further, the sample weights were developed as annual, crosssectional weights, rather than longitudinal weights. Therefore, different weights might be required, depending on the statistical methods employed by the analyst.
One approach to consider in hospitallevel longitudinal analyses is to use repeatedmeasure models that allow hospitals to have missing values for some years. However, the data are not actually missing for some hospitals, such as those that closed during the study period. In any case, the analyses may be more efficient (e.g., produce more precise estimates) if they account for the potential correlation between repeated measures on the same hospital over time, yet incorporate data from all hospitals in the sample during the study period.
When studying trends over time using the NIS, be aware that the sampling frame for the NIS changes almost annually, i.e., more states have been added over time. Estimates from earlier years of the NIS may be subject to more sampling bias than later years of the NIS. In order to facilitate analysis of trends using multiple years of NIS data, an alternate set of NIS discharge and hospital weights for the 19881997 HCUP NIS were developed. These alternative weights were calculated in the same way as the weights for the 1998 and later years of the NIS. The special report, Using the HCUP Nationwide Inpatient Sample to Estimate Trends, includes details regarding the alternate weights and other recommendations for trends analysis. Both the NIS Trends Report and the alternative weights are available on the HCUPUS web site under Methods Series (http://www.hcupus.ahrq.gov/reports/methods/methods_topic.jsp). The NIS Trends Report is also available on the NIS Documentation CDROM.
The two nonoverlapping 10% subsamples of discharges were drawn from the NIS file for each year for several reasons pertinent to data analysis. One reason for creating the subsamples was to reduce processing costs for selected studies that will not require the entire NIS. Another reason was that the two subsamples may be used to validate models and obtain unbiased estimates of standard errors. That is, one subsample may be used to estimate statistical models, while the other subsample may be used to test the fit of those models on new data. This is a very important analytical step, particularly in exploratory studies, where one runs the risk of fitting noise in the data.
For example, it is well known that the percentage of variance explained by a regression, R 2, is generally overestimated by the data used to fit a model. The regression model could be estimated from the first subsample and then applied to the second subsample. The squared correlation between the actual and predicted value in the second subsample is an unbiased estimate of the model's true explanatory power when applied to new data.
In this report, we have described the development and use of the NIS sample and weights and summarized the contents of the 2003 NIS. We have included cumulative information for all previous years to provide a longitudinal view of the database. Once again, the nationwide representation of the sample has been enhanced by incorporating data from additional HCUP State Partners, a total of 37 participants for the year 2003. We have highlighted important considerations for data analysis and have provided references to detailed reports on this subject.
/* FIRST ESTABLISH SHORTTERM
BEDS DEFINITION */
IF BDH NE
. THEN BEDTEMP = BDH ; /* SHORT TERM BEDS */
ELSE IF BDH =. THEN BEDTEMP=BDTOT ; /* TOTAL BEDS PROXY
*/
/*******************************************************/
/* NEXT ESTABLISH
TEACHING STATUS BASED ON FT & PT */
/* RESIDENT/INTERN
STATUS FOR HOSPITALS.
*/
/*******************************************************/
RESINT = (FTRES + .5*PTRES)/BEDTEMP
;
IF (MAPP3 = . AND MAPP8 =
.) THEN DO ;
IF RESINT > .10 THEN ST_TEACH = 1 ;
ELSE ST_TEACH = 0 ;
END ;
IF (MAPP3=1 OR MAPP8=1) THEN
ST_TEACH=1 ; /* 1=TEACHING */
ELSE ST_TEACH=0 ;
/* 0=NONTEACHING */
/*******************************************************/
/* CREATE TEACHING CATEGORY
VARIABLES TO FURTHER */
/* REFINE TEACHING STATUS
DEFINITION.
*/
/*******************************************************/
IF ST_TEACH = 1 THEN DO ;
IF
0 < RESINT < .15 THEN TEACHCAT=0 ; /* MINOR CATEGORY */
ELSE
IF RESINT GE .15 THEN TEACHCAT=1 ; /* MAJOR CATEGORY */
ELSE
ST_TEACH = 0 ;
/* NONTEACH STATUS*/
END ;
/*******************************************************/
/* FIRST ESTABLISH SHORTTERM
BEDS DEFINITION */
/*******************************************************/
IF BDH NE . THEN BEDTEMP =
BDH ; /* SHORT TERM BEDS */
ELSE IF BDH =. THEN BEDTEMP
= BDTOT ; /* TOTAL BEDS PROXY */
/*******************************************************/
/* ESTABLISH IRB NEEDED FOR
TEACHING STATUS
*/
/* BASED ON FT PT RESIDENT
INTERN STATUS
*/
/*******************************************************/
IRB = (FTRES + .5*PTRES) /
BEDTEMP ;
/*******************************************************/
/* CREATE TEACHING STATUS
VARIABLE */
/*******************************************************/
IF (MAPP8 EQ 1) OR (MAPP3
EQ 1) THEN HOSP_TEACH = 1 ;
ELSE IF (IRB GE 0.25) THEN
HOSP_TEACH = 1 ;
ELSE HOSP_TEACH = 0;
Internet Citation: 2003 NIS Design Report. Healthcare Cost and Utilization Project (HCUP). February 2011. Agency for Healthcare Research and Quality, Rockville, MD. www.hcupus.ahrq.gov/db/nation/nis/reports/NIS_2003_Design_Report.jsp. 
Are you having problems viewing or printing pages on this Website? 
If you have comments, suggestions, and/or questions, please contact hcup@ahrq.gov. 
Privacy Notice, Viewers & Players 
Last modified 2/1/11 