Corri modelos binomiais negativos e quase-Poisson com base em uma abordagem de teste de hipóteses. Meus modelos finais, usando os dois métodos, têm covariáveis e interações diferentes. Parece que não há padrões quando planto meus resíduos nos dois casos. Assim, eu queria saber qual teste eu...
Estou tentando aprender algumas estatísticas usando o livro Biometry by Sokal and Rohlf (3e). Este é um exercício no quinto capítulo que aborda a probabilidade, a distribuição binomial e a distribuição de Poisson. Percebo que existe uma fórmula para produzir uma resposta a esta pergunta: No...
Estou tentando fazer testes A / B da maneira bayesiana, como em Programação Probabilística para Hackers e testes A / B bayesianos . Ambos os artigos pressupõem que o tomador de decisão decide qual das variantes é melhor com base apenas na probabilidade de algum critério, por exemplo, ; portanto, é...
Atualmente preso a isso, eu sei que provavelmente deveria usar o desvio médio da distribuição binomial, mas não consigo
Eu me deparei com o problema dos coletores de cupons e estava tentando descobrir uma fórmula para uma generalização. Se houver NNN objetos distintos e que pretende coletar pelo menos kkk cópias de cada de qualquer mmm deles (onde m ≤ Nm≤Nm \le N ), qual é a expectativa de quantos objetos...
No presente trabalho (disponível gratuitamente via PubMed Central), os autores usam regressão binomial negativa para modelar a pontuação em um instrumento de triagem de 10 itens marcou 0-40. Este procedimento assume dados de contagem, o que claramente não é o caso aqui. Gostaria de suas opiniões...
Em uma configuração binomial, a variável aleatória X, que fornece o número de sucessos, é distribuída binomialmente. A proporção da amostra pode então ser calculada como que é o tamanho da sua amostra. Meu livro declara que nXnXn\frac{X}{n}nnn Essa proporção não possui distribuição binomial no...
Eu tenho dados com um pico duplo que estou tentando modelar e há sobreposição suficiente entre os picos que não posso tratá-los independentemente. Um histograma dos dados pode ser algo como isto: Eu criei dois modelos para isso: um usa duas distribuições de Poisson e o outro usa duas...
Meu conjunto de dados é composto pela mortalidade total ou pela sobrevivência de um organismo em três tipos de locais, costeiros, canais intermediários e offshore. Os números na tabela abaixo representam o número de sites. 100% Mortality 100% Survival Inshore 30 31 Midchannel 10 20 Offshore 1...
Tenho quase as mesmas perguntas como esta: Como modelar eficientemente a soma das variáveis aleatórias de Bernoulli? Mas a configuração é bem diferente: P ( X i = 1 ) = p i N p iS=∑i=1,NXiS=∑i=1,NXiS=\sum_{i=1,N}{X_i} , , ~ 20, ~ 0,1P(Xi=1)=piP(Xi=1)=piP(X_{i}=1)=p_iNNNpipip_i Temos os dados...
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...
Os seguintes enxertos são retirados deste artigo . Eu sou novato no bootstrap e estou tentando implementar o bootstrap paramétrico, semiparamétrico e não paramétrico para o modelo misto linear com o R bootpacote. Código R Aqui está o meu
Estou tentando escrever um script R para simular a interpretação de experimentos repetidos de um intervalo de confiança de 95%. Descobri que superestima a proporção de vezes em que o verdadeiro valor populacional de uma proporção está contido no IC de 95% da amostra. Não é uma grande diferença -...
Se são binomiais negativos, então qual é a distribuição de ( x 1 , x 2 , … , x n ) dadax1,x2,…,xnx1,x2,…,xnx_1, x_2, \ldots, x_n(x1,x2,…,xn)(x1,x2,…,xn)(x_1, x_2, \ldots, x_n) x1+x2+…+xn=Nx1+x2+…+xn=Nx_1 + x_2 + \ldots + x_n = N\quad? é fixo.NNN Se são Poisson, então, condicional ao total, ( x...
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 tenho um bloqueio intuitivo com isso. Para um problema binomial, o desvio padrão de uma contagem é . Por outro lado, o desvio padrão da proporção da amostra diminui com o aumento de n e é \ sqrt {\ frac {p (1-p)} {n}} . Eu posso fazer a divisão por n, mas não sinto por que os desvios padrão se...
Eu tenho esses dados: set.seed(1) predictor <- rnorm(20) set.seed(1) counts <- c(sample(1:1000, 20)) df <- data.frame(counts, predictor) Fiz uma regressão de Poisson poisson_counts <- glm(counts ~ predictor, data = df, family = "poisson") E uma regressão binomial...
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...