Diferenças

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

--- roteiros:veror [2013/10/21 16:07] – adalardo
+++ roteiros:veror [2024/01/12 10:40] (atual) – edição externa 127.0.0.1
@@ Linha 98: / Linha 98: @@
 	{
 	p<-seqp[j]
-	resulta[i,j]<- -1*sum(log(veropsi(psi, p, nocor, nt)))
+	resulta[i,j]<- -1*sum(log(verosim(psi, p, nocor, nt)))
 	}
 }
@@ Linha 108: / Linha 108: @@
 psiEst<-dimnames(resulta)[[1]][[indmin[1]]]
 pEst<-dimnames(resulta)[[2]][[indmin[2]]]
-pEst
+psiEst; pEst
 ######################################################
+## O gráfico da superfície de Verossimilhança
+contour(resulta, xlab="Ocorrência (psi)", ylab="Detecção (p)", nlevels=20)
+abline(v=as.numeric(psiEst), col="red", lty=2)
+abline(h=as.numeric(pEst), col="red", lty=2)
+########################################################
+# Um conjunto de dados que não conhecemos os parâmetros
+########################################################
+## dados de Ploeophana em 10 árvores
+####################################
+ploe<-rbind(c(	0,	0,	0,	0,	0,	0),
+c(	0,	0,	0,	0,	0,	1),
+c(	0,	0,	0,	0,	0,	0),
+c(	0,	0,	0,	0,	0,	0),
+c(	0,	0,	0,	0,	0,	1),
+c(	0,	1,	1,	0,	0,	0),
+c(	0,	1,	1,	1,	0,	0),
+c(	0,	0,	0,	0,	0,	0),
+c(	0,	0,	0,	1,	0,	0),
+c(	1,	0,	1,	0,	1,	0),
+c(	1,	1,	1,	1,	1,	1))
+#####################################################################
+# calcule os parâmetros e faça o gráfico da função de verossimilhança
+#####################################################################
 </code>