Não foi possível enviar o arquivo. Será algum problema com as permissões?
Diferenças
Aqui você vê as diferenças entre duas revisões dessa página.
| Ambos lados da revisão anterior Revisão anterior Próxima revisão | Revisão anterior Próxima revisão Ambos lados da revisão seguinte | ||
|
disciplinas:ce071-2014-01 [2014/03/19 09:57] walmes |
disciplinas:ce071-2014-01 [2014/03/19 10:08] walmes |
||
|---|---|---|---|
| Linha 92: | Linha 92: | ||
| * Função para quadro de anova particionado. | * Função para quadro de anova particionado. | ||
| * Função para calcular o valor predito com IC. | * Função para calcular o valor predito com IC. | ||
| + | * Entregar o código impresso das funções programadas no dia 24/03/14. | ||
| <code R> | <code R> | ||
| Linha 118: | Linha 119: | ||
| ... | ... | ||
| } | } | ||
| + | </code> | ||
| + | |||
| + | === Trabalho 2 === | ||
| + | |||
| + | * Fazer estudo de simulação para estudar a distribuição amostral dos estimadores e das estatísticas do testes. | ||
| + | * Verificar que E(\hat\beta) = \beta, var(\hat\beta) = \sigma^2(X'X)^{-1}, e que \hat\betas têm distribuição Normal. | ||
| + | * Verificar que E(\hat\sigma^2) = \sigma^2 e que (n-p)*\hat\sigma/\sigma têm distribuição qui-quadrado. | ||
| + | * Verificar que F = (A\hat\beta-m)'[A(X'X)^{-1}A']^{-1}(A\hat\beta-m)/(r QMRes) têm distribuição F sob H0 que A\betas = m. | ||
| + | * Estudar a distribuição da estatística F = QMReg/QMres e comparar com o F anterior. | ||
| + | * Entregar código impresso com gráficos e tabelas que sobre os resultados solicitados no dia 24/03/14. | ||
| + | |||
| + | <code R> | ||
| + | ## Função que retorna estimativas de parâmetros e estatísticas sob uma | ||
| + | ## amostra aleatória simulada ao ser executada. | ||
| + | mysimula <- function(X, beta, sigma, A, m=beta){ | ||
| + | ... | ||
| + | } | ||
| + | |||
| + | results <- replicate(10000, mysimula) | ||
| </code> | </code> | ||
