data <- read.csv("http://blog.lib.umn.edu/wage0005/echinacea/angustifoliacrossdata_2012_all_locations.csv")

str(data)
data[ ,6]
data[ ,12]
str(data)
levels(data$styleStatus)

data$compat <- data$styleStatus

levels(data$compat)
levels(data$compat)[5] <- "DNC"
levels(data$compat)

levels(data$compat)[4] <- "DNC"
levels(data$compat)

levels(data$compat)[5] <- "DNC"
levels(data$compat)

levels(data$compat)[5] <- "DNC"
levels(data$compat)

levels(data$compat)[6] <- "DNC"
levels(data$compat)

levels(data$compat)[6] <- "DNC"
levels(data$compat)

levels(data$compat)[10] <- "DNC"
levels(data$compat)

levels(data$compat)[2] <- "DC"
levels(data$compat)

levels(data$compat)[14] <- "DC"
levels(data$compat)

levels(data$compat)[14] <- "DC"
levels(data$compat)

levels(data$compat)[13] <- "DC"
levels(data$compat)

levels(data$compat)[10] <- "DC"
levels(data$compat)

levels(data$compat)[8] <- "DC"
levels(data$compat)

levels(data$compat)[5] <- "MNC"
levels(data$compat)

levels(data$compat)[1] <- "US"
levels(data$compat)

levels(data$compat)[7] <- "US"
levels(data$compat)

levels(data$compat)[8] <- "US"
levels(data$compat)

levels(data$compat)[6] <- "MNC"
levels(data$compat)

levels(data$compat)[7] <- "MC"
levels(data$compat)

levels(data$compat)[6] <- "DC"
levels(data$compat)

data[1:14]

str(data)

data$maternal.ID
as.integer(round(data$maternal.ID,0))

data$mat <- as.integer(round(data$maternal.ID,0))
data$mat

str(data)

########### load dataset with distances
load(url("http://blog.lib.umn.edu/wage0005/echinacea/distancesForLydia.RData"))

# merge Lydia's data with distances
dataPlus <- merge(data, distances)
str(dataPlus) # notice one record is missing
setdiff(data$paternalID, distances$paternalID) # what happened here?