require(xlsx) 
dados=read.xlsx("peso.xls",1)
attach(dados)
names(dados)

###Distribuições Amostrais######################################

##Distribuição Normal
##densidade no ponto 1 de uma normal padrão
dnorm(1,2,3)
##probabilidade x<1 de uma normal padrão
pnorm(1,2,3)
##quantil de deixa 95% de probabilidade abaixo dele, de uma normal padrão
qnorm(0.95,2,3)
##gerar uma amostra de 10 elementos, de uma normal padrão
rnorm(10,2,3)
##densidade no ponto 1,  com média 2 e desvio padrão 3
dnorm(1,2,3)
##probabilidade x<1, com média 2 e desvio padrão 3
pnorm(1,2,3)
##quantil de deixa 95% de probabilidade abaixo dele, com média 2 e desvio padrão 3
qnorm(0.95,2,3)
##gerar uma amostra de 10 elementos, com média 2 e desvio padrão 3
rnorm(10,2,3)

##Distribuição t
##densidade no ponto 2, com grau de liberdade 12
dt(1,12)
##probabilidade x<2, com grau de liberdade 12
pt(2,12)
##quantil de deixa 95% de probabilidade abaixo dele, com grau de liberdade 12
qt(0.95,12)
##gerar uma amostra de 10 elementos, com grau de liberdade 12
rt(10,12)

##Distribuição Qui-quadrado
##densidade no ponto  5 com grau de liberdade 15
dchisq(1,12)
##probabilidade x<2 com grau de liberdade 15
pchisq(2,12)
##quantil de deixa 95% de probabilidade abaixo dele, com grau de liberdade 15
qchisq(0.95,15)
##gerar uma amostra de 10 elementos, com grau de liberdade 15
rchisq(10,15)

##Distribuição F
##densidade no ponto 5, com grau de liberdade 12 e 20
df(5,12,20)
##probabilidade x<5, com grau de liberdade 12 e 20
pf(5,12,20)
##quantil de deixa 95% de probabilidade abaixo dele, com grau de liberdade 12 e 20
qf(0.95,12,20)
##gerar uma amostra de 10 elementos,  com grau de liberdade 12 e 20
rchisq(10,12,20)

###Intervalo de Confiança e Teste de Hipótese########################

##Intervalo de confiança para proporção de homens
fa=table(SEXO)
fr=fa/sum(fa)
fr
p=fr[2] ##proporção de homens
p
q=1-p
q
n=sum(fa)
p.ic=p+qnorm(c(0.025,0.975))*sqrt(q*p/n)
p.ic

binom.test(12,29, p = 0.5, conf.level = 0.95, alternative ="less")

##Normalidade
##Testar normalidade das variáveis Idade, Peso, Peso1 e Altura
shapiro.test(IDADE)
shapiro.test(P)
shapiro.test(P1)
shapiro.test(A)


##Média
##Intervalo de 95% confiança para altura
XB=mean(A)
XB
Sigma=0.5
n=length(A)
n
XB.ic=XB+qnorm(c(0.025,0.975))*(Sigma/sqrt(n))
XB.ic

t.test(A,mu=1.72)

wilcox.test(IDADE,mu=36)

##Diferença de médias
##Testa diferença de Pesos
t.test(P,P1,m=0,paired=TRUE)

##Testar valores de peso dentro de sexo
##Obter os valores da idade por sexo
Grupos=tapply(P,SEXO,sort) 
Grupos
##valores de idade para mulheres
x=Grupos$F
x
##valores de idade para homens
y=Grupos$M
y
var.test(x,y) ##Testar igualdade variância dos grupos
t.test(x,y,m=0,var.equal=TRUE)

##Obter os valores da idade por sexo
Grupos=tapply(IDADE,SEXO,sort) 
Grupos
##valores de idade para mulheres
x=Grupos$F
x
##valores de idade para homens
y=Grupos$M
y
wilcox.test(x,y,m=0)