JoFrhwld

\documentclass{beamer}
\usetheme{Singapore}
\usepackage{tipa}

%% Make the r-code small
\ifdefined\knitrout
  \renewenvironment{knitrout}{\begin{footnotesize}}{\end{footnotesize}}
\else
\fi

JoFrhwld

## ggplot2 and mgcv for the plot
library(ggplot2)
library(mgcv)

## scales packages to define revlog
library(scales)

revlog_trans <- function(base = exp(1)){
    ## Define the desired transformation.
    trans <- function(x){

JoFrhwld

library(rgdal) #This may be a pain to get installed. The OGR functions come from it.
## install gdal, available here: http://www.gdal.org/
## install proj.4, available here: http://trac.osgeo.org/proj/
## For some reason, R wasn't expecting for libproj to be installed where it was,
## so I had to tell it to look at /usr/local/bin

library(ggplot2)
library(reshape2)
library(plyr)
gpclibPermit()

JoFrhwld

library(rgdal) #This may be a pain to get installed. The OGR functions come from it.
## install gdal, available here: http://www.gdal.org/
## install proj.4, available here: http://trac.osgeo.org/proj/
## For some reason, R wasn't expecting for libproj to be installed where it was,
## so I had to tell it to look at /usr/local/bin

library(ggplot2)
library(reshape2)
library(plyr)
gpclibPermit()

JoFrhwld

##########
## Block 1: Preparing the data

library(plyr)
library(reshape)

## Create a Subsystem Column
all_philly$Subsystem <- all_philly$VClass
levels(all_philly$Subsystem) <- c("Vhr", "Vw", "V", "misc", "Vh", "Vh", "Vhr", "Vw", "Vy", 
"Vy", "V", "Vy", "Vy", "Vy", "V", "Vw", "Vy", "Vy", "Vhr", 

JoFrhwld

clip.border <- function(deldir.obj, border){
	## deldir.obj = output of deldir::deldir
	## border = list of borders defined in data frames
	##		Importantly, with columns called "x" and "y" 
	require(gpclib)

	## bord will be a gpc polygon of all of the regions given to the argument "border"
	bord <- as(matrix(ncol = 2), "gpc.poly")
	for(i in seq(along = border)){
		b <- as(border[[i]], "gpc.poly")

JoFrhwld

scored_vclass <- function(df,class,vowels, dims = c("F1","F2")){
	require(MASS)
	df <- df[df[ ,class] %in% vowels, ]
	

	df[ ,class] <- as.character(df[ ,class])
	df[,class] <- as.factor(df[,class])
	
	
	for(i in dims){