In [1]:
# rm(list=ls())
options(OutDec = ",")
#==============================================================================
# Exemplo:
#
# Integracao Monte Carlo
#
# Calculando P(Z<1,96), Z ~ N(0,1)
#
#==============================================================================
set.seed(12345) # Semente
In [2]:
#==============================================================================
# Para n = 100 # Cem
#==============================================================================
n <- 100
z <- rnorm(n)
y <- as.integer(z<1.96)
print("===============================")
print(mean(y))
print("===============================")
print(sqrt(var(y)/n)) # Erro de Monte Carlo
[1] "===============================" [1] 0,94 [1] "===============================" [1] 0,02386833
In [3]:
#==============================================================================
# Para n = 1000 # Mil
#==============================================================================
n <- 1000
z <- rnorm(n)
y <- as.integer(z<1.96)
print("===============================")
print(mean(y))
print("===============================")
print(sqrt(var(y)/n)) # Erro de Monte Carlo
[1] "===============================" [1] 0,973 [1] "===============================" [1] 0,005128089
In [4]:
#==============================================================================
# Para n = 10000 # Dez mil
#==============================================================================
n <- 10000
z <- rnorm(n)
y <- as.integer(z<1.96)
print("===============================")
print(mean(y))
print("===============================")
print(sqrt(var(y)/n)) # Erro de Monte Carlo
[1] "===============================" [1] 0,9755 [1] "===============================" [1] 0,001546032
In [5]:
#==============================================================================
# Para n = 100000 # Cem mil
#==============================================================================
n <- 100000
z <- rnorm(n)
y <- as.integer(z<1.96)
print("===============================")
print(mean(y))
print("===============================")
print(sqrt(var(y)/n)) # Erro de Monte Carlo
[1] "===============================" [1] 0,97468 [1] "===============================" [1] 0,0004967811
In [6]:
#==============================================================================
# Para n = 1000000 # Um milhao
#==============================================================================
n <- 1000000
z <- rnorm(n)
y <- as.integer(z<1.96)
print("===============================")
print(mean(y))
print("===============================")
print(sqrt(var(y)/n)) # Erro de Monte Carlo
[1] "===============================" [1] 0,974846 [1] "===============================" [1] 0,0001565928
In [7]:
# rm(list=ls())
#==============================================================================
# graphics.off()
#==============================================================================
# Fim
#==============================================================================