#Import original dataset
muleys <-read.csv("DCmuleysedited.csv", header=T)
str(muleys)
#Remove outlier locations and Mule deer D12 with too few locations
muleys <-subset(muleys, muleys$X > 599000 & muleys$X < 705000 & muleys$Y > 4167000 & muleys$id != "D12")
muleys$id <- factor(muleys$id)#This step must be done to completely remove D12
summary(muleys$id)
