LAB 02 - 7.10.2009

  Pokracujeme v praci s viacrozmernymi datami "cars" z minuleho tyzdna. 
  Cielom je vytvorit vizualizaciu typu "timeline", zobrazujucu niekolko 
  premennych naraz tak aby sa dali vidiet trendy v case. Pripravil
  som pre Vas pociatocne kroky, Vasou ulohou bude ich rozvinut dalej.
  Pripravena je ukazka konstrukcie jednoduchej linearnej casovej osi
  a grafickych prvkov nad nou. Nasleduje idea spiralovitej casovej osi,
  s ktorou budete mat za ulohu dalej pracovat a vytvorit vlastne zobrazenie 
  na nej. Mozete napriklad menit podla dat hrubku a farbu spiraly, umiestnovat 
  nu body roznej velkosti alebo doplnit ju obrazkami ilustrujucimi trendy v
  datach.
  Vysledky prace si porovname na projektore. Za uspesne vytvoreny "timeline"
  dostanete 5 bodov, bud priamo na cviceni alebo do zaciatku buduceho (pocas 
  semestra budu este 2 takto hodnotene cvicenia).

  Ak Vam budu pri praci chybat nejake balicky v R, mali by ste byt schopni
  si ich nainstalovat prikazom "Install packages.." v menu Packages.
  Budu instalovane vo Vasom domovskom adresari.

  Pokial nemate uchovanu historiu z minuleho tyzdna, vykonajte nasledovne
  prikazy (nezabudnite na "Change dir..."):

cardata <- read.table("cars.data")
names(cardata) <- c("mpg","cylinders","displacement","horsepower","weight","acceleration","model.year","origin")
carnames <- read.table("cars.names")
names(carnames) <- c("name")
cars <- data.frame(c(carnames,cardata))

  Naucite sa pouzivat graficke prikazy z balika grid. Pracu z jeho 
  funkciami vykonavame vzdy v ramci nejakej podoblasti zobrazovacej plochy -
  tzv. viewport, ktory ma vlastnu velkost, orientaciu a suradnicovy system.
  Tieto podoblasti sa aktivuju prikazom pushViewport() a deaktivuju 
  popViewport(), po ktorom nastane navrat do predchadzajucej nadradenej
  podoblasti. Vytvorime si oblast s okrajmi:

library(grid)
grid.newpage()
pushViewport(plotViewport())

  V nom si vytvorime podoblast pre data od roku 1969 do 1983

pushViewport(viewport(xscale=c(1969,1983)))
years <- c(1969:1983)

  Vykreslime body pre jednotlive roky, kazdy desiaty vyraznejsie

grid.circle(years,0,0.15,default.units="native",gp=gpar(fill="red"))
syears <- years[which(years/10 == floor(years/10))]
grid.circle(syears,0,0.25,default.units="native",gp=gpar(fill="blue"))

  Popiseme modre body rokmi

grid.text(syears,syears,-0.1,default.units="native",rot=90,gp=gpar(fontsize=14,fontface="bold"))

  Aby sme mohli opakovat rozne vizualizacie s tou istou casovou osou, ulozime
  si predchadzajuce prikazy ako funkciu initTimeline()

initTimeline <- function(years){
  ...
}

  Vytvorime farebnu plochu, ktora bude znazornovat spotrebu aut v danom roku
  Najprv ale potrebujeme funkciu, ktora spocita pocet aut pre dany rok

countCars <- function(year){
  length(which(cars$model.year == year-1900))
}

  Vyskusajte, ci funguje:

countCars(1971)

  vysledok by mal byt 29. Vyborne! 

  podobna funkcia nam spocita priemerny pocet valcov pre dany rok

countCylinders <- function(year){
  mean(cars$cylinders[which(cars$model.year == year-1900)])
}

countCylinders(1971)

  Je vysledok 5.517241?

a priemernu spotrebu

countMpg <- function (year)
{
  if(countCars(year)>0){
    mean(cars$mpg[which(cars$model.year == year-1900)],na.rm=TRUE)
  } else {
    0
  }
}

  Mozeme zobrazit oblast podla poctu aut v danom roku

