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 Ambos lados da revisão seguinte | ||
disciplinas:ce071-2014-01 [2014/04/03 17:06] walmes |
disciplinas:ce071-2014-01 [2014/04/07 18:36] walmes |
||
---|---|---|---|
Linha 163: | Linha 163: | ||
results <- replicate(10000, mysimula) | results <- replicate(10000, mysimula) | ||
+ | </code> | ||
+ | |||
+ | === Trabalho 3 === | ||
+ | |||
+ | * Programar funções para obter: | ||
+ | * Resíduos ordinários, padronizados e studentizados; | ||
+ | * Valores de alavancagem; | ||
+ | * Distância de Cook; | ||
+ | * DFfits, DFbetas; | ||
+ | * As funções devem receber como argumentos as matrizes X e y e retornas as respectivas medidas; | ||
+ | |||
+ | * Resíduos crus | ||
+ | <latex> | ||
+ | oi | ||
+ | </latex> | ||
+ | |||
+ | |||
</code> | </code> | ||
<code R> | <code R> | ||
- | da <- data.frame(cat=gl(3,3), x=10*(1:9)) | + | ##----------------------------------------------------------------------------- |
- | X <- model.matrix(~x, da); X | + | |
- | X <- model.matrix(~cat, da); X | + | require(car) |
- | X <- model.matrix(~cat+x, da); X | + | |
- | X <- model.matrix(~cat*x, da); X | + | ##----------------------------------------------------------------------------- |
- | X <- model.matrix(~cat:x, da); X | + | ## Dados sobre o preço de leitão (price) de relógios antigos (do avô) em |
- | X <- model.matrix(~cat/x, da); X | + | ## função da idade do relógio (age) e do número de potenciais |
- | X <- model.matrix(~-1+cat/x, da); X | + | ## compradores (bidders). |
da <- | da <- | ||
- | read.table("http://www.leg.ufpr.br/~walmes/data/business_economics_dataset/EXAMPLES/REALESTA.DAT", | + | read.table("http://www.leg.ufpr.br/~walmes/data/business_economics_dataset/EXAMPLES/GFCLOCKA.DAT", |
header=FALSE) | header=FALSE) | ||
- | names(da) <- c("number","saleprice","landvalue","improvvalue","area") | + | |
str(da) | str(da) | ||
+ | ## Essa coluna removida é o produto de age*bidders | ||
+ | da <- da[,-4] | ||
- | </code> | + | names(da) <- c("age", "bidders", "price") |
+ | str(da) | ||
- | <code R> | + | ##----------------------------------------------------------------------------- |
- | du <- | + | ## Dados sobre o salário em função dos anos de experiência de uma |
- | read.table("http://www.leg.ufpr.br/~walmes/data/duster_venda_260314.txt", | + | ## amostra de trabalhadores sociais. |
- | header=TRUE, sep="\t") | + | |
- | str(du) | + | db <- |
+ | read.table("http://www.leg.ufpr.br/~walmes/data/business_economics_dataset/EXAMPLES/SOCWORK.DAT", | ||
+ | header=FALSE) | ||
+ | |||
+ | names(db) <- c("yrsexp", "salary") | ||
+ | ##----------------------------------------------------------------------------- | ||
</code> | </code> | ||