Relatório de Teste de Hipóteses

Relatório para a disciplina MAF 160 (Teste de Hipóteses)

1 -Trace uma curva normal e sobreie a área desejada obtendo então a informação.

Área à direita de $z=1$

pnormGC(1, region = "above", mean = 0, sd=1, graph = T)

## [1] 0.1586553

Área à esquerda de $Z=1$

pnormGC(1, region = "below", mean = 0, sd=1, graph = T)

## [1] 0.8413447

Área entre $Z=0$ e $Z=1.5$)

pnormGC(c(0,1.5), region = "between", mean = 0, sd=1, graph = T)

## [1] 0.4331928

Área entre $Z=-0.56$ e $Z=-0.2$

pnormGC(c(-0.56, -0.2), region = "between", mean = 0, sd=1, graph = T)

## [1] 0.1330006

e)Área entre $Z= 0.5$ e $Z= 0,5$

pnormGC(c(0.5, 0.5), region = "between", mean = 0, sd=1, graph = T)

## [1] 0

Área entre $Z= 0$ e $Z = −2.5$

pnormGC(c(0, -2.5), region = "between", mean = 0, sd=1, graph = T)

## [1] 0.4937903

1. Área de $0.0505$ à esquerda de Z

qnorm(0.5005, mean = 0, sd=1, lower.tail = T)

## [1] 0.001253314

pnormGC(qnorm(0.0505, mean = 0, sd=1, lower.tail = T), region = "below", mean = 0, sd=1, graph = T)

## [1] 0.0505

Área de $0.0228$ à direita de Z

qnorm(0.0228, mean = 0, sd=1, lower.tail = T)

## [1] -1.999077

pnormGC(qnorm(0.0228, mean = 0, sd=1, lower.tail = T), region = "above", mean = 0, sd=1, graph = T)

## [1] 0.9772

Área de $0.0228$ à esquerda de Z

qnorm(0.0228)

## [1] -1.999077

pnormGC(qnorm(0.0228, mean = 0, sd=1), region="below", mean=0,
        sd=1, graph=TRUE)

## [1] 0.0228

$0.4772$ entre 0 e z.

qnorm(0.5+0.4772)

## [1] 1.999077

3- Consultando a tabela, determine a probabilidade de certo valor padronizado de Z estar entre Z0 = −1, 20 e Z1 = 2, 00. Desenhe o gráfico.

pnorm(q=2.00, mean = 0, sd=1)-pnorm(q=-1.20, mean = 0, sd=1)

## [1] 0.8621802

pnormGC(c(-1.20,2.00), region="between", mean=0,
        sd=1, graph=TRUE)

## [1] 0.8621802

4- Dado uma variável X com distribuição normal de média 25 e desvio-padrão 2, determine os valores de Z para os seguintes valores (x) : Sabemos que Z=X−μσ, logo:

mu=25
sigma=2
x=23
Z=(x-mu)/sigma
Z

## [1] -1

x<- 23.5
mu<- 25
sigma<- 2

Z<-(x-mu)/sigma
Z

## [1] -0.75

x=24
Z=(x-mu)/sigma
Z

## [1] -0.5

x=25.2
Z=(x-mu)/sigma
Z

## [1] 0.1

x= 25.5
Z=(x-mu)/sigma
Z

## [1] 0.25

5- Determine a probabilidade de certo valor padronizado de Z estar entre Z0 = $−1.30$ e Z1 = $1.5$. Desenhe o gráfico $−1.30$ e Z1 = $1.5$

pnorm(q=1.5, mean = 0, sd=1)-pnorm(q=-1.30, mean = 0, sd=1)

## [1] 0.8363923

pnormGC(c(-1.30,1.5), region="between", mean=0,
        sd=1, graph=TRUE)

## [1] 0.8363923

6- Uma população normal tem média 40 e desvio-padrão 3. Determine os valores da população correspondentes aos seguintes de Z: Sabemos que Z=X−μσ, logo, dado Z, temos que X=Zσ+μ, assim:

0,10

mu=40
sigma=3
Z=0.1
x=mu+sigma*Z
x

## [1] 40.3

2,00

Z=2
x=mu+sigma*Z
x

## [1] 46

0,75

Z=0.75
x=mu+sigma*Z
x

## [1] 42.25

-3,00

Z=-3
x=mu+sigma*Z
x

## [1] 31

-2.53

Z=-2.53
x=mu+sigma*Z
x

## [1] 32.41

7- a) Teste de hipótese; É a estátitica que testa duas hipóteses, permitindo tomar decisões de qual são corretas.

b)Hipótese nula e alternativa;

Nula: afirma um valor hipótetico Alternativa: afirma que mu é diferente do valor da nula

Erros do tipo I e II;

Tipo 1, reita-se a H0, quando verdadeira

Tipo 2, H0 é falsa, e não rejeitamos

nível de signicância; chamado de apha é número que irá determinar o limite para um resultado ser sginificativo ou não

8- Enuncie a hipótese nula e a hipótese alternativa em cada um dos casos a seguir.

Enuncie a hipótese nula e a hipótese alternativa em cada um dos casos a seguir.

H0:μ=40 toneladas por hectare

H1:μ>40 toneladas por hectare

#H0:mu=40
#H1:mu>40

Um sindicato de empregados de certa categoria deseja verificar se a taxa de desemprego em certo município é maior que a taxa de 12% observada seis meses antes.

#H0:mu=12
#H1:mu>12

9- O fabricante de certa marca de suco informa que as embalagens de seu produto têm em média 500 ml, com desvio padrão igual a 10 ml. Tendo sido encontradas no mercado algumas embalagens com menos de 500 ml, suspeita-se que a informação do fabricante seja falsa. Para verificar se isto ocorre, um fiscal analisa uma amostra de 200 embalagens escolhidas aleatoriamente no mercado e constata que as mesmas contêm em média 498 ml. Considerando-se um nível de significância de 5%, pode-se afirmar que o fabricante está mentindo? Calcule o valor da prova para esta amostra.

mu=500
sigma=10
n=200
alpha=0.05
xbarra=498

# unilateral a esquerda
# utilizar o Teste z
# H0:mu=500ml
# H1:mu<500ml
# Estatística do teste

Zcal <- (xbarra-mu)/(sigma/sqrt(n))
Zcal

## [1] -2.828427

Ztab <- qnorm(0.05)
Ztab

## [1] -1.644854

ConclusaoZ <- 
ifelse(abs(Zcal)>abs(Ztab), paste("Como |Zcal|>|Ztab| Rejeita-se H0 ao nível de", alpha,"de significância"), paste
("Como |Zcal|<|Ztab| Não rejeita-se H0 ao nível de", alpha,"de significância"))
ConclusaoZ

## [1] "Como |Zcal|>|Ztab| Rejeita-se H0 ao nível de 0.05 de significância"

pnormGC(Zcal, region = "below", mean = 0, sd=1, graph = T)

## [1] 0.002338867

# Nesse caso o fabricante está mentindo! O real valor da prova deve ser:
(x=mu+(sigma/sqrt(n))*Ztab)

## [1] 498.8369

10- A duração das lâmpadas produzidas por certo fabricante tem distribuição normal com média igual a 1200 horas e desvio padrão igual a 300 horas. O fabricante introduz um novo processo na produção das lâmpadas.Para verificar se o novo processo produz lâmpadas de maior duração, o fabricante observa 100 lâmpadas produzidas pelo novo processo e constata que as mesmas duram em média 1265 horas. Admitindo-se um nível de significância de 5%, pode-se concluir que o novo processo produz lâmpadas com maior duração?

n <- 100
xbarra <- 1265
mu <- 1200
sigma <- 300
alpha <- 0.05
#H0:mu = 1200h
#H1:mu > 1200h (Teste-z unilateral)
#Estatística do Teste
Zcal <- (xbarra-mu)/(sigma/sqrt(n))
Zcal

## [1] 2.166667

pnormGC(Zcal, region="above", mean=0,
        sd=1, graph=TRUE)

## [1] 0.01513014

RRH0 <- ("Como |Zcal|>|Ztab| Rejeita-se H0 ao nível alpha=0.05 de significância")
RNH0 <- ("Como |Zcal|<|Ztab| Não rejeita-se H0 ao nível alpha=0.05 de significância")

(Conclusão <- ifelse(abs(Zcal)>abs(Ztab),RRH0, RNH0))

