############Resoluções dos Exercícios de Estatístic Computacional###############
#####Wagner Hugo Binat - LEG/UFPR 29/02/2008####################################
################################################################################

## Ex 1 Função de probabilidade da Binomial
x <- 0:10
fx <- dbinom(x, 10, 0.35)
plot(x, fx, type='h')
title("Função de Probabilidade da Binomial")

##Ex 2 Função Densidade de Probabilidade da Normal
plot(function(x) dnorm(x, 70, 10),30 , 110, ylab='f(x)')
title("Função Densidade Probabilidade da Normal")

## aula 27/02
## Ex 1 Sequancia dos multiplos de 10 até 200##
resul <- list()
for(i in 1:200){
resul[i] <- ifelse(x[i]%%10 == 0, x[i],0)}
resul <- resul[which(resul !=0)]
resul <- data.frame(resul)
resul

## Ex 2 Criando um vetor
a1 <- c(23,45,21,29,40,22,29,37,44,37,31,33,36)

## Ex 3 maiores que 30
maior.30 <-a1[which(a1 > 30)]

## Ex 4 menores 25 ou maiores de 40
dupla <- a1[which(a1 < 25 | a1 > 40)]
dupla

## Ex 5 maiores de 30 e menores de 40
dupla <- a1[which( a1 > 30 & a1 < 40)]
dupla

## Ex 6 posições onde a1 > 30
posi <-which(a1 < 30)

## Ex 7 posição do maior elemento
posi.max <- which(a1 == max(a1))
posi.max

## Ex 8 posição do minimo elemento
posi.min <- which(a1 ==min(a1))
posi.min

## Ex 9 resto da divisão por 2 é igual a 3
a3 <- a1[which(a1%%3 ==2)]
a3

## Ex 10 extrair os multiplos de 4
a4 <- a1[which(a1%%4 == 0)]
a4

## Ex 11 substituição
a1[which(a1 == 37)] <- 36

## Ex 12 substituição
a1[which(a1 == 40)]  <- "NA"

## Ex 13 posição
which(a1 == "NA")

## Ex 14 vetores
sexo <- c(1,2,2,2,1,2,1,1,1,2,2,1,2)
sexo1 <- factor(sexo,lev=1:2,lab=c("M","F"))
sexo1

## Ex 15 localizando posições
which(sexo1 =="M")

## Ex 16 juntado os vetores, obtendo informações conjuntas
x <- data.frame(sexo1,a1)
x[which(x$sexo1 == "M"),]

##  Ex 17
x[which(x$sexo1 == "F"),]

## Ex 18
sort(a1) #Ordena em ordem crescente
order(a1) # Da posição segundo a ordem de cada valor
a1[order(a1)] # Ordena os valores pelas posições obtidos no comando order
sort(a1,dec=TRUE) # Ordem decrescente

## Ex 19
a1.ord <- sort(a1)
a1.ord

## Ex 20
x[order(x$sexo1),]

## Ex 21
sex <- split(x,x$sexo1)
notas <-rbind(sex[[1]][order(sex[[1]]$a1),],sex[[2]][order(sex[[2]]$a1),])
notas

###############################FIM DO CÓDIGO####################################


















xxxx






xxxxxxxxxx



x





xxxxxxxxxxxxxxxxxxxxxx