#eksponnsial
#dr soal diominta cdf dan pdf
#lambda = 1/miu
#peluang P(X < 5) pdf (tepat)
eks <- function(x, lambda){lambda*exp((-1) * lambda*x)}
eks(15, 0.05)
## [1] 0.02361833
#Untuk P(X <= 15) cdf
cdfexp <- function(x, lambda = 0.05){
pdf <- function(t) {
lambda * exp(-lambda * t)
}
hasil <- integrate(pdf, lower = 0, upper = x)
return(hasil$value)
}
cdfexp(15, 0.05)
## [1] 0.5276334
#Distribusi Uniform
#a = nilai terkecil di interval (batas bwh), b = nilai terbesar di interval (batas atas), x = yang dicari \
#kalau x diluar interval bernilai 0, kalaudidalam pake rumus 1/b-a
uniform <- function (x, a = 2, b = 10) {1/(b-a)}
uniform(6, a = 2, b = 10)
## [1] 0.125
dunif(6, min = 2, max = 10) #karna nyari pdf jadinya pake dunif, kalau nyari cdf pakenya punif
## [1] 0.125
#distribusi normal
#pakai pdf
#mencari nilai tepat
#sigma=sd
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