speciesList.R

#install rglobi (see ropensci.org for more info)
install.packages(c('rglobi'))

# read species list after exporting as csv
krak <- read.csv(file='Krak_list_DELETE.csv', stringsAsFactors=F)
# pick unique families
taxonList <- unique(krak$Family)

# find and aggregate interactions for each species . . . one by one
interactions <- Reduce(function(accum, taxon) {
    message(paste(nrow(accum), '.', sep=''), appendLF=F)
    rbind(accum, rglobi::get_interactions_by_taxa(taxon, interactiontype='interactsWith', returnobservations=T))
  }, 
  taxonList, 
  init=data.frame())

# write interactions to file
write.csv(interactions, file='interactions.csv')