# Resolução da Lista 01 # Número USP: ABCDEFG A=2 B=3 C=4 D=6 E=1 F=1 G=1 # # Definição inicial de valores P1=86+F P2=91+G P3=33+E P4=1+A+G P5=15+B P6=73+C # Amostras Amostra1 <- c(P1,P4,P2,25,28,20,58,5,37,70) Amostra2 <- c(P2,P5,P3,92,40,13,57,80,16,10) Amostra3 <- c(P3,P6,P1,18,95,30,73,65,9,27) # # Leitura dos dados no arquivo Cap01_Corrar_etal_2007.csv dados <- read.table("Cap01_Corrar_etal_2007.csv",header=TRUE,sep=";",dec=",") # # Extrai amostras da base de dados dados1 <- dados[Amostra1,] dados2 <- dados[Amostra2,] dados3 <- dados[Amostra3,] # # 1.1 Medidas de centro para PL # Amostra 1 mean(dados1$PL) median(dados1$PL) (max(dados1$PL)+min(dados1$PL))/2 # Amostra 2 mean(dados2$PL) median(dados2$PL) (max(dados2$PL)+min(dados2$PL))/2 # Amostra 3 mean(dados3$PL) median(dados3$PL) (max(dados3$PL)+min(dados3$PL))/2 # # 1.2.a Resumo dos cinco números: # mínimo, quartil 1, mediana, quartil 3, máximo quantile(dados1$PL, probs = c(0.0, 0.25, 0.50, 0.75, 1)) quantile(dados2$PL, probs = c(0.0, 0.25, 0.50, 0.75, 1)) quantile(dados3$PL, probs = c(0.0, 0.25, 0.50, 0.75, 1)) # # Dados para boxplot dados_bp <- cbind(dados1$PL,dados2$PL,dados3$PL) dados_bp <- as.data.frame(dados_bp) names(dados_bp) <- c("PL1","PL2","PL3") dados_bp # # 1.2.b Boxplot boxplot(dados_bp$PL1,dados_bp$PL2,dados_bp$PL3) # #----------------------------------------------- # Tabela Exercício 2 Casados <- c(15+C,5+E) Solteiros <- c(4+F,20+D) Tab <- cbind(Casados,Solteiros) rownames(Tab) <- c("C.Exper","S.Exper") Tab # 2.1 Escolhendo ao acaso um candidato entrevistado: # 2.1.a A probabilidade de selecionar um candidato # com nenhuma experiência. (Tab[2]+Tab[4])/sum(Tab) # 2.1.b A probabilidade de selecionar um candidato solteiro. (Tab[3]+Tab[4])/sum(Tab) # 2.1.c A probabilidade de selecionar um candidato # com alguma experiência e solteiro. Tab[3]/sum(Tab) # 2.1.d A probabilidade de selecionar um candidato # casado ou com alguma experiência. (sum(Tab)-Tab[4])/sum(Tab) # 2.1.e Dado que o candidato selecionado é casado, # a probabilidade que o mesmo tenha alguma experiência. Tab[1]/(Tab[1]+Tab[2]) # #----------------------------------------------- # Valores para os exercícios 3 e 4 C1 = 300 + E + F + G C2 = 230 + F C3 = 90 + G # # Observe a função normal acumulada: # pnorm(2.0, mean = 0, sd = 1) # # 3 m3 = C1 dp3 = 30 # # 3.1 vmin = 310 vmax = 355 # Probabilidade até o valor pvmin = pnorm(vmin, mean = m3, sd = dp3) pvmax = pnorm(vmax, mean = m3, sd = dp3) # Probabilidade entre os valores P3 = pvmax - pvmin P3 # 3.2 valor = 320 pmaior = 1-pnorm(valor, mean = m3, sd = dp3) pmaior # #----------------------------------------------- # 4 nc = C3/100 alfa = 1-nc a2 = alfa/2 # # 4.1 x = C2 n = C1 p = x/n q = 1 - p z_a2 =qnorm(nc+a2,0,1) E = z_a2 * sqrt(p*q/n) Intervalo = c(p-E,p+E) Intervalo # 4.2 E_max = 0.02 n2 = (z_a2^2)*0.25/(E_max^2) n2 # Arredondar o resultado de n2 para o inteiro acima # # 4.3 N = 1000 E3 = E * sqrt((N-n)/(N-1)) Intervalo3 = c(p-E3,p+E3) Intervalo3