## [1] "Como |Zcal|>|Ztab| Rejeita-se H0 ao nível alpha=0.05 de significância"

#Logo, o fabricante está mentindo!
(pvalor <- pnorm(Zcal, lower.tail = FALSE))

## [1] 0.01513014

11- O custo de produção de certo artigo numa localidade tem distribuição normal com média igual a R$42,00.Desenvolve-se uma política de redução de custos na empresa para melhorar a competitividade do referido produto no mercado. Observando-se os custos de 10 unidades deste produto, obtiveram-se os seguintes valores: 34, 41, 36, 41, 29, 32, 38, 35, 33 e 30. Admitindo-se um nível de significância de 5%, pode-se afirmar que o custo do produto considerado diminuiu?

mu <- 42
n <- 10
x <- c(34, 41, 36, 41, 29, 32, 38, 35, 33, 30)
xbarra <- mean(x)
s <- sd(x)
alpha <- 0.05
#H0: mu = 42 
#H1: mu < 42 (Teste-t unilateral)
gl <- n-1
#Estatística do Teste
Tcal <- (xbarra-mu)/(s/sqrt(n))
Tcal

## [1] -5.377348

Ttab <- qt(alpha, df=gl, lower.tail = TRUE,)
Ttab

## [1] -1.833113

ConclusaoT <- ifelse(abs(Tcal)>abs(Ttab),paste(
"Como |Tcal|>|Ttab| Rejeita-se H0 ao nível de"
, alpha ,
"de significância"
), paste(
"Como |Tcal|<|Ttab| Não Rejeita-se H0 ao nível de"
, alpha ,
"de significância"
))
ConclusaoT

## [1] "Como |Tcal|>|Ttab| Rejeita-se H0 ao nível de 0.05 de significância"

#Logo, o fabricante está mentindo!
(pvalor <- pt(Tcal, df = gl, lower.tail = TRUE))

## [1] 0.0002230215

ConclusaoT <- ifelse(pvalor>alpha,paste(
"Como p-valor>", alpha, " Não Rejeita-se H0"
), paste(
"Como p-valor<", alpha, " Rejeita-se H0")
)
ConclusaoT

## [1] "Como p-valor< 0.05  Rejeita-se H0"

12- O controle de qualidade das peças produzidas por certa fábrica exige que o diâmetro médio das mesmas seja 57 mm. Para verificar se o processo de produção está sob controle, observam-se os diâmetros de 10 peças, constatando-se os seguintes valores em mm: 56,5; 56,6; 57,3; 56,9; 57,1; 56,7; 57,1; 56,8; 57,1; 57,0. Admitindo-se um nível de significância de 5%, pode-se concluir que o processo de produção está sob controle?

mu <- 57
n <- 10
x <- c(56.5, 56.6, 57.3, 56.9, 57.1, 56.7, 57.1, 56.8, 57.1, 57.0)
xbarra <- mean(x)
s <- sd(x)
alpha <- 0.05
#H0: mu = 57 
#H1: mu != 57 (Teste-t unilateral)
gl <- n-1
#Estatística do Teste
Tcal <- (xbarra-mu)/(s/sqrt(n))
Tcal

## [1] -1.112516

Ttab <- qt(alpha, df=gl, lower.tail = TRUE,)
Ttab

## [1] -1.833113

ConclusaoT <- ifelse(abs(Tcal)>abs(Ttab),paste(
"Como |Tcal|>|Ttab| Rejeita-se H0 ao nível de"
, alpha ,
"de significância"
), paste(
"Como |Tcal|<|Ttab| Não Rejeita-se H0 ao nível de"
, alpha ,
"de significância"
))
ConclusaoT

## [1] "Como |Tcal|<|Ttab| Não Rejeita-se H0 ao nível de 0.05 de significância"

#o processo de produção está sob controle
(pvalor <- pt(Tcal, df = gl, lower.tail = TRUE))

## [1] 0.1473741

ConclusaoT <- ifelse(pvalor>alpha,paste(
"Como p-valor>", alpha, " Não Rejeita-se H0"
), paste(
"Como p-valor<", alpha, " Rejeita-se H0")
)
ConclusaoT

## [1] "Como p-valor> 0.05  Não Rejeita-se H0"

15- Suponha que o tempo necessário para que estudantes completem uma prova tenha distribuição normal com média 90 minutos e desvio padrão 15 minutos.

a)Qual é a probabilidade do estudante terminar a prova em menos de 80 minutos?

(pMenor80 <- pnorm(80, mean = 90, sd = 15, lower.tail = TRUE))

## [1] 0.2524925

b)Em mais de 120 minutos?

(pMaior120 <- pnorm(120, mean = 90, sd = 15, lower.tail = FALSE))

## [1] 0.02275013

c)Entre 75 e 85 minutos?

pnorm(q=85, mean = 90, sd=15)-75

## [1] -74.63056

pnormGC(c(75,85), region="between", mean=90,
        sd=15, graph=TRUE)

## [1] 0.2107861

d)Qual é o tempo necessário para que 98% dos estudantes terminem a prova?

qnorm(0.98, mean = 90, sd = 15)

## [1] 120.8062

16- Uma v.a. X tem distribuição normal, com média 100 e desvio padrão 10

Qual a P(90 < X < 110)?

mu=100
sigma=10

#P=P(90<X<110)=P(-1>Z<1)=2*P(0<Z<1)
P <- 2*(pnorm(1)-pnorm(0))
P

## [1] 0.6826895

b)Se X¯ for a média de uma amostra de 16 elementos retirados dessa população, calcule P(90<X¯<110).

#P=P(90<Xbarra<110=P(-4<Z<4)=2*P(0<Z<4)
P <- 2*(pnorm(4)-pnorm(0))
P

## [1] 0.9999367

c)Represente, num único gráfico, as distribuições de X e Xbarra.

n <- 16
x <- seq(-50, 150, length=1000)
mu <- 100
sigma <- 10/sqrt(n)
z <- dnorm(x, mean = 100, sd = 10)
colors <- "blue"
plot(x, z, type="l", lty=2, xlab="x", ylim = c(0,0.2),
    ylab="Densidade", main="Comparação de distribuições normais")
lines(x, dnorm(x,mu,sigma), lwd=2, col=colors)

d)Que tamanho deveria ter a amostra para que P(90 < Xbarra < 110) = 0, 95?

qnorm(0.98, mean = 90, sd = 15)

## [1] 120.8062

17- Nas situações abaixo, escolha como hipótese nula, H0, aquela que para você leva a um erro tipo I mais importante. Descreva quais os dois erros em cada caso:

O trabalho de um operador de radar é detectar aeronaves inimigas. Quando surge alguma coisa estranha na tela, ele deve decidir entre as hipotéses:

está começando um ataque; R- Hipótese nula, pois há um maior indício de que ocorra um erro do tipo 1. Então rejeitariamos a hipótese de que o país está sendo atacado quando na verdade concluimos que está sendo.

tudo bem, apenas uma leve interferência. R-Hipótese alternativa.

b)Num júri, um indivíduo está sendo julgado por um crime. As hipóteses sujeitas ao júri são: 1. o acusado é inocente; R- Hipótese alternativa. 2. o acusado é culpado. R-H0:o acusado é inocente

Um pesquisador acredita que descobriu uma vacina contra resfriado. Ele irá conduzir uma pesquisa de laboratório para verificar a veracidade da afirmação. De acordo com o resultado, ele lançará ou não a vacina no mercado. As hipóteses que pode testar são:
a vacina é eficaz;

a vacina não é eficaz. H0:a vacina é eficaz

Quando se rejeita H0 neste caso você fala que a vacina não é eficaz, assim uma vacina que poderia salvar vidas será desconsiderada, o que acho pior do considerar uma vacina não eficaz em eficaz.

18- Uma fábrica de automóveis anuncia que seus carros consomem, em média, 11 litros por 100 km, com desvio padrão de 0,8 litros. Uma revista resolve testar essa afirmação e analisa 35 automóveis dessa marca, obtendo 11,3 litros por 100 km como consumo médio (considerar distribução normal). O que a revista pode concluir sobre o anúncio da fábrica, no nível de 10%?

mu=11
sigma=0.8
n=35
alpha=0.10
xbarra=11.3

# Vamos utilizar o Teste - z
# H0:mu=11
# H1:mu>11
# Valor calculado do teste
Zcal <- (xbarra-mu)/(sigma/sqrt(n))
Zcal

