[BioC] Standard deviation

Heidi Dvinge heidi at ebi.ac.uk
Mon Aug 24 16:41:37 CEST 2009


Hello,

try having a look at the plotCI function from gplots, for plotting error
bars and confidence intervals.

Cheers
\Heidi

> Dear BioCs,
> I'd like to receive some help of you. I need to plot a growth curve (which
> I
> already did) with standard deviation bars (which I didn't do yet).
>
> I already have the standard deviation values calculated. I believe I can
> read the file containing those values and plot above the points in my
> graphic, or calculate those values using R, and plot then.
>
> Can you help me giving some suggestions?
>
> tabela=read.table("irradiacao2.txt",sep="\t",skip=1)
> tabelan2=tabela[,2:7]
> colnames(tabelan2)=c("Pre-irradiation","0H","4H","24H","48H","96H")
> rownames(tabelan2)=tabela[,1]
> tabela2 = as.matrix(tabelan2)
>
> plot(tabela2[1,],type="l", col=" blue", bty="o", lty=1, font.main=2,
> ylim=c(10000000,120000000), main="Growth curves", ylab="cell number",
> xlab="Hours after treatment", xaxt="n")
> points(tabela2[1,], pch=1, col="blue", cex=1.0)
> axis(1, at=1:6, labels=colnames(tabelan2))
> lines(tabela2[2,],type="l", col="black")
> points(tabela2[2,], pch=1, col="black", cex=1.0)
> lines(tabela2[3,], type="l",col="blue")
> points(tabela2[3,], pch=2, col="blue", cex=1.0)
> lines(tabela2[4,],type="l", col="black")
> points(tabela2[4,], pch=2, col="black", cex=1.0)
> legend("topleft", cex=0.9, legend=c("Control RI", "Treated RI", "Control
> RII", "Treated RII"), lty=c(1,1,1,1), pch=c(1,2,1,2), col=c("blue",
> "blue",
> "black", "black"))
>
>
> --
> Priscila Grynberg, B.Sc., M.Sc.
> Doutoranda em Bioinformática (Bioinformatics D.Sc student)
> Laboratório de Genética Bioquímica
> Universidade Federal de Minas Gerais
> Tel: +55 31 3409-2628
> CV: http://lattes.cnpq.br/8808643075395963
>
> 	[[alternative HTML version deleted]]
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives:
> http://news.gmane.org/gmane.science.biology.informatics.conductor



More information about the Bioconductor mailing list