BUGS é um acrônimo para inferência bayesiana usando a amostragem de Gibbs; O BUGS também é um pacote de software para fazer isso.
Estou prestes a experimentar um ambiente de estilo BUGS para estimar modelos bayesianos. Existem vantagens importantes a serem consideradas na escolha entre OpenBugs ou JAGS? É provável que um substitua o outro em um futuro próximo? Usarei o Gibbs Sampler escolhido com R. Ainda não tenho uma...
Eu encontrei algumas distribuições para as quais BUGS e R têm diferentes parametrizações: Normal, log-Normal e Weibull. Para cada uma delas, entendo que o segundo parâmetro usado por R precisa ser transformado inversamente (1 / parâmetro) antes de ser usado no BUGS (ou JAGS no meu caso). Alguém...
Vou explicar meu problema com um exemplo. Suponha que você queira prever a renda de um indivíduo, com alguns atributos: {Idade, Sexo, País, Região, Cidade}. Você tem um conjunto de dados de treinamento como esse train <- data.frame(CountryID=c(1,1,1,1, 2,2,2,2, 3,3,3,3), RegionID=c(1,1,1,2,...
Olá gurus estatísticos e assistentes de programação R, Estou interessado em modelar capturas de animais em função das condições ambientais e do dia do ano. Como parte de outro estudo, tenho contagens de capturas em ~ 160 dias em três anos. Em cada um desses dias, tenho temperatura, precipitação,...
Esta pergunta é um acompanhamento técnico desta pergunta . Tenho dificuldade em entender e replicar o modelo apresentado em Raftery (1988): Inferência para o binômio NNN parâmetro : uma abordagem hierárquica de Bayes no WinBUGS / OpenBUGS / JAGS. Não se trata apenas de código, portanto, ele deve...
Eu queria saber qual pacote estatístico de software vocês recomendam para executar a inferência bayesiana. Por exemplo, eu sei que você pode executar o openBUGS ou o winBUGS como autônomo ou também pode chamá-los de R. Mas o R também possui vários de seus próprios pacotes (MCMCPack, BACCO), que...
Fechadas. Esta questão está fora de tópico . No momento, não está aceitando respostas. Deseja melhorar esta pergunta? Atualize a pergunta para que ela esteja no tópico de Validação cruzada. Fechado há 11 meses . Estou seguindo um curso de estatística...
Eu tenho um conjunto de dados de séries temporais no qual estou tentando ajustar um Modelo de Markov oculto (HMM) para estimar o número de estados latentes nos dados. Meu pseudo-código para fazer isso é o seguinte: for( i in 2 : max_number_of_states ){ ... calculate HMM with i states ......
Agora, existem várias abordagens diferentes para realizar uma meta-análise de rede ou comparação de tratamentos mistos. Os mais usados e acessíveis são provavelmente os seguintes: em uma estrutura bayesiana : abordagem de interação design por tratamento no WinBUGS (por exemplo, Jackson et al...
Parece que condicionais completos geralmente são bastante difíceis de derivar, mas programas como JAGS e BUGS os derivam automaticamente. Alguém pode explicar como eles geram algoritmos condicionais completos para qualquer especificação arbitrária de
Gelman e Hill (2006) dizem: No Bugs, os resultados ausentes em uma regressão podem ser manipulados facilmente, simplesmente incluindo o vetor de dados, NAs e tudo. Os bugs modelam explicitamente a variável de resultado e, portanto, é trivial usar esse modelo para, de fato, atribuir valores...
Eu sou um novo usuário do WinBUGS e tenho uma pergunta para sua ajuda. Depois de executar o código a seguir, obtive parâmetros de beta0through beta4(estatísticas, densidade), mas não sei como obter a previsão do último valor de h, que defini NAcomo modelo no código. Alguém pode me dar uma dica?...
O que acontece quando você não tem uma idéia da distribuição de parâmetros? Que abordagem devemos usar? Na maioria das vezes, pretendemos subestimar se uma determinada variável tem alguma influência sobre a presença / ausência de uma determinada espécie, e a variável é aceita ou não, de acordo com...
Estou tentando construir um modelo em que a resposta seja proporcional (na verdade, é a parcela de votos que um partido obtém nos círculos eleitorais). Sua distribuição não é normal, então decidi modelá-la com uma distribuição beta. Eu também tenho vários preditores. No entanto, não sei como...
Estou tentando ajustar um modelo de tempo discreto no R, mas não sei como fazê-lo. Eu li que você pode organizar a variável dependente em linhas diferentes, uma para cada observação no tempo, e usar a glmfunção com um link logit ou cloglog. Neste sentido, tem três colunas: ID, Event(1 ou 0, em...
No Rque pudermos "peso antes" uma glmregressão via pesos parâmetro. Por exemplo: glm.D93 <- glm(counts ~ outcome + treatment, family = poisson(), weights=w) Como isso pode ser feito em um JAGSou BUGSmodelo? Encontrei algum artigo discutindo isso, mas nenhum deles fornece um exemplo. Estou...
Em R, c (3,1,0) * c (2,0,1) == c (6,0,0). Este não é um produto escalar e não é um produto cruzado. Primeiro, qual é o nome desse produto e, segundo, ele funciona no WinBUGS, OpenBUGS e / ou
Fechadas. Esta questão está fora de tópico . No momento, não está aceitando respostas. Deseja melhorar esta pergunta? Atualize a pergunta para que ela esteja no tópico de Validação cruzada. Fechado há 2 anos . R parece ser capaz de gerar gráficos de...
Estou interessado em montar uma ANOVA Bayesiana de Dois Fatores em ERROS ou utilizando algum pacote R. Infelizmente, estou tendo dificuldade em encontrar recursos sobre este tópico. Alguma sugestão? Mesmo um artigo descrevendo a abordagem seria