## [1] 2.21853

Ztab <- qnorm(0.10)
Ztab

## [1] -1.281552

ConclusaoZ <- 
ifelse(abs(Zcal)>abs(Ztab), paste("Como |Zcal|>|Ztab| Rejeita-se H0 ao nível de", alpha,"de significância"), paste
("Como |Zcal|<|Ztab| Não rejeita-se H0 ao nível de", alpha,"de significância"))
ConclusaoZ

## [1] "Como |Zcal|>|Ztab| Rejeita-se H0 ao nível de 0.1 de significância"

pnormGC(c(qnorm(0.05), qnorm(0.05, lower.tail = F)), region="between", mean=0, sd=1, graph=TRUE)

## [1] 0.9

ConclusaoZ <- ifelse(pvalor>alpha,paste(
"Como p-valor>", alpha, " Não Rejeita-se H0"
), paste(
"Como p-valor<", alpha, " Rejeita-se H0")
)
ConclusaoZ

## [1] "Como p-valor> 0.1  Não Rejeita-se H0"

19 - Duas máquinas, A e B, são usadas para empacotar pó de café. A experiência passada garante que o desvio padrão para ambas é de 10 g. Porém, suspeita-se que elas têm médias diferentes. Para verificar, sortearam-se duas amostras: uma com 25 pacotes da máquina A e outra com 16 pacotes da máquina B. As médias foram, respectivamente, xbarraA = 502, 74g e xbarraB = 496, 60g. Com esses números, e com o nível de 5%, qual seria a conclusão do teste H0 : µA = µB?

xAbarra <- 502.74
xBbarra <- 496.6
alpha <- 0.05
nA <- 25
nB <- 16
sigma <- 10

# H0:muA=muB
# H1:muA!muB
# Estatística do Teste:

Zcal <- (xAbarra-xBbarra)/sqrt(sigma^2*((1/nA)+(1/nB))) 
Zcal

## [1] 1.917814

Ztab <- qnorm(alpha/2)
Ztab

## [1] -1.959964

ConclusaoZ <- 
ifelse(abs(Zcal)>abs(Ztab), paste("Como |Zcal|>|Ztab| Rejeita-se H0 ao nível de", alpha,"de significância"), paste
("Como |Zcal|<|Ztab| Não rejeita-se H0 ao nível de", alpha,"de significância"))
ConclusaoZ

## [1] "Como |Zcal|<|Ztab| Não rejeita-se H0 ao nível de 0.05 de significância"

pnormGC(c(qnorm((alpha)/2), qnorm((alpha)/2, lower.tail = F)), region="between", mean=0, sd=1, graph=TRUE)

## [1] 0.95

ConclusaoZ <- ifelse(pvalor>alpha,paste(
"Como p-valor>", alpha, " Não Rejeita-se H0"
), paste(
"Como p-valor<", alpha, " Rejeita-se H0")
)
ConclusaoZ

## [1] "Como p-valor> 0.05  Não Rejeita-se H0"

20- Uma fábrica de embalagens para produtos químicos está estudando dois processos para combater a corrosão de suas latas especiais. Para verificar o efeito dos tratamentos, foram usadas amostras cujos resultados estão no quadro abaixo (em porcentagem de corrosão eliminada). Qual seria a conclusão sobre os dois tratamentos?

# Será utilizado o teste t
# H0>muA=muB
# H1:muA!=muB
nA=15
nB=12
xBarraA=48
xBarraB=52
sA=10
sB=15
alpha=0.05

# Será utilizado o teste-F
# H0:SigmaA^2=SigmaB^2
# H1:SigmaA^2<SigmaB^2

Fcal <- (sB^2)/(sA^2)
Fcal

## [1] 2.25

(pvalor <- pf(q=Fcal, df1 = nB-1, df2 = nA-1))

## [1] 0.9224523

#Conclusao: Não rejeita-se H0 para todo alpha < pvalor. Assim, para realizar o teste - t, vamos considerar que as variâncias do método A 
#e método B são iguais (SigmaA^2=SigmaB^2). Dessa forma, procedemos 
# o teste-t para variâncias populacionais iguais.

#H0:muA=muB
#H1:muA!=muB (Teste-t bilateral)

#A <- (sA^2)/nA
#B <- (sB^2)/nB
#
#df <- ((A+B)^2)/(((A^2)/(nA-1))+((B^2)/(nB-1)))
df <- nA+nB-2  
Sc <- sqrt((((sA^2)*(nA-1))+((sB^2)*(nB-1)))/(nA+nB-2))  
#Estatística do Teste
Tcal <- (xAbarra-xBbarra)/(Sc*sqrt((1/nA)+(1/nB)))
Tcal

## [1] 1.273377

(pvalor <- pt(Tcal, df = gl, lower.tail = TRUE))

## [1] 0.8826059

ConclusaoT <- "Não rejeita-se H0 para alpha menor que pvalor"

21- Para investigar a influência da opção profissional sobre o salário inicial de recém-formados, investigaram-se dois grupos de profissionais: um de liberais em geral e outro de formandos em Administração de Empresas. Com os resultados abaixo, expressos em salários mínimos, quais seriam suas conclusões?

Liberais 6,6 10,3 10,8 12,9 9,2 12,3 7,0 Administradores 8,1 9,8 8,7 10,0 10,2 8,2 8,7

Li <- c(6.6, 10.3, 10.8, 12.9, 9.2, 12.3, 7.0)
Ad <- c(8.1, 9.8, 8.7, 10.0, 10.2, 8.2, 8.7, 10.1)
# H0:Sigma2Li=Sigma2Ad
# H1:Sigma2Li!=Sigma2Ad
var.test(Li,Ad,alternative = "two.sided")

## 
##  F test to compare two variances
## 
## data:  Li and Ad
## F = 7.5128, num df = 6, denom df = 7, p-value = 0.01768
## alternative hypothesis: true ratio of variances is not equal to 1
## 95 percent confidence interval:
##   1.467755 42.789180
## sample estimates:
## ratio of variances 
##           7.512844

## 
##  F test to compare two variances
## 
## data:  Li and Ad
## F = 7.5128, num df = 6, denom df = 7, p-value = 0.01768
## alternative hypothesis: true ratio of variances is not equal to 1
## 95 percent confidence interval:
##   1.467755 42.789180
## sample estimates:
## ratio of variances 
##           7.512844

23- Os dados abaixo referem-se a medidas de determinada variável em 19 pessoas antes e depois de uma cirurgia.Verifique se as medidas pré e pós-operatórias apresentam a mesma média. Que suposições você faria para resolver o problema?

Pre <- c(50.0,50.0,50.0,87.5,32.5,35.0,40.0,45.0,62.5,40.0,50.0,75.0,92.5,38.0,46.5,50.0,30.0,35.0,39.4)
Pos <- c(42.0,42.0,78.0,33.0,96.0,82.0,44.0,31.0,87.0,50.0,48.0,52.0,74.0,47.5,49.0,58.0,42.0,60.0,28.0)
Dif <- Pre-Pos
sd(Dif)

## [1] 26.35174

mean(Dif)

## [1] -4.978947

#H0:muDif=0
#H1:muDif!=0

t.test(Dif, alternative = "two.sided")

## 
##  One Sample t-test
## 
## data:  Dif
## t = -0.82358, df = 18, p-value = 0.421
## alternative hypothesis: true mean is not equal to 0
## 95 percent confidence interval:
##  -17.680077   7.722183
## sample estimates:
## mean of x 
## -4.978947

## 
##  One Sample t-test
## 
## data:  Dif
## t = -0.82358, df = 18, p-value = 0.421
## alternative hypothesis: true mean is not equal to 0
## 95 percent confidence interval:
##  -17.680077   7.722183
## sample estimates:
## mean of x 
## -4.978947

print("Não rejeita-se H0 ao nível de 5% de significância, pois pvalor é maior que alpha=0.05")

## [1] "Não rejeita-se H0 ao nível de 5% de significância, pois pvalor é maior que alpha=0.05"

24- Uma empresa deseja estudar o efeito de uma pausa de dez minutos para um cafezinho sobre a produtividade de seus trabalhadores. Para isso, sorteou seis operários, e contou o número de peças produzidas durante uma semana sem intervalo e uma semana com intervalo. Os resultados sugerem se há ou não melhora na produtividade? Caso haja melhora, qual deve ser o acréscimo médio de produção para todos os trabalhadores da fábrica?

