Figure out the age control by year data for Neotoma.

Dating_neotoma_Pollen

library(neotoma)
library(reshape2)

pollen_data <- get_dataset(datasettype= 'pollen')
full_records <- get_download(pollen_data) # this takes a bit of time

pub.dates <- get_publication(pollen.data) # so does this

earliest.year <- do.call(rbind.data.frame, 
                         lapply(pub.dates, function(x)data.frame(id = NA, 
                                                                 year = min(sapply(x, function(y)y$meta$year)))))

earliest.year$id <- sapply(pollen.data, function(x) x$dataset$dataset.meta$dataset.id)

aa <- lapply(full_records, function(x){
  control.table <- try(get_chroncontrol(x$sample.meta$chronology.id[1])$chron.control$control.type)
  if(!class(control.table) == 'try-error'){
    output <- data.frame(control = as.character(control.table), id = x$dataset$dataset.meta$dataset.id)
  } else {
    output <- data.frame(control = NA, id = NA)
  }
  output})

chron.table <- do.call(rbind.data.frame, aa)
chron.table$year <- earliest.year$year[match(chron.table$id, earliest.year$id)]

model.table <- dcast(chron.table, formula = control~year, length)

It's because it should be pulling from full_data, not the pollen_data that's getting pulled. That's my mistake.