[R] Fw: showing ratio in ggplot bar plot with facets

Melanie Zoelck mzoelck at gmx.com
Thu Mar 14 20:27:43 CET 2013


Dear R help members,

I have constructed a bar plot using ggplot, where I am showing the number of particles in 10m depth bins (from 0-150m) by release area (3 different ones) and Model (also 3 different ones). My problem is that for each release area a different number of particles was released (say 35, 45 and 56). I now need to be able to show the number of particles in each depth bin proportional to the number of particles originally released. So that all bars that show release area 1 are divided by 35, those for release area 2 divided by 45 and those for release area 3 by 56. Is it possible to do this in ggplot? I thought that it would work using the following code, but I don't think its doing what I need it to do:

AreaTidalEnd.df<-read.csv("./Analyses/Dispersal_Modelling/Data/AreaTidalEnd.csv", header=TRUE)
AreaTidalEnd.df$Model<-factor(as.character(AreaTidalEnd.df$Model),levels=c("extP","extL","extT"))
levels(AreaTidalEnd.df$Model)[levels(AreaTidalEnd.df$Model)=="extP"] <- "PTS"
levels(AreaTidalEnd.df$Model)[levels(AreaTidalEnd.df$Model)=="extL"] <- "LGS"
levels(AreaTidalEnd.df$Model)[levels(AreaTidalEnd.df$Model)=="extT"] <- "TDS"
AreaTidalEnd.df$Year<-as.factor(as.character(AreaTidalEnd.df$Year))

breaks<-seq(0,150,20)

AreaTidal.df<-subset(AreaTidalEnd.df, Depth>=0)
AreaTidal2003<-subset(AreaTidal.df, Year=="2003")

##Counts per Model and RA
update_stat_defaults("bin",aes(y=..density..))

qplot(Depth, data = AreaTidal2003, binwidth = 10, origin=0,xlim=c(0,150), ylab="Larval Ratio")+ scale_x_continuous(breaks = breaks, labels=breaks)+facet_grid(Release_Area~Model)

I would be hugely grateful for your help!

Melanie

____________________________________________
Melanie Zölck (Zoelck)
PhD Candidate
Galway-Mayo Institute of Technology
Marine and Freshwater Research Centre
Commercial Fisheries Research Group
Department of Life Science
Dublin Road
Galway
Republic of Ireland
Mobile: +353 (0) 85 7246196
Skype: melaniezoelck
E-mail: mzoelck at gmx.com or mzoelck at hotmaill.com



More information about the R-help mailing list