Sem <- c(23,35,29,33,43,32)
Com <- c(28,38,29,37,42,30)
Dif <- Sem-Com
sd(Dif)

## [1] 2.880972

mean(Dif)

## [1] -1.5

#H0:muDif=0
#H1:muDif!=0

t.test(Dif, alternative = "two.sided")

## 
##  One Sample t-test
## 
## data:  Dif
## t = -1.2753, df = 5, p-value = 0.2582
## alternative hypothesis: true mean is not equal to 0
## 95 percent confidence interval:
##  -4.523395  1.523395
## sample estimates:
## mean of x 
##      -1.5

print(“Não rejeita-se H0 ao nível de 5% de significância, pois pvalor é maior que alpha=0.05, portanto não houve acréscimo na produtividade.”)

25- Num levantamento feito com os operários da indústria mecânica, chegou-se aos seguintes números: salário médio = 3,64 salários mínimos e desvio padrão = 0,85 salário mínimo. Suspeita-se que os salários de subclasse formada pelos torneiros mecânicos são diferentes dos salários do conjunto todo, tanto na média como na variância. Que conclusões você obteria se uma amostra de 25 torneiros apresentasse salário médio igual a 4,22 salários mínimos e desvio padrão igual a 1,25 salário mínimo?

mu <- 3.64
sigma <- 0.85
n <- 35
xbarra <- 4.22
alpha <- 0.5
#H0:mu=3,64
#H1:mu>3,64

#Estatística do Teste
Zcal <- (xbarra-mu)/(sigma/sqrt(n))
Zcal

## [1] 4.036854

Ztab <- qnorm(alpha, mean = 0, sd=1, lower.tail = T)
Ztab

## [1] 0

ConclusaoZ <- 
ifelse(abs(Zcal)>abs(Ztab), paste("Como |Zcal|>|Ztab| Rejeita-se H0 ao nível de", alpha,"de significância"), paste
("Como |Zcal|<|Ztab| Não rejeita-se H0 ao nível de", alpha,"de significância"))
ConclusaoZ

## [1] "Como |Zcal|>|Ztab| Rejeita-se H0 ao nível de 0.5 de significância"

#Logo, conclui- se que formados em torneiro mecânico possui de fato o salário maior!
(pvalor <- pnorm(Zcal, lower.tail = FALSE))

## [1] 2.708633e-05

ConclusaoZ <- ifelse(pvalor>alpha,paste(
"Como p-valor>", alpha, " Não Rejeita-se H0"
), paste(
"Como p-valor<", alpha, " Rejeita-se H0")
)
ConclusaoZ

## [1] "Como p-valor< 0.5  Rejeita-se H0"

26- Um partido afirma que a porcentagem de votos masculinos a seu favor será de 10 % a mais do que a porcentagem de votos femininos. Numa pesquisa feita entre 400 homens, 170 votariam no partido, enquanto entre 625 mulheres, 194 lhe seriam favoráveis. A afirmação do partido é verdadeira ou não?

##teste da variância
n = 25
S = 1.25
sigma = 0.85
alpha = 0.05
#H0: Sigma = 0.85
#H1: Sigma =! 0.85
XsqCAL = ((n-1)*S^2)/sigma^2
XsqTAB = 30.144
ConclusãoXsq = ifelse(abs(XsqCAL)>abs(XsqTAB),paste("Como |XsqCAL|>|XsqTAB| Rejeita-se H0 ao nível de", alpha, "de significância"),paste("Como |XsqCAL|<|XsqTAB| Não Rejeita-se H0 ao nível de",alpha,"de significância"))
ConclusãoXsq

## [1] "Como |XsqCAL|>|XsqTAB| Rejeita-se H0 ao nível de 0.05 de significância"

## Os valores do salário da subclasse possuem variâncias diferentes.
#teste z (variância populacional conhecida)
#H0: mu = 3.64
#H1: mu =! 3.64
#Estatística do teste
mu = 3.64
sigma = 0.85
n = 25
xbarra = 4.22
s = 1.25
alpha = 0.05
zcal = (xbarra-mu)/(sigma/sqrt(n))
zcal

## [1] 3.411765

ztab = qnorm(alpha/2)
ztab

## [1] -1.959964

ConclusãoZ = ifelse(abs(zcal)>abs(ztab),paste("Como |zcal|>|ztab| Rejeita-se H0 ao nível de", alpha, "de significância"),paste("Como |zcal|<|ztab| Não Rejeita-se H0 ao nível de",alpha,"de significância"))
ConclusãoZ

## [1] "Como |zcal|>|ztab| Rejeita-se H0 ao nível de 0.05 de significância"

#O salário médio dos torneiros mecânicos é diferente do que o salário da indústria mecânica.

27-

28- Para verificar o grau de adesão de uma nova cola para vidros, preparam-se dois tipos de montagem: cruzado (A), onde a cola é posta em forma de X, e quadrado (B), onde a cola é posta apenas nas quatro bordas. Os resultados da resistência para as duas amostras de 10 cada estão abaixo. Que tipo de conclusão poderia ser tirada? Método A 16 14 19 18 19 20 15 18 17 18 Método B 13 19 14 17 21 24 10 14 13 15

n1=10
n2=10
xbarra1 <- mean(c(16,14,19,18,19,20,15,18,17,18))
xbarra1

## [1] 17.4

xbarra2 <- mean(c(13,19,14,17,21,24,10,14,13,15))
xbarra2

## [1] 16

dp1 <- sd(c(16,14,19,18,19,20,15,18,17,18))
dp1

## [1] 1.897367

dp2 <- sd(c(13,19,14,17,21,24,10,14,13,15))
dp2

## [1] 4.242641

alpha=0.05

# Será utilizado o Teste F
# H0:sigma1^2=sigma2^2
# H1:sigma1^2<sigma2^2

Fcal <- (dp2^2)/(dp1^2)
Fcal

## [1] 5

Ftab <- qf(alpha,n2-1,n1-1)
Ftab

## [1] 0.3145749

ConclusaoF <- 
ifelse((Fcal)>(Ftab), paste("Como |Fcal|>|Ftab| Rejeita-se H0 ao nível de", alpha,"de significância"), paste("Como |Fcal|<|Ftab| Não rejeita-se H0 ao nível de", alpha,"de significância"))
ConclusaoF

## [1] "Como |Fcal|>|Ftab| Rejeita-se H0 ao nível de 0.05 de significância"

# Será utilizado o teste t
# H0:mu1=mu2
# H1:mu1!=mu2

Tcal <- (xbarra1-xbarra2)/sqrt(((dp1^2)/n1)+((dp2^2)/n2))
Tcal

## [1] 0.9525793

A <- (dp1^2)/n1
B <-(dp2^2)/n2 

df <- ((A+B)^2)/(((A^2)/(n1-1))+((B^2)/(n2-1)))
(df <- round(df,digits = 0))

## [1] 12

Ttab <- qt(alpha/2,df)
Ttab

## [1] -2.178813

ConclusaoT <- 
ifelse(abs(Tcal)>abs(Ttab), paste("Como |Tcal|>|Ttab| Rejeita-se H0 ao nível de", alpha ,"de significância"), paste("Como |Tcal|<|Ttab| Não Rejeita-se H0 ao nível de", alpha ,"de significância"))
ConclusaoT

## [1] "Como |Tcal|<|Ttab| Não Rejeita-se H0 ao nível de 0.05 de significância"

## [1] "Como |Tcal|<|Ttab| Não Rejeita-se H0 ao nível de 0.05 de significância"

29- Em um estudo para comparar os efeitos de duas dietas, A e B, sobre o crescimento, 6 ratos foram submetidos à dieta A, e 9 ratos à dieta B. Após 5 semanas, os ganhos em peso foram: | Dieta A: 15; 18; 12; 11; 14; 15 | Dieta B: 11; 11; 12; 16; 12; 13; 8; 10; 13

a)Admitindo que temos duas amostras independentes de populações normais, teste a hipótese de que não há diferença entre as duas dietas, contra a alternativa que a dieta A é mais eficaz, usando o teste t de Student, no nível de α = 0, 01.

