#Exemplo 3.2: teor de gordura #entra com os dados gordura=c( 278, 182, 247, 227, 277, 194, 196, 276, 244, 192, 118, 219, 255, 201, 209, 219, 228, 209, 209, 171, 213, 233, 226, 209, 200, 200, 363, 209, 200, 179, 167, 192, 277, 317, 146, 217, 292, 217, 255, 212, 233, 250, 243, 150, 209, 184, 199, 250, 479, 175, 194, 221, 233, 184, 217, 150, 167, 265, 242, 180, 255, 170, 209, 161, 196, 165, 234, 179, 248, 184, 291, 185, 242, 276, 243, 229, 242, 250) ################################################### # Faz a tabela de classes de frequencias ################################################### # calcula a amplitude dos dados range(gordura) #calcula o numero de classes para o histograma usando o metodo Sturges k=nclass.Sturges(gordura) #Acha o minimo e o maximo dos dados min.gor=min(gordura) max.gor=max(gordura) #Cria as classes para gerar a tabela (e arredonda h) h=(max.gor-min.gor)/k h=45 #Arredonda o minimo e o maximo l.min=100 l.max=500 #Quebra o eixo x e grupos x.axis=seq(l.min, l.max, by=h) x.axis=c(x.axis,max(x.axis)+h) #Gera a tabela de contagem dos dados de gordura tab_gordura=table(cut(gordura, x.axis, right=FALSE)) tab_gordura #Faz a tabela de proporcao dos dados e arredonda em 2 casas decimas tab=round(prop.table(tab_gordura),2) tab #Faz o histograma hist(gordura, breaks=x.axis, right=FALSE, main="Histograma", ylab="Frequência", xaxt = "n") axis(1,x.axis) #Faz a ogiva (grafico das funcao acumulada plot(x.axis,cumsum(c(0,tab)),axes = FALSE, ylab="Freq. relativa acumulada", xlab="x") lines(x.axis,cumsum(c(0,tab))) axis(1,x.axis) axis(2,seq(0,1,by=0.1)) #Faz um resumo dos dados de gordura summary(gordura) #Parametros de locacao #media mean(gordura) #mediana median(gordura) #Calcula estatistica de variabilidade #variancia var(gordura) #desvio padrao sd(gordura) #coeficiente de variacao cv=sd(gordura)/mean(gordura) cv #Calcula quantis P10=quantile(gordura,0.10) P90=quantile(gordura,0.90) P10 P90 #Calcula amplitude amplitude=max(gordura)-min(gordura) amplitude #Faz o diagrama de folhas stem(gordura) #Gera o box-plot passo a passo #acha o os limites do box-plot Q1=quantile(gordura,0.25) Q3=quantile(gordura,0.75) Q1 Q3 #medida de dispercao DI=Q3-Q1 DI #limites LI=Q1-1.5*DI LS=Q3+1.5*DI LI LS #Faz o box-plot boxplot(gordura)