Exercícios recomendados da CE-223 Estatística computacional, 2008

Semana 1

Aula 25/02

Fazer um gráfico da função de probabilidade de uma v.a. <latex>X \sim Bin(n=10, p=0.03)</latex>

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

Fazer um gráfico da função de densidade de probabilidade de uma v.a. <latex>X \sim N(70, 10^2)</latex>
```
curve(dnorm(x, mean=70, sd=10), xlim=c(-3,3))
```

Aula 27/02

Mostrar o comando para obter uma sequência dos múltiplos de 10 até 200.
```
seq(10, 200, by=10)
```
Criar um vetor a1 com os elementos (23, 45, 21, 29, 40, 22, 29, 37, 44, 37, 31, 33, 36)
```
a1 <- c(23, 45, 21, 29, 40, 22, 29, 37, 44, 37, 31, 33, 36)
```
Extrair os elementos de a1 que sejam maiores que 30.
```
a1[a1 > 30]
```
Extrair os elementos de a1 que sejam menores que 25 ou maiores que 40. Guardar estes valores em um vetor a2
```
a2 <- a1[a1 < 25 | a1 > 40]
```
Extrair os elementos de a1 que sejam maiores que 30 e menores que 40.
```
a1[a1 > 30 & a1 < 40]
```
Obter as posições dos elementos de a1 que sejam menores que 30
```
which(a1[a1 < 30])
```
Obter a posição do maior elemento da a1
```
wihch.max(a1)
```
Obter a posição do menor elemento da a1
```
    which.min(a1)
```
Criar um vetor a3 com os elementos de a1 para os quais o resto da divisão por 3 seja igual a 2.
```
a3 <- a1[a1 %% 3 == 2]
```
Extrair os elementos de a1 que sejam múltiplos de 4
```
a1[a1 %% 4 == 0]
```
Substituir em a1 os elementos iguais a 37 pelo valor 36
```
a1[a1 == 37] <- 36
```
Substituir em a1 os elementos maiores que 40 pelo código de valor perdido NA
```
a1[a1 > 40] <- NA   
```
Obter as posições de a1 onde estão os valores perdidos
```
which(is.na(a1))
```

Crie um vetor chamado sexo com os comandos a seguir:

sexo <- c(1, 2, 2, 2, 1, 2, 1, 1, 1, 2, 2, 1, 2)
sexo <- factor(sexo, lev=1:2, lab=c("M","F"))

Obter as posições em sexo que possuem o valor “M“
```
which(sexo == "M")
```

Obter os valores de a1 para os quais o valor correspondente em sexo é “M“

a1[which(sexo == "M")]
<code>
    - Obter os valores de ''a1'' para os quais o valor correspondente em ''sexo'' é ''“F“''<code R>
a1[which(sexo == "F")]

Descrever o resultado de cada um dos comandos a seguir:

sort(a1)             # Ordena o vetor a1 em ordem crescente
order(a1)            # Retorna a ordem de cada elemento no vetor a1
a1[order(a1)]        # Ordena o vetor a1 em ordem crescente
sort(a1, dec = TRUE) # Ordena o vetor a1 em ordem decrescente

Criar um objeto a1.ord com os elementos de a1 em ordem crescente
```
a1.ord <- sort(a1)
```
Ordenar os objetos de a1 de forma a exibir primeiro todos os elementos correpondentes a “M“ e depois os correspondentes a “F“
```
a1[order(sexo)]
```
Criar um objeto chamado notas que possua os elementos de a1 com valores correspondentes de sexo sendo “M“ ordenados de forma crescente, seguidos pelos correspondentes a “F“ também ordenados de forma crescente. Em outras palavras, o objeto notas deverá ter as notas dos homes ordenadas seguidas pelas das mulheres também ordenadas.
```
notas <- a1[order(sexo, a1)]
```