Summary of the Data

The summary table, named “Table_summ_0429.csv”, can be downloaded from google drive.

View Missing Patterns

The following variables have more than 7000 records missing.

The plot below shows the missing pattern among the above variables, where the 16943 rows indicating records in our database, columns indicating variables, and cells are colored according to missingness: black = missing and grey = present. The plot, named “Rplot_miss_pattern.png,” can be downloaded from google drive.

Understand missing STC

How missing is coded in LAD and STC?

                     
                       -99  -98   15 4906610
  Chlorine               1    0    0       0
  Corn Syrup             2    0    0       0
  Not Applicable         0    3    0       0
  Potassium Hydroxide    1    0    0       0
  Unknown             8433    1    0       1
  Unknown but Loaded     0    0  657       0

The R output shows the number of incidents with missing STC along years.


1980 1981 1982 1983 1984 1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 
  71  117   94  121  166  274  259  359  501  408  521  434  383  355  304  270 
1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 
 296  323  282  220  267  265  226  244  271  263  203  242  352  227  214  193 
2012 2013 2014 2015 2016 2017 2018 
  69   91   64   45   33   36   36 

There are 8437 records with STC missing (coded as -99), most of which have INC starting with “AAR” or “DTE”. In the whole database, there are

  • 6242 records with INC starting with “AAR” (among which, 4425 with STC missing),

  • 2454 records with INC starting with “DTE” (among which, 2333 with STC missing).

INC entry that begins “DTE” (for Date To be Established) basically means that the accident’s DOA is unknown.

Mismatch between DOA and DOA.1

There are five records whose DOA and DOA.1 are different.

Mismatch between DOE and DOE.1

There are 1820 records whose DOE and DOE.1 are different (including the 1st record, “5557”, in the table above). Among them, most have DOE latter than DOE.1

These 1820 records, ordered by the difference between the two dates, are in file “subset_data_DOE_mismatch.csv” in google drive.

