Analisis Angka Kematian Bayi dan Angka Harapan Hidup Kab/Kota di JawaTimur dengan Menggunakan Library Gregmisc

Archive for the ‘R Scripts’ Category

Analisis Angka Kematian Bayi dan Angka Harapan Hidup Kab/Kota di JawaTimur dengan Menggunakan Library Gregmisc

June 15th, 2005 by Unung Istopo

Untuk analisis data ini kita menggunakan data di enciety data management.

Berikut ini script yang dapat kita gunakan :

# Mengambil data dari enciety Data Manajemen
read.csv("http://www.enciety.com/edm/manajemen/data/AKB%20dan%20AHH%20Kab%20Kota%201999%20-%202001%20-%202003%20Jawa%20Timur.csv") -> data.AKB.AHH
require(gregmisc)
# Misal untuk analisis tahun 2003
x <- data.frame(data.AKB.AHH[1:37,c(2,5,8)],row.names=1)
as.matrix(x) -> x
rc <- rainbow(nrow(x), start=0, end=.3)
cc <- rainbow(ncol(x), start=0, end=.3)
hv <- heatmap.2(x, col = cm.colors(256), scale="column", RowSideColors = rc, ColSideColors = cc, margin=c(5,10), xlab = "Angka Kematian Bayi        Angka HarapanHidup", ylab= "", main = "Angka Kematian Bayi & Angka Harapan Hidup 2003", tracecol="black", density="density")
str(hv)  # the two re-ordering index vectors

Hasil Analisisnya sebagai berikut :

Salam,

Unung Istopo H

Tags: heatmap
Posted in R Scripts | Comments Off

Analisis enciety CheckClock Report menggunakan R

May 28th, 2005 by Unung Istopo

Analisis enciety CheckClock Report dengan menggunakan R dapat kita lakukan sebagai berikut :

Mis: Analisis untuk divisi IT: Network & Information Service Provider pada bulan Mei 2005.
1. Akses data ke aplikasi enciety CheckClock Report

read.csv("http://10.20.32.209/sugengrawuh/report_division_csv.php?div=6&bulan=05&tahun=2005&order=start") -> IT.NISP
names(IT.NISP)

2. Analisis Barplot untuk menggambarkan Persentase masuk dan keluar.

barplot(t(as.matrix(IT.NISP[,c(3,5)])),beside=T,ylim=c(0,100),axisname=F,legend=T)  -> barplot.IT.NISP
# Membuat Nama Label
for (i in 1:length(IT.NISP[,1])){      mtext(as.character(IT.NISP[i,1]),at=(barplot.IT.NISP[1,i]+barplot.IT.NISP[2,i])/2,side=1,line=1)
}
# Membuat Abline untuk persentase = 50%
abline(h=seq(10,50,10),col="red",lty=3)
abline(h=seq(60,90,10),col="blue",lty=3)
# Menampilkan abline sebagai background
par(new=T)
barplot(t(as.matrix(IT.NISP[,c(3,5)])),beside=T,ylim=c(0,100),axisname=F,legend=T)
box()
title("enciety Checkclock Report Summary - IT NISP Division - May 2005")
# Membuat text untuk besar persentase masuk dan keluar
for (i in 1:length(IT.NISP[,1])) {
   for (j in 1:2){
      text(barplot.IT.NISP[j,i],IT.NISP[,c(3,5)][i,j]+1,IT.NISP[,c(3,5)][i,j])
   }
}

3. Hasilnya adalah sebagai berikut :

Salam,

Unung Istopo H

Tags: chart, descriptive, web data
Posted in R Scripts | Comments Off

Memulai R dan Analisis Sederhana (2)

May 4th, 2005 by Unung Istopo

Oleh: Unung Istopo H

Untuk materi kedua ini kita akan membahas tentang bagaimana melakukan pengolahan data yang disimpan di SPSS (*.sav). Untuk materi ini data dapat didownload dan silahkan simpan di folder tertentu yang mudah diakses.

Analisis sederhana dengan menggunakan data SPSS.

Membaca data SPSS dengan function foreign dan menampilkan atribut-nya

#Memanggil library foreign yang digunakan untuk membaca file SPSS, MINITAB, SAS, Systat, dll.
require(foreign)
read.spss("/home/unung/work/REncietyCourse/DataSPSS.sav") -> data.SPSS
names(data.SPSS)
#Untuk mengetahui atribut yang ada di data.SPSS maka gunakan function attributes
attributes(data.SPSS) -> attr.data.SPSS
names(attr.data.SPSS)

Analisis Statistika

attach(data.SPSS)
#Diagram Dahan Daun
stem(A1)
#Menghitung Frekuensi dan Tabulasi Silang
table(A2,B1) -> data.1
data.1
names(dimnames(data.1)) <- c("Jenis Kelamin", "Langganan Koran")
data.1
total=sum(data.1)
data.1/total
row.total <- apply(data.1,1,sum)
col.total <- apply(data.1,2,sum)
row.percentage <- sweep(data.1,1,row.total,"/")
col.percentage <- sweep(data.1,1,col.total,"/")
round(col.percentage,3)
#Merubah data.frame ke dalam matrix table
data.frame(data.1) -> data.1.frame
attach(data.1.frame)
tapply(Freq,list(Jenis.Kelamin,Langganan.Koran),c)

Membuat Plot

attach(data.SPSS)
layout(matrix(c(1,3,6,1,4,6,2,4,7,2,5,7),3,4))
par(mar = c(2,2,2,2))
dotchart(table(B2,A2),main="Hobby & Jenis Kelamin")
par(mar = c(2,2,2,2))
dotchart(table(B3,A2),main="Pengeluaran & Jenis Kelamin")
par(mar = c(2,2,2,2))
pie(c(1,2),c("Perempuan","Laki-laki"),col=c("black","gray"),radius=sum(c(1,2))/20+0.65,main=paste("Tidak Berlangganan Koran n = ",sum(c(1,2))))
par(mar = c(2,2,2,2))
barplot(table(A2,B1),col=c("black","gray"),main="Langganan Koran & Jenis Kelamin",beside=T)
pie(c(4,3),c("Perempuan","Laki-laki"),col=c("black","gray"),radius=sum(c(4,3))/20+0.65,main=paste("Berlangganan Koran n =",sum(c(4,3))))
par(mar = c(2,2,2,2))
hist(A1,main="Histogram Usia Responden",col="darkgray")
par(mar = c(2,2,2,2))
boxplot(A1,horizontal=T,main="Boxplot Usia Responden",col="darkgray")

Hasil Grafiknya adalah sebagai berikut :

(*) Dokumentasi R enciety Internal Course, 2 dan 3 Mei 2005

Tags: descriptive
Posted in R Scripts | Comments Off

Statistical Graphics & Data Visualization

Innovation and New Experimentation Model with R

Archive for the ‘R Scripts’ Category