Pie Charts:

DATA LINK: https://drive.google.com/file/d/0BwoPt0jyGdzORkM3cVA0WjJodVk/view?usp=sharing
Mydata<-read.csv(file="final_analysis_candy_analysis.csv",head=TRUE,sep=",") dhfr.Arg <- table(Mydata$dhfr.Arg.59.163.137.) dhfr.Ile <- table(Mydata$dhfr.Ile.51.214.65.) dhfr.Asn108 <- table(Mydata$dhfr.Asn108.328.372.) glu.540 <- table(Mydata$glu.540.538.326.200.) gly.437 <- table(Mydata$gly.437.848.300.) library(plotrix) par(op) op <-par(mfrow=c(2,3),mar=c(0,0,1,0)) pct <- round(dhfr.Arg/sum(dhfr.Arg)*100) lbls <- paste(names(dhfr.Arg), pct) # add percents to labels lbls <- paste(lbls,"%",sep="") # ad % to labels lp<-pie3D(dhfr.Arg,radius=0.8,labels=lbls,explode=0.1, labelrad=1.4,main="dhfr Arg 59(163,137)") pct <- round(dhfr.Ile/sum(dhfr.Ile)*100) lbls <- paste(names(dhfr.Ile), pct) # add percents to labels lbls <- paste(lbls,"%",sep="") # ad % to labels lp<-pie3D(dhfr.Ile,radius=0.8,labels=lbls,explode=0.1, labelrad=1.4,main="dhfr Ile 51(214,65)") pct <- round(dhfr.Asn108/sum(dhfr.Asn108)*100) lbls <- paste(names(dhfr.Asn108), pct) # add percents to labels lbls <- paste(lbls,"%",sep="") # ad % to labels lp <- pie3D(dhfr.Asn108,radius=0.8,labels=lbls,explode=0.1, labelrad=1.4,main="dhfr Asn108(328,372)") pct <- round(glu.540/sum(glu.540)*100) lbls <- paste(names(glu.540), pct) # add percents to labels lbls <- paste(lbls,"%",sep="") # ad % to labels lp<-pie3D(glu.540,radius=0.8,labels=lbls,explode=0.1, labelrad=1.4,main="glu 540(538,326,200)") pct <- round(gly.437/sum(gly.437)*100) lbls <- paste(names(gly.437), pct) # add percents to labels lbls <- paste(lbls,"%",sep="") # ad % to labels lp <- pie3D(gly.437,radius=0.8,labels=lbls,explode=0.1, labelrad=1.4,main="gly 437(848,300)") par(op)
