roteiros:veror
Diferenças
Aqui você vê as diferenças entre duas revisões dessa página.
Próxima revisão | Revisão anterior | ||
roteiros:veror [2013/10/21 16:04] – criada adalardo | roteiros:veror [2024/01/12 10:40] (atual) – edição externa 127.0.0.1 | ||
---|---|---|---|
Linha 10: | Linha 10: | ||
ii) Representar graficamente a superfície de verosimilhança.\\ | ii) Representar graficamente a superfície de verosimilhança.\\ | ||
- | 1) **Construa uma base de dados fictícia**. __Dicas de comandos__: '' | + | 1) **Construa uma base de dados fictícia**. __Dicas de comandos__: |
* Esta base de dados pode ser compreendida como uma matriz [i,j], com número i de linhas e j de colunas.\\ | * Esta base de dados pode ser compreendida como uma matriz [i,j], com número i de linhas e j de colunas.\\ | ||
* Como tratam-se de dados de presença e ausência, vamos representar as ausências por zeros, e as presenças por um.\\ | * Como tratam-se de dados de presença e ausência, vamos representar as ausências por zeros, e as presenças por um.\\ | ||
Linha 40: | Linha 40: | ||
8) Adicione no gráfico o **ponto de coordenadas (__psi__, | 8) Adicione no gráfico o **ponto de coordenadas (__psi__, | ||
+ | |||
+ | ===== Uma solução ===== | ||
+ | |||
+ | < | ||
+ | ################################# | ||
+ | ## criando os dados #### | ||
+ | ## com parâmetros conhecidos #### | ||
+ | ################################# | ||
+ | ## Alexandre Adalardo de Oliveira | ||
+ | ## 21 de outubro de 2013 | ||
+ | # processo de interesse: ocorrência | ||
+ | # prob de ocorrenciam psi=0.7 | ||
+ | ocorre< | ||
+ | #ou | ||
+ | # | ||
+ | ########################## | ||
+ | ## processo de observação | ||
+ | ######################### | ||
+ | # prob de detecção p=0.4 | ||
+ | ######################## | ||
+ | # 10 sessões de amostra em cada localidade | ||
+ | observa= sample(c(1, | ||
+ | detecta< | ||
+ | ## conectando os dois processos | ||
+ | detecta[ocorre==0, | ||
+ | amostra< | ||
+ | ############################ | ||
+ | ## nome das colunas e linhas | ||
+ | rownames(amostra)< | ||
+ | colnames(amostra)< | ||
+ | ############################### | ||
+ | ## Calculando as probabilidade | ||
+ | ## para diferentes parâmetros | ||
+ | ############################### | ||
+ | nocor< | ||
+ | nt<- dim(amostra)[2] ## número de visitas a cada localidade | ||
+ | ## sequencia de valores de psi e p | ||
+ | seqpsi< | ||
+ | ## criando uma matrix para guardar os resultados | ||
+ | resulta< | ||
+ | ## um função para o cálculo de verossimilhança para valores de p e psi | ||
+ | verosim< | ||
+ | { | ||
+ | n0 = nocor==0 | ||
+ | res< | ||
+ | res[n0]< | ||
+ | return(res) | ||
+ | } | ||
+ | ############################################################# | ||
+ | ## calculando a Verossimilhança para a sequência de valores | ||
+ | ## de p e psi | ||
+ | ############################################################# | ||
+ | for (i in 1: | ||
+ | { | ||
+ | psi< | ||
+ | for (j in 1: | ||
+ | { | ||
+ | p< | ||
+ | resulta[i, | ||
+ | } | ||
+ | } | ||
+ | ####################################################### | ||
+ | ## encontrando a estimativa de Maxima Verossimilhança | ||
+ | ## ou minima -LogVerossimilhança | ||
+ | ####################################################### | ||
+ | indmin< | ||
+ | psiEst< | ||
+ | pEst< | ||
+ | psiEst; pEst | ||
+ | ###################################################### | ||
+ | ## O gráfico da superfície de Verossimilhança | ||
+ | contour(resulta, | ||
+ | abline(v=as.numeric(psiEst), | ||
+ | abline(h=as.numeric(pEst), | ||
+ | ######################################################## | ||
+ | # Um conjunto de dados que não conhecemos os parâmetros | ||
+ | ######################################################## | ||
+ | ## dados de Ploeophana em 10 árvores | ||
+ | #################################### | ||
+ | ploe< | ||
+ | c( 0, | ||
+ | c( 0, | ||
+ | c( 0, | ||
+ | c( 0, | ||
+ | c( 0, | ||
+ | c( 0, | ||
+ | c( 0, | ||
+ | c( 0, | ||
+ | c( 1, | ||
+ | c( 1, | ||
+ | ##################################################################### | ||
+ | # calcule os parâmetros e faça o gráfico da função de verossimilhança | ||
+ | ##################################################################### | ||
+ | </ |
roteiros/veror.1382371478.txt.gz · Última modificação: 2024/01/12 10:39 (edição externa)