aqui vai um exemplinho especifico da normal considerando a verossimilhanca
para ambos parametros:
dados <- rnorm(50)
mu.vals <- seq(-3,3, len=100)
sigma.vals <- seq(0.5, 5, l=100)
logV <- function(pars, y)
sum(dnorm(y, mean=pars[1] , sd=pars[2], log=TRUE))
ll <- apply(expand. grid(mu.vals, sigma.vals), 1, logV, y=dados)
range(ll)
contour(mu.vals, sigma.vals, matrix(ll, nc=length(sigma. vals)),
levels=seq(- 200, -70, by=10))
Paulo Justiniano Ribeiro Jr
LEG (Laboratorio de Estatistica e Geoinformacao)
Universidade Federal do Parana
Caixa Postal 19.081
CEP 81.531-990
Curitiba, PR - Brasil
Tel: (+55) 41 3361 3573
Fax: (+55) 41 3361 3141
e-mail: paulojus AT ufpr br
http://www.leg. ufpr.br/~ paulojus
On Thu, 18 Dec 2008, Daniela Recchia wrote:
> Olá Gledson,
>
> muito obrigada pela ajuda!
>
> Daniela
>
>
>
>
> ____________ _________ _________ __
> De: Gledson Picharski <gledson.picharski@ yahoo.com. br>
> Para: R_STAT@yahoogrupos. com.br
> Enviadas: Quinta-feira, 18 de Dezembro de 2008 0:43:08
> Assunto: Res: [R_STAT] gráfico da função de verossimilhanç a no R
>
>
> Olá,
> veja na curso do Professor Paulo Justiniano.
> http://www.leg. ufpr.br/~ paulojus/ embrapa/Rembrapa /Rembrapase17. html#x18- 10500017
>
>
>
> ____________ _________ __
> Gledson Luiz Picharski
> PET Estatística - UFPR
> http://www.leg. ufpr.br/~ gledson
> http://www.leg. ufpr.br/pet
>
>
>
>
>
> ____________ _________ _________ __
> De: Daniela Recchia <daniela_recchia@ yahoo.com. br>
> Para: stat-math <stat-math@yahoogru p os.com.br>; R_STAT@yahoogrupos. com.br
> Enviadas: Quarta-feira, 17 de Dezembro de 2008 22:22:31
> Assunto: [R_STAT] gráfico da função de verossimilhanç a no R
>
>
> Olá,
>
> estou tentando fazer um gráfico da função de verossimilhanç a, mas não para modelos, um caso mais simples, para distribuições.
>
> Por exemplo, como crio um gráfico para a função de verossimilhanç a da distribuição normal?
>
> Vi que o comando boxcox do MASS gera um gráfico da verossimilhanç a para modelos, porém não para as distribuições. Assim como não mostra qual o valor do parâmetro, só um intervalo.
>
> Alguém saberia como fazer?
>
> Grata,
> Daniela
> ____________ _________ _________ __
> Veja quais são os assuntos do momento no Yahoo! + Buscados: Top 10 - Celebridades - Música - Esportes
> ____________ _________ _________ __
> Veja quais são os assuntos do momento no Yahoo! + Buscados: Top 10 - Celebridades - Música - Esportes
>
>
> Veja quais são os assuntos do momento no Yahoo! +Buscados
> http://br.maisbusca dos.yahoo. com