x = c(16,14,19,18,19,20,15,18,17,18)
y = c(13,19,14,17,21,24,10,14,13,15)
nA=10
nB=10
xAbarra = mean(x)
xBbarra = mean(y)
sA = sd(x)
sB = sd(y)
alpha = 0.05
#teste-F (desvio populacional desconhecido)
#H0:SigmaA^2 = SigmaB^2
#H1:SigmaA^2 < SigmaB^2
Fcal = (sB^2)/(sA^2)
Fcal

## [1] 5

Ftab = qf(alpha,nB-1,nA-1)
Ftab

## [1] 0.3145749

ConclusaoF = ifelse(Fcal>Ftab,paste("Como Fcal>Ftab Rejeita-se H0 ao nível de", alpha ,"de significância"), paste("Como Fcal<Ftab Não Rejeita-se H0 ao nível de", alpha ,"de significância"))
ConclusaoF

## [1] "Como Fcal>Ftab Rejeita-se H0 ao nível de 0.05 de significância"

# A variância das populações são diferentes.
#teste-t
#H0:muA = muB
#H1:muA!= muB
Tcal = (xAbarra-xBbarra)/sqrt(((sA^2)/nA)+((sB^2)/nB))
Tcal

## [1] 0.9525793

A = (sA^2)/nA
B = (sB^2)/nB
df = ((A+B)^2)/(((A^2)/(nA-1))+((B^2)/(nB-1)))
(df = round(df, digits = 0))

## [1] 12

df

## [1] 12

Ttab = qt(alpha/2, df)
Ttab

## [1] -2.178813

ConclusaoT = ifelse(abs(Tcal)>abs(Ttab),paste("Como |Tcal|>|Ttab| Rejeita-se H0 ao nível de", alpha ,"de significância"), paste("Como |Tcal|<|Ttab| Não Rejeita-se H0 ao nível de", alpha,"de significância"))
ConclusaoT

## [1] "Como |Tcal|<|Ttab| Não Rejeita-se H0 ao nível de 0.05 de significância"

#Não existe diferença entre os tratamentos analisados.

29-1 Admitindo que temos duas amostras independentes de populações normais, teste a hipótese de que não há diferença entre as duas dietas, contra a alternativa que a dieta A é mais eficaz, usando o teste t de Student, no nível de α=0,01.

#H0: muA = muB
#H1: muA > muB
x = c(15, 18, 12, 11, 14, 15)
y = c(11, 11, 12, 16, 12, 13, 8, 10, 13)
xAbarra = mean(x)
xBbarra = mean(y)
nA = 6
nB = 9
sA = sd(x)
sB = sd(y)
alpha = 0.01
#teste-F (desvio populacional desconhecido)
#H0:SigmaA^2 = SigmaB^2
#H1:SigmaA^2 < SigmaB^2
Fcal = (sB^2)/(sA^2)
Fcal

## [1] 0.8018018

#H0: muA = muB
#H1: muA > muB
x = c(15, 18, 12, 11, 14, 15)
y = c(11, 11, 12, 16, 12, 13, 8, 10, 13)
xAbarra = mean(x)
xBbarra = mean(y)
nA = 6
nB = 9
sA = sd(x)
sB = sd(y)
alpha = 0.01
#teste-F (desvio populacional desconhecido)
#H0:SigmaA^2 = SigmaB^2
#H1:SigmaA^2 < SigmaB^2
Fcal = (sB^2)/(sA^2)
Fcal

## [1] 0.8018018

ConclusaoF <- ifelse(Fcal>Ftab,paste("Como Fcal>Ftab Rejeita-se H0 ao nível de", alpha ,"de significância"), paste("Como Fcal<Ftab Não Rejeita-se H0 ao nível de", alpha ,"de significância"))
ConclusaoF

## [1] "Como Fcal>Ftab Rejeita-se H0 ao nível de 0.01 de significância"

# A variância das populações são diferentes.

#teste-t
#H0:muA = muB
#H1:muA > muB
Tcal = (xAbarra-xBbarra)/sqrt(((sA^2)/nA)+((sB^2)/nB))
Tcal

## [1] 1.902208

A <- (sA^2)/nA
B <- (sB^2)/nB
df = ((A+B)^2)/(((A^2)/(nA-1))+((B^2)/(nB-1)))
(df = round(df, digits = 0))

## [1] 10

df

## [1] 10

Ttab <- qt(alpha/2, df)
Ttab

## [1] -3.169273

ConclusaoT <- ifelse(abs(Tcal)>abs(Ttab),paste("Como |Tcal|>|Ttab| Rejeita-se H0 ao nível de", alpha ,"de significância"), paste("Como |Tcal|<|Ttab| Não Rejeita-se H0 ao nível de", alpha,"de significância"))
ConclusaoT

## [1] "Como |Tcal|<|Ttab| Não Rejeita-se H0 ao nível de 0.01 de significância"

#Não existe diferença entre os tratamentos analisados

30- Suponha que o tempo necessário para atendimento de clientes em uma central de atendimento telefônico siga uma distribuição normal de média de 8 minutos e desvio padrão de 2 minutos. a) Qual é a probabilidade de que um atendimento dure menos de 5 minutos?

mu=8
sigma=2
pnorm(q=5, mean = mu, sd=sigma)

## [1] 0.0668072

pnormGC(5, region="below", mean=mu,
        sd=sigma, graph=TRUE)

## [1] 0.0668072

E mais do que 9,5 minutos?

mu=8
sigma=2
pnorm(9.5, mean = mu, sd=sigma, lower.tail = FALSE)

## [1] 0.2266274

pnormGC(9.5, region="above", mean=mu,
        sd=sigma,graph=TRUE)

## [1] 0.2266274

E entre 7 e 10 minutos?

mu=8
sigma=2
pnorm(10, mean = mu, sd=sigma)-pnorm(7, mean = mu, sd=sigma)

## [1] 0.5328072

pnormGC(c(7,10), region="between", mean=mu, sd=sigma,graph=TRUE)

## [1] 0.5328072

75% das chamadas telefônicas requerem pelo menos quanto tempo de atendimento?

qnorm(0.75, mean = 8, sd = 2, lower.tail = FALSE)

## [1] 6.65102

qnormGC(0.75, region = "above", mean = 8, sd = 2, graph = T )

## [1] 6.65102

31- A distribuição dos pesos de coelhos criados numa granja pode muito bem ser representada por uma distribuição Normal, com média 5 kg e desvio padrão 0,9 kg. Um abatedouro comprará 5000 coelhos e pretende classificá-los de acordo com o peso do seguinte modo: 15% dos mais leves como pequenos, os 50% seguintes como médios, os 20% seguintes como grandes e os 15% mais pesados como extras. Quais os limites de peso para cada classificação? - Pequenos

qnorm(0.15, mean = 5, sd = 0.9, lower.tail = TRUE)

## [1] 4.06721

qnormGC(0.15, region = "below", mean = 5, sd = 0.9, graph = TRUE )

## [1] 4.06721

Médios

qnorm(0.15, mean = 5, sd = 0.9, lower.tail = FALSE)- qnorm(0.65, mean = 5, sd = 0.9, lower.tail = T)

## [1] 0.5860016

qnormGC(0.5860016, region = "between", mean = 5, sd = 0.9, graph = T )

## [1] 4.26481 5.73519

Grandes

qnorm(0.85, mean = 5, sd = 0.9, lower.tail = TRUE)-qnorm(0.85, mean = 5, sd = 0.9, lower.tail = FALSE)

## [1] 1.86558

qnormGC( 0.86558, region = "between", mean = 5, sd = 0.9, graph = T )

## [1] 3.652792 6.347208

qnorm(0.15, mean = 5, sd = 0.9, lower.tail = TRUE)-1

## [1] 3.06721

qnormGC( 0.86558, region = "between", mean = 5, sd = 0.9, graph = T )

## [1] 3.652792 6.347208

32- Uma enchedora automática de refrigerantes está regulada para que o volume médio de líquido em cada garrafa seja de 1000cm3 e desvio padrão de 10m3. Admita que o volume siga uma distribuição normal.

Qual é a porcentagem de garrafas em que o volume de líquido é menor que 990cm3?

mu=1000
sigma=10*1000000
pnorm(q=990, mean = mu, sd=sigma)

## [1] 0.4999996

pnormGC(990, region="below", mean=mu, sd=sigma, graph=T)

## [1] 0.4999996

Qual é a porcentagem de garrafas em que o volume de líquido não se desvia da média em mais do que dois desvios padrões?

