Question La variance d'une distribution binomiale négative (NB) est toujours supérieure à sa moyenne. Lorsque la moyenne d'un échantillon est supérieure à sa variance, essayer d'ajuster les paramètres d'un NB avec une probabilité maximale ou avec une estimation de moment échouera (il n'y a pas de...
J'ai l'histogramme suivant des données de comptage. Et je voudrais y adapter une distribution discrète. Je ne sais pas comment je dois procéder. Dois-je d'abord superposer une distribution discrète, disons distribution binomiale négative, sur l'histogramme afin d'obtenir les paramètres de la...
Existe-t-il un tel ensemble qui prévoit une estimation du modèle binomial à effets mixtes négatifs gonflés zéro dans R? J'entends par là: Zero-inflation où vous pouvez spécifier le modèle binomial pour zéro inflation, comme dans la fonction zeroinfl dans le package pscl: zeroinfl (y ~ X | Z, dist =...
Quelle est la stratégie appropriée pour décider quel modèle utiliser avec les données de comptage? J'ai des données de comptage que j'ai besoin de modéliser en tant que modèle à plusieurs niveaux et il m'a été recommandé (sur ce site) que la meilleure façon de le faire soit par le biais de bogues...
J'ai souvent vu les conseils pour vérifier si un ajustement du modèle de Poisson est trop dispersé, ce qui implique de diviser la déviance résiduelle par les degrés de liberté. Le rapport résultant doit être "environ 1". La question est de savoir de quelle plage parlons-nous pour "approximative" -...
La distribution binomiale négative est devenue un modèle populaire pour les données de comptage (en particulier le nombre attendu de lectures de séquençage dans une région donnée du génome d'une expérience donnée) en bioinformatique. Les explications varient: Certains l'expliquent comme quelque...
Une distribution de Poisson peut mesurer des événements par unité de temps et le paramètre est . La distribution exponentielle mesure le temps jusqu'au prochain événement, avec le paramètre . On peut convertir une distribution dans l'autre, selon qu'il est plus facile de modéliser des événements ou...
Quel est le nombre de fois que vous devez lancer un dé jusqu'à ce que chaque camp apparaisse 3 fois? Cette question a été posée à l'école primaire en Nouvelle-Zélande et elle a été résolue à l'aide de simulations. Quelle est la solution analytique à ce
Je viens de lancer un GLM binomial négatif et voici la sortie: Call: glm.nb(formula = small ~ method + site + depth, data = size.dat, init.theta = 1.080668549, link = log) Deviance Residuals: Min 1Q Median 3Q Max -2.2452 -0.9973 -0.3028 0.3864 1.8727 Coefficients: Estimate Std. Error z value...
J'ai dirigé un glm.nb par glm1<-glm.nb(x~factor(group)) avec groupe étant une catégorie et x une variable métrique. Lorsque j'essaie d'obtenir le résumé des résultats, j'obtiens des résultats légèrement différents, selon que j'utilise summary()ou summary.glm. summary(glm1)Donne moi ......
Donc, je veux adapter un modèle binomial à effets aléatoires négatifs. Pour un tel modèle, STATA peut produire des coefficients exponentiels. Selon le fichier d'aide, ces coefficients peuvent être interprétés comme des rapports de taux d'incidence. Malheureusement, je ne suis pas un anglophone...
J'ai un très petit ensemble de données sur l'abondance des abeilles solitaires que j'ai du mal à analyser. Ce sont des données de comptage, et presque tous les comptages sont dans un traitement avec la plupart des zéros dans l'autre traitement. Il existe également quelques valeurs très élevées (une...
J'ai rencontré un certain nombre de questions pratiques lors de la modélisation des données de comptage issues de la recherche expérimentale à l'aide d'une expérience intra-sujet. Je décris brièvement l'expérience, les données et ce que j'ai fait jusqu'à présent, suivi de mes questions. Quatre...
J'ai rencontré trois propositions pour traiter la surdispersion dans une variable de réponse de Poisson et un modèle de départ à effets fixes: Utilisez un quasi modèle; Utiliser un GLM binomial négatif; Utilisez un modèle mixte avec un effet aléatoire au niveau du sujet. Mais lequel choisir...
J'ai un ensemble de données de série chronologique auquel j'essaie d'adapter un modèle de Markov caché (HMM) afin d'estimer le nombre d'états latents dans les données. Mon pseudo-code pour ce faire est le suivant: for( i in 2 : max_number_of_states ){ ... calculate HMM with i states ......
J'essaie de modéliser les intensités moyennes des parasites affectant un hôte dans R en utilisant un modèle binomial négatif. Je reçois toujours 50 avertissements ou plus qui disent: In dpois(y, mu, log = TRUE) : non-integer x = 251.529000 Comment puis-je gérer cela? Mon code ressemble à ceci:...
La question est la suivante: Un échantillon aléatoire de n valeurs est collecté à partir d'une distribution binomiale négative avec le paramètre k = 3. Trouvez l'estimateur du maximum de vraisemblance du paramètre π. Trouvez une formule asymptotique pour l'erreur standard de cet estimateur....
Dans cet article (disponible gratuitement via PubMed Central), les auteurs utilisent une régression binomiale négative pour modéliser le score sur un instrument de dépistage à 10 éléments noté 0-40. Cette procédure suppose des données de comptage, ce qui n'est clairement pas le cas ici. J'aimerais...
J'ai exécuté des modèles binomiaux et quasi-poissons négatifs basés sur une approche de test d'hypothèse. Mes modèles finaux utilisant les deux méthodes ont différentes covariables et interactions. Il semble qu'il n'y ait aucun motif lorsque je trace mes résidus dans les deux cas. Ainsi, je me...
J'ai un GLMM du formulaire: lmer(present? ~ factor1 + factor2 + continuous + factor1*continuous + (1 | factor3), family=binomial) Lorsque j'utilise drop1(model, test="Chi"), j'obtiens des résultats différents de ceux que j'utilise à Anova(model, type="III")partir du package de voiture ou...