HowTo: Primer uporabe barplot A. Blejec 23. april 2009 Povzetek Opisan je naˇcinprikaza dveh porazdelitev, ki omogoˇcalažjo primerjavo razlik. Kazalo 1 Uvod1 2 Spremljanje dnevnih ˇcasopisov2 3 Spremljanje televizijskih programov3 4 Priljubljena alternativa :)4 5 Delno prekriti trakovi5 Dodatek6 A R ukazi in funkcije6 1 Uvod V soboto 11. aprila 2009 so bili v Veˇceru objavljeni grafi o spremljanju medijev pri študentih in poslancih. Grafi so bili urejeni po padajoˇcipogostosti uporabe, loˇcenoza študente in poslance. Primerjava uporabe ni enostavna, ker je treba za primerjavo v grafih posebej poiskati stolpec za isti medij. V prikazu so tudi grafi za poslance urejeni glede na pogostost uporabe pri študentih. 1 2 Spremljanje dnevnih ˇcasopisov studenti poslanci Žurnal 24 36.7 1 Delo 18.2 25 Slovenske novice 12.3 4 Dnevnik 11.7 19 Vecerˇ 8.2 16 Primorske novice 4.5 7 Finance 2.3 13 drugo 1.4 8 ni navedel 0.1 1 Ekipa 0.1 0 Indirekt 0.0 3 študenti poslanci Žurnal 24 Delo Delo Dnevnik Slovenske novice Ve..er Dnevnik Finance Ve..er drugo Primorske novice Primorske novice Finance Slovenske novice drugo Indirekt Ekipa Žurnal 24 ni navedel ni navedel Indirekt Ekipa 0 5 15 25 35 0 5 10 15 20 25 Žurnal 24 Žurnal 24 Delo Delo Slovenske novice Slovenske novice Dnevnik Dnevnik Ve..er Ve..er Primorske novice Primorske novice Finance Finance drugo drugo Ekipa Ekipa ni navedel ni navedel Indirekt Indirekt 0 5 15 25 35 0 5 15 25 35 Slika 1: Spremljanje dnevnih ˇcasopisov Na zgornjih dveh panelih slike1 je primerjava razliˇcnostuporabe dnevnih ˇca- sopisov vidna s prebiranjem imen ˇcasopisov, na spodnjih dveh pa iz trakov grafa, ki se za poslance ne spreminjajo enako kot za študente. Primerjaj tudi sliko2. 2 3 Spremljanje televizijskih programov studenti poslanci Pop TV 68.3 25 TV Slovenija 1 31.0 32 Kanal A 18.3 10 brez odgovora 14.0 0 TV Slovenija 2 11.3 18 drugo 9.5 2 CNN 3.5 0 TV 3 2.7 3 Discovery 1.6 0 Info TV 1.6 7 BBC World 1.2 0 lokalna TV 0.2 4 študenti poslanci Pop TV TV Slovenija 1 TV Slovenija 1 Pop TV Kanal A TV Slovenija 2 brez odgovora Kanal A TV Slovenija 2 Info TV drugo lokalna TV CNN TV 3 TV 3 drugo Info TV brez odgovora Discovery CNN BBC World Discovery lokalna TV BBC World 0 20 40 60 0 5 10 20 30 Pop TV Pop TV TV Slovenija 1 TV Slovenija 1 Kanal A Kanal A brez odgovora brez odgovora TV Slovenija 2 TV Slovenija 2 drugo drugo CNN CNN TV 3 TV 3 Info TV Info TV Discovery Discovery BBC World BBC World lokalna TV lokalna TV 0 20 40 60 0 20 40 60 Slika 2: Spremljanje televizijskih programov 3 4 Priljubljena alternativa :) Indirekt ni navedel Ekipa drugo Finance Primorske novice studenti 0 0.1 0.1 1.4 2.3 4.5 poslanci 3 1.0 0.0 8.0 13.0 7.0 Vecerˇ Dnevnik Slovenske novice Delo Žurnal 24 studenti 8.2 11.7 12.3 18.2 36.7 poslanci 16.0 19.0 4.0 25.0 1.0 Žurnal 24 Delo Slovenske novice Dnevnik Ve..er Primorske novice Finance drugo Ekipa ni navedel Indirekt 0 5 10 15 20 25 30 35 Slika 3: Primerjalni trakovi: priljubljena alternativa prikaza podatkov 4 5 Delno prekriti trakovi Z uporabo argumenta space=c(-0.4,0.5) lahko stolpce delno prekrijemo: Žurnal 24 Delo Slovenske novice Dnevnik Ve..er Primorske novice Finance drugo Ekipa ni navedel Indirekt 0 5 10 15 20 25 30 35 Slika 4: Primerjalni trakovi: referenˇcnaporazdelitev (študentje, oranžno) je lahko delno skrita za primerjano porazdelitvijo. 5 A R ukazi in funkcije ################################################### ### chunk number 1: initialize ################################################### options(width=70) library(Hmisc) #library(xlsReadWrite) ################################################### ### chunk number 2: function barplot ################################################### barplots <- function(X,i1=1,i2=2){ par(mfrow=c(2,2),las=1,mar=c(4,8,1,1)) X <- X[order(X[,i1]),] labels <- dimnames(X)[[1]] x <- X[,i1] y <- X[,i2] names(x) <- labels names(y) <- labels barplot(x,horiz=TRUE,col="orange",main="študenti") barplot(y[order(y)],horiz=TRUE,col="brown",main="poslanci") barplot(x,horiz=TRUE,col="orange",xlim=c(0,max(X))) barplot(y,horiz=TRUE,col="brown",xlim=c(0,max(X))) invisible(NULL) } ################################################### ### chunk number 3: read casopis data ################################################### casopis <- read.delim("../data/dnevniCasopis.txt",row.names=1) casopis ################################################### ### chunk number 4: barplot-casopis ################################################### barplots(casopis) ################################################### ### chunk number 5: read TV data ################################################### TV <- read.delim("../data/TV.txt",row.names=1) TV ################################################### ### chunk number 6: barplot-TV ################################################### 6 barplots(TV) ################################################### ### chunk number 7: barplot-beside ################################################### par(mar=c(4,8,2,1)) X <- as.matrix(t(casopis[order(casopis[,1]),])) X barplot(X, beside=TRUE, horiz=TRUE, col=c("orange","brown"), las=1 ) ################################################### ### chunk number 8: barplot-overlapped ################################################### par(mar=c(4,8,2,1)) X <- as.matrix(t(casopis[order(casopis[,1]),])) barplot(X, beside=TRUE, horiz=TRUE, col=c("orange","brown"), las=1, space=c(-.4,.5) ) ################################################### ### chunk number 9: ################################################### rCode <- readLines("../doc/HowTo-barplot.R") rCode <- paste(rCode,collapse="\n") cat(rCode) ################################################### ### chunk number 10: sessionInfo ################################################### cat(win.version(),"\n") toLatex(sessionInfo()) 7 SessionInfo Windows XP (build 2600) Service Pack 3 • R version 2.8.0 (2008-10-20), i386-pc-mingw32 • Locale: LC_COLLATE=Slovenian_Slovenia.1250;LC_CTYPE=Slovenian_Slovenia.1250;LC_MONETARY=Slovenian_Slovenia.1250;LC_NUMERIC=C;LC_TIME=Slovenian_Slovenia.1250 • Base packages: base, datasets, graphics, grDevices, methods, stats, utils • Other packages: Hmisc 3.5-2 • Loaded via a namespace (and not attached): cluster 1.11.13, grid 2.8.0, lat- tice 0.17-22 Revision Id: HowTo-barplot.Rnw 11 2009-04-23 07:15:22Z ablejec (c) A. Blejec 8.