#Distribus exponensial
# Untuk P(X = 4)
eks <- function(x, lambda) {
lambda * exp(-lambda * x)
}
eks(4, 0.16)
## [1] 0.08436679
# Untuk P(X <= 4)
cdfexp <- function(x, lambda = 0.16) {
pdf <- function(t) {
lambda * exp(-lambda * t)
}
hasil <- integrate(pdf, lower = 0, upper = x)
return(hasil$value)
}
cdfexp(4, 0.16)
## [1] 0.4727076
#Distrbusi Uniform
uniform<-function (x, a=0,b=1){1/(b-a)}
uniform(6,a=2,b=10)
## [1] 0.125
dunif(6,min=2,max=10)
## [1] 0.125
# Fungsi Distribusi Normal
norm <- function(x, miu, sigma) {(1 / (sigma * sqrt(2 * pi))) *exp(-1/2 * ((x - miu) / sigma)^2)}
norm(5, 4, 1)
## [1] 0.2419707
# Fungsi Bawaan R
dnorm(5, mean = 4, sd = 1)
## [1] 0.2419707