Determine whether crushing tablet medication creates an aerosol exposure risk for healthcare workers, and whether the risk differs by pill crusher and crushing method.
The three pill crushers used included MAXCRUSH, Silent Knight, and SafeCrush. Two different pill containers were used with the MAXCRUSH pill crusher labeled as “cup” and “packaged”.
Tablet medications were crushed in a controlled laboratory setting. Medication preparation was divided into three phases. “Crushing” was when the pill was actively being crushed. “Transitioning” included the time immediately following crushing until the medication was poured from the container. “Pouring” was when the crushed tablet was poured from the container into a cup of water.
Breathing zone measurements were taken under controlled and uncontrolled conditions. “Controlled” involved carefully crushing and pouring the crushed medication out of the pill container and “uncontrolled” involved vigorously crushing and pouring the crushed medication out of the pill container. Twenty trials were performed for each condition.
## # A tibble: 4 × 2
## container p_value
## <chr> <dbl>
## 1 cup 0.263
## 2 packaged 0.651
## 3 safecrush 0.215
## 4 silent knight 0.578
Indicates data follows a normal distribution
## # A tibble: 4 × 2
## container p_value
## <chr> <dbl>
## 1 cup 0.646
## 2 packaged 0.0518
## 3 safecrush 0.369
## 4 silent knight 0.0872
Indicates data follows a normal distribution
##
## Bartlett test of homogeneity of variances
##
## data: log_particle_count by container
## Bartlett's K-squared = 1.597, df = 3, p-value = 0.6601
Indicates that variance of groups are equal
##
## Bartlett test of homogeneity of variances
##
## data: log_particle_count by container
## Bartlett's K-squared = 1.1122, df = 3, p-value = 0.7741
Indicates that variance of groups are equal
## Levene's Test for Homogeneity of Variance (center = median)
## Df F value Pr(>F)
## group 3 0.3387 0.7974
## 76
Indicates that variance of groups are equal
## Levene's Test for Homogeneity of Variance (center = median)
## Df F value Pr(>F)
## group 3 0.3425 0.7947
## 75
Indicates that variance of groups are equal
## Df Sum Sq Mean Sq F value Pr(>F)
## container 3 33.13 11.044 13.34 4.4e-07 ***
## Residuals 76 62.91 0.828
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Indicates that at least one group is significantly different from the others
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = log_particle_count ~ container, data = .)
##
## $container
## diff lwr upr p adj
## packaged-cup 1.6574293 0.90168833 2.4131702 0.0000010
## safecrush-cup 0.1833568 -0.57238415 0.9390977 0.9196153
## silent knight-cup 0.6988412 -0.05689971 1.4545821 0.0801888
## safecrush-packaged -1.4740725 -2.22981340 -0.7183316 0.0000130
## silent knight-packaged -0.9585880 -1.71432896 -0.2028471 0.0071538
## silent knight-safecrush 0.5154844 -0.24025649 1.2712254 0.2852687
Indicates significant difference for packaged-cup, safecrush-packaged, and silent knight-packaged
## Df Sum Sq Mean Sq F value Pr(>F)
## container 3 112.43 37.48 41.86 5.3e-16 ***
## Residuals 75 67.14 0.90
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Indicates that at least one group is significantly different from the others
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = log_particle_count ~ container, data = .)
##
## $container
## diff lwr upr p adj
## packaged-cup 3.0651508 2.2789751 3.85132647 0.0000000
## safecrush-cup 0.9719537 0.1857780 1.75812940 0.0092200
## silent knight-cup 2.3406182 1.5441653 3.13707115 0.0000000
## safecrush-packaged -2.0931971 -2.8793728 -1.30702137 0.0000000
## silent knight-packaged -0.7245326 -1.5209855 0.07192037 0.0877120
## silent knight-safecrush 1.3686645 0.5722116 2.16511744 0.0001336
Indicates significant difference for packaged-cup, safecrush-cup, silent knight-cup, safecrush-packaged, and silent knight-safecrush
## # A tibble: 2 × 2
## scenario p_value
## <chr> <dbl>
## 1 controlled 0.263
## 2 uncontrolled 0.646
Indicates data follows a normal distribution
##
## Bartlett test of homogeneity of variances
##
## data: log_particle_count by scenario
## Bartlett's K-squared = 9.3168e-06, df = 1, p-value = 0.9976
Indicates that variance of groups are equal
## Levene's Test for Homogeneity of Variance (center = median)
## Df F value Pr(>F)
## group 1 0.0103 0.9198
## 38
Indicates that variance of groups are equal
##
## Two Sample t-test
##
## data: . and cup %>% filter(scenario == "uncontrolled") %>% pull(log_particle_count)
## t = -1.4811, df = 38, p-value = 0.1468
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## -0.9009372 0.1396240
## sample estimates:
## mean of x mean of y
## 5.183078 5.563735
Indicates no significant difference between groups
## # A tibble: 2 × 2
## scenario p_value
## <chr> <dbl>
## 1 controlled 0.651
## 2 uncontrolled 0.0518
Indicates data follows a normal distribution
##
## Bartlett test of homogeneity of variances
##
## data: log_particle_count by scenario
## Bartlett's K-squared = 0.58253, df = 1, p-value = 0.4453
Indicates that variance of groups are equal
## Levene's Test for Homogeneity of Variance (center = median)
## Df F value Pr(>F)
## group 1 0.4265 0.5176
## 38
Indicates that variance of groups are equal
##
## Two Sample t-test
##
## data: . and packaged %>% filter(scenario == "uncontrolled") %>% pull(log_particle_count)
## t = -6.1104, df = 38, p-value = 4.016e-07
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## -2.380868 -1.195888
## sample estimates:
## mean of x mean of y
## 6.840507 8.628885
Indicates significant difference between groups
## # A tibble: 2 × 2
## scenario p_value
## <chr> <dbl>
## 1 controlled 0.578
## 2 uncontrolled 0.0872
Indicates data follows a normal distribution
##
## Bartlett test of homogeneity of variances
##
## data: log_particle_count by scenario
## Bartlett's K-squared = 0.027252, df = 1, p-value = 0.8689
Indicates that variance of groups are equal
## Levene's Test for Homogeneity of Variance (center = median)
## Df F value Pr(>F)
## group 1 0.025 0.8752
## 37
Indicates that variance of groups are equal
##
## Two Sample t-test
##
## data: . and silent_knight %>% filter(scenario == "uncontrolled") %>% pull(log_particle_count)
## t = -6.08, df = 37, p-value = 4.893e-07
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## -2.696417 -1.348450
## sample estimates:
## mean of x mean of y
## 5.881919 7.904353
Indicates significant difference between groups
## # A tibble: 2 × 2
## scenario p_value
## <chr> <dbl>
## 1 controlled 0.215
## 2 uncontrolled 0.369
Indicates data follows a normal distribution
##
## Bartlett test of homogeneity of variances
##
## data: log_particle_count by scenario
## Bartlett's K-squared = 0.023154, df = 1, p-value = 0.8791
Indicates that variance of groups are equal
## Levene's Test for Homogeneity of Variance (center = median)
## Df F value Pr(>F)
## group 1 0.0876 0.7689
## 38
Indicates that variance of groups are equal
##
## Two Sample t-test
##
## data: . and safe_crush %>% filter(scenario == "uncontrolled") %>% pull(log_particle_count)
## t = -3.9975, df = 38, p-value = 0.0002845
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## -1.7613738 -0.5771331
## sample estimates:
## mean of x mean of y
## 5.366435 6.535688
Indicates significant difference between groups
## # A tibble: 3 × 2
## phase p_value
## <chr> <dbl>
## 1 crushing 0.134
## 2 pouring 0.405
## 3 transitioning 0.828
Indicates data follows a normal distribution
## # A tibble: 3 × 2
## phase p_value
## <chr> <dbl>
## 1 crushing 0.388
## 2 pouring 0.417
## 3 transitioning 0.421
Indicates data follows a normal distribution
##
## Bartlett test of homogeneity of variances
##
## data: log_particle_count by phase
## Bartlett's K-squared = 15.142, df = 2, p-value = 0.0005152
Indicates that variance of groups are unequal
##
## Bartlett test of homogeneity of variances
##
## data: log_particle_count by phase
## Bartlett's K-squared = 1.9155, df = 2, p-value = 0.3837
Indicates that variance of groups are equal
## Levene's Test for Homogeneity of Variance (center = median)
## Df F value Pr(>F)
## group 2 6.5417 0.002769 **
## 57
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Indicates that variance of groups are unequal
## Levene's Test for Homogeneity of Variance (center = median)
## Df F value Pr(>F)
## group 2 1.559 0.2192
## 57
Indicates that variance of groups are equal
## [1] 3.35799e-05
Indicates that at least one group is significantly different from the others
## Comparison Z P.unadj P.adj
## 1 crushing - pouring -3.3464118 8.186472e-04 2.455942e-03
## 2 crushing - transitioning 0.9826405 3.257844e-01 9.773532e-01
## 3 pouring - transitioning 4.3290523 1.497524e-05 4.492571e-05
Indicates significant difference for crushing-pouring and pouring-transitioning
## Df Sum Sq Mean Sq F value Pr(>F)
## phase 2 25.08 12.542 8.759 0.000482 ***
## Residuals 57 81.62 1.432
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Indicates that at least one group is significantly different from the others
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = log_particle_count ~ phase, data = .)
##
## $phase
## diff lwr upr p adj
## pouring-crushing 1.0148760 0.1042766 1.9254754 0.0255023
## transitioning-crushing -0.5455612 -1.4561606 0.3650382 0.3266477
## transitioning-pouring -1.5604372 -2.4710366 -0.6498378 0.0003559
Indicates significant difference for pouring-crushing and transitioning-pouring
## # A tibble: 3 × 2
## phase p_value
## <chr> <dbl>
## 1 crushing 0.000102
## 2 pouring 0.991
## 3 transitioning 0.0540
Indicates data does not follow a normal distribution
## # A tibble: 3 × 2
## phase p_value
## <chr> <dbl>
## 1 crushing 0.0455
## 2 pouring 0.0567
## 3 transitioning 0.0120
Indicates data does not follow a normal distribution
## [1] 1.586744e-08
Indicates that at least one group is significantly different from the others
## Comparison Z P.unadj P.adj
## 1 crushing - pouring -5.804954 6.438349e-09 1.931505e-08
## 2 crushing - transitioning -1.611984 1.069654e-01 3.208961e-01
## 3 pouring - transitioning 4.192970 2.753256e-05 8.259767e-05
Indicates significant difference for crushing-pouring and pouring-transitioning
## [1] 9.122313e-08
Indicates that at least one group is significantly different from the others
## Comparison Z P.unadj P.adj
## 1 crushing - pouring -5.3011625 1.150676e-07 3.452028e-07
## 2 crushing - transitioning -0.8510833 3.947231e-01 1.000000e+00
## 3 pouring - transitioning 4.4500792 8.583862e-06 2.575159e-05
Indicates significant difference for crushing-pouring and pouring-transitioning
## # A tibble: 3 × 2
## phase p_value
## <chr> <dbl>
## 1 crushing 0.201
## 2 pouring 0.109
## 3 transitioning 0.0923
Indicates data follows a normal distribution
## # A tibble: 3 × 2
## phase p_value
## <chr> <dbl>
## 1 crushing 0.00196
## 2 pouring 0.230
## 3 transitioning 0.0596
Indicates data does not follow a normal distribution
##
## Bartlett test of homogeneity of variances
##
## data: log_particle_count by phase
## Bartlett's K-squared = 0.030444, df = 2, p-value = 0.9849
Indicates that variance of groups are equal
## Levene's Test for Homogeneity of Variance (center = median)
## Df F value Pr(>F)
## group 2 0.4389 0.6469
## 57
Indicates that variance of groups are equal
## Df Sum Sq Mean Sq F value Pr(>F)
## phase 2 32.1 16.050 9.026 0.000393 ***
## Residuals 57 101.3 1.778
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Indicates that at least one group is significantly different from the others
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = log_particle_count ~ phase, data = .)
##
## $phase
## diff lwr upr p adj
## pouring-crushing 1.0986051 0.08388705 2.1133232 0.0309214
## transitioning-crushing -0.6763645 -1.69108260 0.3383535 0.2522578
## transitioning-pouring -1.7749696 -2.78968772 -0.7602516 0.0002676
Indicates significant difference for pouring-crushing and transitioning-pouring
## [1] 6.450195e-08
Indicates that at least one group is significantly different from the others
## Comparison Z P.unadj P.adj
## 1 crushing - pouring -4.197828 2.694871e-05 8.084614e-05
## 2 crushing - transitioning 1.309683 1.903030e-01 5.709090e-01
## 3 pouring - transitioning 5.507511 3.639421e-08 1.091826e-07
Indicates significant difference for crushing-pouring and pouring-transitioning
## # A tibble: 3 × 2
## phase p_value
## <chr> <dbl>
## 1 crushing 0.710
## 2 pouring 0.0388
## 3 transitioning 0.0512
Indicates data does not follow a normal distribution
## # A tibble: 3 × 2
## phase p_value
## <chr> <dbl>
## 1 crushing 0.482
## 2 pouring 0.0675
## 3 transitioning 0.639
Indicates data follows a normal distribution
## [1] 2.433138e-07
Indicates that at least one group is significantly different from the others
## Comparison Z P.unadj P.adj
## 1 crushing - pouring -4.020905 5.797510e-05 1.739253e-04
## 2 crushing - transitioning 1.263325 2.064725e-01 6.194175e-01
## 3 pouring - transitioning 5.284229 1.262351e-07 3.787053e-07
Indicates significant difference for crushing-pouring and pouring-transitioning
##
## Bartlett test of homogeneity of variances
##
## data: log_particle_count by phase
## Bartlett's K-squared = 6.2474, df = 2, p-value = 0.04399
Indicates that variance of groups are equal
## Levene's Test for Homogeneity of Variance (center = median)
## Df F value Pr(>F)
## group 2 1.6982 0.1921
## 57
Indicates that variance of groups are unequal
## [1] 2.420254e-08
Indicates that at least one group is significantly different from the others
## Comparison Z P.unadj P.adj
## 1 crushing - pouring -4.545968 5.468328e-06 1.640499e-05
## 2 crushing - transitioning 1.014240 3.104683e-01 9.314048e-01
## 3 pouring - transitioning 5.560208 2.694531e-08 8.083593e-08
Indicates significant difference for crushing-pouring and pouring-transitioning
## # A tibble: 4 × 2
## container p_value
## <chr> <dbl>
## 1 cup 0.0710
## 2 packaged 0.000563
## 3 safecrush 0.00000126
## 4 silent knight 0.00000228
## # A tibble: 4 × 2
## container p_value
## <chr> <dbl>
## 1 cup 0.00709
## 2 packaged 0.197
## 3 safecrush 0.0446
## 4 silent knight 0.000848
## [1] 5.001664e-06
Indicates that at least one group is significantly different from the others
## Comparison Z P.unadj P.adj
## 1 cup - packaged -4.7391656 2.146001e-06 1.287601e-05
## 2 cup - safecrush -0.4660917 6.411499e-01 1.000000e+00
## 3 packaged - safecrush 4.2730739 1.927965e-05 1.156779e-04
## 4 cup - silent knight -1.9119964 5.587665e-02 3.352599e-01
## 5 packaged - silent knight 2.8271691 4.696152e-03 2.817691e-02
## 6 safecrush - silent knight -1.4459048 1.482039e-01 8.892235e-01
Indicates significant difference for packaged-cup, safecrush-packaged, and silent knight-packaged
## [1] 8.444565e-11
Indicates that at least one group is significantly different from the others
## Comparison Z P.unadj P.adj
## 1 cup - packaged -6.459159 1.052861e-10 6.317167e-10
## 2 cup - safecrush -2.232286 2.559610e-02 1.535766e-01
## 3 packaged - safecrush 4.226874 2.369603e-05 1.421762e-04
## 4 cup - silent knight -5.008838 5.475967e-07 3.285580e-06
## 5 packaged - silent knight 1.366974 1.716334e-01 1.000000e+00
## 6 safecrush - silent knight -2.805357 5.026083e-03 3.015650e-02
Indicates significant difference for packaged-cup, silent knight-cup, safecrush-packaged, and silent knight-safecrush. Tukey showed difference between safecrush-cup
## # A tibble: 2 × 2
## scenario p_value
## <chr> <dbl>
## 1 controlled 0.0710
## 2 uncontrolled 0.00709
##
## Wilcoxon rank sum test with continuity correction
##
## data: particle_count by scenario
## W = 145.5, p-value = 0.1441
## alternative hypothesis: true location shift is not equal to 0
Indicates no significant difference between groups
## # A tibble: 2 × 2
## scenario p_value
## <chr> <dbl>
## 1 controlled 0.000563
## 2 uncontrolled 0.197
##
## Wilcoxon rank sum exact test
##
## data: particle_count by scenario
## W = 37, p-value = 1.654e-06
## alternative hypothesis: true location shift is not equal to 0
Indicates significant difference between groups
## # A tibble: 2 × 2
## scenario p_value
## <chr> <dbl>
## 1 controlled 0.00000228
## 2 uncontrolled 0.000848
##
## Wilcoxon rank sum exact test
##
## data: particle_count by scenario
## W = 33, p-value = 1.503e-06
## alternative hypothesis: true location shift is not equal to 0
Indicates significant difference between groups
## # A tibble: 2 × 2
## scenario p_value
## <chr> <dbl>
## 1 controlled 0.00000126
## 2 uncontrolled 0.0446
##
## Wilcoxon rank sum exact test
##
## data: particle_count by scenario
## W = 76, p-value = 0.0005305
## alternative hypothesis: true location shift is not equal to 0
Indicates significant difference between groups
## # A tibble: 3 × 2
## phase p_value
## <chr> <dbl>
## 1 crushing 0.0000243
## 2 pouring 0.00550
## 3 transitioning 0.0320
## [1] 3.35799e-05
Indicates that at least one group is significantly different from the others
## Comparison Z P.unadj P.adj
## 1 crushing - pouring -3.3464118 8.186472e-04 2.455942e-03
## 2 crushing - transitioning 0.9826405 3.257844e-01 9.773532e-01
## 3 pouring - transitioning 4.3290523 1.497524e-05 4.492571e-05
Indicates significant difference for crushing-pouring and pouring-transitioning
## # A tibble: 3 × 2
## phase p_value
## <chr> <dbl>
## 1 crushing 0.00101
## 2 pouring 0.000514
## 3 transitioning 0.0000131
## [1] 0.001893469
Indicates that at least one group is significantly different from the others
## Comparison Z P.unadj P.adj
## 1 crushing - pouring -1.924419 0.0543020363 0.162906109
## 2 crushing - transitioning 1.611984 0.1069653693 0.320896108
## 3 pouring - transitioning 3.536404 0.0004056145 0.001216843
Indicates significant difference for transitioning-pouring. Tukey indicated significant difference for pouring-crushing
## # A tibble: 3 × 2
## phase p_value
## <chr> <dbl>
## 1 crushing 0.0000000216
## 2 pouring 0.000101
## 3 transitioning 0.0000000133
## [1] 1.586744e-08
## Comparison Z P.unadj P.adj
## 1 crushing - pouring -5.804954 6.438349e-09 1.931505e-08
## 2 crushing - transitioning -1.611984 1.069654e-01 3.208961e-01
## 3 pouring - transitioning 4.192970 2.753256e-05 8.259767e-05
Indicates significant difference for crushing-pouring and pouring-transitioning
## # A tibble: 3 × 2
## phase p_value
## <chr> <dbl>
## 1 crushing 0.000000399
## 2 pouring 0.0315
## 3 transitioning 0.0000000836
## [1] 9.122313e-08
Indicates that at least one group is significantly different from the others
## Comparison Z P.unadj P.adj
## 1 crushing - pouring -5.804954 6.438349e-09 1.931505e-08
## 2 crushing - transitioning -1.611984 1.069654e-01 3.208961e-01
## 3 pouring - transitioning 4.192970 2.753256e-05 8.259767e-05
Indicates significant difference for crushing-pouring and pouring-transitioning
## # A tibble: 3 × 2
## phase p_value
## <chr> <dbl>
## 1 crushing 0.0000368
## 2 pouring 0.000000180
## 3 transitioning 0.000000654
## # A tibble: 3 × 2
## phase p_value
## <chr> <dbl>
## 1 crushing 0.00000000882
## 2 pouring 0.000599
## 3 transitioning 0.0000000104
## [1] 0.001129727
Indicates that at least one group is significantly different from the others
## Comparison Z P.unadj P.adj
## 1 crushing - pouring -2.259495 0.0238526257 0.0715578771
## 2 crushing - transitioning 1.390110 0.1644954758 0.4934864274
## 3 pouring - transitioning 3.649605 0.0002626441 0.0007879323
Indicates significant difference for transitioning-pouring
## [1] 6.450195e-08
Indicates that at least one group is significantly different from the others
## Comparison Z P.unadj P.adj
## 1 crushing - pouring -4.197828 2.694871e-05 8.084614e-05
## 2 crushing - transitioning 1.309683 1.903030e-01 5.709090e-01
## 3 pouring - transitioning 5.507511 3.639421e-08 1.091826e-07
Indicates significant difference for crushing-pouring and pouring-transitioning
## # A tibble: 3 × 2
## phase p_value
## <chr> <dbl>
## 1 crushing 0.00173
## 2 pouring 0.000000117
## 3 transitioning 0.0000000243
## # A tibble: 3 × 2
## phase p_value
## <chr> <dbl>
## 1 crushing 0.00283
## 2 pouring 0.0370
## 3 transitioning 0.587
## [1] 2.433138e-07
Indicates that at least one group is significantly different from the others
## Comparison Z P.unadj P.adj
## 1 crushing - pouring -4.020905 5.797510e-05 1.739253e-04
## 2 crushing - transitioning 1.263325 2.064725e-01 6.194175e-01
## 3 pouring - transitioning 5.284229 1.262351e-07 3.787053e-07
Indicates significant difference for crushing-pouring and pouring-transitioning
## [1] 2.420254e-08
Indicates that at least one group is significantly different from the others
## Comparison Z P.unadj P.adj
## 1 crushing - pouring -4.545968 5.468328e-06 1.640499e-05
## 2 crushing - transitioning 1.014240 3.104683e-01 9.314048e-01
## 3 pouring - transitioning 5.560208 2.694531e-08 8.083593e-08
Indicates significant difference for crushing-pouring and pouring-transitioning
## # A tibble: 3 × 2
## medication p_value
## <chr> <dbl>
## 1 calcium 0.104
## 2 senokot 0.0204
## 3 tylenol 0.276
##
## Kruskal-Wallis rank sum test
##
## data: particle_per_gram by medication
## Kruskal-Wallis chi-squared = 2.2924, df = 2, p-value = 0.3178