Topik : membuat fungsi sederhana
Kasus : Plot Min – Max Temperature Jawa Timur 2000 – 2003
Berikut ini adalah script sederhana untuk melakukan plot time series temperatur udara di Propinsi Jawa Timur berdasarkan data dari Stasiun Metereologi Juanda Surabaya yang dimuat dalam publikasi BPS – Jawa Timur dalam angka tahun 2000 – 2003. Di mana struktur data yang disediakan adalah data minimum dan maksimum, sehingga bentuk visualisasi datanya hendaknya mencerminkan kondisi tersebut.
Menggunakan fungsi segments
bartemp <- function(at,min,max,delta,col="black") {
segments(at-delta,min,at+delta,min,col=col)
segments(at,min,at,max,col=col)
segments(at-delta,max,at+delta,max,col=col)
}temp < - read.csv("/home/unung/work/temp.csv")
plot(c(1:48),ylim=c(min(temp[2])-5,max(temp[1]+5)),type="n",axes=F,xlab="",
ylab="temperatur (C)")
axis(1,seq(1,48,6),c("Jan-2000","Jul-2000","Jan-2001","Jul-2001",
"Jan-2002","Jul-2002","Jan-2003","Jul-2003"),cex=0.6)
box()
axis(2,cex=0.4)
bartemp(at=c(1:48),temp[[2]],temp[[1]],0.2,col="blue")
title(main="Plot Series Temperatur Udara Jawa Timur")
Hasil grafiknya adalah sebagai berikut :
Nach, rekan-rekan semua dapat memodifikasi fungsi di atas dan menambahkan fungsi-fungsi lain untuk mempercantik plot. Plot ini dapat anda gunakan jika menemui kasus data min – max.
Oke, selamat mencoba ya 
Salam,
Unung Istopo
