n <- 8
p <- 1/2
probabilidade8 <- dbinom(3, size = n, prob = p)
print(paste("Probabilidade de 3 caras: ", probabilidade8))
## [1] "Probabilidade de 3 caras: 0.21875"
n <- 10
p <- 1/6
probabilidade5 <- dbinom(2, size = n, prob = p)
print(paste("Probabilidade de 2 vezes o número 5: ", probabilidade5))
## [1] "Probabilidade de 2 vezes o número 5: 0.290710049201722"
n <- 15
p <- 9/10
probabilidade12 <- dbinom(12, size = n, prob = p)
print(paste("Probabilidade de 12 serem de boa qualidade: ", probabilidade12))
## [1] "Probabilidade de 12 serem de boa qualidade: 0.128505439098855"
n <- 20
p <- 1/4
probabilidade15 <- sum(dbinom(15:20, size = n, prob = p))
print(paste("Probabilidade de acertar pelo menos 15 perguntas: ", probabilidade15))
## [1] "Probabilidade de acertar pelo menos 15 perguntas: 3.81302743335255e-06"
probabilidade2 <- dhyper(2, m = 8, n = 5, k = 3)
print(paste("Probabilidade de 2 serem vermelhas: ", probabilidade2))
## [1] "Probabilidade de 2 serem vermelhas: 0.489510489510489"
print(paste("Em percentual temos: ", round(probabilidade2 * 100, 2), "%"))
## [1] "Em percentual temos: 48.95 %"
n <- 5
p <- 1/4
probabilidade3 <- dbinom(3, size = n, prob = p)
print(paste("Probabilidade de acertar exatamente 3 questões: ", probabilidade3))
## [1] "Probabilidade de acertar exatamente 3 questões: 0.087890625"
n <- 6
p <- 0.4
probabilidade2 <- dbinom(2, size = n, prob = p)
print(paste("Probabilidade de 2 ímpares: ", probabilidade2))
## [1] "Probabilidade de 2 ímpares: 0.31104"
n <- 20
p <- 0.3
probabilidade6 <- dbinom(6, size = n, prob = p)
print(paste("Probabilidade de 6 sucessos: ", probabilidade6))
## [1] "Probabilidade de 6 sucessos: 0.191638982753443"
n <- 3
p <- 4/12
probabilidade2oumais <- dbinom(2, size = n, prob = p) + dbinom(3, size = n, prob = p)
print(paste("Probabilidade de pelo menos 2 serem defeituosas: ", probabilidade2oumais))
## [1] "Probabilidade de pelo menos 2 serem defeituosas: 0.259259259259259"
print(paste("Em percentual temos: ", round(probabilidade2oumais * 100, 2), "%"))
## [1] "Em percentual temos: 25.93 %"
n <- 5
p <- 0.9
probabilidade_4oumais <- dbinom(4, size = n, prob = p) + dbinom(5, size = n, prob = p)
print(paste("Probabilidade de pelo menos 4 delas funcionarem: ", probabilidade_4oumais))
## [1] "Probabilidade de pelo menos 4 delas funcionarem: 0.91854"
print(paste("Em percentual temos: ", round(probabilidade_4oumais * 100, 2), "%"))
## [1] "Em percentual temos: 91.85 %"
lambda <- 2
probabilidade_3_defeitos <- dpois(3, lambda)
print(paste("A probabilidade de haver exatamente 3 defeitos em um lote: ", probabilidade_3_defeitos))
## [1] "A probabilidade de haver exatamente 3 defeitos em um lote: 0.180447044315484"
print(paste("Em percentual temos: ", round(probabilidade_3_defeitos*100,2), "%"))
## [1] "Em percentual temos: 18.04 %"
lambda <- 4
probabilidade_6_rem <- 1 - ppois(5, lambda)
print(paste("A probabilidade de haver pelo menos 6 reclamações: ", probabilidade_6_rem))
## [1] "A probabilidade de haver pelo menos 6 reclamações: 0.214869612969595"
print(paste("Em percentual temos: ", round(probabilidade_6_rem*100,2), "%"))
## [1] "Em percentual temos: 21.49 %"
lambda <- 8
probabilidade_10_cli <- 1 - ppois(9, lambda)
print(paste("A probabilidade de haver pelo menos 10 clientes: ", probabilidade_10_cli))
## [1] "A probabilidade de haver pelo menos 10 clientes: 0.283375741272989"
print(paste("Em percentual temos: ", round(probabilidade_10_cli*100,2), "%"))
## [1] "Em percentual temos: 28.34 %"
lambda <- 0.5
probabilidade_1_dis <- dpois(1, lambda)
print(paste("A probabilidade de haver exatamente 1 disparo: ", probabilidade_1_dis))
## [1] "A probabilidade de haver exatamente 1 disparo: 0.303265329856317"
print(paste("Em percentual temos: ", round(probabilidade_1_dis*100,2), "%"))
## [1] "Em percentual temos: 30.33 %"
lambda <- 3
probabilidade_5_at <- 1- ppois(4, lambda)
print(paste("A probabilidade de haver pelo menos 5 atrasos: ", probabilidade_5_at))
## [1] "A probabilidade de haver pelo menos 5 atrasos: 0.184736755476228"
print(paste("Em percentual temos: ", round(probabilidade_5_at*100,2), "%"))
## [1] "Em percentual temos: 18.47 %"
lambda <- 12
probabilidade_10_ped <- dpois(10, lambda)
print(paste("A probabilidade de receber exatamente 10 pedidos: ", probabilidade_10_ped))
## [1] "A probabilidade de receber exatamente 10 pedidos: 0.104837255883659"
print(paste("Em percentual temos: ", round(probabilidade_10_ped*100,2), "%"))
## [1] "Em percentual temos: 10.48 %"
lambda <- 1.5
probabilidade_2_ent <- dpois(2, lambda)
print(paste("A probabilidade de realizar exatamente 2 entregas: ", probabilidade_2_ent))
## [1] "A probabilidade de realizar exatamente 2 entregas: 0.251021430166984"
print(paste("Em percentual temos: ", round(probabilidade_2_ent*100,2), "%"))
## [1] "Em percentual temos: 25.1 %"
lambda <- 5
probabilidade_8_def <- 1 - ppois(7, lambda)
print(paste("A probabilidade de haver pelo menos 8 carros: ", probabilidade_8_def))
## [1] "A probabilidade de haver pelo menos 8 carros: 0.133371674070007"
print(paste("Em percentual temos: ", round(probabilidade_8_def*100,2), "%"))
## [1] "Em percentual temos: 13.34 %"
lambda <- 0.2
probabilidade_0_ev <- dpois(0, lambda)
print(paste("A probabilidade de não ocorrer nenhum evento: ", probabilidade_0_ev))
## [1] "A probabilidade de não ocorrer nenhum evento: 0.818730753077982"
print(paste("Em percentual temos: ", round(probabilidade_0_ev*100,2), "%"))
## [1] "Em percentual temos: 81.87 %"
lambda <- 7
probabilidade_6_nas <- dpois(6, lambda)
print(paste("A probabilidade de haver exatamente 6 nascimentos: ", probabilidade_6_nas))
## [1] "A probabilidade de haver exatamente 6 nascimentos: 0.149002779674338"
print(paste("Em percentual temos: ", round(probabilidade_6_nas*100,2), "%"))
## [1] "Em percentual temos: 14.9 %"
.
mu <- 170
sigma <- 10
altura <- 185
probabilidade_maior_185 <- pnorm(altura, mean = mu, sd = sigma, lower.tail = FALSE)
# Exibindo resultado
print(paste("Probabilidade de uma pessoa aleatória ter altura superior a 185 cm: ",
probabilidade_maior_185))
## [1] "Probabilidade de uma pessoa aleatória ter altura superior a 185 cm: 0.0668072012688581"
print(paste("Em percentual:", round(probabilidade_maior_185 * 100, 2), "%"))
## [1] "Em percentual: 6.68 %"
mu <- 800
sigma <- 50
dias <- 750
probabilidade_750 <- pnorm(dias, mean = mu, sd = sigma, lower.tail = FALSE)
# Exibindo resultado
print(paste("Probabilidade de uma bateria durar pelo menos 750 dias: ",
probabilidade_750))
## [1] "Probabilidade de uma bateria durar pelo menos 750 dias: 0.841344746068543"
print(paste("Em percentual:", round(probabilidade_750 * 100, 2), "%"))
## [1] "Em percentual: 84.13 %"
mu <- 100
sigma <- 15
ponto <- 120
probabilidade_sup_120 <- pnorm(ponto, mean = mu, sd = sigma, lower.tail = FALSE)
# Exibindo resultado
print(paste("Probabilidade de um aluno ter uma pontuação superior a 120: ",
probabilidade_sup_120))
## [1] "Probabilidade de um aluno ter uma pontuação superior a 120: 0.0912112197258678"
print(paste("Em percentual:", round(probabilidade_sup_120 * 100, 2), "%"))
## [1] "Em percentual: 9.12 %"
mu <- 5
sigma <- 0.5
peso <- 4.2
probabilidade_saco <- pnorm(peso, mean = mu, sd = sigma )
# Exibindo resultado
print(paste("Probabilidade de um saco ter peso inferior a 4,2 kg: ",
probabilidade_saco))
## [1] "Probabilidade de um saco ter peso inferior a 4,2 kg: 0.054799291699558"
print(paste("Em percentual:", round(probabilidade_saco * 100, 2), "%"))
## [1] "Em percentual: 5.48 %"
mu <- 25
sigma <- 3
dia <- 30
probabilidade_temp <- pnorm(dia, mean = mu, sd = sigma, lower.tail = FALSE)
# Exibindo resultado
print(paste("Probabilidade de um dia ter temperatura superior a 30°C: ",
probabilidade_temp))
## [1] "Probabilidade de um dia ter temperatura superior a 30°C: 0.0477903522728147"
print(paste("Em percentual:", round(probabilidade_temp * 100, 2), "%"))
## [1] "Em percentual: 4.78 %"
mu <- 50
sigma <- 8
vel <- 40
probabilidade_conex <- pnorm(vel, mean = mu, sd = sigma)
# Exibindo resultado
print(paste("Probabilidade de probabilidade de uma conexão ter velocidade inferior a 40 Mbps: ",
probabilidade_conex))
## [1] "Probabilidade de probabilidade de uma conexão ter velocidade inferior a 40 Mbps: 0.105649773666855"
print(paste("Em percentual:", round(probabilidade_conex * 100, 2), "%"))
## [1] "Em percentual: 10.56 %"
media <- 70
desvio <- 10
prob_ate_80 <- pnorm(80, mean = media, sd = desvio)
prob_ate_60 <- pnorm(60, mean = media, sd = desvio)
prob_entre <- prob_ate_80 - prob_ate_60
print(paste("A probabilidade de um aluno tirar nota entre 60 e 80 é:", round(prob_entre, 4)))
## [1] "A probabilidade de um aluno tirar nota entre 60 e 80 é: 0.6827"
print(paste("Em percentual:", round(prob_entre * 100, 2), "%"))
## [1] "Em percentual: 68.27 %"
mu <- 2000
sigma <- 300
cal <- 2500
probabilidade_consumo <- pnorm(cal, mean = mu, sd = sigma, lower.tail = FALSE)
# Exibindo resultado
print(paste("A probabilidade de uma pessoa consumir mais de 2500 calorias por dia: ",
probabilidade_consumo))
## [1] "A probabilidade de uma pessoa consumir mais de 2500 calorias por dia: 0.0477903522728147"
print(paste("Em percentual:", round(probabilidade_consumo * 100, 2), "%"))
## [1] "Em percentual: 4.78 %"
mu <- 120
sigma <- 10
press <- 130
probabilidade_pressao <- pnorm(press, mean = mu, sd = sigma, lower.tail = FALSE)
# Exibindo resultado
print(paste("A probabilidade de uma pessoa ter pressão superior a 130 mmHg: ",
probabilidade_pressao))
## [1] "A probabilidade de uma pessoa ter pressão superior a 130 mmHg: 0.158655253931457"
print(paste("Em percentual:", round(probabilidade_pressao * 100, 2), "%"))
## [1] "Em percentual: 15.87 %"
mu <- 8
sigma <- 1
comp <- 6.5
probabilidade_elet <- pnorm(comp, mean = mu, sd = sigma)
# Exibindo resultado
print(paste("A probabilidade de um componente ter medida inferior a 6,5 cm: ",
probabilidade_elet))
## [1] "A probabilidade de um componente ter medida inferior a 6,5 cm: 0.0668072012688581"
print(paste("Em percentual:", round(probabilidade_elet * 100, 2), "%"))
## [1] "Em percentual: 6.68 %"
R: Amostragem aleatória simples
R: Amostragem Aleatória Simples
R: Amostragem Estratificada
R: Amostragem Estratificada
R: Amostragem por Conglomerados
R: Amostragem por Conglomerados
R: Amostragem Aleatória Simples
R: Amostragem Estratificada
R: Amostragem Estratificada
R: Amostragem por Conglomerados
R: Amostragem por conglomerados
n <- 100
media_amostral <- 175
desvio_padrao_populacional <- 8
nivel_confianca <- 0.95
# Cálculo do intervalo de confiança
erro_padrao <- desvio_padrao_populacional / sqrt(n)
margem_erro <- qnorm((1 + nivel_confianca) / 2) * erro_padrao
intervalo_confianca <- c(media_amostral - margem_erro,
media_amostral + margem_erro)
print("Intervalo de confiança: ")
## [1] "Intervalo de confiança: "
print(intervalo_confianca)
## [1] 173.432 176.568
n <- 200
satisfeitos <- 150
proporcao_amostral <- satisfeitos / n
nivel_confianca <- 0.90
erro_padrao <- sqrt(proporcao_amostral * (1 - proporcao_amostral) / n)
margem_erro <- qnorm((1 + nivel_confianca) / 2) * erro_padrao
intervalo_confianca <- c(proporcao_amostral - margem_erro, proporcao_amostral + margem_erro)
print("Intervalo de confiança:")
## [1] "Intervalo de confiança:"
print(intervalo_confianca)
## [1] 0.6996368 0.8003632
n <- 30
media_amostral <- 50
desvio_padrao_amostral <- 6
nivel_confianca <- 0.99
erro_padrao <- desvio_padrao_amostral / sqrt(n)
margem_erro <- qt((1 + nivel_confianca) / 2, df = n - 1) * erro_padrao
intervalo_confianca_maca <- c(media_amostral - margem_erro, media_amostral + margem_erro)
print("Intervalo de confiança:")
## [1] "Intervalo de confiança:"
print(intervalo_confianca_maca)
## [1] 46.98053 53.01947
n <- 50
media_amostral <- 1200
desvio_padrao_populacional <- 100
nivel_confianca <- 0.95
# Cálculo do intervalo de confiança
erro_padrao <- desvio_padrao_populacional / sqrt(n)
margem_erro <- qnorm((1 + nivel_confianca) / 2) * erro_padrao
intervalo_confianca_lamp <- c(media_amostral - margem_erro,
media_amostral + margem_erro)
print("Intervalo de confinaça: ")
## [1] "Intervalo de confinaça: "
print(intervalo_confianca_lamp)
## [1] 1172.282 1227.718
n <- 500
media_amostral <- 0.04
desvio_padrao_populacional <- 0.001 # 0.1% em decimal
nivel_confianca <- 0.98
erro_padrao <- desvio_padrao_populacional / sqrt(n)
margem_erro <- qnorm((1 + nivel_confianca) / 2) * erro_padrao
intervalo_confianca_epi <- c(media_amostral - margem_erro, media_amostral + margem_erro)
print("Intervalo de confiança:")
## [1] "Intervalo de confiança:"
print(intervalo_confianca_epi)
## [1] 0.03989596 0.04010404
n_46 <- 20
media_amostral_46 <- 25
desvio_padrao_amostral_46 <- 3
nivel_confianca_46 <- 0.90
erro_padrao_46 <- desvio_padrao_amostral_46 / sqrt(n_46)
margem_erro_46 <- qt((1 + nivel_confianca_46) / 2, df = n_46 - 1) * erro_padrao_46
intervalo_confianca_46 <- c(media_amostral_46 - margem_erro_46,
media_amostral_46 + margem_erro_46)
print("Intervalo de confiança :")
## [1] "Intervalo de confiança :"
print(intervalo_confianca_46)
## [1] 23.84006 26.15994
n <- 100
media_amostral <- 2000
desvio_padrao_populacional <- 300
nivel_confianca <- 0.95
# Cálculo do intervalo de confiança
erro_padrao <- desvio_padrao_populacional / sqrt(n)
margem_erro <- qnorm((1 + nivel_confianca) / 2) * erro_padrao
intervalo_confianca <- c(media_amostral - margem_erro,
media_amostral + margem_erro)
print(intervalo_confianca)
## [1] 1941.201 2058.799
n <- 50
media_amostral <- 500
desvio_padrao_populacional <- 100
nivel_confianca <- 0.99
erro_padrao <- desvio_padrao_populacional / sqrt(n)
margem_erro <- qnorm((1 + nivel_confianca) / 2) * erro_padrao
intervalo_confianca_sal <- c(media_amostral - margem_erro, media_amostral + margem_erro)
print("Intervalo de confiança:")
## [1] "Intervalo de confiança:"
print(intervalo_confianca_sal)
## [1] 463.5723 536.4277
n_49 <- 25
media_amostral_49 <- 12
desvio_padrao_amostral_49 <- 2
nivel_confianca_49 <- 0.96
# Cálculo do intervalo de confiança usando a distribuição t
erro_padrao_49 <- desvio_padrao_amostral_49 / sqrt(n_49)
margem_erro_49 <- qt((1 + nivel_confianca_49) / 2, df = n_49 - 1) * erro_padrao_49
intervalo_confianca_49 <- c(media_amostral_49 - margem_erro_49,
media_amostral_49 + margem_erro_49)
print("Intervalo de confiança:")
## [1] "Intervalo de confiança:"
print(intervalo_confianca_49)
## [1] 11.13138 12.86862
n_50 <- 150
satisfeitos_50 <- 120
proporcao_amostral_50 <- satisfeitos_50 / n_50
nivel_confianca_50 <- 0.99
# Cálculo do intervalo de confiança para proporção
erro_padrao_50 <- sqrt(proporcao_amostral_50 * (1 - proporcao_amostral_50) / n_50)
margem_erro_50 <- qnorm((1 + nivel_confianca_50) / 2) * erro_padrao_50
intervalo_confianca_50 <- c(proporcao_amostral_50 - margem_erro_50,
proporcao_amostral_50 + margem_erro_50)
print("Intervalo de confiança:")
## [1] "Intervalo de confiança:"
print(intervalo_confianca_50)
## [1] 0.7158738 0.8841262
z <- qnorm(0.975)
erro <- 100
desvio <- 500
n_51 <- (z * desvio / erro)^2
ceiling(n_51)
## [1] 97
z <- qnorm(0.95)
erro <- 0.05
p <- 0.5
N <- 5000
n_bruto <- (z^2 * p * (1 - p)) / erro^2
n_52 <- n_bruto / (1 + ((n_bruto - 1) / N))
ceiling(n_52)
## [1] 257
z <- qnorm(0.975)
erro <- 0.03
p <- 0.5
N <- 8000
n_bruto <- (z^2 * p * (1 - p)) / erro^2
n_54 <- n_bruto / (1 + ((n_bruto - 1) / N))
ceiling(n_54)
## [1] 942
z <- qnorm(0.95)
erro <- 1
# desvio de 6 meses = 0.5 ano
desvio <- 0.5
n_55 <- (z * desvio / erro)^2
ceiling(n_55)
## [1] 1
z <- qnorm(0.98)
erro <- 0.02
p <- 0.5
N <- 10000
n_bruto <- (z^2 * p * (1 - p)) / erro^2
n_56 <- n_bruto / (1 + ((n_bruto - 1) / N))
ceiling(n_56)
## [1] 2087
z <- qnorm(0.99)
erro <- 10
desvio <- 50
n_57 <- (z * desvio / erro)^2
ceiling(n_57)
## [1] 136
z <- qnorm(0.97)
erro <- 0.04
p <- 0.5
N <- 6000
n_bruto <- (z^2 * p * (1 - p)) / erro^2
n_58 <- n_bruto / (1 + ((n_bruto - 1) / N))
ceiling(n_58)
## [1] 507
z <- qnorm(0.995)
erro <- 5
desvio <- 10
n_59 <- (z * desvio / erro)^2
ceiling(n_59)
## [1] 27
z <- qnorm(0.96)
erro <- 0.02
p <- 0.5
N <- 12000
n_bruto <- (z^2 * p * (1 - p)) / erro^2
n_60 <- n_bruto / (1 + ((n_bruto - 1) / N))
ceiling(n_60)
## [1] 1652