Hi, May be: dat1<- read.table(text="city year sex obs 1 1990 M 25 1 1990 F 32 1 1991 M 15 1 1991 F 22 2 1990 M 42 2 1990 F 36 2 1991 M 12 2 1991 F 16",sep="",header=TRUE,stringsAsFactors=FALSE)
library(plyr) #by city d1 <- ddply(dat1,.(city),summarize,Tot=sum(obs)) d2 <- ddply(dat1,.(city,sex),summarize,Tot=sum(obs)) res <- merge(d1,d2,by="city") res$percent <- round((res$Tot.y/res$Tot.x)*100,1) library(plotrix) plot(percent~city,data=res,type="p") thigmophobe.labels(res$city,res$percent,labels=res$sex) #by year d1 <- ddply(dat1,.(year),summarize,Tot=sum(obs)) d2 <- ddply(dat1,.(year,sex),summarize,Tot=sum(obs)) res2 <- merge(d1,d2,by="year") res2$percent <- round((res2$Tot.y/res2$Tot.x)*100,1) plot(percent~year,data=res2,type="p",xaxt="n") axis(1,at=res2$year,labels=res2$year) thigmophobe.labels(res2$year,res2$percent,labels=res2$sex) A.K. On Thursday, October 10, 2013 10:37 PM, Val <valkr...@gmail.com> wrote: Hi All, I have a huge data set with the following type; city year sex obs 1 1990 M 25 1 1990 F 32 1 1991 M 15 1 1991 F 22 2 1990 M 42 2 1990 F 36 2 1991 M 12 2 1991 F 16 I want to calculate the percentage of M and F by city, year and year within city and also plot. city 1 total Male= 40; total female= 54; %M= 40/(40+54)=42.6 %F= 54/(40+54)=57.4 and so on. Can any body help me out? Thanks in advance [[alternative HTML version deleted]] ______________________________________________ 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. ______________________________________________ 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.
