# ==============================================================================
# Simulacao: erro de medida na veri?vel dependente
# Update: 13/09/2022
# ==============================================================================

rm(list=ls())

# Define seed
set.seed(1234567)

# Define os valores verdadeiro dos paramentros: intercepto e inclinacao
b0 <- 1; b1 <- 0.5

# Define vetores para guardar os resultados
b1hat <- numeric(10000)
b1hat.me <-numeric(10000)

# Sorteia uma amostra para X
x <- rnorm(1000,4,1)

# Estima o modelo r vezes

for (j in 1:10000) {
  
  # sorteia uma amostra de u
  u <- rnorm(1000)
  
  # calcula ystar:
  ystar <- b0 + b1*x + u
  
  # Estima o modelo e guarda os resultados 
  bhat <- coef(lm(ystar~x))
  b1hat[j] <- bhat["x"]
  
  # Define o erro de medida de y:
  e0 <- rnorm(1000)
  y <- ystar+e0
  
  #Estima o modelo com erro de medida e guarda os resultados 
  bhat.me <- coef( lm(y~x))
  b1hat.me[j] <- bhat.me["x"]

}

hist(b1hat, breaks = 30, col = 'blue', xlim = c(0.2,.8))
hist(b1hat.me, breaks = 30, col = 'transparent', add = TRUE, xlim = c(0.2,.8))

# M?dia com e sem o erro de medida
c(mean(b1hat), mean(b1hat.me))

# Vari?ncia com e sem o erro de medida
c(var(b1hat), var(b1hat.me))