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.