aihand: Industrial Hygiene Air Monitoring Data

In STAND: Statistical Analysis of Non-Detects

Description

Data from Mulhausen and Damiano Appendix V is used to illustrate data with "less-than" values (non-detects).

Usage

data(aihand)

Format

A data frame with 15 observations on the following 3 variables

xnd

air monitoring data with 3 non-detects

det

0 if non-detect, 1 if detect

x

air monitoring data (mg/m^3)

Details

The data in column 1 was obtained from the data in column 3 by assuming a limit of detection of 1.9. The original data in column 3 is used as an example in Appendices V, VI, and VII ( see Tables V.1, V.5, V.6, VI.2 ) to illustrate methods of analysis when there are no non-detects.

Source

Table V.2 page 244 in Mulhausen and Damiano (1998) and Table IV.3 page 349 in Ignacio and Bullock (2006)

References

Ignacio, J. S. and W. H. Bullock (2006), A Strategy for Assessing and Managing Occupational Exposures, Third Edition, AIHA Press, Fairfax, VA.

Mulhausen, J. R. and J. Damiano (1998), A Strategy for Assessing and Managing Occupational Exposures, Second Edition, AIHA Press, Fairfax, VA.

Examples

data(aihand) # Table V.2 Mulhausen and Damiano (1998)

# calculate summary statistics for non-detect data
ndt <- IH.summary(aihand,L=5)
x <- aihand$x ;  det <- rep(1,length(x))
aiha<-data.frame(x,det) #  complete data
# calculate summary statistics for complete data
cdt <- IH.summary(aiha,L=5)
# output results in a table
round(cbind(cdt,ndt),3)

#  rm(aiha,aihand,cdt,det,ndt)
# add results for exact method for complete data case
# for Xp and the exceedance fraction

Example output

Loading required package: survival
           aiha aihand
mu        0.908  0.925
se.mu     0.102  0.099
sigma     0.393  0.370
se.sigma  0.072  0.079
GM        2.479  2.522
GSD       1.482  1.447
EX        2.678  2.700
EX.LCL    2.225  2.260
EX.UCL    3.225  3.226
KM.mean   2.680  2.773
KM.LCL    2.162  2.290
KM.UCL    3.198  3.257
KM.se     0.294  0.269
Xp.obs    4.750  4.750
Xp        4.734  4.633
Xp.LCL    3.598  3.521
Xp.UCL    6.228  6.096
NpUTL        NA     NA
Maximum   5.500  5.500
NonDet%   0.000 20.000
n        15.000 15.000
Rsq       0.971  0.969
m        15.000 12.000
f         3.723  3.208
f.LCL     0.594  0.396
f.UCL    14.645 14.767
fnp       6.667  6.667
fnp.LCL   0.341  0.341
fnp.UCL  27.940 27.940
m2logL   41.807 41.304
L         5.000  5.000
p         0.950  0.950
gamma     0.950  0.950

STAND documentation built on May 2, 2019, 3:39 p.m.