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: mzoe...@gmx.com or mzoe...@hotmaill.com
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
