gistfile1.txt

### Recode all character to factor data ## 
arc[sapply(arc, is.character)]<-lapply(arc[sapply(arc, is.character)],as.factor)
str(arc)
## Recode Generation (A/S)
levels(arc$`generation(a/s)`)[levels(arc$`generation(a/s)`)=="05SB"]<-NA
### LOOK AT RAW DATA
dim(arc)
str(arc)
skim(arc)
## tabulate the data ###
table(arc$region)
table(arc$region, arc$slaughteryear)
table(arc$slaughteryear)
heatmap(table(arc$weightclasshalv, arc$slaughterq))
heatmap(table(arc$slaughteryear,arc$regi))