# Leitura dos dados ob<-read.csv("http://www.leg.ufpr.br/~joel/dados/obesidade.csv") # Nomes das variaveis names(ob) # Resumo das variaveis summary(ob) # Analise da variavel altura (HEIGHT) summary(ob$HEIGHT) # alternativamente with(ob,summary(HEIGHT)) # distribuicao normal explica valores de altura hist(ob$HEIGHT,prob=T) rug(ob$HEIGHT) curve(dnorm(x,mean(ob$HEIGHT,na.rm=T),sd(ob$HEIGHT,na.rm=T)),100,200,add=T,col='red') lines(density(ob$HEIGHT),col='blue') # Analise Bivariada- altura em funcao do sexo #media with(ob,by(HEIGHT,SEX,mean)) #mediana with(ob,by(HEIGHT,SEX,median)) #desvio padrao with(ob,by(HEIGHT,SEX,sd)) #boxplot de altura por sexo with(ob,boxplot(HEIGHT~SEX)) ############################### # ILUSTRACAO TESTE DE HIPOTESES ############################### # Criando duas subpopulacoes com variavel altura # meninos subM<-subset(ob,SEX=="M",select="HEIGHT") # Tamanho populacao meninos N1<-dim(subM)[1] # Media de altura meninos m1<-mean(subM$HEIGHT) # Desvio padrao de altura meninos sigma1<-sd(subM$HEIGHT) # meninas subF<-subset(ob,SEX=="F",select="HEIGHT") # Tamanho da populacao meninas N2<-dim(subF)[1] # Media de altura meninas m2<-mean(subF$HEIGHT) # Desvio padrao de altura meninas sigma2<-sd(subF$HEIGHT) # Escolhendo uma subpopulacao set.seed(123) sexo1<-ifelse(rbinom(1,1,0.5)==1,"M","F") subpop1<-subset(ob,SEX==sexo1,select="HEIGHT") # Retirando uma amostra da subpopulacao amostra1<-sample(subpop1[,1],30) # Resumindo a informacao da amostra summary(amostra1) # Tamanho da amostra 1 n1<-length(amostra1) # Media da amostra 1 xbarra1<-mean(amostra1) # Devio padrão da amostra 1 s1<-sd(amostra1) # Esta amostra veio dos meninos ou meninas ? # Confira resposta correta sexo1 # Outra amostra set.seed(321) sexo2<-ifelse(rbinom(1,1,0.5)==1,"M","F") subpop2<-subset(ob,SEX==sexo2,select="HEIGHT") amostra2<-sample(subpop2[,1],30) # Resumo summary(amostra2) # Tamanho da amostra 2 n2<-length(amostra2) # Media da amostra 2 xbarra2<-mean(amostra2) # Devio padrão da amostra 2 s2<-sd(amostra2) ############################################# ### TESTES PARA UMA AMOSTRA ############################################# #HO:m=m1 (meninos) H1: m