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 | |||
|
cursos:rbelem:dia3 [2010/05/28 09:16] paulojus |
cursos:rbelem:dia3 [2010/05/29 11:13] (atual) paulojus |
||
|---|---|---|---|
| Linha 1: | Linha 1: | ||
| + | ===== Dia 3 ===== | ||
| + | |||
| <code R> | <code R> | ||
| ## listas | ## listas | ||
| Linha 37: | Linha 39: | ||
| ## ou acessada do objeto de summary!! | ## ou acessada do objeto de summary!! | ||
| regS$sigma | regS$sigma | ||
| + | |||
| + | ## apagando **tudo** (todos os objetos) da área de trabalho | ||
| + | rm(list(ls()) | ||
| ## Exemplo de Verossimilhança | ## Exemplo de Verossimilhança | ||
| Linha 56: | Linha 61: | ||
| fp <- function(par) sapply(par, function(p) sum(log(p) + dg*log(1-p))) | fp <- function(par) sapply(par, function(p) sum(log(p) + dg*log(1-p))) | ||
| plot(ps, fp(ps), ty="l") | plot(ps, fp(ps), ty="l") | ||
| - | ## reescrevando novamente usando a função já disponível no R | + | ## reescrevendo novamente usando a função já disponível no R |
| - | fp <- fp(par) sapply(par, function(p) dgeom(dg, prob=p, log=T)) | + | fp <- function(par) sapply(par, function(p) sum(dgeom(dg, prob=p, log=T))) |
| + | plot(ps, fp(ps), ty="l") | ||
| + | |||
| + | ## a tornando a função mais genérica para receber qualquer vetor de dados | ||
| + | fp <- function(par, dados) sapply(par, function(p) sum(dgeom(dados, prob=p, log=T))) | ||
| + | plot(ps, fp(ps, dados=dg), ty="l") | ||
| + | |||
| + | ## acrescentando elementos no gráfico | ||
| + | plot(ps, fp(ps, dados=dg), ty="l", xlab="p", ylab="l(p)") | ||
| + | title(eval(substitute(expression(paste("log-verossimilhança do parâmetro da Geométrica\n", hat(p) == a)), list(a=est$max)))) | ||
| </code> | </code> | ||
