#Amostragem set.seed(5712) sexo <- c('Homem', 'Mulher') distancia = runif(n = 100000, min = 100, max = 10000) tabela <- data.frame(glicemia = rnorm(n = 100000, mean = 90, sd = 5), renda = rlnorm(n = 100000, meanlog = 7, sdlog = 0.7) * (1.2 - (distancia)/10000), distancia, sexo = sample(x = sexo, size = 100000, replace = T)) plot(y = tabela$renda, x = tabela$distancia) media_verdadeira_renda <- mean(tabela$renda) hist(tabela$renda) #adicionando linha com a média verdadeira da renda abline(v = media_verdadeira_renda, col = 'blue') set.seed(3705676) #substituir pelo seu número USP numero_amostras <- 1000 tamanho_amostras <- 30 medias_amostras_glicemia <- NA medias_amostras_renda <- NA medias_amostras_distancia <- NA for (amostra in 1:numero_amostras) { ind_amostra <- sample(x = 1:100000, size = tamanho_amostras, prob = 1/distancia ) medias_amostras_glicemia[amostra] <- mean( tabela$glicemia[ ind_amostra] ) medias_amostras_renda[amostra] <- mean( tabela$renda[ ind_amostra] ) medias_amostras_distancia[amostra] <- mean( tabela$distancia[ ind_amostra] ) } ###RENDA #calculando a média das médias amostrais de renda media_medias_renda <- mean(medias_amostras_renda) #histograma das medias amostrais hist(medias_amostras_renda, xlim = c(500, 1700)) #adicionando linha com a média das médias amostrais de renda abline(v = media_medias_renda, col = 'red') #adicionando linha com a média verdadeira da renda abline(v = media_verdadeira_renda, col = 'blue') ###GLICEMIA #média verdadeira da glicemia media_verdadeira_glicemia <- mean(tabela$glicemia) #calculando a média das médias amostrais de renda media_medias_glicemia <- mean(medias_amostras_glicemia) #histograma das medias amostrais hist(medias_amostras_glicemia) #adicionando linha com a média das médias amostrais de renda abline(v = media_medias_glicemia, col = 'red') #adicionando linha com a média verdadeira da renda abline(v = media_verdadeira_glicemia, col = 'blue') #DISTANCIA hist(medias_amostras_distancia) ##amostragem e intervalo de confiança amostra <- sample(x = tabela$renda, size = tamanho_amostras) #estimativa do intervalo de confiança t.test(x = amostra) media_verdadeira_renda