rm(list=ls()) blaa <- c('#084594','#2171b5','#4292c6','#6baed6','#c6dbef') # Fem blaatoner graa <- c('#4D4D4D','#737373','#A6A6A6','#DADADA') # Fire graatoner kontrast <- '#FF7260'; moerkeblaa <- '#000059' # Spesialfarger #### SoYLE ######################################################################################### #### DATA #### dummy <- as.matrix(rbind( 'Sverige totalt'= 46, 'Norge totalt'= 41, 'Helse Midt'= 47, 'Helse Vest'= 45, 'Helse Nord'= 39, 'Helse Sor-ost'= 38) ) KI.opp <- dummy + c(1.5,2,4,5,5,6); KI.ned <- dummy - c(1.5,2,4,5,5,6) # Konfidensintervall antall <- dummy; antall[1:6,1] <- c(3000,1000,1000,1000,6000,8000) # Ant. (N) for soylene overskrift <- c('Søylediagram','med sorterte søyler og 95%-KI') # MERK: Overskrifter skal ikke med på grafikkfiler, kun benyttet som illustrasjon ant.grp <- dim(dummy)[1] # Ant. grupper x.lab <- 'Prosent' # X-akse x.axs <- seq(0,52,10) x.lim <- c(0,52) # Plottets X-akse #### PLOTT #### windows( width=5, height=2.8 ) # Plottdimensjon par( mar=c(2.2,5.5,2,1.5) ) # Marger temp <- barplot( t(dummy)[,ant.grp:1], beside=T, las=1, main = '', xlim=x.lim, col=0, # Dummyplott horiz=T, names.arg='', axes=F, border=NA) title( overskrift[1], line=1.2, cex.main=1, font.main=2, col.main=moerkeblaa) # Overskrifter mtext( overskrift[2], line=0.1, col=moerkeblaa, cex=0.9 ) mtext( x.lab, 1, cex=0.8, line=1, col=moerkeblaa) # Aksemerkinger axis( 1, x.axs, x.axs, cex.axis=0.75, lwd=0, tck=-0.02, lwd.ticks=1, # Akser col=moerkeblaa, col.axis=moerkeblaa, padj=-2) temp <- barplot( t(dummy)[,ant.grp:1], beside=T, las=1, main = '', xlim=x.lim, # Selve plottet col=c(rep(blaa[2],4),rep(blaa[4],2)), horiz=T, names.arg=rep('',ant.grp), axes=F, border=NA, add=T) text( x=0, y=temp, labels=paste0('N = ', antall), col='white', cex=0.7, pos=4, font=2) # N for soylene text( x=rev(dummy), y=temp+0.3, labels=rev(dummy), col=moerkeblaa, # Soyletall font=2, cex=0.7, pos=4, offset=0.15) mtext( rev(row.names(dummy)), side=2, cex=0.9, line=0.3, las=1, # Soylenavnene at=temp, col=moerkeblaa) arrows( x0=KI.opp, y0=rev(temp), x1=KI.ned, y1=rev(temp), code=0, col=blaa[1] ) # Konf.intervall savePlot('Soyle2.pdf', type='pdf') # Lagrer fil #### KUTTET ######################################################################################## #### DATA #### x.lim <- c(30,50) # ellers samme data x.axs <- seq(30,50,5) # som tidligere #### PLOTT #### windows( width=5, height=2.5 ) # Plottdimensjon par( mar=c(2.2,5.5,2,1.5) ) # Marger temp <- barplot( t(dummy)[,ant.grp:1], beside=T, las=1, main='', xlim=x.lim, # Dummyplott col=0, horiz=T, names.arg='', axes=F, border=NA) title( 'Kuttet soylediagram', line=1.2, cex.main=1.1, font.main=2, col.main=moerkeblaa) # MERK: Overskrifter skal ikke med på grafikkfiler mtext( 'sorterte punkter og 95%-KI', line=0.1, col=moerkeblaa, cex=0.9 ) mtext( x.lab, 1, font=1, cex=0.8, line=1, col=moerkeblaa) # Aksemerking abline( v=c(39,41,43), col=graa[c(3,2,2)], lty=3) # Landsnitt axis( 1, x.axs, x.axs, cex.axis=0.75, lwd=0, tck=-0.02, lwd.ticks=1, # Akser col=moerkeblaa, col.axis=moerkeblaa, padj=-2) points( x=dummy, y=rev(temp), pch=19, col=c(rep(blaa[4],2),rep(blaa[1],4)), cex=1.4 ) # Selve plottet text( x=rev(dummy), y=temp+0.3, labels=rev(dummy), col=moerkeblaa, font=2, cex=0.7, # Punktverdier pos=4, offset=0.3) mtext( rev(row.names(dummy)), side=2, cex=0.9, line=0.3, las=1, # Soylenavn at=temp, col=moerkeblaa) arrows( x0=KI.opp, y0=rev(temp), x1=KI.ned, y1=rev(temp), code=0, # Konf.intervall lwd=1.5, col=c(rep(blaa[4],2),rep(blaa[1],4))) #par(xpd=T); legend( locator(1), c('Beskriv.', '95%-KI'), pch=c(19,-1), lty=c(-1,1), # Merkelapp # pt.cex=1.3, cex=0.8, bty='n', col=blaa[1], lwd=c(0,1.5), # text.col=moerkeblaa ); par(xpd=F) savePlot('Kuttet.pdf', type='pdf') # Lagrer fil #### GRUPPERT ###################################################################################### #### DATA #### dummy <- cbind( 'Soyle 2' = c(65,64,58), # Dummydata 'Soyle 1' = c(74,68,60) ) rownames(dummy) <- c('Indikator 1','Indikator 2','Indikator 3') # Gruppenavn maalverdier <- cbind(c(70,80,90),c(70,80,90)) # Malsetninger ant.grp <- dim(dummy)[1] # Ant. grupper ant.var <- dim(dummy)[2] # Ant. variabler x.lab <- 'Prosent' # X-akse x.axs <- seq(0,100,20) x.lim <- c(0,100) # Plottets X-aks #### PLOTT #### windows( width=5, height=2.4 ) # Plottdimensjon par( mar=c(2.2,5.5,2.2,1.5) ) # Marger temp <- barplot( t(maalverdier), beside=T, las=1, main = '', xlim=x.lim, col=0, # Dummyplott horiz=T, names.arg=rep('',ant.grp), axes=F, border=NA) # (evt. col='gray90') title( 'Gruppert soylediagram', line=1.4, cex.main=1.1, font.main=2, col.main=moerkeblaa) # MERK: Overskrifter skal ikke med på grafikkfiler mtext( 'med målsetninger', line=0.3, col=moerkeblaa ) mtext( x.lab, 1, font=1, cex=0.8, line=1, col=moerkeblaa ) # Aksemerking axis( 1, x.axs, x.axs, cex.axis=0.75, lwd=0, tck=-0.02, lwd.ticks=1, # Akser col=moerkeblaa, col.axis=moerkeblaa, padj=-2) for(grp in 1:ant.grp) { # Maalsetninger lines( x=maalverdier[grp,], y=temp[,grp]+c(-0.5,0.5),col=moerkeblaa, lwd=1.5 ) par(xpd=T); text( x=maalverdier[grp,],y=temp[ant.var,grp]+0.1, label=paste0(maalverdier[grp,1],'%'), col=moerkeblaa, cex=0.6, pos=3 ); par(xpd=F) } temp <- barplot( t(dummy[ant.grp:1,]), beside=T, las=1, main = '', xlim=x.lim, # Selve plottet col=blaa[3:4], horiz=T, names.arg=rep('',ant.grp), axes=F, border=NA, add=T, space=c(0.03,1)) text( x=t(dummy)[,ant.grp:1], y=temp+0.1, labels=t(dummy)[,ant.grp:1], # Soyleverdier col=moerkeblaa, font=2, cex=0.7, pos=4, offset=0.2) mtext( rev(row.names(dummy)), side=2, line=0.3, las=1, # Soylenavn at=apply(temp,2,mean), col=moerkeblaa) text( x=-2, y=temp[1:2,3], labels=c('Søyle 2','Søyle 1'), col='white', cex=0.7, pos=4, font=2) par(xpd=T); legend( locator(1), 'Målsetting', cex=0.7, bty='n', text.col=moerkeblaa); par(xpd=F) savePlot('Gruppert.pdf', type='pdf') # Lagrer fil #### LINJE ######################################################################################### #### DATA #### dummy <- as.matrix(rbind( 'Variabel 1' = c(18,20,17,18,18,20), # Dummydata 'Variabel 2' = c( 7,11,10,13,11,14)) ) KI.opp <- dummy + c(2,2.5); KI.ned <- dummy - c(2,2.5) # Konfidensintervall overskrift <- c('Linjediagram', 'Utvikling over tid') # MERK: Overskrifter skal ikke med på grafikkfiler ant.grp <- dim(dummy)[1] # Ant. linjer ant.pkt <- dim(dummy)[2] # Ant. punkter x.lab <- 'Tidspunkt'; y.lab <- 'Verdi' # X- og Y-akse x.lim <- c(1,6.1) ; y.lim <- c(0,25.5) x.axs <- 1:6 ; y.axs <- seq(0,25,5) #### PLOTT #### windows( width=5, height=2.7 ) # Plottdimensjon par( mar=c(2.5,2,2.5,0.6) ) # Marger temp <- plot( dummy[1,], main='', xlab='', ylab='', xlim=x.lim, ylim=y.lim, col=0, axes=F ) # Dummyplott title( overskrift[1], line=1.3, cex.main=1.2, font.main=2, col.main=moerkeblaa) # Overskrift mtext( overskrift[2], line=0.2, col=moerkeblaa ) mtext( x.lab, 1, cex=0.8, line=1.3, col=moerkeblaa) # Aksemerking mtext( y.lab, 2, at=par('usr')[4], padj=-0.7, cex=0.9, col=moerkeblaa, las=1 ) for(grp in 1:ant.grp) { # Konf.int. polygon( x = c(1:ant.pkt, ant.pkt:1), y = c(KI.opp[grp,], rev(KI.ned[grp,])), col = 'gray91', border=0) } axis( 1, x.axs, x.axs, cex.axis=0.9, lwd=0, tck=-0.02, lwd.ticks=1, # Akser col=moerkeblaa, col.axis=moerkeblaa, padj=-1) axis( 2, y.axs, y.axs, cex.axis=0.9, lwd=0, tck=+0.02, lwd.ticks=1, col=moerkeblaa, col.axis=moerkeblaa, hadj=0.3, las=1 ) lines( dummy[1,], type='b', cex=1, pch=16, lty=3, col=blaa[1]) # Linjer/punkt lines( dummy[2,], type='b', cex=0.9, pch=17, lty=3, col=blaa[4]) par(xpd=T) legend( locator(1), 'HF 1', bty='n', text.col=moerkeblaa, text.font=2, cex=0.7) legend( locator(1), 'HF 2', bty='n', text.col=moerkeblaa, text.font=2, cex=0.7) legend( locator(1), '95%-KI',pch=15, pt.cex=2, cex=0.8, bty='n',col='gray87',text.col=moerkeblaa ) par(xpd=F) savePlot('Linje.pdf', type='pdf') # Lagrer fil for(grp in 1:ant.grp) { # Alternativt konf.int. arrows( y0 = KI.opp[grp,], x0 = 1:ant.pkt, y1 = KI.ned[grp,], x1 = 1:ant.pkt, length=0, col=blaa[c(1,4)][grp]) } #### LINJE ALTERNATIV############################################################################### #### PLOTT #### windows( width=5, height=2.7 ) # Plottdimensjon par( mar=c(2.5,2,2.5,0.6) ) # Marger temp <- plot( dummy[1,], main='', xlab='', ylab='', xlim=x.lim, ylim=y.lim, col=0, axes=F ) # Dummyplott title( overskrift[1], line=1.3, cex.main=1.2, font.main=2, col.main=moerkeblaa) # Overskrift mtext( overskrift[2], line=0.2, col=moerkeblaa ) mtext( x.lab, 1, cex=0.8, line=1.3, col=moerkeblaa) # Aksemerking mtext( y.lab, 2, at=par('usr')[4], padj=-0.7, cex=0.9, col=moerkeblaa, las=1 ) axis( 1, x.axs, x.axs, cex.axis=0.9, lwd=0, tck=-0.02, lwd.ticks=1, # Akser col=moerkeblaa, col.axis=moerkeblaa, padj=-1) axis( 2, y.axs, y.axs, cex.axis=0.9, lwd=0, tck=+0.02, lwd.ticks=1, col=moerkeblaa, col.axis=moerkeblaa, hadj=0.3, las=1 ) lines( dummy[1,], type='l', cex=1, pch=16, lty=1, col=blaa[1]) # Linjer/punkt points( y=dummy[1,],x=1:6, cex=0.8, pch=16, lty=1, col=blaa[1]) lines( dummy[2,], type='l', cex=1, pch=17, lty=1, col=blaa[4]) points( y=dummy[2,],x=1:6, cex=0.9, pch=17, lty=1, col=blaa[4]) par(xpd=T) legend( locator(1), 'HF 1', bty='n', text.col=moerkeblaa, text.font=2, cex=0.7) legend( locator(1), 'HF 2', bty='n', text.col=moerkeblaa, text.font=2, cex=0.7) par(xpd=F) savePlot('Linje2.pdf', type='pdf')