Skip to content

Instantly share code, notes, and snippets.

@geneva
geneva / Bimini_plot.R
Created May 29, 2019 23:14
Plotting Bimini
library(dplyr)
library(ggmap)
library(ggsn)
#Download map data for South Bimini
pngMAP_df <- get_map(location = c(lon = -79.275, lat = 25.71),
source = "google",
zoom = 13,
color="bw",
@geneva
geneva / rwty_runthrough.R
Created June 13, 2018 00:08
example execution of R package - rwty
# aRe We There Yet (RWTY) Tutorial
# Bayesian Phylogenetic Inference using Markov Chain Monte Carlo (mcmc)
### Introduced in late 90s, now wildly popular
### Implemented in many software packages: MrBayes, BEAST, PhlyoBayes
# Phylogenetic inference computationally expensive, number of possible
# rooted trees for 60 taxa > number of protons in known universe
# Exhausitive tree search via optimality criteria such as maximum parsimony
# or maximum liklihood is NP-Hard
@geneva
geneva / AFLP_functions.R
Last active October 30, 2019 03:06
A collection of R Functions for analysis of AFLP data
###########################################
# Functions for the analysis of AFLP data #
# Anthony J Geneva #
# 9 April 2013 #
###########################################
###########################################
# This function converts the binary table
# output from RawGeno to a format that can
# be read by AFLPScore
@geneva
geneva / RawGeno_2_AFLPscore
Last active October 30, 2019 03:07
A simple R function that converts the binary table output from RawGeno to a format that can be read by AFLPScore
# This function converts the binary table output from
# RawGeno to a format that can be read by AFLPScore
##################################
# RawGeno_2_AFLPscore function #
##################################
RawGeno_2_AFLPscore <- function(outfilename)
{
@geneva
geneva / plot_morphospecies.R
Created February 5, 2011 21:41
Reads TA data and plots piechart and chao1/rarefaction line graphs3
library(vegan)
class <- "Anthony"
dat <- read.table(paste(class,".dat", sep=""), sep = "\t", header=TRUE)
#######PieChart############
piedata <- matrix(ncol=1, nrow=length(unique(dat[,1])))
names <- levels(unique(dat[,1]))
rownames(piedata) <- names