#Untuk P(X=4)
eks <- function(x, lambda) {
lambda * exp(-lambda*x)
}
eks(4, 1/6)
## [1] 0.08556952
#Untuk P(X≤4)
cdfexp <- function(x, lambda = 1/6) {
pdf <- function(t) {
lambda * exp(-lambda*t)
}
hasil <- integrate(pdf, lower = 0, upper = x)
return(hasil$value)
}
cdfexp(4, 1/6)
## [1] 0.4865829
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
norm <- function(x, miu, sigma) {
1 / (sigma * sqrt(2 * pi)) * exp(-1 / 2 * ((x - miu) / sigma)^2)
}
norm(4, 5, 1)
## [1] 0.2419707
dnorm(4, mean = 5, sd = 1)
## [1] 0.2419707