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/04/07 19:27] walmes |
disciplinas:ce071-2014-01 [2014/04/17 11:34] walmes |
||
---|---|---|---|
Linha 194: | Linha 194: | ||
* Resíduos studentizados (ou externamente studentizados) | * Resíduos studentizados (ou externamente studentizados) | ||
<latex> | <latex> | ||
- | r_i = \dfrac{e_i}{s(e_i)} = \dfrac{e_i}{\hat{\sigma}_{-i}\sqrt{1-h_{i}}}\\ | + | t_i = \dfrac{e_i}{s(e_i)} = \dfrac{e_i}{\hat{\sigma}_{-i}\sqrt{1-h_{i}}}\\ |
- | \hat{\sigma}_{-i}^2 = \dfrac{(n-p)\hat{\sigma}^2-\frac{e_i}{1-h_{i}}}{(n-1)-p} | + | \hat{\sigma}_{-i}^2 = \dfrac{(n-p)\hat{\sigma}^2-\frac{e_i^2}{1-h_{i}}}{(n-1)-p} |
</latex> | </latex> | ||
Linha 206: | Linha 206: | ||
* DFfits | * DFfits | ||
<latex> | <latex> | ||
- | dffits_i = \dfrac{\hat{y}_i-\hat{y}_{i(-i))}}{hat{\sigma}_{-i}\sqrt{1-h_i}} = \left( \dfrac{p\cdot D_i \cdot\hat{\sigma}^2}{\hat{\sigma}^2_{-i}} \right )^{1/2} | + | dffits_i = \dfrac{\hat{y}_i-\hat{y}_{i(-i))}}{\hat{\sigma}_{-i}\sqrt{h_i}} = t_i\left( \dfrac{h_i}{1-h_i} \right )^{1/2} |
</latex> | </latex> | ||
Linha 214: | Linha 214: | ||
\hat{\beta}_{-i} = \hat{\beta}-\dfrac{e_i}{1-h_i}\cdot (X^\top X)^{-1} x_i | \hat{\beta}_{-i} = \hat{\beta}-\dfrac{e_i}{1-h_i}\cdot (X^\top X)^{-1} x_i | ||
</latex> | </latex> | ||
- | |||
- | |||
- | </code> | ||
<code R> | <code R> | ||
##----------------------------------------------------------------------------- | ##----------------------------------------------------------------------------- | ||
+ | ## Definições da sessão. | ||
+ | require(lattice) | ||
+ | require(latticeExtra) | ||
require(car) | require(car) | ||
+ | |||
+ | ## Função vif da página da Professora Dra Sueli Giolo. | ||
+ | source("http://people.ufpr.br/~giolo/CE071/Exemplos/vif.R") | ||
##----------------------------------------------------------------------------- | ##----------------------------------------------------------------------------- | ||
- | ## Dados sobre o preço de leitão (price) de relógios antigos (do avô) em | + | ## Dados de gasto com alimentação (foodcons) em função da receita da |
- | ## função da idade do relógio (age) e do número de potenciais | + | ## família (income) e do número de membros (size). |
- | ## compradores (bidders). | + | ## |
+ | ## Ajuste um modelo com foodcons~income+size. Faça uma análise | ||
+ | ## completa do modelo e verifique se há necessidade de modificações. Se | ||
+ | ## sim, proceda e justifique. | ||
da <- | da <- | ||
- | read.table("http://www.leg.ufpr.br/~walmes/data/business_economics_dataset/EXAMPLES/GFCLOCKA.DAT", | + | read.table("http://www.leg.ufpr.br/~walmes/data/business_economics_dataset/EXERCISE/DCFOOD.DAT") |
- | header=FALSE) | + | names(da) <- c("house","foodcons","income","size") |
- | + | ||
- | str(da) | + | |
- | + | ||
- | ## Essa coluna removida é o produto de age*bidders | + | |
- | da <- da[,-4] | + | |
- | + | ||
- | names(da) <- c("age", "bidders", "price") | + | |
str(da) | str(da) | ||
##----------------------------------------------------------------------------- | ##----------------------------------------------------------------------------- | ||
- | ## Dados sobre o salário em função dos anos de experiência de uma | + | ## Engenheiros mediram o comprimento (length, cm) e o peso (weight, g) e |
- | ## amostra de trabalhadores sociais. | + | ## o nível de DDT (ddt, ppm) para 144 peixes capturados. Além do mais, a |
+ | ## distância de captura rio acima (mile), o rio (river) e a espécie do | ||
+ | ## peixe (species) também foram registrados. | ||
+ | ## | ||
+ | ## Ajuste um modelo com ddt~mile+length+weight. Faça uma análise | ||
+ | ## completa do modelo e verifique se há necessidade de modificações. Se | ||
+ | ## sim, proceda e justifique. | ||
db <- | db <- | ||
- | read.table("http://www.leg.ufpr.br/~walmes/data/business_economics_dataset/EXAMPLES/SOCWORK.DAT", | + | read.table("http://www.leg.ufpr.br/~walmes/data/business_economics_dataset/EXERCISE/DDT.DAT") |
- | header=FALSE) | + | names(db) <- c("river","mile","species","length","weight","ddt") |
- | + | str(db) | |
- | names(db) <- c("yrsexp", "salary") | + | |
##----------------------------------------------------------------------------- | ##----------------------------------------------------------------------------- |