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

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

      curve(pbinom(x, size=10, prob=0.03), xlim=c(0,3))
 

• 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))
 

• 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)]