Contoh Soal Harapan Hidup Lengkap

Contoh 2.3.1a

contoh231a <- function(age=30, n=50) {
exp.kontinu <- function (t) {
      (1-(0.01*(age+t))^2)/(1-(0.01*age)^2) }
exp2.kontinu <- function(t) {
  t*(1-(0.01*(age+t))^2)/(1-(0.01*age)^2) }
output.exp.kontinu <- integrate (exp.kontinu,
                                 0, n)$value
output.exp2.kontinu <- 2*integrate(exp2.kontinu,
                                   0, n)$value
output.var.kontinu <- output.exp2.kontinu-
                      (output.exp.kontinu^2)
cat("Nilai Ekspetasi kontinu     =",
    output.exp.kontinu, "\n")
cat("Nilai (Ekspetasi^2) kontinu =",
    output.exp2.kontinu, "\n")
cat("Nilai Variansi kontinu      =",
    output.var.kontinu, "\n") }

contoh231a()

## Nilai Ekspetasi kontinu     = 37.17949 
## Nilai (Ekspetasi^2) kontinu = 1607.143 
## Nilai Variansi kontinu      = 224.8286

Contoh 2.3.1b

contoh231b <- function(age=30, n=50) {
k = 1 : n
exp.diskret <- (1-(0.01*(age+k))^2)/(1-(0.01*age)^2)
exp2.diskret <- (2*k-1)*(1-(0.01*(age+k))^2)/
                (1-(0.01*age)^2)
output.exp.diskret <- sum(exp.diskret)
output.exp2.diskret <- sum(exp2.diskret)
output.var.diskret <- output.exp2.diskret-
                      (output.exp.diskret^2)
cat("Nilai Ekspetasi diskret        =",
    output.exp.diskret, "\n")
cat("Nilai (Ekspetasi^2) diskret    =",
    output.exp2.diskret,"\n")
cat("Nilai Variansi diskret         =",
    output.var.diskret, "\n") }

contoh231b()

## Nilai Ekspetasi diskret        = 36.87637 
## Nilai (Ekspetasi^2) diskret    = 1589.799 
## Nilai Variansi diskret         = 229.9325