P.S.Pudyastuti's Lovely Rooms…

Learning is a long-life process, so it's never late to learn…

R untuk Analisis Frekuensi (1)

Leave a comment

R menyediakan banyak packages yang dibuat oleh para ahli di bidangnya. Salah satu package yang bermanfaat untuk analisis frekuensi dari data hidrologi adalah lmom package, yang dapat di-download dari CRAN.

Untuk meng-install package yang kita butuhkan, klik menu Packages pada jendela R GUI dan klik sub-menu untuk meng-install.

Untuk contoh di bawah ini, saya menggunakan data debit maksimum tahunan di stasiun pengukuran X. Data saya simpan di file csv (comma-separated value), dengan nama file StasiunX.csv. Dalam contoh ini, saya ingin mendapatkan growth curve untuk mengetahui debit dengan kala ulang tertentu. Di sini saya melakukan standardisasi debit dengan menggunakan nilai median data dengan tujuan untuk mengurangi site-specific factor yang kadang membuat kurva yang kita inginkan tidak sesuai dengan harapan. Distribusi yang saya gunakan adalah distribusi GEV (Generalized Extreme Value).

Berikut ini adalah script yang saya gunakan :

> #–Menggunakan package lmom untuk analisis nilai ekstrem dengan pendekatan L-moments
> library (lmom)

> #–Mengatur folder tempat program dijalankan
> setwd(“D:/tuti/Prog”)

> #–Membaca data debit di stasiun X yang disimpan di file StasiunX.csv
> AM <- read.csv(“StasiunX.csv”,header=TRUE)
> AM

Tahun Tanggal         Q
1 1975 04-Apr-75 1894.96
2 1976 12-Mar-76 1380.64
3 1977 28-Mar-77 1417.32
4 1978 24-Feb-78 1465.72
5 1979 09-May-79 1748.50
6 1980 19-Apr-80 1588.32
7 1981 30-Jan-81 1202.07
8 1982 16-Jan-82 2207.09
9 1983 09-Jan-83 1411.27
10 1984 14-Feb-84 2154.17
11 1985 12-Mar-85 1668.99
12 1986 29-Mar-86 1852.88
13 1987 19-Feb-87 1845.18
14 1988 30-Jan-88 1601.35
15 1989 25-Jan-89 1433.29
16 1990 20-Feb-90 1570.43
17 1991 17-Feb-91 1572.90
18 1992 19-Mar-92 1431.37
19 1993 06-Feb-93 1774.56
20 1994 27-Mar-94 1931.77
21 1995 07-Feb-95 1642.47
22 1996 17-Mar-96 1473.13
23 1997 23-Feb-97 1445.22
24 1998 18-Mar-98 1755.38
25 1999 27-Jan-99 1717.32
26 2000 23-Mar-00 1501.30
27 2001 27-Mar-01 1627.70
28 2002 26-Jan-02 1564.47
29 2003 03-Feb-03 1464.96
30 2004 15-Mar-04 1503.66
31 2005 07-Apr-05 1424.47
32 2006 04-May-06 1449.86

> #–Membaca nama judul tiap kolom dari data di file StasiunX.csv
> names(AM)

[1] “Tahun”   “Tanggal” “Q” 

> #–Menggambar histogram frekuensi dari data debit di stasiun X
> hist(AM$Q,col=”lightblue”,main=”Histogram Data Debit di Stasiun X”,cex.main=1.10,
+          col.main=”blue”,ylim=c(0,20),ylab=”Frekuensi”,xlab=”Debit (m3)”)

> #–Mencari nilai median dari data debit di stasiun X
> Qmed <- median(AM$Q)
> #–Mendefinisikan Q adalah data debit Q di stasiun X
> Q <- c(AM$Q)
> #–Menstandardisasi debit dengan nilai median dari data
> Qstd <- Q/Qmed

> # plot distribusi Generalized Extreme Value (GEV) untuk data debit di stasiun X

> evplot(Qstd,ylim=c(0,2),ylab=”Q/QMed”,pch=17,col=”blue”,main=”Discharge Growth Curve (Standardized by Qmed) \n at X Gauging Station Using GEV Distribution”,col.main=”darkblue”,cex.main=1.00)
> legend(“topleft”, legend=”Stasiun X”,cex=1.0,pch=17,col=”blue”, bg=”white”)

> # menggambar garis fitting distribusi GEV distribution pada plot

> evdistq(quagev, pelgev(samlmu(Qstd)), col=”red”)

Referensi : lmom

Author: P.S.Pudyastuti

I am an academic staff at Civil Engineering Department, Faculty of Engineering, Universitas Muhammadiyah Surakarta - Indonesia. I got my Bachelor Degree (ST) from Gadjah Mada University, Indonesia, and Master Degree (M.Sc) from Loughborough University, England. I am married and have got 3 children.

Leave a Reply

Please log in using one of these methods to post your comment:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s