#slide 6 #lendo os dados no R peso=read.table('birthwt.dat',header=TRUE,sep='') peso[1,] #slide 9 table(peso$race) #freq absoluta em cada categoria table(peso$race)/length(peso$race) #freq relativa em cada categoria #slide 10 h=hist(peso$age,xlab='Idade (anos)',ylab='Frequência absoluta',main='',col=grey(.5)) h$counts h$counts/sum(h$counts)*100 cumsum(round(h$counts/sum(h$counts)*100,2)) #slide 14 barplot(table(peso$race),names.arg=c("Branca","Negra","Outra"),ylab="Frequências absolutas",xlab="Raça") pie(table(peso$race),labels=c("Branca","Negra","Outra"),clockwise=TRUE,col=grey.colors(3)) #slide 16 hist(peso$age,xlab='Idade',ylab='Frequência absoluta',main='',col=grey(.5)) #slide 21 round(table(peso$smoke)/length(peso$smoke)*100,1) round(table(peso$hyper)/length(peso$hyper)*100,1) #round(table(peso$nprem)/length(peso$nprem)*100,1) #slide 28 summary(peso$age) boxplot(peso$age, range=0,xlab='Idade (anos)',horizontal=TRUE) rug(jitter(peso$age,amount=0.5), col='blue') #slide 33 fra=cumsum(h$counts/length(peso$age)) yy=rep(c(0,fra),rep(2,8)) xx=c(0,rep(h$breaks[-8],rep(2,7)),h$breaks[8]) plot(xx,yy,type='s',xlab='Idade (anos)', ylab='Frequência relativa acumulada') lines(h$breaks,c(0,fra),lty=3,lwd=3,col=2) abline(h=0.5,lty=3) abline(v=22,lty=3) text(22,0,'22') abline(h=0.25,lty=3) abline(v=18.3,lty=3) text(18.3,0,'18.3') abline(h=0.75,lty=3) abline(v=26,lty=3) text(26,0,'26') #slide 37 summary(peso$bwt) #slide 38 par(mfrow=c(1,2)) hist(peso$bwt,xlab='peso (g)',ylab='Frequência absoluta',main='',col=grey(.5)) boxplot(peso$bwt,range=0,ylab='peso (g)',col=grey(.5)) #slide 54 m=round(by(peso$bwt,peso$smoke,mean),1) s=round(by(peso$bwt,peso$smoke,sd),1) cv=round(s/m*100) #slide 61 z1=(7-5)/2 z2=(8-6)/4