PostestKomstat_Tiara Adistira

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

# 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