mu=1000
sigma=10*1000000
pnorm(10*1000000*2, mean = mu, sd=sigma, lower.tail = FALSE)

## [1] 0.02275553

pnormGC(10*1000000*2, region="above", mean=mu, sd=sigma,graph=T)

## [1] 0.02275553

33- Uma empresa produz televisores de 2 tipos, tipo A (comum) e tipo B (luxo), e garante a restituição da quantia paga se qualquer televisor apresentar defeito grave no prazo de seis meses. O tempo para ocorrência de algum defeito grave nos televisores tem distribuição normal sendo que, no tipo A, com média de 10 meses e desvio padrão de 2 meses e no tipo B, com média de 11 meses e desvio padrão de 3 meses. Os televisores de tipo A e B são produzidos com lucro de 1200 u.m. e 2100 u.m. respectivamente e, caso haja restituição, com prejuízo de 2500 u.m. e 7000 u.m. Respectivamente.

Calcule as probabilidades de haver restituição nos televisores do tipo A e do tipo B.

#Tipo A
mu=10
sigma=2
proba <- pnorm(q=6, mean = mu, sd=sigma)
pnorm(q=6, mean = mu, sd=sigma)

## [1] 0.02275013

#Tipo B
mu=11
sigma=3
probb <- pnorm(q=6, mean = mu, sd=sigma)
pnorm(q=6, mean = mu, sd=sigma)

## [1] 0.04779035

Calcule o lucro médio para os televisores do tipo A e para os televisores do tipo B.

# Lucro médio tipo A

#Lucro médio Tipo B
1-probb

## [1] 0.9522096

2100*1-probb - 7000*probb

## [1] 1765.42

Baseando-se nos lucros médios, a empresa deveria incentivar as vendas dos aparelhos do tipo A ou do tipo B?

print("A empresa deveria incentivar as vendas dos aparelhos do tipo B, pois o lucro B é maior que o lucro médio de A.")

## [1] "A empresa deveria incentivar as vendas dos aparelhos do tipo B, pois o lucro B é maior que o lucro médio de A."

## [1] "A empresa deveria incentivar as vendas dos aparelhos do tipo B, pois o lucro B é maior que o lucro médio de A."

34- Um estudo comparou dois métodos (A e B) para ensinar matemática a alunos do primeiro grau. Após 10 semanas, o desempenho dos alunos foi avaliado em um teste. Teste a hipótese de que o método A resulta num melhor desempenho médio, ao nível α=5%, com base nos resultados da tabela a seguir:

nA <- 10
nB <- 8
xAbarra <- 8.15
xBbarra <- 7.31
sA <- 1.15
sB <- 1.94
#H0:muA=muB
#H1:muA!=muB (Teste-t bilateral)
#Precisamos proceder antes do teste-T o teste-F.
#H0:SigmaA^2=SigmaB^2
#H1:SigmaA^2<SigmaB^2 (Teste unilateral)

(Fcal <- (sB^2)/(sA^2))

## [1] 2.845822

(pvalor <- pf(q=Fcal, df1 = nB-1, df2 = nA-1))

## [1] 0.9267392

#Conclusao: Não rejeita-se H0 para todo alpha < pvalor. Assim, para realizar o teste - t, vamos considerar que as variâncias do método A 
#e método B são iguais (SigmaA^2=SigmaB^2). Dessa forma, procedemos 
# o teste-t para variâncias populacionais iguais.

#H0:muA=muB
#H1:muA!=muB (Teste-t bilateral)

#A <- (sA^2)/nA
#B <- (sB^2)/nB
#
#df <- ((A+B)^2)/(((A^2)/(nA-1))+((B^2)/(nB-1)))
df <- nA+nB-2  
Sc <- sqrt((((sA^2)*(nA-1))+((sB^2)*(nB-1)))/(nA+nB-2))  
#Estatística do Teste

ConclusaoT <- "Não rejeita-se H0 para alpha menor que pvalor, portanto não a diferença entre os testes"

35- A lei trabalhista estabelece que o pagamento diário mínimo deve ser de 13, 20 U.M. (unidades monetárias). Assuma distribuição normal com desvio padrão igual a 2,0 U.M. Uma amostra aleatória de 40 trabalhadores de uma firma revelou média diária de 12,20 U.M .Esta firma deve ser acusada de estar infringindo a lei? Conclua a 1% de probabilidade.

u=13.20
sigma=2.0
n=40
alpha=0.01
xbarra=12.20

# Será utilizado o teste z
# H0:mu=13.20
# H1:mu<13.20
# Estatística do teste:

Zcal <- (xbarra-mu)/(sigma/sqrt(n))
Zcal

## [1] 3.794733

Ztab <- qnorm(alpha, mean = 0, sd=1, lower.tail = FALSE)
Ztab

## [1] 2.326348

ConclusaoZ <- ifelse(abs(Zcal)>abs(Ztab),paste(
"Como |Zcal|>|Ztab| Rejeita-se H0 ao nível de"
, alpha ,
"de significância"
), paste(
"Como |Zcal|<|Ztab| Não Rejeita-se H0 ao nível de"
, alpha ,
"de significância"
))
ConclusaoZ

## [1] "Como |Zcal|>|Ztab| Rejeita-se H0 ao nível de 0.01 de significância"

#Logo, conclui- se que essa firma deve ser acusada de infrigir a lei!
(pvalor <- pnorm(Zcal, lower.tail = TRUE))

## [1] 0.9999261

ConclusaoZ <- ifelse(pvalor>alpha,paste(
"Como p-valor>", alpha, " Não Rejeita-se H0"
), paste(
"Como p-valor<", alpha, " Rejeita-se H0")
)
ConclusaoZ

## [1] "Como p-valor> 0.01  Não Rejeita-se H0"

36- A tabela a seguir mostra a frequência de acidentes automobilísticos por ano, de acordo som a faixa etária (idade) do motorista, para motoristas com idade inferior a 25 anos. Teste a hipótese de que o número de acidentes independe da idade, a 5% de probabilidade. Isto é, teste a hipótese de que o número anual de acidentes se distribui proporcionalmente nas faixas etárias. A tabela abaixo apresenta o percentual de motoristas em cada faixa etária.

% de motoristas 10 20 20 25 25 idade (anos) 15-16 17-18 19-20 21-22 23-24 número de acidentes 8 15 13 11 8

alpha <- 0.05
(TotalAcidentes <- 8+15+13+11+8)

## [1] 55

Obs <- c(8,15,13,11,8)
#H0:O percentual de acidentes independe da idade
#H1:Não H0
XiQuad <- chisq.test(Obs, p = c(0.1,0.2,0.2,0.25,0.25))
XiQuad

## 
##  Chi-squared test for given probabilities
## 
## data:  Obs
## X-squared = 5.9091, df = 4, p-value = 0.206

## 
##  Chi-squared test for given probabilities
## 
## data:  Obs
## X-squared = 5.9091, df = 4, p-value = 0.206

37- Uma indústria farmacêutica conduziu um estudo para avaliar o tempo médio em dias para recuperação dos efeitos da gripe. O estudo comparou o tempo de indivíduos que tomaram 500 mg diárias de vitamina C, contra indivíduos que não tomaram vitamina C (nenhum suplemento). Com base nos dados a seguir, conclua e interprete a 5% de probabilidade.

Nenhum suplemento 500mg Vit. C Tamanho da amostra 12 12 Tempo médio 7,4 5,8 Variâncias 2,9 2,4

nA <- 12
nB <- 12
xAbarra <- 7.4
xBbarra <- 5.8
sA <- sqrt(2.9)
sB <- sqrt(2.4)
#H0:muA=muB
#H1:muA!=muB (Teste-t bilateral)
#Precisamos proceder antes do teste-T o teste-F.
#H0:SigmaA^2=SigmaB^2
#H1:SigmaA^2<SigmaB^2 (Teste unilateral)

(Fcal <- (sB^2)/(sA^2))

## [1] 0.8275862

(pvalor <- pf(q=Fcal, df1 = nB-1, df2 = nA-1))

## [1] 0.3795903

#Conclusao: Não rejeita-se H0 para todo alpha < pvalor. Assim, para realizar o teste - t, vamos considerar que as variâncias do método A 
#e método B são iguais (SigmaA^2=SigmaB^2). Dessa forma, procedemos 
# o teste-t para variâncias populacionais iguais.

