Research Article Volume 2 Issue 1
1Department of Health Policy and Management, Johns Hopkins Bloomberg School of Public Health, Johns Hopkins University, USA
2Data Analytics center, Kangwon National University, Korea
3Leadership Development Division, The Catholic Education Foundation, Korea
4Performance Improvement Team, Seoul St. Mary’s Hospital, The Catholic University of Korea, Korea
5Center for Patient Education, The Care Quality Research Group, USA
6Annenberg School for Communication, University of Pennsylvania, USA
Correspondence: Heon-Jae Jeong, Department of Health Policy and Management, Bloomberg School of Public Health, Johns Hopkins University, 624 North Broadway, Rm. 455, Baltimore, 21205, USA, Tel 410-955-5315, Fax 410-955-6959
Received: January 20, 2015 | Published: February 6, 2015
Citation: Jeong HJ, Jung SM, An EA, et al. Development of the Safety Attitudes Questionnaire - Korean Version (SAQ-K) and its Novel Analysis Methods for Safety Managers. Biom Biostat Int J. 2015;2(1):8-17. DOI: 10.15406/bbij.2015.02.00020
Patient safety has emerged as the key quality indicator of health care, and safety culture among health care workers (HCWs) is also highlighted as it serves a foundational role in safety improvement endeavors. This study developed and validated the Korean version of the Safety Attitudes Questionnaire (SAQ-K) and administered it to HCWs in a large hospital. Analysis results showed a satisfactory level of internal consistency for all six domains of SAQ-K and a good fit based on confirmatory factor analysis. We have also provided novel analysis methods that can be applied to analyze SAQ data even from a large hospital with multiple clinical areas. First, variance components models were developed to examine whether SAQ-K scores differed by clinical area. Likelihood-ratio tests were conducted, and a statistically significant difference was observed for each of the six domains. Then, we applied the empirical Bayes (EB) method to obtain more accurate estimates of SAQ scores for clinical areas with too much variance. Confidence intervals were significantly reduced for those areas, offering improved accuracy, in return for EB means shrinking toward the overall mean. In sum, the current study provides a valuable instrument with which to measure safety culture in Korea, as well as solid methodologies that can be readily applied in health care organizations not only in Korea but also in other countries.
Keywords: safety culture, safety attitude questionnaire, patient safety, variance components model, random effect, empirical bayes method
As the importance of patient safety has been widely acknowledged, numerous improvement programs have been designed and implemented internationally.1 These programs aimed not only to solve specific problems such as central line-associated bloodstream infection and patient falls, but also to nurture organizational and individual health care workers’ (HCWs’) safety culture.2 Although safety culture has been defined in various ways, the core idea involves “the collective beliefs and perceptions of workers regarding the organization and safety of their workplace operations”3 Safety culture has served as the foundation of safety improvement endeavors.4 Thus, there is high demand for instruments to measure safety culture among HCWs in hospitals.
In response to the demand, various tools have been developed and tested,5-7 and one of the most embraced and thoroughly validated instruments is the Safety Attitudes Questionnaire (SAQ), which is composed of six domains8 (Table 1).
SAQ Domain |
Definition |
Teamwork Climate (TC) |
Perceived quality of collaboration between personnel |
Safety Climate (SC) |
Perception of a strong and proactive organizational commitment to safety |
Job Satisfaction (JS) |
Positivity about the work experience |
Stress Recognition (SR) |
Acknowledgment of how performance is influenced by stressors |
Perception of Management (PM) |
Approval of managerial action |
Working Conditions (WC) |
Perceived quality of the work environment and logistical support |
Table 1 SAQ domain definitions
In addition to the dire need for a validated nationwide tool to measure safety culture, more advanced analytical methods of the collected safety culture data are also required. To date, most analyses of SAQ surveys have been conducted by simply calculating the mean and standard deviation of SAQ domain scores14,18,19 or the percentage of respondents reporting positive perceptions of SAQ domains for each cluster, most frequently a clinical area,15,20 and then examining patterns across clinical areas.19 From a statistical standpoint, this traditional approach leaves huge room for improvement, although its simplicity has attracted hospital safety managers with limited knowledge of statistics.
The most frequently omitted step is testing whether SAQ scores are statistically different across clusters (i.e., testing the effect of clusters in explaining the variance in SAQ before calculating summary statistics like the mean for each cluster). Only if there is a difference in SAQ scores across clusters should cluster-specific statistics be estimated; otherwise, the overall SAQ score (pooled estimate) for the entire hospital is sufficient. For ordinary cases, analysis of variance (ANOVA) can be used, but for a highly unbalanced dataset with unequal variances, such as SAQ scores over clinical areas, the assumptions of ANOVA are easily violated,21 which requires seeking a different approach.
Another issue from the above-mentioned traditional approach arises from many work areas having too small a number of HCWs to provide meaningful information about the precision of clinical area-specific statistics. Compared to hospitals in the US, hospitals in Korea are understaffed, and a lot of clinical areas even have fewer than five HCWs (e.g., outpatient clinics in a hospital). For these small areas, most dispersion measures, such as standard deviation, inter quartile range, or median absolute deviation, cannot be easily calculated, and even if calculated they do not provide much information.
Therefore, we examined methodologies that can be applied to evade the aforementioned issues. First, to test for differences in SAQ scores by cluster, we used the variance components model, which is designed to estimate such between-cluster variances (i.e., within-cluster correlations).22 Meanwhile, the empirical Bayesian (EB) method was used to obtain a more accurate estimate for each cluster by allowing small clusters to borrow information from large clusters.22-25
This study was conducted with three goals:
Development and validation of SAQ-K
With permission of the original developer of SAQ, a team composed of physicians, nurses, patient safety managers, and communication experts translated SAQ into the Korean language. Not only the language difference but also structural and cultural differences in the health care sector were addressed based on the cross-cultural survey adaptation guidelines from.26 Another researcher who did not participate in the translation, back-translated the Korean version into English to check for unintended distortions in meaning. The Korean version was then pilot tested with various health care professionals, including but not limited to physicians, nurses, pharmacists, paramedics, and managerial personnel. Two negatively worded items (item 2: “In your clinical area, it is difficult to speak up if you perceive a problem with patient care” and item 11: “In this clinical area, it is difficult to discuss errors”) of the original survey were pointed out as difficult to understand and answer properly. This might be attributed to the different language structures of English and Korean; Koreans are not accustomed to answering to negatively worded questionnaire items. Since there was no other way to properly translate these items, they were excluded, and a questionnaire composed of 34 items was developed (TC: 5 items, SC: 6 items, JS: 5 items, SR: 4 items, PM, 10 items, and WC: 4 items).
The questionnaire was administered anonymously to health care professionals working at a large metropolitan hospital in Seoul from October 2013 through November 2013. Participants were allowed to choose the modality of survey between the intranet or paper version. SAQ-K responses were measured on a 5-point Likert scale (1 = Disagree Strongly, 2 = Disagree Slightly, 3 = Neutral, 4 = Agree Slightly, 5 = Agree Strongly) and then converted into a scale of 0 to 100 with the interval of 25 as the original SAQ rubric. For each respondent, six domain scores were obtained by calculating the mean of the items in each domain. Cronbach’s alpha was calculated to examine internal consistency, and confirmatory factor analysis (CFA) was conducted to test construct validity. These analyses were performed with statistical software packages, IBM SPSS Statistics 21 and AMOS 21 (SPSS Inc., Chicago Illinois).
Testing differences in SAQ-K domain scores over clinical areas
The following steps apply to each of the six SAQ-K domains. We began by modeling the SAQ-K score of a certain domain for a person
(model 0). This model did not include random effects and is denoted as
Where is the overall hospital mean of the SAQ-K scores of all survey participants from a certain domain and is the random deviation of the SAQ-K score of a respondentfrom the overall mean . is assumed to follow a normal distribution with a mean of 0 and variance . Then we built a variance components model addressing the variance introduced by clinical area (model 1), the model that describes the score of a person in clinical area of a certain domain as
Whereis the overall hospital mean of SAQ scores for a certain domain and is a random deviation of the mean SAQ score for clinical area from. , the random effect of clinical areas, is assumed to follow a normal distribution with a mean of 0 and variance and to be independent over clinical areas. Thus, is the residual within clinical area and random deviation of from the mean of clinical area . is assumed to be normally distributed with a mean of 0 and variance and independent from both clinical areas and survey participants.
Then we tested the null hypothesis, , with the likelihood-ratio (LR) test to compare the above variance components model (model 1) to the model without the random effects (model 0). Because variance, , cannot be negative, the p-value from the LR test should be divided by 2 [22], which is usually done automatically by statistical software; henceforth, p-values in this article denote the divided p-value. If the result of the LR test rejects the null hypothesis, , then we can conclude that clinical area explains the variance in SAQ-K score with statistical significance, and therefore it is worth calculating clinical area-specific SAQ-K scores. If the LR test does not reject the null hypothesis, then the overall mean of SAQ-K for the entire organization is sufficient. For this and the rest of the sections, Stata 13.1 (Stata Corp, College Station, Texas) was used.
Obtaining more accurate estimates of cluster-level SAQ-K scores with the empirical Bayes method
The EB method was used in this study to obtain more precise estimates of the SAQ-K domain score for each clinical area. Because this method calculates the cluster estimates by pulling the ordinary mean of unit scores to the overall mean, it is also called empirical Bayes shrinkage. The amount of shrinkage for clinical area is determined by the shrinkage factor
(3)
Where
is the within-cluster variance of clinical area
. Then, the EB mean of clinical area
,
is denoted as
(4)
Where
is the mean of clinical area
and
is the overall mean score of a certain domain of SAQ-K. Therefore, if clinical area
has a wide variance and the shrinkage factor is large, then the EB estimate of the area moves closer to the overall mean accordingly. This shrinkage is prominent usually for the clusters with a small sample size.
Characteristics of respondents
A total of 1,381 questionnaires were returned. After excluding those with too many missing values, especially in clinical area variables, the main interest of this study, 1,142 questionnaires were analyzed. As depicted in (Table 2), most of the survey respondents (73.7%) were female, which is understandable in that most nurses in Korea are female and the number of nurses is much larger than the number of physicians. Health care professionals with 5 to 10 years of experience (25.4%) formed the largest group and were followed by those with 3 to 4 years (21.8%). The tails on both sides of work years accounted for the smallest proportion at 6.7% for respondents with less than 6 months of experience and 5.3% for those with 21 years or more. For job type, nurses constituted the majority at 53.3%, followed by physicians at 33.1%. Supporting staff was 11.6% of the respondents and there was a small number of pharmacists, administrators, and uncategorized HCWs.
Characteristics |
N |
% |
Gender |
||
Male |
300 |
26.3 |
Female |
842 |
73.7 |
Work years |
||
Less than 6 months |
77 |
6.7 |
7 ~ 11 months |
122 |
10.7 |
1 ~ 2 years |
193 |
16.9 |
3 ~ 4 years |
249 |
21.8 |
5 ~ 10 years |
290 |
25.4 |
11 ~ 20 years |
150 |
13.1 |
More than 21 years |
61 |
5.3 |
Job type |
||
Physician |
378 |
33.1 |
Nurse |
609 |
53.3 |
Pharmacist |
10 |
0.9 |
Supporting Staff |
132 |
11.6 |
Administration |
9 |
0.8 |
Other |
4 |
0.4 |
Total |
1142 |
100.0 |
Table 2 Characteristics of respondents
0* |
222234 |
|
|
0. |
55666666667778888899999 |
|
|
1* |
001222222344 |
|
|
1. |
5788889999 |
2* |
01233 |
|
|
2. |
55567888 |
3* |
133 |
|
|
3. |
7 |
|
|
4* |
03 |
|
|
4. |
8 |
|
|
5* |
3 |
Figure 1 Stem-and-Leaf Plots of the number of respondents by clinical area.
Domain |
Items |
β |
B |
SE |
TC |
Nurse input is well received in this clinical area |
.68 |
1.00 |
|
Disagreements in this clinical area are resolved appropriately (i.e., not who is right, but what is best for the patient) |
.76 |
1.20 |
.05 |
|
I have the support I need from other personnel to care for patients |
.74 |
1.17 |
.05 |
|
It is easy for personnel here to ask questions when there is something that they do not understand |
.70 |
1.09 |
.05 |
|
The physicians and nurses here work together as a well-coordinated team |
.68 |
1.13 |
.06 |
|
SC |
I would feel safe being treated here as a patient |
.71 |
1.00 |
|
Medical errors are handled appropriately in this clinical area |
.73 |
.99 |
.04 |
|
I know the proper channels to direct questions regarding patient safety in this clinical area |
.66 |
1.04 |
.05 |
|
I receive appropriate feedback about my performance |
.72 |
.99 |
.04 |
|
I am encouraged by my colleagues to report any patient safety concerns I may have |
.68 |
1.03 |
.05 |
|
The culture in this clinical area makes it easy to learn from the errors of others |
.62 |
.84 |
.04 |
|
JS |
I like my job |
.78 |
1.00 |
|
Working here is like being part of a family |
.82 |
1.09 |
.04 |
|
This is a good place to work |
.87 |
1.09 |
.03 |
|
I am proud to work in this clinical area |
.84 |
1.01 |
.03 |
|
Morale in this clinical area is high |
.79 |
.93 |
.04 |
|
SR |
When my workload becomes excessive, my performance is impaired |
.56 |
1.00 |
|
I am less effective at work when fatigued |
.57 |
1.10 |
.08 |
|
I am more likely to make errors in tense or hostile situations |
.71 |
1.41 |
.09 |
|
Fatigue impairs my performance during emergency situations (e.g., emergency resuscitation, seizure) |
.70 |
1.39 |
.09 |
|
PM |
Unit management supports my daily efforts |
.79 |
1.00 |
|
Hospital management supports my daily efforts |
.79 |
1.00 |
.04 |
|
Unit management doesn't knowingly compromise patient safety |
.58 |
.73 |
.04 |
|
Hospital management doesn't knowingly compromise patient safety |
.63 |
.77 |
.04 |
|
Unit management is doing a good job |
.80 |
.99 |
.03 |
|
Hospital management is doing a good job |
.77 |
.97 |
.04 |
|
Problem personnel are dealt with constructively by our unit management |
.79 |
.99 |
.03 |
|
Problem personnel are dealt with constructively by our hospital management |
.75 |
.94 |
.03 |
|
I get adequate, timely info about events that might affect my work from unit management |
.80 |
.94 |
.03 |
|
I get adequate, timely info about events that might affect my work from hospital management |
.74 |
.92 |
.03 |
|
WC |
The levels of staffing in this clinical area are sufficient to handle the number of patients |
.59 |
1.00 |
|
This hospital does a good job of training new personnel |
.73 |
1.28 |
.08 |
|
All the necessary information for diagnostic and therapeutic decisions is routinely available to me |
.73 |
1.13 |
.07 |
|
Trainees in my discipline are adequately supervised |
.78 |
1.18 |
.07 |
Table 3 Standardized and unstandardized coefficients for CFA
Β, Standardized coefficients; B, Standardized coefficient; SE, Standard error of B
Goodness-of-Fit Index |
Value |
Chi-square Test of Model Fit (df, p-value) |
1782.13 (497, p< 0.01) |
Goodness-of-Fit Index (GFI) |
0.911 |
Adjusted Goodness-of-Fit Index (AGFI) |
0.894 |
Normed Fit Index (NFI) |
0.924 |
Comparative Fit Index (CFI) |
0.944 |
Root Mean Square Error of Approximation (RMSEA) |
0.044 |
Table 4 Confirmatory factor analysis results
|
TC |
SC |
JS |
SR |
PM |
WC |
Model 0 |
||||||
64.86 |
65.32 |
58.55 |
68.66 |
61.39 |
55.23 |
|
17.60 |
16.47 |
21.26 |
18.51 |
17.23 |
16.70 |
|
Model 1 |
||||||
65.60 |
65.65 |
59.52 |
68.84 |
61.78 |
55.65 |
|
5.35* |
3.74* |
7.84* |
3.29* |
3.79* |
3.88* |
|
16.85 |
16.03 |
20.00 |
18.19 |
16.78 |
16.27 |
Table 5 SAQ-K domain scores and variance components model parameters
The results and fit indices of CFA are provided in (Table 3 & 4). To help readers better understand the CFA results, items in (Table 3) are shown in English as the original SAQ, although actual survey administration and analysis used SAQ-K. GFI, NFI, and CFI were higher than 0.9, AGFI was 0.894, and RMSEA was less than 0.05, indicating a good model fit based on commonly used criteria.29,30.
Effects of clinical area on the score of SAQ-K domains
(Table 5) shows the overall means
and standard deviations
of the six SAQ-K domain scores from model 0; the overall means
and standard deviations
of the clinical area random effect around the overall mean, and standard deviations
of residuals of SAQ-K scores around the clinical area-specific mean for each of the six domains from (model 1).
The means of both models were similar, ranging from approximately 55 (WC) to 69 (SR). Among the six domains, JS showed the largest standard deviation of random effects, at 7.84, and SR showed the smallest at 3.29. Though appeared smaller than the standard deviation of the residuals within clinical areas for all six domains, the LR test results showed that random effects were statistically significant for all six domains, suggesting that SAQ-K scores differed by clinical area and thus clinical area-specific SAQ-K scores are worth estimating.
Empirical Bayesian estimation of SAQ-K domain scores by clinical area
To obtain more accurate estimates of clinical area-specific means, we applied the empirical Bayesian method. (Figure 2) shows the means estimated for the TC domain by ordinary least square (OLS) regression (upper panel) and the EB method (lower panel). Each bar stands for the mean and 95% confidence interval (CI) of a clinical area. The numbers over the CI bars are the clinical area ID numbered by the rank of OLS means; ten clinical areas with the widest CIs by OLS were marked, and these area IDs were marked on the EB means and CIs to clearly show the changes.
For those ten clinical areas, CIs hugely shrunk and their ranks changed, indicating that area mean estimates from the EB method differed from OLS results. Note that any changes in rank were toward the center, the overall TC score mean that was around 65, as denoted in (Table 5). These findings exemplify the shrinkage idea of the EB method, where clusters with large variance borrow information from other clusters to yield more accurate estimates; in other words, CIs become narrower. In return, the means from the EB method are pulled toward the overall mean even to the level at which the rank orders change significantly.
Figure 3 shows the estimates from OLS and the EB method side by side for all six domains, helping the reader understand the effects of the EB method in obtaining more accurate estimates. Due to space limits, IDs of only seven clinical areas are displayed. Although the TC domain was shown in (Figure 2), we show it again to better compare the results.
The primary aim of this study was to develop and validate SAQ-K. The results showed that SAQ-K had acceptable internal consistency and construct validity, suggesting that it is ready for use in health care organizations across Korea. Note that two items from the original SAQ, one in the TC domain and the other in PM, were excluded in the final version of SAQ-K. These two items were negatively phrased and their Korean translation was hard for respondents to answer, which emphasizes the difficulty of cross-cultural adaptation of an instrument and the need for exceptional caution.26
One reason we want to measure safety culture is to predict future safety incidents. Indeed, many studies have shown that safety improvement programs with a component that nurtures safety culture have significantly reduced the incidence rates of harmful events.2,31 Therefore, hospital safety managers naturally ask, “Which areas of my hospital have a low safety culture?” The answer to this question guides resource allocation for safety improvement efforts across various areas of the hospital. In this sense, the more accurate the SAQ scores we obtain, the better we can allocate resources, which will eventually enable a hospital to get the most out of its safety programs.
Traditional approaches have simply calculated the mean and standard deviation of SAQ domain scores by clinical area, and more likely they only calculated the percentage of respondents with positive perceptions. Both of these approaches are limited in use by overly small sample sizes in many clinical areas. As shown in Figure 1, several clinical areas have only a few HCWs working in them and, therefore, the clinical area mean of such units cannot be estimated accurately due to large standard deviations, as shown in (Figures 2 & 3). However, in the most cases, those simple traditional approaches have prevailed, probably because safety managers more often than not have limited statistical background, especially in advanced statistics like the Bayesian approach.
This study, therefore, applied the empirical Bayes method to provide more accurate mean SAQ domain scores for each clinical area. The empirical Bayes predictor is the best linear unbiased predictor (BLUP), and in this study the EB method improves the accuracy of area-specific estimates by borrowing strengths from other clinical areas.22 In this process, the mean of each area shrinks toward the overall mean, yielding results as seen in the lower panel of (Figure 2). The advantage of the EB method is that it can be applied even when there are clusters contain very small sample sizes, like one or two,22-24,32 which is frequently the case in hospitals in Korea.
Despite the above-mentioned statistical strengths of the EB method, one might point out overlaps in the CIs of clinical areas and therefore this EB approach might not be necessary. Though seemingly logical, this doubt is easily resolved if one understands how SAQ scores are used in hospitals. Due to the scarcity of resources, a hospital can run safety improvement programs like teamwork training only for a certain number of clinical areas (e.g., ten units in a given year). Selection of those ten units is based on SAQ scores (e.g., bottom ten areas), regardless of whether CIs overlap or not. Therefore, estimating mean SAQ scores of clinical areas as accurately as possible is essential to guaranteeing the just allocation of resources. In addition, SAQ is usually measured annually, leading safety managers to examine changes in SAQ scores by clinical area over time. One might consider building a longitudinal model with the SAQ score of each HCW as the unit of analysis. However, in most cases, this can be challenging because new HCWs are joining and also many HCWs are leaving the area. This naturally leads us to use the mean SAQ score of a clinical area as the unit of analysis, where a sound longitudinal model can only be built with accurate mean SAQ scores of clinical areas.
Before we applied the EB method, we tested whether SAQ scores were significantly different among clinical areas. This step was particularly important in this study in that it guided us in judging whether EB estimates were worth applying. For a dataset with not many clusters, one-way ANOVA can be used, but this study covered 72 clinical areas which were heavily unbalanced in sample sizes and had unequal variance, as seen in (Figures 1 & 2). Therefore, we built variance components models with clinical area-level random effect and tested the difference by observing whether the random effect was statistically significant. In this study, we used the likelihood-ratio test because we were only interested in whether the random effect model was necessary, as seen in (Table 5). If the estimates of CIs for random effect are of interest, score tests (Lagrange multiplier tests) can be used instead.22,33 Employing variance components models is useful, especially when multiple clustering factors are of interest, such as clinical areas and job type of HCW. In an earlier section, how resources can be allocated over clinical areas was introduced, and the same logic can be applied to job type instead of clinical area. The problem arises when a safety manager needs to decide which clustering factors between clinical area and job type should be used for resource allocation. To illustrate, if variance in TC scores over job type is much bigger than that over clinical area, then the manager can develop a program for the job type with the lowest TC score to match other job types with higher TC score.34
The variance components model allows us to perform such comparison between variances. For example, by expanding the above variance component model, a crossed random effects model can be built on the assumption that the SAQ score of an HCW is correlated with both clinical area and job type, without assuming hierarchy between clinical area and job type. We can test whether each of the clustering criteria is significant and compare their variances. Also, we can test whether the effect of job type is nested within a clinical area, or the other way around, with hierarchical models. If there are only a few different job types in a hospital and just marginal differences in SAQ scores among job types are of interest, a generalized estimating equation (GEE) with clinical area as the clustering variable35 can be used. Though building a GEE model is relatively simple, it is not discussed in this study because it is not within the scope of the study.
Although we developed and validated SAQ-K successfully and provided methodology for its analysis, some issues remain to be addressed. First, many HCWs work in multiple areas of the hospital. For example, a physician in the Department of Anesthesiology and Clinical Care Medicine (ACCM) works in both operating rooms and intensive care units. Which clinical area should the physician choose? A more fundamental question might be which area one perceives as the culture in which one is embedded. SAQ defines one’s clinical area as the area where one ‘typically’ spends time, but this is not clear for many HCWs. Also, sufficient time is required for an HCW to be assimilated into an area’s culture. If the period from when an HCW begins to work in an area to SAQ administration is only one week, it might not be wise to think that the HCW’s responses accurately reflect the area. In hospitals where HCWs frequently change their working areas or where new trainees join regularly, this issue is amplified. These issues probably cannot be solved perfectly. Rather, what we can do is maintain consistency across the hospital. For example, a SAQ participant inclusion criterion in the Johns Hopkins Hospital is “HWCs with a 50% commitment to one clinical area for at least 4 consecutive weeks before the survey administration”17 Though still arbitrary, this offers better internal validity when a comparison of SAQ scores within a hospital is of interest.
Culture is complex issue. Measuring culture, therefore, is rightly difficult, and safety culture in hospitals is no exception. However, such difficulty cannot justify our negligence in measuring it. Indeed, without measuring safety culture, we will never be able to improve it. Thus, instruments to measure safety culture are crucial. Through this study, we developed the Korean version of such measurement tool, and we also developed methodologies to better analyze it; these methodologies can be used in other countries as well. We hope that this study helps Korean hospitals provide patients in Korea with safer care, the care that patients well deserve.
Contact the corresponding author at hjeong1@jhu.edu for SAQ-K and Stata commands for the analyses in this article.
©2015 Jeong, et al. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.
2 7