# Hitung Z score & probability
# Sunu Pradana
# Samarinda, Indonesia, Agustus 2014
#
# Creative Commons Attribution-ShareAlike 4.0
#
# Usage (contoh penggunaan):
# mungkin(1, 0, 1, -2, 2)
# mungkin(5, 200, 31, 187.4, 221.2)
#
# sd1 = Population standard deviation
# xbar1 = Population mean
# n1 = Jumlah sample (sample size)
# value1 = left value
# value2 = right value
# Untuk nilai ekstrim (for extrme value), gunakan minimal:
# kiri (left) = -1E99
# kanan (right) = 1E99
mungkin=function(sd1,xbar1,n1,value1,value2)
{
es1 = sd1/(sqrt(n1))
z1 = (value1-xbar1)/es1
z2 = (value2-xbar1)/es1
kem=pnorm(z2)-pnorm(z1)
cat("\n")
cat( sprintf("%08s % 20.4f
\r%08s % 20.4f
\r%08s % 20.4g
\r\n%08s % 20.4E
\r%08s % 20.4E
\n" ,
"sd1:",sd1,"xbar1:",xbar1,"n1:",n1,"value1:",value1,"value2:",value2), fill = TRUE);
#cat(" \n ");
cat( sprintf("%-31s %1s % 10.4E %10s % 5.2f ","Jumlah sample", ":", n1, "[trunc2:]", round(n1,2)), fill = TRUE )
cat( sprintf("%-31s %1s % 10.4E %10s % 5.2f ","Population mean", ":", xbar1, "[trunc2:]", round(xbar1,2)), fill = TRUE )
cat( sprintf("%-31s %1s % 10.4E %10s % 5.2f ","Population standard deviation", ":", sd1, "[trunc2:]", round(sd1,2)), fill = TRUE )
cat("\n")
cat( sprintf("%-31s %1s % 10.4E %10s % 5.2E ","Nilai batas kiri (left)", ":", value1, "[trunc2:]", round(value1,2)), fill = TRUE )
cat( sprintf("%-31s %1s % 10.4E %10s % 5.2E ","Nilai batas kanan (right)", ":", value2, "[trunc2:]", round(value2,2)), fill = TRUE )
cat("\n")
cat( sprintf("%-31s %1s % 10.4E %10s % 5.2E ","Standard Error", ":", es1, "[trunc2:]", round(es1,2)), fill = TRUE )
cat("\n")
cat( sprintf("%-31s %1s % 10.4g %10s % 5.2E ","Nilai Z kritis kiri (left)", ":", z1, "[trunc2:]", round(z1,2)), fill = TRUE )
cat( sprintf("%-31s %1s % 10.4E %10s % 5.2E ","Nilai Z kritis kanan (right)", ":", z2, "[trunc2:]", round(z2,2)), fill = TRUE )
cat("\n")
cat( sprintf("%-31s %1s % 10.4E %10s % 5.2f ","Nilai kemungkinan (probability)", ":", kem, "[trunc2:]", round(kem,2)), fill = TRUE )
cat( sprintf("%-31s %1s % 10.4E %10s % 5.2f ","Nilai kemungkinan dalam persen", ":", kem*100,"[trunc2:]", round(kem*100,2)),fill = TRUE )
cat("\n")
}
Comments