rm(list=ls()) pdf(file="peso-altura-idade-sexo-criancas.pdf",width=12,height=10) #################################################################### # Alunos de econometria 2014-1 # Altura x peso x sexo #################################################################### altura = c(182,184,194,178,179,196,163,163,172,170,168,180,195,164,169, 172,179,160,164,158,171,176,163,181,179,172,163,164) peso = c(75,80,120,62,72,96,56,50,71,65,54,53,140,60,64,59,82,57,63,60, 61,78,54,75,82,72,52,63) sexo = c(1,1,1,1,1,1,0,0,0,0,0,0,1,0,1,0,1,0,0,0,0,1,0,1,1,1,0,0) cols = c(2,4) par(mfrow=c(1,2)) reg1 = lm(peso~altura) reg2 = lm(peso~altura+sexo) reg3 = lm(peso~altura+sexo+altura*sexo) plot(altura,peso,col=cols,ylab="Peso (kg)",xlab="Altura (cm)", pch=16,ylim=c(30,140),xlim=range(altura)) abline(reg1$coef,lwd=2) abline(reg2$coef[1],reg2$coef[2],col=2,lwd=2) abline(reg2$coef[1]+reg2$coef[3],reg2$coef[2],col=4,lwd=2) abline(reg3$coef[1],reg3$coef[2],col=2,lty=2,lwd=2) abline(reg3$coef[1]+reg3$coef[3],reg3$coef[2]+reg3$coef[4],col=4,lty=2,lwd=2) title("Com outliers") peso1 = peso[-c(3,6,12,13)] altura1 = altura[-c(3,6,12,13)] sexo1 = sexo[-c(3,6,12,13)] reg1 = lm(peso1~altura1) reg2 = lm(peso1~altura1+sexo1) reg3 = lm(peso1~altura1+sexo1+altura1*sexo1) plot(altura1,peso1,col=cols,ylab="Peso (kg)",xlab="Altura (cm)", pch=16,ylim=c(30,140),xlim=range(altura)) abline(reg1$coef,lwd=2) abline(reg2$coef[1],reg2$coef[2],col=2,lwd=2) abline(reg2$coef[1]+reg2$coef[3],reg2$coef[2],col=4,lwd=2) abline(reg3$coef[1],reg3$coef[2],col=2,lty=2,lwd=2) abline(reg3$coef[1]+reg3$coef[3],reg3$coef[2]+reg3$coef[4],col=4,lty=2,lwd=2) title("Sem outliers") par(mfrow=c(1,1)) plot(altura1,peso1,col=cols,ylab="Peso (kg)",xlab="Altura (cm)", pch=16,ylim=range(peso1),xlim=range(altura1)) abline(reg1$coef,lwd=2) abline(reg2$coef[1],reg2$coef[2],col=2,lwd=2) abline(reg2$coef[1]+reg2$coef[3],reg2$coef[2],col=4,lwd=2) abline(reg3$coef[1],reg3$coef[2],col=2,lty=2,lwd=2) abline(reg3$coef[1]+reg3$coef[3],reg3$coef[2]+reg3$coef[4],col=4,lty=2,lwd=2) #################################################################### # Alunos de econometria 2014-1 (Insper) # Alunos de Business Statistics (Chicago Booth) # Peso x Altura #################################################################### data = matrix(c(144.8,70.5,152.4,84.1,154.9,52.3,160,56.8,165.1,56.8,165.1,56.8,167.6, 53.6,167.6,61.4,170.2,61.4,172.7,72.7,172.7,72.7,172.7,77.3,175.3,68.2,175.3,75,175.3, 80.5,177.8,71.8,177.8,77.3,177.8,79.5,177.8,81.8,177.8,84.1,177.8,86.4,177.8,95.5,180.3, 47.7,180.3,59.1,180.3,79.5,180.3,84.1,182.9,54.5,182.9,63.6,182.9,72.7,182.9,77.3,182.9, 86.4,182.9,88.6,182.9,88.6,182.9,93.2,182.9,98.6,185.4,63.6,185.4,106.8,188,77.7,188,81.8, 188,84.1,188,104.5,190.5,57.3,190.5,59.1,190.5,59.1,193,66.4,193,104.5,195.6,68.2,195.6, 71.4,195.6,88.6,200.7,56.8,203.2,70.5,203.2,93.2,203.2,100.9,205.7,85,205.7,102.3),55,2,byrow=TRUE) height = data[,1] weight = data[,2] n = nrow(data) n1 = length(altura) altura1 = c(height,altura) peso1 = c(weight,peso) escola1 = c(rep(1,n),rep(0,n1)) n2 = length(altura1) plot(altura1+rnorm(n2),peso1+rnorm(n2),pch=16, col=2+escola1,ylab="Peso (kg)",xlab="Altura (cm)") legend(150,140,legend=c("Alunos de econometria 2014-1 (Insper)", "Alunos do Full-MBA (Chicago Booth)"),col=2:3,pch=16,bty="n") abline(lm(peso1~altura1),col=1,lwd=2) abline(lm(peso~altura),col=2,lwd=2) abline(lm(weight~height),col=3,lwd=2) #################################################################### # # Peso, altura e idade de criancas # #################################################################### rm(altura) rm(peso) rm(sexo) data = read.table("peso-altura-idade-sexo-criancas.txt",header=TRUE) n = nrow(data) attach(data) par(mfrow=c(2,3)) boxplot(peso[sexo==0],peso[sexo==1],names=c("Feminino","Masculino"), ylab="Peso (kg)") boxplot(altura[sexo==0],altura[sexo==1],names=c("Feminino","Masculino"), ylab="Altura (cm)") boxplot(idade[sexo==0],idade[sexo==1],names=c("Feminino","Masculino"), ylab="Idade (anos)") plot(altura,peso,xlab="Altura (cm)",ylab="Peso (kg)",pch=16,col=sexo+1) abline(lm(peso[sexo==0]~altura[sexo==0])$coef,lwd=2) abline(lm(peso[sexo==1]~altura[sexo==1])$coef,lwd=2,col=2) legend(1.3,79,legend=c("Feminino","Masculino"),col=1:2,pch=16,bty="n") plot(idade,peso,xlab="Idade (anos)",ylab="Peso (kg)",pch=16,col=sexo+1) abline(lm(peso[sexo==0]~idade[sexo==0])$coef,lwd=2) abline(lm(peso[sexo==1]~idade[sexo==1])$coef,lwd=2,col=2) plot(idade,altura,xlab="Idade (anos)",ylab="Altura (cm)",pch=16,col=sexo+1) abline(lm(altura[sexo==0]~idade[sexo==0])$coef,lwd=2) abline(lm(altura[sexo==1]~idade[sexo==1])$coef,lwd=2,col=2) dev.off()