LS0tCnRpdGxlOiAiVENBRCBEYXRhIEVEQSIKYXV0aG9yOiAiRmVuZyBMaWFuZyIKZGF0ZTogIkFwcmlsIDI5LCAyMDIwIgpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sKLS0tCgpgYGB7ciwgZWNobz1GQUxTRSwgaW5jbHVkZT1GQUxTRX0Kb3B0aW9ucyh3aWR0aCA9IDQwMCkKc2V0d2QoIn4vRGVza3RvcC9UQ0FEX0RhdGEvIikKbGlicmFyeShsdWJyaWRhdGUpCmxpYnJhcnkoZ2dwbG90MikgICAgICAgICAgICAgICAgICAgICAgIApsaWJyYXJ5KG5hbmlhcikKZGF0YSA9IHJlYWQuY3N2KCJNYXN0ZXJfVENBRF9GaW5hbF8yMDIwMDQyNl9GaWx0ZXJlZC5jc3YiLCBoZWFkZXI9VFJVRSkKZGF0YSA9IGRhdGFbLCAxOjkxXQp2YXJJbmZvID0gcmVhZC5jc3YoIlRDQURfVmFyX0luZm8uY3N2IiwgaGVhZGVyPVRSVUUsIAogICAgICAgICAgICAgICAgICAgY29sQ2xhc3NlcyA9ICJjaGFyYWN0ZXIiKQp2YXJJbmZvJFR5cGVbdmFySW5mbyRUeXBlPT0iQ2F0ZWdvcmljYWwiXSA9ICJDYXQiCnZhckluZm8kVHlwZVt2YXJJbmZvJFR5cGU9PSJJbnRlZ2VyIl0gPSAiSW50Igp2YXJJbmZvJFR5cGVbdmFySW5mbyRUeXBlPT0iTnVtZXJpY2FsIl0gPSAiTnVtIgp2YXJJbmZvJFR5cGVbdmFySW5mbyRUeXBlPT0iSW5kZXgiXSA9ICJJbmQiCgp2YXJJbmZvJFZhck5hbWVbOTldID0gIktpbGxlZCAoTGFkaW5nIExvc3Mgb25seSkgSGF6bWF0IENhcnMgb25seSIKdmFySW5mbyRWYXJOYW1lWzEwMF0gPSAiSW5qdXJlZCAoTGFkaW5nIExvc3Mgb25seSkgSGF6bWF0IENhcnMgb25seSIKdmFySW5mbyRWYXJOYW1lWzEwMV0gPSAiRXZhY3VhdGVkIChMYWRpbmcgTG9zcyBvbmx5KSBIYXptYXQgQ2FycyBvbmx5IgoKdmFyTmFtZXMgPSBtZXJnZShkYXRhLmZyYW1lKFZhckNvZGUgPSBuYW1lcyhkYXRhKSksIHZhckluZm8sIGJ5PSJWYXJDb2RlIikKCgojIFRpbWUgdmFyaWFibGVzCiMgY2hhbmdlIHRoZSBmb3JtYXQgb2YgbW9udGgvZGF5L3llYXIKZGF0YSRET0VbZGF0YSRET0UgPT0gJyddID0gTkEgIApkYXRhJERPRS4xW2RhdGEkRE9FLjEgPT0gJzEvMC8xOTAwJ10gPSBOQQoKZGF0YSRET0EgPSBtZHkoZGF0YSRET0EpCmRhdGEkRE9FID0gbWR5KGRhdGEkRE9FKSAKZGF0YSRET0EuMSA9IG1keShkYXRhJERPQS4xKQpkYXRhJERPRS4xID0gbWR5KGRhdGEkRE9FLjEpCmBgYAoKCiMjIFN1bW1hcnkgb2YgdGhlIERhdGEKCgotIEhvdyBtYW55IG1pc3NpbmcgKHZhbHVlID0gLTk5IG9yIGJsYW5rIGVudHJ5KTogYG1pc3NpbmdfI2AgKGNvdW50cykKLSBIb3cgbWFueSB1bmlxdWUgdmFsdWVzOiBgdW5pcXVlXyNgCi0gTm90ZSBgTEFEYCBjYW4gdGFrZSB2YWx1ZXMgbGlrZSAiVW5rbm93biIgb3IgIlVua25vd24gYnV0IGxvYWRlZCIsIGFsdGhvdWdoIG5vIG1pc3NpbmcgaXMgcmVwb3J0ZWQgaW4gdGhlIHRhYmxlIGJlbG93LiAKLSBOb3RlIGBDU1QgPSAyYCBmb3IgYWxsIGluY2lkZW50cyBpbiB0aGUgZGF0YQoKVGhlIHN1bW1hcnkgdGFibGUsIG5hbWVkICJUYWJsZV9zdW1tXzA0MjkuY3N2IiwgY2FuIGJlIGRvd25sb2FkZWQgZnJvbSBnb29nbGUgZHJpdmUuIAoKYGBge3IsIGVjaG89RkFMU0V9Cm4gPSBkaW0oZGF0YSlbMV0KcCA9IGRpbShkYXRhKVsyXQpzdW1tID0gbWF0cml4KDAsIHAsIDUpCnN1bW0gPSBkYXRhLmZyYW1lKHN1bW0pCnJvd25hbWVzKHN1bW0pID0gTlVMTApuYW1lcyhzdW1tKSA9IGMoIlZhckNvZGUiLCAibWlzc2luZ18jIiwgICJ1bmlxdWVfIyIsICJWYXJOYW1lIiwgIlR5cGUiKQpzdW1tJFZhckNvZGUgPSB2YXJOYW1lcyRWYXJDb2RlCnN1bW0kVmFyTmFtZSA9IHZhck5hbWVzJFZhck5hbWUKc3VtbSRUeXBlID0gdmFyTmFtZXMkVHlwZQpmb3IoaSBpbiAxOnApewogIAogIHRtcC5kYXRhID0gZGF0YVssIHdoaWNoKG5hbWVzKGRhdGEpPT1zdW1tJFZhckNvZGVbaV0pXQogIHN1bW1baSwgIm1pc3NpbmdfIyJdID0gc3VtKHRtcC5kYXRhICVpbiUgYygtOTkpKSArIAogICAgc3VtKGlzLm5hKHRtcC5kYXRhKSkKICB0bXAgPSBzb3J0KHRhYmxlKHRtcC5kYXRhKSwgIGRlY3JlYXNpbmcgPSBUUlVFKQogIHN1bW1baSwgInVuaXF1ZV8jIl0gPSBsZW5ndGgodG1wKQp9Cgpyb3duYW1lcyhzdW1tKSA9IHN1bW1bLDFdCnByaW50KHN1bW1bLC0xXSwgaW5kaWNlcz1GQUxTRSkKd3JpdGUuY3N2KHN1bW0sIGZpbGU9IlRhYmxlX3N1bW1fMDQyOS5jc3YiKQpgYGAKCgojIyBWaWV3IE1pc3NpbmcgUGF0dGVybnMKCgpUaGUgZm9sbG93aW5nIHZhcmlhYmxlcyBoYXZlIG1vcmUgdGhhbiA3MDAwIHJlY29yZHMgbWlzc2luZy4gCgoKYGBge3IsIGVjaG89RkFMU0V9CnRtcCA9IHN1bW1bLDJdCnRtcFt0bXA9PScnXSA9IDAKdG1wID0gYXMubnVtZXJpYyh0bXApCmlkID0gd2hpY2godG1wPjcwMDApCiN0bXAgPSBkYXRhWywgYXMuY2hhcmFjdGVyKHN1bW1baWQsMV0pXQojZm9yKGkgaW4gMTpkaW0odG1wKVsyXSl7CiMgIHRtcFt3aGljaCh0bXBbLGldICVpbiUgYygtOTksIC05OCwgLTk3LCAtOTYpKSxpXSA9IE5BCiMgfQojbWlzcy5wYXR0ZXJuID0gdmlzX21pc3ModG1wLCBjbHVzdGVyPVRSVUUsIHNvcnRfbWlzcyA9IFRSVUUpCiNzYXZlKG1pc3MucGF0dGVybiwgZmlsZT0icGxvdF9taXNzX3BhdHRlcm4uUmRhdGEiKQpzdW1tW2lkLCBjKCJWYXJDb2RlIiwgIlZhck5hbWUiKV0KYGBgCgpUaGUgcGxvdCBiZWxvdyBzaG93cyB0aGUgbWlzc2luZyBwYXR0ZXJuIGFtb25nIHRoZSBhYm92ZSB2YXJpYWJsZXMsIHdoZXJlIHRoZSAxNjk0MyByb3dzIGluZGljYXRpbmcgcmVjb3JkcyBpbiBvdXIgZGF0YWJhc2UsIGNvbHVtbnMgaW5kaWNhdGluZyB2YXJpYWJsZXMsIGFuZCBjZWxscyBhcmUgY29sb3JlZCBhY2NvcmRpbmcgdG8gbWlzc2luZ25lc3M6IGJsYWNrID0gbWlzc2luZyAgYW5kIGdyZXkgPSBwcmVzZW50LiBUaGUgcGxvdCwgbmFtZWQgIlJwbG90X21pc3NfcGF0dGVybi5wbmcsIiBjYW4gYmUgZG93bmxvYWRlZCBmcm9tIGdvb2dsZSBkcml2ZS4KCgpgYGB7ciwgZWNobz1GQUxTRSwgZmlnLndpZHRoID0gNCwgZmlnLmhlaWdodD00fQpsb2FkKCJwbG90X21pc3NfcGF0dGVybi5SZGF0YSIpCm1pc3MucGF0dGVybgpgYGAKCiMjIFVuZGVyc3RhbmQgbWlzc2luZyBTVEMKCkhvdyBtaXNzaW5nIGlzIGNvZGVkIGluIExBRCBhbmQgU1RDPyAKCmBgYHtyLCBlY2hvPUZBTFNFfQpkYXRhJExBRCA9IGFzLmNoYXJhY3RlcihkYXRhJExBRCkKaWQgPSB3aGljaChkYXRhJExBRCAlaW4lIGMoIlVua25vd24iLCAiTm90IEFwcGxpY2FibGUiLCAiVW5rbm93biBidXQgTG9hZGVkIikpCmlkID0gYyhpZCwgd2hpY2goZGF0YSRTVEMgJWluJSBjKC05OSwgLTk4LCAxNSkpKQppZCA9IHVuaXF1ZShpZCkKdG1wID0gZGF0YVtpZCwgYygiTEFEIiwgIlNUQyIpXQp0YWJsZSh0bXAkTEFELCB0bXAkU1RDKQpgYGAKCi0gVGhyZWUgb2JzIHdpdGggU1RDIHRhZyAiLTk5IiBidXQga25vd24gTEFEIHRhZ3M6ICJDaGxvcmluZSIgKG9icyA3Mjk1KSBhbmQgIkNvcm4gU3lydXAiIChvYnMgMTA2NDEsIDEwNjQyKS4gCgotIFR3byBvYnMgd2l0aCBMQUQgdGFnICJVbmtub3duIiBidXQga25vd24gU1RDIHRhZ3MuIExpa2VseSB0aGUgdHdvIGluY2lkZW50cyBzaG91bGQgaGF2ZSBMQUQgdGFncyAiRXRoeWxlbmUgT3hpZGUiIGFuZCAiU3VsZnVyaWMgQWNpZCIsIHJlc3BlY3RpdmVseS4gCgojIyMjIFRoZSBSIG91dHB1dCBzaG93cyB0aGUgbnVtYmVyIG9mIGluY2lkZW50cyB3aXRoIG1pc3NpbmcgU1RDIGFsb25nIHllYXJzLiAKCmBgYHtyLCBlY2hvPUZBTFNFfQp0YWJsZSh5ZWFyKGRhdGEkRE9BW2lkXSkpCmBgYApUaGVyZSBhcmUgODQzNyByZWNvcmRzIHdpdGggU1RDIG1pc3NpbmcgKGNvZGVkIGFzIC05OSksIG1vc3Qgb2Ygd2hpY2ggaGF2ZSBgSU5DYCBzdGFydGluZyB3aXRoICJBQVIiIG9yICJEVEUiLiBJbiB0aGUgd2hvbGUgZGF0YWJhc2UsIHRoZXJlIGFyZSAKCi0gNjI0MiByZWNvcmRzIHdpdGggYElOQ2Agc3RhcnRpbmcgd2l0aCAiQUFSIiAoYW1vbmcgd2hpY2gsIDQ0MjUgd2l0aCBTVEMgbWlzc2luZyksIAoKLSAyNDU0IHJlY29yZHMgd2l0aCBgSU5DYCBzdGFydGluZyB3aXRoICJEVEUiIChhbW9uZyB3aGljaCwgMjMzMyB3aXRoIFNUQyBtaXNzaW5nKS4gCgpJTkMgZW50cnkgdGhhdCBiZWdpbnMgIkRURSIgKGZvciBEYXRlIFRvIGJlIEVzdGFibGlzaGVkKSBiYXNpY2FsbHkgbWVhbnMgdGhhdCB0aGUgYWNjaWRlbnQncyBET0EgaXMgdW5rbm93bi4gCgojIyBNaXNtYXRjaCBiZXR3ZWVuIERPQSBhbmQgRE9BLjEKClRoZXJlIGFyZSBmaXZlIHJlY29yZHMgd2hvc2UgYERPQWAgYW5kIGBET0EuMWAgYXJlIGRpZmZlcmVudC4gCgpgYGB7ciwgZWNobz1GQUxTRX0KaWQgPSB3aGljaChkYXRhJERPQSAhPSBkYXRhJERPQS4xKQpkYXRhW2lkLCBjKCJJTkMiLCAiRE9BIiwgIkRPQS4xIiwgIlNUQyIsICJMQUQiKV0KYGBgCgojIyBNaXNtYXRjaCBiZXR3ZWVuIERPRSBhbmQgRE9FLjEKClRoZXJlIGFyZSAxODIwIHJlY29yZHMgd2hvc2UgYERPRWAgYW5kIGBET0UuMWAgYXJlIGRpZmZlcmVudCAoaW5jbHVkaW5nIHRoZSAxc3QgcmVjb3JkLCAiNTU1NyIsICBpbiB0aGUgdGFibGUgYWJvdmUpLiBBbW9uZyB0aGVtLCBtb3N0IGhhdmUgYERPRWAgbGF0dGVyIHRoYW4gYERPRS4xYCAKCi0gVGhyZWUgKHdoaWNoIGFjdHVhbGx5IGNvcnJlc3BvbmQgdG8gdGhlIHNhbWUgaW5jaWRlbnQpIGhhdmUgRE9FIDI1Ni1kYXkgPGI+ZWFybGllcjwvYj4gdGhhbiBET0UuMTsgCgotIE9uZSBoYXMgRE9FIDE4Mi1kYXkgPGI+ZWFybGllcjwvYj4gdGhhbiBET0UuMTsKCi0gT25lIGhhcyBET0UgNS1kYXkgPGI+ZWFybGllcjwvYj4gdGhhbiBET0UuMTsKCi0gNDk0IGhhdmUgdGhlIHR3byBkYXRlcyAxLWRheSBhcGFydDsgCgotIDE3NyBoYXZlICBET0UgMi1kYXkgbGF0dGVyIHRoYW4gRE9FLjE7IAoKLSAyNTggaGF2ZSBET0UgMy1kYXkgbGF0dGVyIHRoYW4gRE9FLjE7IAoKLSA4ODYgaGF2ZSBET0UgbW9yZSB0aGFuIDMtZGF5IGxhdHRlciB0aGFuIERPRS4xOyAKClRoZXNlIDE4MjAgcmVjb3Jkcywgb3JkZXJlZCBieSB0aGUgZGlmZmVyZW5jZSBiZXR3ZWVuIHRoZSB0d28gZGF0ZXMsIGFyZSBpbiBmaWxlICJzdWJzZXRfZGF0YV9ET0VfbWlzbWF0Y2guY3N2IiBpbiBnb29nbGUgZHJpdmUuCgpgYGB7ciwgZWNobz1GQUxTRX0KaWQgPSB3aGljaChkYXRhJERPRSAhPSBkYXRhJERPRS4xKQp0bXAub3JkZXIgPSBvcmRlcihkYXRhJERPRVtpZF0gLSBkYXRhJERPRS4xW2lkXSkKd3JpdGUuY3N2KGRhdGFbaWRbdG1wLm9yZGVyXSwgYygiRE9FIiwgIkRPRS4xIiwgc2V0ZGlmZihuYW1lcyhkYXRhKSwgYygiRE9FIiwgIkRPRS4xIikpKV0sIGZpbGU9InN1YnNldF9kYXRhX0RPRV9taXNtYXRjaC5jc3YiKQpgYGA=