Pacotes Estatísticos
Universidade Federal de Minas Gerais
Instituto de Ciências Exatas
Departamento de Estatística

Disciplina: Pacotes Estatísticos (EST-008)
Prof.: Frederico R. B. Cruz
Período: 1º Semestre de 2022
Sala: 4069 (ICEx-Pampulha)
Divulgação: 23/06/2022 Entrega: 23/06/2022 (em PDF via Moodle)

R - Roteiro de Estudos 4 (Lista 8)

O acompanhamento deste roteiro de estudos pressupõe acesso contínuo ao relatório técnico Noções Básicas de S-PLUS for Windows® [1], e às Notas de Aula Cap. 2. O R [2].

8. escrevendo funções R (continuação)

  • implemente uma função para estimar o vício, E(ˆλ)-λ e a variância, V(ˆλ), do estimador de máxima verossimilhança (EMV) do parâmetro λ da distribuição de Poisson truncada em zero, implementado em roteiro anterior, para um dado tamanho de amostra, via 1000 simulações Monte Carlo (maiores detalhes serão vistos em roteiros futuros); em outras palavras, vocâ deverá implementar o seguinte algoritmo:

    função MonteCarlo(tam)
      para i de 1 até 1000 faça
        gerar dados
        tirar contagens nulas
        ajustar modelo e guardar ajuste
      fim para
      retornar o vício e a variância das estimativas
    fim função

  • teste a função MonteCarlo, para alguns valores de tamanhos de amostra:
    MonteCarlo(10)
    MonteCarlo(100)
    MonteCarlo(1000)

Lista de Exercícios

Valendo como lista de exercícios, envie pelo Moodle um arquivo PDF com:

  • a função MonteCarlo;
  • resultados de testes com a sua função.

Dica:
Para traçar gráficos use os comandos:

#dividir a tela gráfica em duas colunas
par(mfrow=c(1,2))

#fazer as simulações e guardar resultados
res <- matrix(nrow = 3, ncol = 2)
res[1,] <- MonteCarlo(10)
res[2,] <- MonteCarlo(100)
res[3,] <- MonteCarlo(1000)

#criar eixo x e fazer gráficos
x <- c(10,100,1000)

#gráficos da média:
plot(x, res[,1], type="b")

#gráficos da variância:
plot(x, res[,2], type="b")

Referências:

[1] E. A. Reis. Noções Básicas de S-PLUS for Windows®. RTE-03/1997, EST-ICEx-UFMG, Belo Horizonte, 1997 (disponível através da homepage de Relatórios Técnicos - Série Ensino, Departamento de Estatística, UFMG, ou em RTE-03/97).
[2] F. R. B. Cruz. Notas de Aula - O R. EST-ICEx-UFMG, Belo Horizonte, 2015. (disponível em O R).
Endereços:
Atenção:
O material aqui contido deve ser tomado como um guia bastante confiável do conteúdo do curso. Entretanto, você não poderá reivindicar quaisquer direitos baseado neste material. Em particular, reservo-me o direito de alterar datas, programação ou critério de avaliação. Avisos oficiais serão sempre aqueles feitos em sala de aula.

Saturday, 18-May-2024 08:54:55 -03
Last updated: March 04, 2024
Copyright © 2000-2024, F. R. B. Cruz. All Righs Reserved.