#H0:muA=muB
#H1:muA!=muB (Teste-t bilateral)

#A <- (sA^2)/nA
#B <- (sB^2)/nB
#
#df <- ((A+B)^2)/(((A^2)/(nA-1))+((B^2)/(nB-1)))
df <- nA+nB-2  
Sc <- sqrt((((sA^2)*(nA-1))+((sB^2)*(nB-1)))/(nA+nB-2))  
#Estatística do Teste
Tcal <- (xAbarra-xBbarra)/(Sc*sqrt((1/nA)+(1/nB)))
Tcal

## [1] 2.407535

ConclusaoT <- "Não rejeita-se H0 para alpha menor que pvalor"

38- Um pesquisa de opinião entrevistou 50 pessoas em dois distritos. O objetivo era verificar se a distribuição das opiniões era homogênea nos dois distritos. Com base nos dados da tabela, teste a hipótese de homogeneidade de opiniões usando α=5%.

Sim Indeciso Não Total Distrito A 20 9 21 50 Distrito B 26 3 21 50 Total 46 12 42 100

Fobs <- data.frame(Distrito=c("Distrito A", "Distrito B"), Sim=c(20,26),
                   Indeciso=c(9,3), Nao=c(21,21), row.names = TRUE)

#H0: As opiniões são homogêneas
#H1: Não H0

chisq.test(Fobs, correct = TRUE)

## 
##  Pearson's Chi-squared test
## 
## data:  Fobs
## X-squared = 3.7826, df = 2, p-value = 0.1509

## 
##  Pearson's Chi-squared test
## 
## data:  Fobs
## X-squared = 3.7826, df = 2, p-value = 0.1509

39- Uma associação comercial afirma que o número médio de dias de trabalho perdidos anualmente, devido a problemas de saúde, é igual a 60. Uma extensa campanha educacional visando a conscientizar os trabalhadores quanto a importância de uma alimentação balanceada, higiene pessoal, prática de esportes etc, foi conduzida com o intuito de melhorar este quadro. Um ano após esta campanha, um estudo com 30 trabalhadores forneceu média igual a 55 dias. Assuma que o número de dias de trabalho perdidos anualmente é normalmente distribuído com variância σ2=275. Pede-se:

a)Pode-se afirmar que a campanha foi eficaz ao nível de α=1% de probabilidade?

mu <- 60
sigma <- 275
n <- 30
xbarra <- 55
alpha <- 0.01
#H0:mu=60
#H1:mu<60

#Estatística do Teste
Zcal <- (xbarra-mu)/(sigma/sqrt(n))
Zcal

## [1] -0.09958592

Ztab <- qnorm(alpha, mean = 0, sd=1, lower.tail = FALSE)
Ztab

## [1] 2.326348

ConclusaoZ <- ifelse(abs(Zcal)>abs(Ztab),paste(
"Como |Zcal|>|Ztab| Rejeita-se H0 ao nível de"
, alpha ,
"de significância"
), paste(
"Como |Zcal|<|Ztab| Não Rejeita-se H0 ao nível de"
, alpha ,
"de significância"
))
ConclusaoZ

## [1] "Como |Zcal|<|Ztab| Não Rejeita-se H0 ao nível de 0.01 de significância"

#Logo, conclui- se que essa firma deve ser acusada de infrigir a lei!
(pvalor <- pnorm(Zcal, lower.tail = TRUE))

## [1] 0.4603365

ConclusaoZ <- ifelse(pvalor>alpha,paste(
"Como p-valor>", alpha, " Não Rejeita-se H0"
), paste(
"Como p-valor<", alpha, " Rejeita-se H0")
)
ConclusaoZ

## [1] "Como p-valor> 0.01  Não Rejeita-se H0"

Para qual nível de significância se pode afirmar que a campanha educacional foi eficaz?

print("Podemos afirmar que a campanha educacional foi eficaz para valores abaixo de 46% de sgnificancia")

## [1] "Podemos afirmar que a campanha educacional foi eficaz para valores abaixo de 46% de sgnificancia"

40- Um gerente comercial acredita que um número excessivo de horas estejam sendo desperdiçadas em contatos comerciais, via telefone, entre os seus vendedores e os clientes em potencial. Ele deseja no máximo quinze horas por semana por vendedor. Este gerente comercial contratou uma empresa especializada para treinar seus vendedores. Após este treinamento, uma amostra de 36 vendedores revelou média igual a 17h por semana por vendedor. O que pode ser concluído quanto a eficácia do treinamento? Assuma σ2=9 e utilize α=5%.

mu <- 15
sigma <- 9
n <- 36
xbarra <- 17
alpha <- 0.05
#H0:mu=15
#H1:mu>15

#Estatística do Teste
Zcal <- (xbarra-mu)/(sigma/sqrt(n))
Zcal

## [1] 1.333333

Ztab <- qnorm(alpha, mean = 0, sd=1, lower.tail = TRUE)
Ztab

## [1] -1.644854

ConclusaoZ <- ifelse(abs(Zcal)>abs(Ztab),paste(
"Como |Zcal|>|Ztab| Rejeita-se H0 ao nível de"
, alpha ,
"de significância"
), paste(
"Como |Zcal|<|Ztab| Não Rejeita-se H0 ao nível de"
, alpha ,
"de significância"
))
ConclusaoZ

## [1] "Como |Zcal|<|Ztab| Não Rejeita-se H0 ao nível de 0.05 de significância"

#Logo, conclui- se que o treinamento foi eficaz!
(pvalor <- pnorm(Zcal, lower.tail = FALSE))

## [1] 0.09121122

ConclusaoZ <- ifelse(pvalor>alpha,paste(
"Como p-valor>", alpha, " Não Rejeita-se H0"
), paste(
"Como p-valor<", alpha, " Rejeita-se H0")
)
ConclusaoZ

## [1] "Como p-valor> 0.05  Não Rejeita-se H0"

41- Com base em dados obtidos de 400 mulheres, apresentados na tabela abaixo, pode-se concluir que o nível educacional e a adaptação à vida conjugal são independentes? Conclua a 5% de probabilidade.

Fobs <- data.frame(Distrito=c("Universidade", "2° grau", "3° grau"), Ruim=c(18,17,11), Razoável=c(29,28,10), boa=c(70,30,11 ), Muitoboa=c(115,41,20),
row.names = TRUE)

#H0: Nível educacional e a adaptação à vida conjugal são independentes
#H1: Não H0

chisq.test(Fobs, correct = T)

## 
##  Pearson's Chi-squared test
## 
## data:  Fobs
## X-squared = 19.943, df = 6, p-value = 0.002835

## 
##  Pearson's Chi-squared test
## 
## data:  Fobs
## X-squared = 19.943, df = 6, p-value = 0.002835

print("Não são independentes uma vez que se rejeita H0 pelo p valor muito baixo")

## [1] "Não são independentes uma vez que se rejeita H0 pelo p valor muito baixo"

## [1] "Não são independentes uma vez que se rejeita H0 pelo p valor muito baixo"

42- Uma cooperativa de produtores possui uma máquina de encher vasilhame com um litro de leite. Para assegurar que em média cada vasilhame não terá leite a mais e nem a menos, o responsável pelo controle de qualidade amostra, semanalmente, 75 vasilhames enchidos pela máquina. Se uma amostra fornecer 63, 97 litros e desvio padrão s=0,25 litros, deve-se parar a máquina para regulagem ou continuar a produção? Qual deve ser o procedimento adotado a α=5% de probabilidade?

mu <- 75
sigma <- 0.25
n <- 75
xbarra <- 67.97
alpha <- 0.05
#H0:mu=75
#H1:mu<75

#Estatística do Teste
Zcal <- (xbarra-mu)/(sigma/sqrt(n))
Zcal

## [1] -243.5263

## [1] -243.5263

Ztab <- qnorm(alpha, mean = 0, sd=1, lower.tail = FALSE)
Ztab

## [1] 1.644854

ConclusaoZ <- ifelse(abs(Zcal)>abs(Ztab),paste(
"Como |Zcal|>|Ztab| Rejeita-se H0 ao nível de"
, alpha ,
"de significância"
), paste(
"Como |Zcal|<|Ztab| Não Rejeita-se H0 ao nível de"
, alpha ,
"de significância"
))
ConclusaoZ

