library(lme4) out <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd), data = cbpp, family = binomial, contrasts = list(period = "contr.sum")) summary(out) Fixed effects: Estimate Std. Error z value Pr(>|z|) (Intercept) -2.32337 0.22129 -10.499 < 2e-16 *** period1...
Eu tinha uma tarefa de casa para expressar a distribuição binomial negativa como uma família exponencial de distribuições, uma vez que o parâmetro de dispersão era uma constante conhecida. Isso foi bastante fácil, mas eu me perguntava por que eles exigiriam que mantivéssemos esse parâmetro fixo. Eu...
Estou simulando ensaios de Bernoulli com um aleatório entre grupos e, em seguida, encaixo o modelo correspondente em o pacote:logitθ∼N(logitθ0,12)logitθ∼N(logitθ0,12)\text{logit}\, \theta \sim {\cal N}(\text{logit}\, \theta_0, 1^2)lme4 library(lme4) library(data.table) I <- 30 # number of...
Espero que alguém possa ajudar com o que considero uma pergunta relativamente simples, e acho que sei a resposta, mas sem confirmação ela se tornou algo que eu simplesmente não posso ter certeza. Eu tenho alguns dados de contagem como uma variável de resposta e quero medir como essa variável muda...
É possível ajustar um GAMM (Modelo Misto de Aditivo Generalizado) para dados inflados com zero em R? Caso contrário, é possível ajustar um GAM (Modelo Aditivo Generalizado) para dados inflados com zero com uma distribuição binomial negativa ou quase Poisson em R? (Encontrei funções COZIGAM ::...
Eu tenho ~ 1 milhão de pontos de dados. Aqui está o link para o arquivo data.txt. Cada um deles pode ter um valor entre 0 e 145. É um conjunto de dados discreto. Abaixo está o histograma do conjunto de dados. No eixo x é a contagem (0-145) e no eixo y é a densidade. fonte de dados : eu tenho...
Quase todos os livros de texto que discute a aproximação normal para a distribuição binomial menciona a regra de ouro que a aproximação pode ser usado se e . Alguns livros sugerem . A mesma constante geralmente aparece nas discussões sobre quando mesclar células no . Nenhum dos textos que encontrei...
Suponha que eu ajuste uma regressão binomial e obtenha as estimativas pontuais e a matriz de variância-covariância dos coeficientes de regressão. Isso me permitirá obter um IC para a proporção esperada de sucessos em um experimento futuro, , mas preciso de um IC para a proporção observada. Foram...
Eu calculei a proporção de pintos nascidos do número de ovos nascidos em cada ano usando prop.test()em R. Eu vejo que isso me dá a proporção nascida, mas também o intervalo de confiança de 95%, que é o que eu estou procurando. Depois de ler as excelentes informações de outra pergunta neste site...
Para quais distribuições x, exceto beta, a distribuição x-binomial é boa? As distribuições beta e binomial são famosamente conjugadas, mas estou curioso para saber se outras distribuições não conjugadas fornecerão pmfs compostos comparativamente simples. Por legal, quero dizer que é agradável...
Tentando replicar os resultados do artigo publicado recentemente, Aghion, Philippe, John Van Reenen e Luigi Zingales. 2013. "Inovação e Propriedade Institucional". American Economic Review, 103 (1): 277-304. (O código de dados e dados está disponível em
O seguinte é um problema de Bayesian Data Analysis 2nd ed , p. 97. Andrew Gelman não incluiu sua solução no guia em seu site e isso me deixou louco o dia todo. Literalmente o dia todo. yyyNNNθθ\thetaNNNPr(N|μ)=Poisson(μ)Pr(N|μ)=Poisson(μ)\Pr(N|\mu) = Poisson(\mu)μμ\mu(N,θ)(N,θ)(N,...
Estamos testando uma campanha de marketing por email. Em nosso teste inicial, enviamos dois tipos de email diferentes e tivemos um terceiro grupo de controle que não recebeu um email. Agora estamos recuperando "resultados" como proporção de usuários que retornaram ao nosso aplicativo. Aqui estão os...
Situação: Digamos que eu tenha um processo de Poisson, como decaimento radioativo, produzindo partículas R por segundo. Eu meço com um detector. Existe uma probabilidade P de que uma partícula seja detectada pelo detector. Coisas que acho que sei: O tempo entre a chegada da emissão de...
Alguns de vocês podem ter lido este belo artigo: O'Hara RB, Kotze DJ (2010) Não faça log de dados de contagem de transformação. Métodos em Ecologia e Evolução 1: 118–122. Klick . Atualmente, estou comparando modelos binomiais negativos com modelos gaussianos em dados transformados. Ao contrário...
Antes de mais, gostaria de precisar que não sou especialista no assunto. Suponha que duas variáveis aleatórias e sejam binomiais, respectivamente e observe aqui que é o mesmo. Eu sei queXXXYYYX∼B(n1,p)X∼B(n1,p)X\sim B(n_1,p)Y∼B(n2,p),Y∼B(n2,p),Y\sim
Eu quero ajustar um modelo linear por R com family=binomial(link="identity"), no entanto, a família binomial não tem vínculo de identidade. O que devo
Eu tenho uma amostra de uma distribuição discreta, como tal: Type: 0 1 2 3 4 5 Occurrences: 88 12 52 43 21 5 Minha tarefa é testar se uma Distribuição Binomial (n = 5, p) se encaixa nesses dados. Entendo que devo usar o teste de hipóteses e que o teste do qui-quadrado é um grampo para esse...
Suponha que eu saiba gerar variáveis aleatórias binomiais independentes. Como posso gerar duas variáveis aleatórias e modo que X \ sim \ text {Bin} (8, \ dfrac {2} {3}), \ quad Y \ sim \ text {Bin} (18, \ dfrac {2 } {3}) \ \ text {e} \ \ text {Corr} (X, Y) =
Imagine uma eleição em que pessoas façam uma escolha binária: votam a favor ou contra. O resultado é que m pessoas votam em A e, portanto, o resultado de A é p = m / n .nnnmmmp=m/np=m/np=m/n Se eu quiser modelar estas eleições, eu posso assumir que cada pessoa vota para Um independentemente com...