grid.polygon(years,mapply(countMpg,years)/100+0.05,default.units="native",gp=gpar(fill="grey"))

  Nad nou zobrazime prislusny pocet valcov tak, aby sa zdoraznilo, ze jedna 
  klesa zatial co druha stupa. Potrebujeme funkciu drawCyl()

drawCyl <- function(){
  grid.rect(x=0.5,y=0.7,width=0.8,height=0.6,gp=gpar(fill="brown"))
  grid.rect(x=0.5,y=0.9,width=0.8,height=0.08,gp=gpar(fill="brown"))
  grid.rect(x=0.5,y=0.3,width=0.1,height=0.2,gp=gpar(fill="brown"))
  grid.circle(0.5,0.1,0.2,gp=gpar(fill="grey"))
}

  ... a mozeme zacat vykreslovat, viewport vzdy nastavime na spravne miesto 
  pomocou x= a y=

drawSymbols <- function(year){
  n <- round(countCylinders(year))
  for(i in 1:n){
    pushViewport(viewport(x=year,y=1.0-i*0.075,width=0.5,height=0.06,default.units="native"))
    drawCyl()
    popViewport()
  }
}

  Funkcia mapply aplikuje funkciu drawSymbols na vsetky prvky vektora (roky) 

mapply(drawSymbols,c(1970:1982))

  Pokuste sa dokoncit nasledovne vyobrazenie casovych trendov, kde casova os
  je vytvorena vo forme spiraly:

  Prevod roka na polohu na spirale v suradnicovej sustave:

time2pos <- function (time,size=50,spacing=30,density=7.3)
{
  dif <- (time-1950)/density
  angle <- 2*pi*dif
  n <- floor(dif)
  x <- sin(angle)
  y <- cos(angle)
  as.table(trunc(size*c(x,y) * angle*spacing/2/pi))
}

  Suradnice pre jednotlive roky ulozime do d:

d <- data.frame(years,t(mapply(time2pos,years))

  Funkcia spiralplot vykresli roky v spirale na plochu:

spiralplot <- function (years)
{
  d <- data.frame(years,t(mapply(time2pos,years)))
  plot(d[,2:3],type="b", col=c("yellow","orange","red","brown","black"))
  text(d[,2:3]+500,labels=d[,1],
col=c("yellow","orange","red","brown","black"))
}

  S pouzitim grafickych primitiv balicka grid mame k dispozicii krajsie 
  graficke prvky:

newspiralplot <- function (years)
{
  d <- data.frame(years,t(mapply(time2pos,years)))
  grid.newpage()
  pushViewport(plotViewport(xscale=c(-3000,3000),yscale=c(-3000,3000)))
  pushViewport(dataViewport(d[,2],d[,3]))
  grid.circle(d[,2],d[,3],100,default.units="native",gp=gpar(fill="red"))
  grid.text(d[,1],d[,2],d[,3],default.units="native",just=c("left","bottom"))
  grid.xspline(d[,2],d[,3],default.units="native",shape=-0.8,open=TRUE)
}

  Funkcia, ktora sa moze hodit pri zobrazovani vykresli obrazok vo formate
  GIF na lubovolne miesto suradnicoveho systemu 

plotImage <- function (image,xpos,ypos,size)
{
  require(pixmap)
  require(RGraphics)
  require(colorspace)
  imagemat <- getChannels(image,colors="all")
  pushViewport(viewport(x=xpos,y=ypos,width=size,height=size,default.units="native"))
  grid.imageFun(130,200,byrow=FALSE,rev(t(hex(RGB(as.vector(imagemat[,,1]),as.vector(imagemat[,,2]),as.vector(imagemat[,,3]))))))
  popViewport()
}

f1 <- read.pnm("car_f1.pnm")

  Cvicny obrazok najdete medzi datami v ISe. Ine obrazky je potrebne do
formatu pnm konvertovat pomocou programov balika ImageMagick

plotImage(f1,1000,2000,3000)