<!DOCTYPE article
PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD with MathML3 v1.3 20210610//EN" "JATS-archivearticle1-3-mathml3.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="1.3" xml:lang="en" article-type="research-article"><?properties manuscript?><processing-meta base-tagset="archiving" mathml-version="3.0" table-model="xhtml" tagset-family="jats"><restricted-by>pmc</restricted-by></processing-meta><front><journal-meta><journal-id journal-id-type="nlm-journal-id">101530522</journal-id><journal-id journal-id-type="pubmed-jr-id">39042</journal-id><journal-id journal-id-type="nlm-ta">Proc Am Stat Assoc</journal-id><journal-id journal-id-type="iso-abbrev">Proc Am Stat Assoc</journal-id><journal-title-group><journal-title>Proceedings. American Statistical Association. Annual Meeting</journal-title></journal-title-group><issn pub-type="epub">1543-3218</issn></journal-meta><article-meta><article-id pub-id-type="pmid">32336961</article-id><article-id pub-id-type="pmc">7183738</article-id><article-id pub-id-type="manuscript">HHSPA1045618</article-id><article-categories><subj-group subj-group-type="heading"><subject>Article</subject></subj-group></article-categories><title-group><article-title>Negative Binomials Regression Model in Analysis of Wait Time at Hospital Emergency Department</article-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Cai</surname><given-names>Bill</given-names></name><xref rid="A1" ref-type="aff">1</xref></contrib><contrib contrib-type="author"><name><surname>Shimizu</surname><given-names>Iris</given-names></name><xref rid="A1" ref-type="aff">1</xref></contrib></contrib-group><aff id="A1"><label>1</label>National Center for Health Statistic, 3311 Toledo Road, Hyattsville, MD 20782</aff><pub-date pub-type="nihms-submitted"><day>12</day><month>8</month><year>2019</year></pub-date><pub-date pub-type="ppub"><year>2014</year></pub-date><pub-date pub-type="pmc-release"><day>26</day><month>4</month><year>2020</year></pub-date><elocation-id>4262</elocation-id><abstract id="ABS1"><p id="P1">Wait time is the differences between the time a patient arrives in the emergency department (ED) and the time an ED provider examines that patient. This study focuses on the development of a negative binomial model to examine factors associated with ED wait time using the National Hospital Ambulatory Medical Care Survey (NHAMCS). Conducted by National Center for Health Statistics (NCHS), NHAMCS has been gathering, analyzing, and disseminating information annually about visits made for medical care to hospital outpatient department and EDs since 1992. To analyze ED wait times, a negative binomial model was fit to the ED visit data using publically released micro data from the 2009 NHAMCS. In this model, the wait time is the dependent variable while hospital, patient, and visit characteristics are the independent variables. Wait time was collapsed into discrete values representing 15 minutes intervals. The findings are presented.</p></abstract><kwd-group><kwd>Wait Time</kwd><kwd>ED</kwd><kwd>NHAMCS</kwd><kwd>Poisson Distribution</kwd><kwd>Poisson Regression</kwd><kwd>Negative Binomial Regression</kwd><kwd>Clustered Data</kwd></kwd-group></article-meta></front><body><sec id="S1"><label>1.</label><title>Introduction</title><p id="P2">Wait time is the difference between the time a patient arrives in the emergency department (ED) and the time an ED provider examines that patient. Wait time is a problem faced by both patients and hospitals. The mean wait time in EDs increased 25%, from 46.5 minutes in 2003 to 58.1 minutes in 2009. [<xref rid="R1" ref-type="bibr">1</xref>] The time it takes a patient to see a doctor (wait time) can be critical to the patient&#x02019;s health and to a hospital&#x02019;s service quality. Clinically, prolonged ED wait times may result in protracted pain and suffering and in delays in diagnosis and treatment. Many factors, such as ED overcrowding, can affect the wait time. Overcrowding in EDs may also place patients at greater risk for medical errors. (refer to American College of Emergency Physicians web site &#x0201c;<ext-link xlink:href="http://newsroom.acep.org/index.php?s=20301&#x00026;item=29937" ext-link-type="uri">Fact Sheets</ext-link>&#x0201d;)</p><p id="P3">Many studies examining ED wait times have used data from the National Hospital Ambulatory Medical Care Survey (NHAMCS, see <ext-link xlink:href="http://www.cdc.gov/nchs/ahcd/about_ahcd.htm" ext-link-type="uri">http://www.cdc.gov/nchs/ahcd/about_ahcd.htm</ext-link>). Conducted by the National Center for Health Statistics (NCHS) annually since 1992, NHAMCS collects data about sample visits made for medical care to EDs and outpatient departments (OPDs) of a national probability sample of non-Federal, general and short-stay hospitals. The survey is a component of the National Health Care Survey, which measures health care utilization across a variety of health care providers. There are two micro-data files produced from NHAMCS and released for public use, one for OPD visits and one for ED visits. [<xref rid="R2" ref-type="bibr">2</xref>]</p><p id="P4">Wait time is a topic of many publications by staff at NCHS and researchers from different organizations. A majority of the publications use descriptive statistics or quote results from annual Emergency Department Summaries using NHAMCS, such as that released by NCHS for 2007 [<xref rid="R3" ref-type="bibr">3</xref>]. Those summaries report total numbers of visits in the following intervals: fewer than 15 minutes, 15 to 59 minutes, 1 hour but less than 2 hours, 2 hours but less than 3 hours, 3 hours but less than 4 hours, 4 hours but less than 6 hours, 6 hours or more, not seen by a physician, and missing blank. They also give a median wait time (e.g., 33 minutes in 2009 [<xref rid="R1" ref-type="bibr">1</xref>]).</p><p id="P5">To further analyze wait time using NHAMCS data, some statistical model-based methodologies were also introduced. [<xref rid="R4" ref-type="bibr">4</xref>], [<xref rid="R5" ref-type="bibr">5</xref>], [<xref rid="R6" ref-type="bibr">6</xref>]</p><p id="P6">In [<xref rid="R4" ref-type="bibr">4</xref>], multivariate linear regression methodology was applied to the NHAMCS 1997 to 2004 ED data. The change in wait time to see ED physicians was evaluated for all adults, patients diagnosed with acute myocardial infarction (AMI), and patients whose triage personnel designated as needing &#x0201c;emergent&#x0201d; attention. SURVEYFREQ and SURVEYREG procedures were used in SAS to better adjust for NHAMCS&#x02019;s complex sample design.</p><p id="P7">In [<xref rid="R5" ref-type="bibr">5</xref>], multivariate logistic regression was used to estimate the association between wait time and the patient&#x02019;s age, sex, payment status, and geography of that visit using NHAMCS 1997 to 2006 data. Wait times were analyzed for visits with diagnoses of acute pancreatitis, appendicitis, cholecystitis, and upper gastrointestinal hemorrhage (UGIH). For these analyses, the association between the patient&#x02019;s race/ethnicity and frequency of delay relative to triage assignment was evaluated by year.</p><p id="P8">In [<xref rid="R6" ref-type="bibr">6</xref>], both descriptive statistics and linear regression analysis were conducted. Data for ED visits by patients 15 years of age or less from 1997 to 2000 were examined. Sample weights were applied to the identified patient records to yield national estimates. For the purposes of that study, ED wait time was analyzed for 3 major groups, i.e., non-Hispanic white (NHW), non-Hispanic black (NHB), and Hispanic white (HW). Using a linear regression analysis with logarithmically transformed wait time as a dependent variable and with adjustment for potential confounders, including hospital location, geographic region, and payer status, the study found that NHB and HW patients waited longer than NHW patients..</p><p id="P9">This paper will examine the 2009 NHAMCS micro data on ED visits that were publicly released by NCHS. The data include hospital, patient, and visit characteristics for each visit. The sample ED visits from a hospital are treated as a cluster and hospitals are mutually independent. Wait time was collapsed into discrete values based on 15-minute intervals. The resulting frequency plots showed a Poisson distribution for the wait time. A Poisson regression model was fitted to the data with wait time as the dependent variable; hospitals as subject variables (all visits to one hospital as a cluster); patient&#x02019;s arrival time, arrival by ambulance, age, race/ethnicity, hospital location, etc. as covariates. Because of an overdispersion problem associated with the Poisson model for the data set, a negative binomial regression model was used to fit the data instead.</p><p id="P10">The authors have found that neither generalized linear models (to which both Poisson regression model and negative binomial regression model belong) nor the treatment of ED visits to each hospital as a cluster have been examined in the literature on ED wait time. The objective of this study is to find a fitting model to examine characteristics associated with ED wait time.&#x0201d;</p></sec><sec id="S2"><label>2.</label><title>Statement of the problem</title><p id="P11">Scientific questions concerning wait time for ED patient visits are broad. Most of them can be answered with descriptive statistics, but some of them require statistical inference based on a statistical model.</p><p id="P12">Examples of descriptive statistics are: what&#x02019;s the average wait time, what&#x02019;s the median wait time, what&#x02019;s the average wait time for females, at what time of the day is the wait time longest, etc. Those questions could be answered by applying SAS, R or other statistical software directly to the ED visit data. For a complex multistage survey like NHAMCS, the PROC SURVEYFREQ procedure in SAS is quite useful to produce those results. Also, the wait times published in NCHS annual Emergency Department Summaries report frequently used statistics on wait time by time intervals.</p><p id="P13">When one investigates the attribution of some independent factors on ED wait time, one needs to adjust for confounding factors to avoid the possibility of spurious relationship. For example, according to many publications [<xref rid="R4" ref-type="bibr">4</xref>], [<xref rid="R6" ref-type="bibr">6</xref>], racial minority patients have longer wait times than their white counterparts. However, a majority of racial minority patients live in metropolitan areas, and hospitals in metropolitan areas have on average longer wait times than hospitals in non-metropolitan areas [<xref rid="R1" ref-type="bibr">1</xref>], [<xref rid="R4" ref-type="bibr">4</xref>]. One should adjust for hospital location when investigating the differences in mean wait times between racial groups. Statistical modeling and statistical inferences are required to answer those kinds of questions.</p><p id="P14">To compare differences in ED wait time between population subgroups, such as those defined by age, sex, race, etc., one needs to test for significance of the difference. If the p-value of the test statistic is less than the predetermined significant level &#x003b1;, then the null-hypothesis of no difference is rejected and one concludes that there is a difference in wait time between the groups. The probability of a Type I (&#x003b1;) error for a statistical test describes the probability of incorrectly rejecting the null-hypothesis and concluding there is a difference in ED visit wait time. The smaller the &#x003b1;, the less frequently one would reject the null hypothesis. It is common practice to set the probability of a Type I error at 5%.</p><p id="P15">The statistical models and the hypothesis testing described above are used to examine the wait time in hospital EDs among different population subgroups.</p></sec><sec id="S3"><label>3.</label><title>Methodologies</title><sec id="S4"><label>I.</label><title>NHAMCS ED data</title><p id="P16">The sample for the 2009 NHAMCS consisted of 489 hospitals of which 389 were in scope and had eligible EDs. Of these, 356 EDs participated by providing data for a sample of their ED visits for a four-week period. The basic sampling unit for NHAMCS is the patient visit or encounter. During the 2009 NHAMCS, data were collected on 34,942 ED visits.</p><p id="P17">In the data file, the wait time variable is WAITTIME, which is minutes to be seen by a physician, physician assistant, or nurse practitioner after the patient&#x02019;s arrival at the ED, with values from 0 to 1440. Due to the fact that the initial 15 minutes are critical (emergent visits) in the treatment of the patients with serious illness, wait time was collapsed into discrete values based on 15-minute intervals. WAITTIME values of &#x02212;7 (&#x02018;Not applicable&#x02019;) or &#x02212;9 (&#x02018;Blank&#x02019;) were treated as missing (.). Wait time was not reported for about one quarter of the sample ED visits. The current study analyzed only the reported data and did not impute missing wait times. An analysis of potential study bias due to the missing data is needed but was not done in the current study.</p><p id="P18">The purpose of this study is to examine broadly risk factors associated with ED wait time, instead of focusing on one particular variable, such as arrival mode or race. Based on publications on ED visit wait time and their most frequently analyzed variables, the following variables were selected for this investigation:</p><table-wrap position="anchor" id="T1"><table frame="void" rules="none"><colgroup span="1"><col align="left" valign="middle" span="1"/><col align="left" valign="middle" span="1"/></colgroup><tbody><tr><td align="left" valign="top" rowspan="1" colspan="1">&#x025cb; ARRTIME:</td><td align="left" valign="top" rowspan="1" colspan="1">arrival time (military time), converted to minutes.</td></tr><tr><td rowspan="2" align="left" valign="top" colspan="1">&#x025cb; ARREMS:</td><td align="left" valign="top" rowspan="1" colspan="1">arrival by ambulance, with 1='Yes';</td></tr><tr><td align="left" valign="top" rowspan="1" colspan="1">2='No'; &#x02212;8='Unknown'; &#x02212;9='Blank'. For the current study, the categories for &#x02212;8 and &#x02212;9 were replaced by 3=&#x0201d;unknown or blank&#x02019;</td></tr><tr><td align="left" valign="top" rowspan="1" colspan="1">&#x025cb; HOSPCODE:</td><td align="left" valign="top" rowspan="1" colspan="1">hospital number (identifier), starting from 1.</td></tr><tr><td align="left" valign="top" rowspan="1" colspan="1">&#x025cb; SEX:</td><td align="left" valign="top" rowspan="1" colspan="1">patient sex, 1='Female'; 2='Male'.</td></tr><tr><td align="left" valign="top" rowspan="1" colspan="1">&#x025cb; RACER:</td><td align="left" valign="top" rowspan="1" colspan="1">patient race.</td></tr><tr><td rowspan="5" align="left" valign="top" colspan="1">&#x025cb; ETHIM:</td><td align="left" valign="top" rowspan="1" colspan="1">patient ethnicity. For the current study, the variables RACER and ETHIM were combined into</td></tr><tr><td align="left" valign="top" rowspan="1" colspan="1">1='non-Hispanic white (NHW) ',</td></tr><tr><td align="left" valign="top" rowspan="1" colspan="1">2='non-Hispanic black (NHB) ',</td></tr><tr><td align="left" valign="top" rowspan="1" colspan="1">3='Hispanic white (HW) ', and</td></tr><tr><td align="left" valign="top" rowspan="1" colspan="1">4='others'.</td></tr><tr><td align="left" valign="top" rowspan="1" colspan="1">&#x025cb; AGER:</td><td align="left" valign="top" rowspan="1" colspan="1">patient age, 1='less than 15 years old'; 2='15 &#x02013; 24 years old'; 3='25 &#x02013; 44 years old'; 4='45 &#x02013; 64 years old'; 5='65 &#x02013; 74 years old'; 6='75 years or older'.</td></tr></tbody></table></table-wrap></sec><sec id="S5"><label>II.</label><title>Model</title><p id="P19">First, all the patient visits to one hospital could be treated as repeated measurements; therefore, each hospital could be treated as a cluster. Secondly, as seen in the frequency plots for wait time shown in <xref rid="F1" ref-type="fig">Figure 1</xref> (for visits in some individual hospitals) and <xref rid="F2" ref-type="fig">Figure 2</xref> (over all visits and over all hospitals), the assumption of Poisson distribution for wait time is justified as a starting point of the analysis.</p><p id="P20">In statistics, Poisson regression is a form of regression analysis commonly used to model count data in longitudinal and clustered data analysis. Poisson regression assumes the response variable has a Poisson distribution, and assumes the logarithm of its expected value can be modeled by a linear combination of unknown parameters.</p><p id="P21">For a fixed Poisson regression model for ED wait time, with linear covariates; and no interactions, let
<disp-formula id="FD1">
<mml:math id="M1" display="block" overflow="scroll"><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo>&#x000a0;</mml:mo><mml:mo>~</mml:mo><mml:mtext>&#x000a0;Independent&#x000a0;Poisson&#x000a0;</mml:mtext><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>&#x003bb;</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mspace linebreak="newline"/><mml:mrow><mml:mtext>log</mml:mtext><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>&#x003bb;</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:msubsup><mml:mi>x</mml:mi><mml:mi>i</mml:mi><mml:mi>T</mml:mi></mml:msubsup><mml:mi>&#x003b2;</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mi>&#x003b2;</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x003b2;</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:msub><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:msub><mml:mrow/><mml:mn>1</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x003b2;</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:msub><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:msub><mml:mrow/><mml:mn>2</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x003b2;</mml:mi><mml:mn>3</mml:mn></mml:msub><mml:msub><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:msub><mml:mrow/><mml:mn>3</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x003b2;</mml:mi><mml:mn>4</mml:mn></mml:msub><mml:msub><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:msub><mml:mrow/><mml:mn>4</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x003b2;</mml:mi><mml:mn>5</mml:mn></mml:msub><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mn>5</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:math>
</disp-formula>
where,
<list list-type="simple" id="L2"><list-item><p id="P22"><italic>y</italic><sub><italic>ij</italic></sub>: wait time (1, 2, &#x02026;,) of <italic>j</italic><sup><italic>th</italic></sup> patient visit to <italic>i</italic><sup><italic>th</italic></sup> hospital ED, with E(<italic>y</italic><sub><italic>ij</italic></sub>) = <italic>&#x003bb;</italic><sub><italic>i</italic></sub></p></list-item><list-item><p id="P23"><italic>x</italic><sub><italic>i</italic>1</sub> represents patient visit to ED arrival time effects</p></list-item><list-item><p id="P24"><italic>x</italic><sub><italic>i</italic>2</sub> represents patient visit to ED arrival mode effects</p></list-item><list-item><p id="P25"><italic>x</italic><sub><italic>i</italic>3</sub> represents patient sex effects</p></list-item><list-item><p id="P26"><italic>x</italic><sub><italic>i</italic>4</sub> represents patient race/ethnic effects</p></list-item><list-item><p id="P27"><italic>x</italic><sub><italic>i</italic>5</sub> represents patient age effects</p></list-item></list></p><p id="P28">However, for a Poisson distributed variable <italic>Y</italic><sub><italic>i</italic></sub>, if Var(<italic>Y</italic><sub><italic>i</italic></sub>) &#x0003e; E(<italic>Y</italic><sub><italic>i</italic></sub>), then overdispersion is present, which is the presence of greater variability (statistical dispersion) in a data set than would be expected based on a given simple statistical model.</p><p id="P29">For the wait time data,
<disp-formula id="FD2">
<mml:math id="M2" display="block" overflow="scroll"><mml:mrow><mml:mtext>Var</mml:mtext><mml:mrow><mml:mo>(</mml:mo><mml:mi>y</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mtext>&#x000a0;</mml:mtext><mml:mo>=</mml:mo><mml:mtext>&#x000a0;</mml:mtext><mml:mn>19.44</mml:mn><mml:mtext>&#x000a0;</mml:mtext><mml:mo>&#x0003e;</mml:mo><mml:mtext>&#x000a0;</mml:mtext><mml:mn>4.21</mml:mn><mml:mtext>&#x000a0;</mml:mtext><mml:mo>=</mml:mo><mml:mtext>&#x000a0;E</mml:mtext><mml:mrow><mml:mo>(</mml:mo><mml:mi>y</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mo>&#x000a0;</mml:mo><mml:mtext>with</mml:mtext><mml:mspace width="1pt"/><mml:mi>y</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn>3</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:mtext>&#x000a0;</mml:mtext><mml:mo>&#x02026;</mml:mo><mml:mtext>and</mml:mtext><mml:mo>,</mml:mo><mml:mtext>&#x000a0;Var</mml:mtext><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mtext>&#x000a0;</mml:mtext><mml:mo>&#x0003e;</mml:mo><mml:mtext>&#x000a0;E</mml:mtext><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mtext>&#x000a0;for&#x000a0;vast&#x000a0;majority&#x000a0;of</mml:mtext><mml:mo>&#x000a0;</mml:mo><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mtext>&#x000a0;</mml:mtext><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mtext>&#x000a0;</mml:mtext><mml:mn>2</mml:mn><mml:mo>,</mml:mo><mml:mtext>&#x000a0;</mml:mtext><mml:mn>3</mml:mn><mml:mtext>&#x000a0;</mml:mtext><mml:mo>&#x02026;</mml:mo></mml:mrow></mml:math>
</disp-formula></p><p id="P30">To address the overdispersion issue, a negative binomial regression (NB) model was used instead.</p><p id="P31">If <italic>Y</italic><sub><italic>i</italic></sub> is a negative binomial distributed count with mean &#x003bc;<sub><italic>i</italic></sub> and dispersion parameter &#x003b1;, a general form of the probability mass function (pmf) of <italic>Y</italic><sub><italic>i</italic></sub> ~ NB(&#x003bc;<sub><italic>i</italic></sub>, &#x003b1;) is given by
<disp-formula id="FD3">
<mml:math id="M3" display="block" overflow="scroll"><mml:mi>f</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>,</mml:mo><mml:mo>&#x000a0;</mml:mo><mml:msub><mml:mi mathvariant="normal">&#x003bc;</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>,</mml:mo><mml:mtext>&#x000a0;</mml:mtext><mml:mi>&#x003b1;</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mtext>&#x000a0;</mml:mtext><mml:mo>=</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:mfrac><mml:mrow><mml:mi>&#x00393;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:mo>&#x000a0;</mml:mo><mml:msup><mml:mi>&#x003b1;</mml:mi><mml:mrow><mml:mo>&#x02212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:msubsup><mml:mi mathvariant="normal">&#x003bc;</mml:mi><mml:mi>i</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mo>&#x02212;</mml:mo><mml:mi>k</mml:mi></mml:mrow></mml:msubsup></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x000a0;</mml:mo><mml:msup><mml:mi>&#x003b1;</mml:mi><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:msup><mml:msubsup><mml:mi mathvariant="normal">&#x003bc;</mml:mi><mml:mi>i</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:msub><mml:mi>y</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:msubsup></mml:mrow><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>!</mml:mo><mml:mo>&#x000a0;</mml:mo><mml:mi>&#x00393;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msup><mml:mi>&#x003b1;</mml:mi><mml:mrow><mml:mo>&#x02212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:msubsup><mml:mi mathvariant="normal">&#x003bc;</mml:mi><mml:mi>i</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mo>&#x02212;</mml:mo><mml:mi>k</mml:mi></mml:mrow></mml:msubsup></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:mo>&#x000a0;</mml:mo><mml:mi>&#x003b1;</mml:mi><mml:msubsup><mml:mi mathvariant="normal">&#x003bc;</mml:mi><mml:mi>i</mml:mi><mml:mi>k</mml:mi></mml:msubsup></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:mo>&#x000a0;</mml:mo><mml:msup><mml:mi>&#x003b1;</mml:mi><mml:mrow><mml:mo>&#x02212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:msubsup><mml:mi mathvariant="normal">&#x003bc;</mml:mi><mml:mi>i</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mo>&#x02212;</mml:mo><mml:mi>k</mml:mi></mml:mrow></mml:msubsup></mml:mrow></mml:msup></mml:mrow></mml:mfrac></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mn>0</mml:mn><mml:mo>&#x000a0;</mml:mo><mml:mo>&#x000a0;</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mrow><mml:mo>,</mml:mo><mml:mtable columnalign="left"><mml:mtr columnalign="left"><mml:mtd columnalign="left"><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>&#x000a0;</mml:mo><mml:mo>=</mml:mo><mml:mo>&#x000a0;</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mo>&#x000a0;</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mo>&#x000a0;</mml:mo><mml:mo>&#x02026;</mml:mo><mml:mo>,</mml:mo><mml:mo>&#x000a0;</mml:mo><mml:mi>&#x003b1;</mml:mi><mml:mo>&#x000a0;</mml:mo><mml:mo>&#x0003e;</mml:mo><mml:mo>&#x000a0;</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:mtd></mml:mtr><mml:mtr columnalign="left"><mml:mtd columnalign="left"><mml:mrow><mml:mi>o</mml:mi><mml:mi>t</mml:mi><mml:mi>h</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi><mml:mi>w</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi><mml:mi>e</mml:mi></mml:mrow></mml:mtd></mml:mtr></mml:mtable><mml:mspace linebreak="newline"/><mml:mtext>with</mml:mtext><mml:mspace width="1pt"/><mml:mtext>E</mml:mtext><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>Y</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:msub><mml:mtext>&#x000a0;&#x003bc;</mml:mtext><mml:mi>i</mml:mi></mml:msub><mml:mspace width="1pt"/><mml:mtext>and</mml:mtext><mml:mspace width="1pt"/><mml:mtext>Var</mml:mtext><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>Y</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:msub><mml:mtext>&#x003bc;</mml:mtext><mml:mi>i</mml:mi></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:mi>&#x003b1;</mml:mi><mml:msubsup><mml:mtext>&#x003bc;</mml:mtext><mml:mi>i</mml:mi><mml:mi>k</mml:mi></mml:msubsup><mml:mo stretchy="false">)</mml:mo></mml:math>
</disp-formula></p><p id="P32">Here &#x003b1; is assumed to be a constant. The index k identifies various forms of the NB distribution, but two well-known models are given by k = 0 and 1. For k = 0 we have a linear-variance NB regression, or NB1 model, with Var(<italic>Y</italic><sub><italic>i</italic></sub>) = &#x003bc;<sub><italic>i</italic></sub>(1+&#x003b1;) [this is often approximated by fitting the constant overdispersion quasi-likelihood model with Var(<italic>Y</italic><sub><italic>i</italic></sub>) = &#x003c6;&#x003bc;<sub><italic>i</italic></sub>, where &#x003c6; is a constant]. Taking k = 1 gives the more usual quadratic-variance NB, or NB2 model, with Var(<italic>Y</italic><sub><italic>i</italic></sub>) = &#x003bc;<sub><italic>i</italic></sub>(1+&#x003b1;&#x003bc;<sub><italic>i</italic></sub>). As &#x003b1; &#x02192; 0, the NB model reduces to the Poisson model. For both models, some specific regression model is assumed for the mean, i.e. <inline-formula><mml:math id="M4" display="inline" overflow="scroll"><mml:mrow><mml:mtext>log</mml:mtext><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mtext>&#x003bc;</mml:mtext><mml:mi>i</mml:mi></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:msubsup><mml:mi>x</mml:mi><mml:mi>i</mml:mi><mml:mi>T</mml:mi></mml:msubsup><mml:mi>&#x003b2;</mml:mi></mml:mrow></mml:math></inline-formula>.</p></sec><sec id="S6"><label>III.</label><title>Model selection</title><p id="P33">The final negative binomial regression model adopted for the wait time of patient ED visits in the 2009 NHAMCS is:
<disp-formula id="FD4">
<mml:math id="M5" display="block" overflow="scroll"><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo>~</mml:mo><mml:mtext>&#x000a0;independent&#x000a0;NB</mml:mtext><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mtext>&#x003bc;</mml:mtext><mml:mi>i</mml:mi></mml:msub><mml:mo>,</mml:mo><mml:mi>&#x003b1;</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mspace linebreak="newline"/><mml:mrow><mml:mtext>log</mml:mtext><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mtext>&#x003bc;</mml:mtext><mml:mi>i</mml:mi></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:msubsup><mml:mi>x</mml:mi><mml:mi>i</mml:mi><mml:mi>T</mml:mi></mml:msubsup><mml:mi>&#x003b2;</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mi>&#x003b2;</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x003b2;</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:msub><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:msub><mml:mrow/><mml:mn>1</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x003b2;</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:msub><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:msub><mml:mrow/><mml:mn>2</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x003b2;</mml:mi><mml:mn>3</mml:mn></mml:msub><mml:msub><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:msub><mml:mrow/><mml:mn>3</mml:mn></mml:msub></mml:mrow></mml:math>
</disp-formula>
where,
<list list-type="simple" id="L4"><list-item><p id="P34">&#x003b1; = 0.55, &#x003c3;(standard error) = 0.0061;</p></list-item><list-item><p id="P35"><italic>&#x003b2;</italic><sub>0</sub> = 0.73, &#x003c3; =0.0657;</p></list-item><list-item><p id="P36"><italic>&#x003b2;</italic><sub>1</sub> = 0.18, &#x003c3; =0.0211;</p></list-item><list-item><p id="P37"><italic>&#x003b2;</italic><sub>2</sub> = 0.05, &#x003c3; =0.0139;</p></list-item><list-item><p id="P38"><italic>&#x003b2;</italic><sub>3</sub> = 0.13, &#x003c3; =0.0180.</p></list-item><list-item><p id="P39"><italic>y</italic><sub><italic>ij</italic></sub>: wait time (1, 2, &#x02026;) of <italic>j</italic><sup><italic>th</italic></sup> patient visit to <italic>i</italic><sup><italic>th</italic></sup> hospital ED, with E(<italic>y</italic><sub><italic>ij</italic></sub>) = &#x003bc;<sub><italic>i</italic></sub></p></list-item><list-item><p id="P40"><italic>x</italic><sub><italic>i</italic>1</sub> represents patient ED arrival mode effects</p></list-item><list-item><p id="P41"><italic>x</italic><sub><italic>i</italic>2</sub> represents patient sex effects</p></list-item><list-item><p id="P42"><italic>x</italic><sub><italic>i</italic>3</sub> represents patient race/ethnic effects</p></list-item></list></p></sec></sec><sec id="S7"><label>4.</label><title>Results</title><p id="P43">The negative binomial model applied in this research accounts for the properties of the NHAMCS data &#x02013; properties such as the Poisson distribution of wait time and correlated visits to each ED, which have not been investigated before to examine ED visit wait time with NHAMCS data. <xref rid="T2" ref-type="table">Table 1</xref> shows that collapsing wait time values into 15-minute units increased the mean minutes of wait time from 58.8 to 63 (=4.2*15) minutes. That is, collapsing wait times into 15 minute-based discrete values introduced imprecision into the data and, hence, the analysis results.</p><p id="P44"><xref rid="T3" ref-type="table">Table 2</xref> presents results of comparisons between visit groups defined by patient race and ethnicity and groups defined by arrival mode. While the wait time difference between blacks and Hispanic whites is not significant, the ED wait time differences between non-Hispanic whites and Hispanic whites and between non-Hispanic whites and non-Hispanic blacks are both significant. Non-Hispanic whites have the shortest mean wait time compared to non-Hispanic blacks and Hispanic whites.</p><p id="P45"><xref rid="T3" ref-type="table">Table 2</xref> also shows that the ED wait time difference between arrivals by ambulance and arrivals by non-ambulance methods is significant.</p><p id="P46">Also, from the model selection process, it cannot be concluded that patient age significantly effects wait times (not shown).</p></sec><sec id="S8"><label>5.</label><title>Summary</title><p id="P47">This study focuses on the development of a fitting model to examine factors associated with wait time experienced by patients in hospital EDs, while some previously well-known results are also investigated. The study analyzed publicly released micro data on ED visits from the 2009 NHAMCS, and the negative binomial model was chosen.</p><p id="P48">Previously known results that non-Hispanic black patients and Hispanic white patients waited longer than non-Hispanic white patients were independently confirmed in this analysis. It was also shown that wait time was less for patients who arrived by ambulance than for those who arrived by other means.</p><p id="P49">There are limitations to this study. First, ED wait time was not reported for about one quarter of the sample ED visits; therefore, some bias analysis is needed. Second, the main focus of this study is on the development of a fitting model that is used to examine factors associated with ED wait time. Not all relevant patient- and ED-related factors have been included in the study. Other research questions concerning ED wait time can be studied using the model derived in this study. For example, will the difference remain between whites and non-whites in ED wait time if the model adjusts for hospital location?</p><p id="P50">As discussed in this paper, the negative binomial regression model with associated clustered data methods is an appropriate statistical method for use in studying wait time for ED visits using NHAMCS data. The model reflects characteristics of ED wait time data, such as Poisson distribution of wait time and clustered data property of the visits while models used previously in literature, such as linear regression model and logistic model, did not.</p></sec></body><back><fn-group><fn id="FN1"><p id="P51">The findings and conclusions in this report are those of the authors and do not necessarily represent the views of the National Center for Health Statistics or the Centers for Disease Control and Prevention.</p></fn></fn-group><ref-list><title>References</title><ref id="R1"><label>[1]</label><mixed-citation publication-type="journal"><name><surname>Hing</surname><given-names>E</given-names></name>, <name><surname>Bhuiya</surname><given-names>F</given-names></name>, <year>2012</year>. <article-title>Wait Time for Treatment in Hospital Emergency Departments: 2009</article-title>. <source>NCHS Data Brief, No. 102</source>, <month>August</month>
<comment>2012.</comment></mixed-citation></ref><ref id="R2"><label>[2]</label><mixed-citation publication-type="book"><collab>National Center for Health Statistics (NCHS)</collab>. <year>2012</year>. <source>2009 National Hospital Ambulatory Medical Care Survey Micro-data File Documentation</source>. <publisher-loc>Hyattsville, Md.</publisher-loc>:<publisher-name>NCHS</publisher-name>.</mixed-citation></ref><ref id="R3"><label>[3]</label><mixed-citation publication-type="journal"><name><surname>Niska</surname><given-names>R</given-names></name>, <name><surname>Bhuiya</surname><given-names>F</given-names></name>, <name><surname>Xu</surname><given-names>J</given-names></name>, <year>2010</year>. <article-title>National Hospital Ambulatory Medical Care Survey: 2007 Emergency Department Summary</article-title>. <source>National Health Statisticals Reports, No. 26</source>.</mixed-citation></ref><ref id="R4"><label>[4]</label><mixed-citation publication-type="journal"><name><surname>Wilper</surname><given-names>A</given-names></name>, <name><surname>Woolhandler</surname><given-names>S</given-names></name>, <name><surname>Lasser</surname><given-names>K</given-names></name>, <name><surname>Mccormick</surname><given-names>D</given-names></name>, <name><surname>Cutrona</surname><given-names>S</given-names></name>, <name><surname>Bor</surname><given-names>D</given-names></name>, <name><surname>Himmelstein</surname><given-names>D</given-names></name>, <year>2008</year>. <article-title>Wait to See an Emergency Department Physician: U.S</article-title>. <source>Trends and Predictions, 1997&#x02013;2004. Health Affairs</source> (<comment>Published online January 15, 2008</comment>).</mixed-citation></ref><ref id="R5"><label>[5]</label><mixed-citation publication-type="journal"><name><surname>Wu</surname><given-names>B</given-names></name>, <name><surname>Banks</surname><given-names>P</given-names></name>, <name><surname>Conwell</surname><given-names>D</given-names></name>, <year>2009</year>. <article-title>Disparities in emergency department wait times for acute gastrointestinal illnesses: results from the National Hospital Ambulatory Medical Care Survey, 1997&#x02013;2006</article-title>. <source>Am J Gastroenterol</source>. <comment>2009</comment>
<month>May</month>
<day>12</day>.</mixed-citation></ref><ref id="R6"><label>[6]</label><mixed-citation publication-type="journal"><name><surname>James</surname><given-names>CA</given-names></name>, <name><surname>Bourgeois</surname><given-names>FT</given-names></name>, <name><surname>Shannon</surname><given-names>MW</given-names></name>. <article-title>Association of race/ethnicity with emergency department wait times</article-title>. <source>Pediatrics</source>. <year>2005</year>
<month>Mar</month>.</mixed-citation></ref></ref-list></back><floats-group><fig position="float" id="F1"><label>Figure 1.</label><caption><p id="P52">Wait Time Frequency Plot &#x02013; Visits to Hospital 1,2,3,4</p></caption><graphic xlink:href="nihms-1045618-f0001" position="float"/></fig><fig position="float" id="F2"><label>Figure 2.</label><caption><p id="P53">Wait Time Frequency Plot &#x02013; Visits to All Hospitals</p></caption><graphic xlink:href="nihms-1045618-f0002" position="float"/></fig><table-wrap position="float" id="T2"><label>Table 1.</label><caption><p id="P54">Means and standard deviation of ED wait time before and after collapsing</p></caption><table frame="void" rules="none"><colgroup span="1"><col align="left" valign="middle" span="1"/><col align="left" valign="middle" span="1"/><col align="left" valign="middle" span="1"/></colgroup><thead><tr><th align="left" valign="middle" rowspan="1" colspan="1">Variable</th><th align="left" valign="middle" rowspan="1" colspan="1">Mean</th><th align="left" valign="middle" rowspan="1" colspan="1">Std Dev</th></tr></thead><tbody><tr><td align="left" valign="middle" rowspan="1" colspan="1">ED Wait Time (unit: 1 minute)</td><td align="left" valign="middle" rowspan="1" colspan="1">58.8</td><td align="left" valign="middle" rowspan="1" colspan="1">83</td></tr><tr><td align="left" valign="middle" rowspan="1" colspan="1">ED Wait Time (unit: 15 minutes)</td><td align="left" valign="middle" rowspan="1" colspan="1">4.2 (=63min)</td><td align="left" valign="middle" rowspan="1" colspan="1">4.4 (=66min)</td></tr></tbody></table></table-wrap><table-wrap position="float" id="T3"><label>Table 2.</label><caption><p id="P55">Test of mean differences in ED wait times between selected groups of patients</p></caption><table frame="void" rules="none"><colgroup span="1"><col align="left" valign="middle" span="1"/><col align="left" valign="middle" span="1"/><col align="left" valign="middle" span="1"/><col align="left" valign="middle" span="1"/><col align="left" valign="middle" span="1"/><col align="left" valign="middle" span="1"/></colgroup><thead><tr><th rowspan="2" align="center" valign="middle" colspan="1">Label</th><th rowspan="2" align="center" valign="middle" colspan="1">EstimatedMean Difference</th><th colspan="2" align="center" valign="middle" rowspan="1">Mean</th><th rowspan="2" align="center" valign="middle" colspan="1">Chi-Square</th><th rowspan="2" align="center" valign="middle" colspan="1">P-value</th></tr><tr><th colspan="2" align="center" valign="middle" rowspan="1">Confidence Limits</th></tr></thead><tbody><tr><td align="center" valign="top" rowspan="1" colspan="1">NHB - NHW</td><td align="center" valign="middle" rowspan="1" colspan="1">0.7165</td><td align="center" valign="middle" rowspan="1" colspan="1">0.6537</td><td align="center" valign="middle" rowspan="1" colspan="1">0.7853</td><td align="center" valign="middle" rowspan="1" colspan="1">50.77</td><td align="center" valign="middle" rowspan="1" colspan="1">&#x0003c;.0001</td></tr><tr><td align="center" valign="top" rowspan="1" colspan="1">HW - NHW</td><td align="center" valign="middle" rowspan="1" colspan="1">0.7727</td><td align="center" valign="middle" rowspan="1" colspan="1">0.6941</td><td align="center" valign="middle" rowspan="1" colspan="1">0.8602</td><td align="center" valign="middle" rowspan="1" colspan="1">22.2</td><td align="center" valign="middle" rowspan="1" colspan="1">&#x0003c;.0001</td></tr><tr><td align="center" valign="top" rowspan="1" colspan="1">HW - NHB</td><td align="center" valign="middle" rowspan="1" colspan="1">1.0784</td><td align="center" valign="middle" rowspan="1" colspan="1">0.9579</td><td align="center" valign="middle" rowspan="1" colspan="1">1.2141</td><td align="center" valign="middle" rowspan="1" colspan="1">1.56</td><td align="center" valign="middle" rowspan="1" colspan="1">0.2116</td></tr><tr><td align="center" valign="top" rowspan="1" colspan="1">(arrival by non-ambulance)<break/>- (by ambulance)</td><td align="center" valign="middle" rowspan="1" colspan="1">0.8039</td><td align="center" valign="middle" rowspan="1" colspan="1">0.7562</td><td align="center" valign="middle" rowspan="1" colspan="1">0.8545</td><td align="center" valign="middle" rowspan="1" colspan="1">48.97</td><td align="center" valign="middle" rowspan="1" colspan="1">&#x0003c;.0001</td></tr></tbody></table></table-wrap></floats-group></article>