Duvida - Ellipse
Preciso construir elipses de confiança para um conjunto de dados bidimensionais. Mas não quero assumir normalidade. Conheço as posteriores condicionais dos parâmetros, mas meu interesse e em uma variável defendida como a multiplicação da raiz quadrada dos valores singulares pelas coordenadas dos autovetores. Se eu fizer esta multiplicação nas cadeias consigo gerar uma amostra a posteriori, mas não conheço a distribuição conjunta destas variáveis. Desta forma tenho que construir elipses baseadas em agrupamento de pontos. E não sei como construí- las. Logo se alguém souber como posso fazer peço que me ajude. Carlos Pereira da Silva Bacharel em Estatística – UEPB Mestrando em Estatística e Experimentação Agropecuária - UFLA Lavras - MGEmail: ccpsilva81@hotmail.com __________________________________________ Tenha Um Dia Abençoado Por DEUS!!! Sl – 34, 8
Se tivesse a posteriori de um parâmetro dava para fazer um histograma/densityplot e pegar intervalos com a função quantile. Como o caso é bidimensional, tem-se que fazer uma densidade kernel em duas dimensões. Para isso veja o help da função::kde2d() ( http://stat.ethz.ch/R-manual/R-devel/library/MASS/html/kde2d.html). Outra opção é usar KernSmooth::bkde2D(). Inclusive, dá para fazer os gráficos com contornos partindo desse código parcial da mensagem. https://stat.ethz.ch/pipermail/r-help/2008-August/170377.html À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 skype: walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
Ou, para considerar correcao de borda, usar kernel2d() do pacote 'splacs'. Particularmente util para espaco parametrico limitado com estimativas proximas 'a borda. On 11/04/2013 05:19 PM, walmes . wrote:
Se tivesse a posteriori de um parâmetro dava para fazer um histograma/densityplot e pegar intervalos com a função quantile. Como o caso é bidimensional, tem-se que fazer uma densidade kernel em duas dimensões. Para isso veja o help da função::kde2d() (http://stat.ethz.ch/R-manual/R-devel/library/MASS/html/kde2d.html). Outra opção é usar KernSmooth::bkde2D(). Inclusive, dá para fazer os gráficos com contornos partindo desse código parcial da mensagem.
https://stat.ethz.ch/pipermail/r-help/2008-August/170377.html
À disposição. Walmes.
========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 skype: walmeszeviani homepage: http://www.leg.ufpr.br/~walmes <http://www.leg.ufpr.br/%7Ewalmes> linux user number: 531218 ==========================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
participantes (3)
-
Carlos Pereira -
Elias T Krainski -
walmes .