top of page

Put Sd bars on Barplot on R


Data in the file "Saison.txt":

("diversité.génétique" means genetics diversity and "diversité.haplotypique" means haplotypic diversity).

Mars Septembre

diversité.génétique 0.2342 0.311125

sd.dg 0.124675 0.097975

diversité.haplotypique 0.0006455 0.0001065

sd.dh 0.00100925 0.00072425

Script:

as.matrix(read.table("Saison.txt")) -> saison

saison

par(mfrow=c(1,2))

barplot(saison[1,],ylab="Diversité génétique",col=c("lightblue", "lightblue", "lightblue", "lightblue"), border = "blue")-> marks

segments(marks,saison[1,]-saison[2,], marks,saison[1,]+saison[2,])

segments(marks-0.1,saison[1,]-saison[2,], marks+0.1,saison[1,]-saison[2,]);segments(marks-0.1,saison[1,]+saison[2,], marks+0.1,saison[1,]+saison[2,])

barplot(saison[3,],ylab="Diversité haplotypique",col=c("lightblue", "lightblue", "lightblue", "lightblue"), border = "blue")-> marks

segments(marks,saison[3,]-saison[4,], marks,saison[3,]+saison[4,])

segments(marks-0.1,saison[3,]-saison[4,], marks+0.1,saison[3,]-saison[4,]);segments(marks-0.1,saison[3,]+saison[4,], marks+0.1,saison[3,]+saison[4,])


bottom of page