Diferenças

Aqui você vê as diferenças entre duas revisões dessa página.

--- 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>
 ## listas
@@ Linha 37: / Linha 39: @@
 ## ou acessada do objeto de summary!!
 regS$sigma
+## apagando **tudo** (todos os objetos) da área de trabalho
+rm(list(ls())
 ## Exemplo de Verossimilhança
@@ Linha 56: / Linha 61: @@
 fp <- function(par) sapply(par, function(p) sum(log(p) + dg*log(1-p)))
 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>