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 0.92498...
J'ai eu un devoir pour exprimer la distribution binomiale négative comme une famille exponentielle de distributions étant donné que le paramètre de dispersion était une constante connue. C'était assez facile, mais je me demandais pourquoi ils exigeraient que nous maintenions ce paramètre fixe. J'ai...
Je simule des essais de Bernoulli avec un aléatoire entre les groupes, puis j'adapte le modèle correspondant avec le paquet: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...
Est-il possible d'adapter un GAMM (Generalized Additive Mixed Model) pour les données gonflées à zéro dans R? Sinon, est-il possible d'ajuster un GAM (Generalized Additive Model) pour des données gonflées à zéro avec une distribution binomiale négative ou quasi Poisson dans R? (J'ai trouvé COZIGAM...
Presque tous les manuels qui traitent de l'approximation normale de la distribution binomiale mentionnent la règle générale selon laquelle l'approximation peut être utilisée si n p ≥ 5np≥5np\geq5 et n ( 1 - p ) ≥ 5n(1−p)≥5n(1-p)\geq 5. Certains livres suggèrentn p ( 1 - p ) ≥ 5np(1−p)≥5np(1-p)\geq...
Supposons que j'adapte une régression binomiale et que j'obtienne les estimations ponctuelles et la matrice de variance-covariance des coefficients de régression. Cela me permettra d'obtenir un IC pour la proportion attendue de succès dans une future expérience, , mais j'ai besoin d'un IC pour la...
Aujourd'hui, j'ai une question sur la régression binomiale / logistique, basée sur une analyse qu'un groupe de mon département a fait et sur laquelle je sollicitais des commentaires. J'ai composé l'exemple ci-dessous pour protéger leur anonymat, mais ils étaient impatients de voir les réponses....
Je travaille sur un projet où nous observons le comportement sur une tâche (par exemple le temps de réponse) et modélisons ce comportement en fonction de plusieurs variables expérimentalement manipulées ainsi que de plusieurs variables observées (sexe du participant, QI du participant, réponses sur...
J'ai calculé la proportion de poussins à l'envol par rapport au nombre d'œufs éclos chaque année en utilisant prop.test()R. Je vois que cela me donne la proportion à l'envol, mais aussi l'intervalle de confiance à 95%, ce que je recherche. Après avoir lu les excellentes informations d'une autre...
J'ai récemment commencé à lire Gelman et Hill's, "Data Analysis Using Regression and Multilevel / Hierarchical Models" et la question est basée sur cela: L'échantillon contient 6 observations sur les proportions: p1,p2,…,p6p1,p2,…,p6p_{1}, p_{2}, \dots, p_{6} Chaque pipip_{i} a signifie πiπi\pi_{i}...
Pour quelles distributions x, autre que bêta, la distribution x-binomiale est-elle agréable? Les distributions bêta et binomiale sont réputées conjuguées mais je suis curieux de savoir si d'autres distributions non conjuguées donneront des pmfs composés relativement simples. Par gentil, je veux...
Tenter de reproduire les résultats de l'article récemment publié, Aghion, Philippe, John Van Reenen et Luigi Zingales. 2013. «Innovation and Institutional Ownership». American Economic Review, 103 (1): 277-304. (Les données et le code stata sont disponibles sur
Ce qui suit est un problème de Bayesian Data Analysis 2nd ed , p. 97. Andrew Gelman n'a pas inclus sa solution dans le guide sur son site Web et cela m'a rendu fou toute la journée. Littéralement toute la journée. Pour certaines données , modélisées comme une distribution binomiale avec la...
Nous testons une campagne de marketing par e-mail. Lors de notre test initial, nous avons envoyé deux types d'e-mails différents et un troisième groupe de contrôle n'a pas reçu d'e-mail. Maintenant, nous récupérons les "résultats" en proportion des utilisateurs qui sont revenus sur notre...
Situation: Disons que j'ai un processus de Poisson, comme la désintégration radioactive, produisant des particules R par seconde. Je mesure avec un détecteur. Il y a une probabilité P qu'une particule soit détectée par le détecteur. Choses que je pense savoir: Le temps inter-arrivée de l'émission...
Certains d'entre vous ont peut-être lu ce bel article: O'Hara RB, Kotze DJ (2010) Ne pas enregistrer les données de comptage de transformation. Méthodes en écologie et évolution 1: 118-122. Klick . Actuellement, je compare des modèles binomiaux négatifs avec des modèles gaussiens sur des données...
Je voudrais tout d'abord préciser que je ne suis pas un expert du sujet. Supposons que deux variables aléatoires et soient binomiales, respectivement X \ sim B (n_1, p) et Y \ sim B (n_2, p), notez ici que p est le même. Je sais que Z = X + Y \ sim B (n_1 + n_2, p).XXXYYYX∼B(n1,p)X∼B(n1,p)X\sim...
Je veux adapter un modèle linéaire de R avec family=binomial(link="identity"), cependant, la famille binomiale n'a pas de lien d'identité. Que devrais-je
J'ai un échantillon d'une distribution discrète en tant que telle: Type: 0 1 2 3 4 5 Occurrences: 88 12 52 43 21 5 Ma tâche consiste à tester si une distribution binomiale (n = 5, p) correspond ou non à ces données. Je comprends que je suis censé utiliser des tests d'hypothèse et que le test du khi...
Supposons que je sache comment générer des variables aléatoires binomiales indépendantes. Comment puis-je générer deux variables aléatoiresXXX et YYY tel que X∼Bin(8,23),Y∼Bin(18,23) and Corr(X,Y)=0.5X∼Bin(8,23),Y∼Bin(18,23) and Corr(X,Y)=0.5X\sim \text{Bin}(8,\dfrac{2}{3}),\quad Y\sim...