setwd("/Users/jameshanley/Desktop/bios601-2008")

ds=read.table("sexratio.csv",header=T,sep=",")

ds

ds[,"p.male.c"]=ds[,"mc"]/(ds[,"mc"]+ds[,"fc"])
ds[,"sexratio.c"]=ds[,"mc"]/ds[,"fc"]

ds[,"p.male.pei"]=ds[,"mpei"]/(ds[,"mpei"]+ds[,"fpei"])
ds[,"sexratio.pei"]=ds[,"mpei"]/ds[,"fpei"]

plot(ds$year,ds$p.male.c)
plot(ds$year,ds$sexratio.c)
plot(ds$year,(ds$fc+ds$mc))

plot(ds$year,ds$p.male.pei)
plot(ds$year,ds$sexratio.pei)
plot(ds$year,(ds$fpei+ds$mpei))