Summary of the Data
- How many missing (value = -99 or blank entry):
missing_# (counts)
- How many unique values:
unique_#
- Note
LAD can take values like “Unknown” or “Unknown but loaded”, although no missing is reported in the table below.
- Note
CST = 2 for all incidents in 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
Three obs with STC tag “-99” but known LAD tags: “Chlorine” (obs 7295) and “Corn Syrup” (obs 10641, 10642).
Two obs with LAD tag “Unknown” but known STC tags. Likely the two incidents should have LAD tags “Ethylene Oxide” and “Sulfuric Acid”, respectively.
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
Three (which actually correspond to the same incident) have DOE 256-day earlier than DOE.1;
One has DOE 182-day earlier than DOE.1;
One has DOE 5-day earlier than DOE.1;
494 have the two dates 1-day apart;
177 have DOE 2-day latter than DOE.1;
258 have DOE 3-day latter than DOE.1;
886 have DOE more than 3-day 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=