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 | Última revisão Ambos lados da revisão seguinte | ||
cursos:mct:sessao0 [2011/05/05 07:57] paulojus |
cursos:mct:sessao0 [2011/05/09 10:58] paulojus |
||
---|---|---|---|
Linha 277: | Linha 277: | ||
plot(1:13, 1:13, col=gray(seq(1,0,len=13)), pch=19, cex=2) | plot(1:13, 1:13, col=gray(seq(1,0,len=13)), pch=19, cex=2) | ||
par(bg="white") | par(bg="white") | ||
+ | </code> | ||
+ | <fs medium>Um breve exemplo de operações matriciais - coeficientes de regressao(modelo) linear | ||
+ | </fs><code R> | ||
+ | require(MASS) | ||
+ | data(hills) | ||
+ | names(hills) | ||
+ | X <- cbind(1, hills$dist) | ||
+ | y <- hills$time | ||
+ | ## Há formas e formas de fazer operações... | ||
+ | ## uma forma não muito adequada: | ||
+ | solve(t(X) %*% X) %*% t(X) %*% y | ||
+ | ## outra melhor: | ||
+ | solve(crossprod(X), crossprod(X, y)) | ||
+ | </code> | ||
+ | |||
+ | Há funções que já implementam modelos estatísticos usuais. | ||
+ | <code R> | ||
with(hills, plot(time~dist)) | with(hills, plot(time~dist)) | ||