#Mudar Diretório
setwd("C:/Users/anderson/Documents")

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

summary(dados)

###Distribuição de Frequencia###################################################
##Distribuição de frequencia Grau de Instrução.
fa=table(GI) ##frequencia absoluta
fa
fr=fa/sum(fa) ##frequencia relativa
fr
fp=100*fr ##fp=frequencia percentual
fp
dist=cbind(fa,fr,fp) ##distribuição de frequencias
dist

fa=table(SEXO,AT) ##frequencia absoluta
fa
fr=fa/sum(fa) ##frequencia relativ
fr
fp=100*fr ##fp=frequencia percentual
fp

fa=table(SEXO,AT,EC,GI) ##frequencia absoluta
fa
fr=fa/sum(fa) ##frequencia relativ
fr
fp=100*fr ##fp=frequencia percentual
fp


fa=table(NFD) ##frequencia absoluta
fa
fr=fa/sum(fa) ##frequencia relativa
fr
fp=100*fr ##fp=frequencia percentual
fp
dist=cbind(fa,fr,fp) ##distribuição de frequencias
dist


require(fdth)
x=fdt(IDADE)
x
x1=fdt(IDADE,k=4)
x1
x2=fdt(IDADE,start=20,end=50,h=5)
x2
x3=fdt(IDADE,breaks="Scott")
x3

###Graficos#################################################################
##Grafico de barras
fa=table(SEXO) #frequência absoluta
barplot(fa) #plota gráfico.
plot(SEXO) ##plota gráfico

##Grafico de barras
fa=table(SEXO)
barplot(fa,col="red")

plot(SEXO,main="Funcionários",xlab="Sexo",ylab="frequencia absoluta",xlim=c(0,20)
,col="red",horiz=TRUE)

plot(EC,horiz=TRUE)

##Grafico de pizza
fa=table(EC)
fr=fa/sum(fa)
n=length(fr) ##contar quantos classes
pie(fr,main="Estado Civil",col=rainbow(n))

##Gráfico Dispersao
plot(A,P,xlab="Altura",ylab="Peso",col="red")
plot(A,P,xlab="Altura",ylab="Peso",col="red",pch=2)

#histogram para variavel peso.
hist(P,main="Histograma de Peso",xlab="Peso",ylab="Frequencia")

h=hist(P,main="Histograma de Peso",xlab="Peso",ylab="Frequencia")
##poligono de frequência com histograma
lines(c(min(h$breaks), h$mids, max(h$breaks)), c(0,h$counts, 0), type = "l")
##poligono de frequência
plot(c(min(h$breaks), h$mids, max(h$breaks)), c(0,h$counts, 0),
type = "l",main="Polígono de frequencia de Peso",xlab="Peso",ylab="Frequencia")

##histograma pela densidade de frequência
hist(P,freq=FALSE)
h1=density(P) ##obtendo a densidade dos dados
lines(h1) ##adicinar um gráfico de linhas

#Boxplot para variável peso.
boxplot(P,main="Peso")

#Boxplot para variável altura dentro de cada sexo.
boxplot(A~SEXO,main="Altura",xlab="Sexo")

##Divir janele para gráficos
par(mfrow=c(1,2)) ##uma linha e duas colunas
plot(SEXO,main="Sexo")
plot(GI,main="Grau de Instrução")


###Medidas de Posição###################################################################
##Média de Peso
mean(P) 
##Média de Altura por Sexo
tapply(A, SEXO, median)

##Mediana de Peso
mean(P) 
##Mediana de Altura por Sexo
tapply(A, SEXO, median)

require(modeest)
##Moda de Idade
mfv(IDADE)
##Moda de Idade por Sexo
tapply(IDADE, SEXO,mfv)

##moda de Peso
mlv(P)
mlv(P, method = "lientz", bw = 0.2)
mlv(P, method = "naive", bw = 1/3)
mlv(P, method = "venter")
mlv(P, method = "grenander", p=4)
mlv(P, method = "hrm", bw = 0.3)
mlv(P, method = "hsm")
mlv(P, method = "parzen")
mlv(P, method = "tsybakov")
mlv(P, method = "asselin")
mlv(P, method = "vieu")


###Medida de Dispersão########################################################################
##Amplitude de Altura
range(A)
diff(range(A))
X=tapply(A, SEXO,range)
X
diff(X$F) ##Ampliture Altura Mulheres
diff(X$M) ##Ampliture Altura Homens

##Variância de Altura
var(A) 
##Variância de Altura por Sexo
tapply(A, SEXO,var)
##Desvio Padrão de Altura
sd(A) 
##Desvio Padrão de Altura por Sexo
tapply(A, SEXO,sd)