library(psych)
library(ggplot2)
AllAbuseData <- read.csv(file="C:/MSDA/D607/Project2/DrugFatalitiesCT/Accidental_Drug_Related_Deaths__2012-2016.csv")Two questions - 1) Are people of age less than 40 more likely than those above 40 to be a drug abuse fatality? 2) Are people who live in rural areas more likely than those in urban areas to be a drug abuse fatality?
What are the cases, and how many are there?
Each case represents a drug abuse fatality in the state of Connecticut during the period of 2012-2015. There 3045 observations in the given data set.
Describe the method of data collection.
Data are derived from an investigation by the Office of the Chief Medical Examiner which includes the toxicity report, death certificate, as well as a scene investigation.
What type of study is this (observational/experiment)?
This is an observational study.
If you collected the data, state self-collected. If not, provide a citation/link.
The site for the drug death data is https://catalog.data.gov/dataset/accidental-drug-related-deaths-january-2012-sept-2015. For the categorization of towns or cities as rural, the data is taken from http://www.ruralhealthct.org/towns.htm, which is maintained by CT Office of Rural Health.
What is the response variable, and what type is it (numerical/categorical)?
There are 2 response variables - Age and Location, which are numerical and categorical variables respectively.
What is the explanatory variable, and what type is it (numerical/categorival)?
The explanatory variables are 1) number of incidences for age groups 2) number of incidences for death locations classified as rural
All variables which have been captured in study.
names(AllAbuseData)## [1] "CaseNumber" "Date"
## [3] "Sex" "Race"
## [5] "Age" "Residence.City"
## [7] "Residence.State" "Residence.County"
## [9] "Death.City" "Death.State"
## [11] "Death.County" "Location"
## [13] "DescriptionofInjury" "InjuryPlace"
## [15] "ImmediateCauseA" "Heroin"
## [17] "Cocaine" "Fentanyl"
## [19] "Oxycodone" "Oxymorphone"
## [21] "EtOH" "Hydro.codeine"
## [23] "Benzodiazepine" "Methadone"
## [25] "Amphet" "Tramad"
## [27] "Morphine..not.heroin." "Other"
## [29] "Any.Opioid" "MannerofDeath"
## [31] "AmendedMannerofDeath" "DeathLoc"
Description of Age
describe(AllAbuseData$Age)## vars n mean sd median trimmed mad min max range skew kurtosis
## X1 1 3043 41.77 12.34 42 41.65 14.83 14 87 73 0.07 -0.84
## se
## X1 0.22
A graph showing distrubution of cases by age.
ggplot(AllAbuseData, aes(x=Age)) + geom_histogram()Listing the names of counties in drug data.
unique(AllAbuseData$Residence.County)## [1] WINDHAM LITCHFIELD HARTFORD
## [4] NEW HAVEN MIDDLESEX FAIRFIELD
## [7] TOLLAND DUTCHESS
## [10] NEW LONDON PUTNAM MUSCOGEE
## [13] USA NORFOLK HILLSBOROUGH
## [16] MIAMI-DADE ORANGE WESTCHESTER
## [19] QUEENS NEW CASTLE COOK
## [22] HAMPDEN NEW YORK WINNEBAGO
## [25] ESSEX COLLIER SARATOGA
## [28] LOS ANGELES DAKOTA CHESAPEAKE CITY
## [31] SANTA CRUZ WORCESTER MADISON
## [34] ST JOHN THE BAPTIST BERKSHIRE CUMBERLAND
## [37] CUYAHOGA HARRIS SUFFOLK
## [40] NASSAU LAKE PASSAIC
## [43] KINGS GLOUCESTER Hartford
## [46] New Haven AROOSTOOK Fairfield
## [49] NEWPORT Windham Tolland
## [52] Middlesex Litchfield PROVIDENCE
## [55] KNOX WESTMORELAND ROGERS
## [58] OCEAN ULSTER WASHINGTON
## [61] WELD
## 61 Levels: AROOSTOOK BERKSHIRE CHESAPEAKE CITY COLLIER ... WORCESTER
unique(AllAbuseData$Residence.City)## [1] ONECO WINCHESTER AVON
## [4] NEW HAVEN MIDDLETOWN NEWTOWN
## [7] HAMDEN MERIDEN HARTFORD
## [10] NORTH CANAAN WETHERSFIELD STAMFORD
## [13] ENFIELD WATERBURY CHESHIRE
## [16] SOUTHINGTON NEW BRITAIN MANCHESTER
## [19] VINE GROVE NEWINGTON BROOKLYN
## [22] WALLINGFORD VERNON
## [25] WILTON CLINTON GALES FERRY
## [28] SOUTH WINDSOR PLAINVILLE STANFORDVILLE
## [31] WEST HARTFORD WOLCOTT GUILFORD
## [34] ANDOVER EAST HAMPTON DAYVILLE
## [37] WEST HAVEN LEBANON NEW LONDON
## [40] WILLIMANTIC BRIDGEPORT RIDGEFIELD
## [43] COVENTRY SEYMOUR ELLINGTON
## [46] MIDDLEBURY NORWICH NAUGATUCK
## [49] SHELTON TORRINGTON BREWSTER
## [52] OLD LYME THOMASTON MILFORD
## [55] EAST HAVEN STRATFORD WATERTOWN
## [58] BRISTOL DANBURY COLUMBUS
## [61] LISBON WINDSOR LOCKS PROSPECT
## [64] ANSONIA NEW MILFORD SALEM
## [67] GREENWICH DEEP RIVER NORTH HAVEN
## [70] KILLINGWORTH SALISBURY PITTSBURGH
## [73] EAST WINDSOR VOLUNTOWN FARMINGTON
## [76] STAFFORD HARWINTON NIANTIC
## [79] UNCASVILLE EAST HARTFORD GLASTONBURY
## [82] GROTON WELLESLEY MARLBOROUGH
## [85] BRANFORD HANCOCK BOLTON
## [88] WINSTED FRANKLIN LEDYARD
## [91] WILLACOOCHEE JEWETT CITY BROOKFIELD
## [94] WATERFORD NEW BRIT SOUTH GLASTONBURY
## [97] DOVER PLAINS DANIELSON AVENTURA
## [100] TRUMBULL NORWALK NO HAVEN
## [103] ORANGE OAKDALE TERRYVILLE
## [106] OAKVILLE CANTON NORTH SALEM
## [109] MANSFIELD NEW FAIRFIELD COLCHESTER
## [112] DERBY SUFFIELD PRESTON
## [115] EAST CANAAN MOOSUP FLUSHING
## [118] MONROE TAFTVILLE MADISON
## [121] EAST LYME NEW CASTLE AMSTON
## [124] ROCKY HILL UNIONVILLE ROCKVILLE
## [127] TOLLAND ARLINGTON HEIGHTS MYSTIC
## [130] WASHINGTON CROMWELL PAWCATUCK
## [133] MONTVILLE BERLIN DURHAM
## [136] BALTIC NORTH GROSVENORDALE PORTCHESTER
## [139] FAIRFIELD PUTNAM WINDSOR
## [142] PORTLAND PALMER TARIFFVILLE
## [145] IVORYTON THOMPSON WOODBURY
## [148] S GLASTONBURY NEW YORK DARIEN
## [151] WINGDALE WOODBRIDGE PORT CHESTER
## [154] ESSEX BURLINGTON OSHKOSH
## [157] ALFRED STATION STAFFORD SPRINGS REDDING
## [160] SAUGUS POMFRET MORRIS
## [163] STAFFORD SPGS BOZRAH SUDBURY
## [166] PLAINFIELD N HAVEN SOMERS
## [169] QUAKER HILL JASPER NORTH BRANFORD
## [172] NORTHFORD OXFORD NEW BEDFORD
## [175] RYE VERNON ROCKVL HANOVER
## [178] SHERMAN EAST GRANBY WINDHAM
## [181] KILLINGLY SOUTHBURY NEW CANAAN
## [184] SAKEM BROAD BROOK CHESTER
## [187] EASTFORD HEBRON DEER LODGE
## [190] BEACON FALLS EAST HADDAM WESTPORT
## [193] WATERBRUY HURLEY SIMSBURY
## [196] ALLSTON BETHEL OLD SAYBROOK
## [199] MILLERTON STERLING BASS HARBOR
## [202] BETHANY SHARON CANTERBURY
## [205] ZIONSVILLE EAST WOODSTOCK COLUMBIA
## [208] WESTBROOK NAPLES SARATOGA SPGS
## [211] EAST BERLIN CHARLESTOWN BLOOMFIELD
## [214] PLANTSVILLE HAWAIIAN GARDENS INVER GROVE HEIGHTS
## [217] MOODUS HADDAM CHESAPEAKE
## [220] OLD GREENWICH WEST HAVEN WILLINGTON
## [223] ROGERS GOSHEN SANTA CRUZ
## [226] EAST HARTLAND MAPLEWOOD GRISWOLD
## [229] GRANBY GROTON LONG POINT WEBSTER
## [232] SANDY HOOK PAWTUCKET W HAVEN
## [235] NORTH STONINGTON QUINCY ORLANDO
## [238] TACONIC ASTORIA HAMPTON BAYS
## [241] PLYMOUTH EAST MEADOW KINGSTON
## [244] CANAAN WEATOGUE SPRAGUE
## [247] QUINEBAUG COLLINSVILLE MAHOPAC
## [250] SOUTH LYME HAMILTON LITCHFIELD
## [253] NUTLEY FALLS VILLAGE KENSINGTON
## [256] CHICOPEE LA PLACE GRANVILLE
## [259] MELROSE WESTON LEOMINSTER
## [262] PITTSFIELD MT VERNON CHAPLIN
## [265] YORKTOWN HEIGHTS STONINGTON SOUTHPORT
## [268] NOANK VINELAND MILLDALE
## [271] PARMA KINGWOOD NEW HARTFORD
## [274] CHELSEA WILLISTON PARK CLERMONT
## [277] MIDDLEFIELD PASSAIC BANTAM
## [280] RED HOOK NEWFIELD HODGDON
## [283] NEWPORT WOODSTOCK CHEPACHET
## [286] WHITE PLAINS HIGGANUM ASHFORD
## [289] KNOXVILLE VANDERGRIFT CLAREMORE
## [292] HOUSTON JACKSON EASTON
## [295] HIGHLAND RICHMOND COS COB
## [298] MILLIKEN WEST SUFFIELD
## 299 Levels: ALFRED STATION ALLSTON AMSTON ANDOVER ... ZIONSVILLE