## [1] "Como |Zcal|>|Ztab| Rejeita-se H0 ao nível de 0.05 de significância"

#Logo, conclui- se que deve-se paarar a máquina para regulagem!
(pvalor <- pnorm(Zcal, lower.tail = TRUE))

## [1] 0

ConclusaoZ <- ifelse(pvalor>alpha,paste(
"Como p-valor>", alpha, " Não Rejeita-se H0"
), paste(
"Como p-valor<", alpha, " Rejeita-se H0")
)
ConclusaoZ

## [1] "Como p-valor< 0.05  Rejeita-se H0"

43-

44- Assuma que o consumo mensal per capita de determinado produto tem distribuição normal com desvio padrão igual a 5 kg. Com a atual crise (do dólar, do apagão, do futebol…várias opções!) o departamento de vendas da fábrica decidiu que irá retirar o produto do mercado, caso o consumo médio (μ) per capita seja inferior a 10kg. Se uma pesquisa de mercado, com uma amostra de 100 indivíduos, revelar consumo médio mensal per capita de 9 kg, pede-se: Qual deve ser a afirmação, ao nível de significância de 1,5%?

mu <- 10
sigma <- 5
n <- 100
xbarra <- 9
alpha <- 0.015
#H0:mu=10
#H1:mu<10

#Estatística do Teste
Zcal <- (xbarra-mu)/(sigma/sqrt(n))
Zcal

## [1] -2

Ztab <- qnorm(alpha, mean = 0, sd=1, lower.tail = FALSE)
Ztab

## [1] 2.17009

ConclusaoZ <- ifelse(abs(Zcal)>abs(Ztab),paste(
"Como |Zcal|>|Ztab| Rejeita-se H0 ao nível de"
, alpha ,
"de significância"
), paste(
"Como |Zcal|<|Ztab| Não Rejeita-se H0 ao nível de"
, alpha ,
"de significância"
))
ConclusaoZ

## [1] "Como |Zcal|<|Ztab| Não Rejeita-se H0 ao nível de 0.015 de significância"

#Logo, conclui- se que a fábrica não deve retirar o produto do mercado!
(pvalor <- pnorm(Zcal, lower.tail = T))

## [1] 0.02275013

ConclusaoZ <- ifelse(pvalor>alpha,paste(
"Como p-valor>", alpha, " Não Rejeita-se H0"
), paste(
"Como p-valor<", alpha, " Rejeita-se H0")
)
ConclusaoZ

## [1] "Como p-valor> 0.015  Não Rejeita-se H0"

45- No quadro abaixo estão as opiniões, com respeito ao desempenho e a potência do motor, de proprietários de veículos de um determinado fabricante. As opiniões foram classificadas pela idade do proprietário.

Idade Ruim Bom Jovem 30 20 Experiente 20 30

O que pode ser afirmado quanto à seguinte hipótese de nulidade? H0 : Idade e opinião são independentes.

Fobs <- data.frame(Distrito=c("Jovem", "Experiente"), Ruim=c(30,20),
                   Bom=c(20,30),  row.names = TRUE)

#H0: Idade e opinião são independentes
#H1: Não H0

chisq.test(Fobs, correct = TRUE)

## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  Fobs
## X-squared = 3.24, df = 1, p-value = 0.07186

## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  Fobs
## X-squared = 3.24, df = 1, p-value = 0.07186

46- Para comparar duas marcas de pará-choques, montaram-se seis de cada marca em 12 carros compactos, fazendo-se cada carro colidir com um muro de concreto, a uma velocidade de 40 km Registraram-se os seguintes custos de reparo:

Teste (α=5%) a hipótese de igualdade entre os custos médios de reparo dos pará-choques.

A <- c(320,310,380,360,320,345)
B <- c(305,290,340,315,280,305)
sd(A)

## [1] 27.27942

sd(B)

## [1] 20.83667

nA <- length(A)
nB <- length(B)
#Teste-F
#H0:SigmaA^2=SigmaB^2
#H1:SigmaA^2!=SigmaB^2

var.test(A, B, alternative = "two.sided")

## 
##  F test to compare two variances
## 
## data:  A and B
## F = 1.714, num df = 5, denom df = 5, p-value = 0.5687
## alternative hypothesis: true ratio of variances is not equal to 1
## 95 percent confidence interval:
##   0.2398433 12.2489808
## sample estimates:
## ratio of variances 
##           1.714012

print("Para alpha > p-value=0.5687 rejeita-se H0, logo, devemos proceder o teste-t para variâncias desiguais.")

## [1] "Para alpha > p-value=0.5687 rejeita-se H0, logo, devemos proceder o teste-t para variâncias desiguais."

#H0:muLi=muAd
#H1:muLi!=muAd

t.test(A, B, alternative = "two.sided", var.equal = FALSE)

## 
##  Welch Two Sample t-test
## 
## data:  A and B
## t = 2.3786, df = 9.3527, p-value = 0.04033
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
##   1.813021 64.853645
## sample estimates:
## mean of x mean of y 
##  339.1667  305.8333

print("Para alpha menor do que p-value=0.04033, não rejeita-se H0. Ou seja nesse caso, os custos médios de reparo são desiguais.")

## [1] "Para alpha menor do que p-value=0.04033, não rejeita-se H0. Ou seja nesse caso, os custos médios de reparo são desiguais."

47- Se um dado não é viciado cada uma das seis faces ocorre com igual probabilidade. Um determinado dado foi lançado 720 vezes, obtendo-se: Face 1 2 3 4 5 6 Total Frequência observada 129 107 98 132 136 118 720 O dado será considerado viciado para qual nível de significância? Explique sua resposta.

(TotalAcidentes <- 720)

## [1] 720

Obs <- c(129,107,98,132,136,118)
#H0:Dado viciado
#H1:Não H0
XiQuad <- chisq.test(Obs, p = c(1/6,1/6,1/6,1/6,1/6,1/6))
XiQuad

## 
##  Chi-squared test for given probabilities
## 
## data:  Obs
## X-squared = 9.4833, df = 5, p-value = 0.09127

print("o dado será considerado um dado vicido para um nivel de significância maior que 0.09127, que foi o p-valor calculado.")

## [1] "o dado será considerado um dado vicido para um nivel de significância maior que 0.09127, que foi o p-valor calculado."

## [1] "o dado será considerado um dado vicido para um nivel de significância maior que 0.09127, que foi o p-valor calculado."

48- O tempo médio, por operário, para executar uma tarefa, tem sido 100 minutos. Introduziu-se uma modificação para diminuir esse tempo, e, após certo período, sorteou-se uma amostra de 16 operários, medindo-se o tempo de execução de cada um. O tempo médio da amostra foi 85 minutos, e o desvio padrão foi 12 minutos. Estes resultados trazem evidências estatísticas da melhora desejada, considerando α=5%? Apresente as suposições teóricas usadas para resolver problema.

mu <- 100
sigma <- 12
n <- 16
xbarra <- 85
alpha <- 0.05
#H0:mu=100
#H1:mu<100

#Estatística do Teste
Zcal <- (xbarra-mu)/(sigma/sqrt(n))
Zcal

## [1] -5

Ztab <- qnorm(alpha, mean = 0, sd=1, lower.tail = FALSE)
Ztab

## [1] 1.644854

ConclusaoZ <- ifelse(abs(Zcal)>abs(Ztab),paste(
"Como |Zcal|>|Ztab| Rejeita-se H0 ao nível de"
, alpha ,
"de significância"
), paste(
"Como |Zcal|<|Ztab| Não Rejeita-se H0 ao nível de"
, alpha ,
"de significância"
))
ConclusaoZ

## [1] "Como |Zcal|>|Ztab| Rejeita-se H0 ao nível de 0.05 de significância"

#Logo, conclui- se que o tempo para um operário realizar uma tarefa não melhorou!
(pvalor <- pnorm(Zcal, lower.tail = T))

## [1] 2.866516e-07

ConclusaoZ <- ifelse(pvalor>alpha,paste(
"Como p-valor>", alpha, " Não Rejeita-se H0"
), paste(
"Como p-valor<", alpha, " Rejeita-se H0")
)
ConclusaoZ

## [1] "Como p-valor< 0.05  Rejeita-se H0"

Relatório de Teste de Hipóteses

Diego José Aparecido Paulino

2022-07-26

Relatório para a disciplina MAF 160 (Teste de Hipóteses)