Pardonnez-moi si j'ai raté quelque chose d'assez évident. Je suis physicien avec ce qui est essentiellement une distribution (histogramme) centrée sur une valeur moyenne qui se rapproche d'une distribution normale. La valeur importante pour moi est l'écart type de cette variable aléatoire...
Supposons que nous ayons l'ensemble de données suivant: Men Women Dieting 10 30 Non-dieting 5 60 Si je lance le test exact de Fisher dans R, qu'est-ce que cela alternative = greaterimplique (ou moins)? Par exemple: mat = matrix(c(10,5,30,60), 2,2) fisher.test(mat, alternative="greater") Je reçois...
Nous essayons de créer des valeurs aléatoires auto-corrélées qui seront utilisées comme série temporelle. Nous n'avons aucune donnée existante à laquelle nous nous référons et souhaitons simplement créer le vecteur à partir de zéro. D'une part, nous avons bien sûr besoin d'un processus aléatoire...
Je fais des recherches sur le domaine de la réponse fonctionnelle des acariens. Je voudrais faire une régression pour estimer les paramètres (taux d'attaque et temps de manipulation) de la fonction Rogers type II. J'ai un ensemble de données de mesures. Comment puis-je déterminer au mieux les...
J'ai exécuté un glm dans R, et près du bas de la summary()sortie, il indique (Dispersion parameter for gaussian family taken to be 28.35031) J'ai fait quelques recherches sur Google et j'ai appris que le paramètre de dispersion est utilisé pour ajuster les erreurs standard. J'espère que quelqu'un...
Je veux faire une régression logistique ordinale dans R sans l'hypothèse de cotes de proportionnalité. Je sais que cela peut être fait directement en utilisant la vglm()fonction en Rdéfinissant parallel=FALSE. Mais mon problème est de savoir comment corriger un ensemble particulier de coefficients...
Nous avons demandé à 60 personnes de répertorier autant de franchises de restaurants à Atlanta que possible. La liste globale comprenait plus de 70 restaurants, mais nous avons éliminé ceux qui étaient mentionnés par moins de 10% des gens, nous en laissant 45. Pour ces 45, nous avons calculé la...
J'essaie d'appliquer une série chronologique aux données échantillonnées trimestriellement (biomasse animale) sur une période de 10 ans avec 3 répétitions par trimestre. Donc 40 dates mais 120 observations au total. J'ai lu SARIMA'a dans Shumway et Stoffer's Time Series Analysis et ses...
Je souhaite visualiser les résultats d'un clustering (produit avec protoclust{protoclust}) en créant des graphiques de scater pour chaque paire de variables utilisées pour classer mes données, en coloriant par classes et en chevauchant les ellipses pour l'intervalle de confiance à 95% pour chacune...
Salut, je prends un cours d'études supérieures en statistique et nous avons couvert les statistiques de test et d'autres concepts. Cependant, je suis souvent en mesure d'appliquer les formules et de développer une sorte d'intuition sur le fonctionnement des choses, mais j'ai souvent le sentiment...
Je voudrais tester dans quelle régression correspond le mieux à mes données. Ma variable dépendante est un nombre et a beaucoup de zéros. Et j'aurais besoin d'aide pour déterminer le modèle et la famille à utiliser (poisson ou quasipoisson, ou régression poisson zéro gonflé), et comment tester les...
Je voudrais utiliser lme4pour ajuster une régression à effets mixtes et multcomppour calculer les comparaisons par paire. J'ai un ensemble de données complexe avec plusieurs prédicteurs continus et catégoriques, mais ma question peut être démontrée en utilisant l' ChickWeightensemble de données...
La documentation R pour l'un ou l'autre n'apporte pas beaucoup de lumière. Tout ce que je peux obtenir de ce lien est que l'utilisation de l'un ou l'autre devrait être correcte. Ce que je ne comprends pas, c'est pourquoi ils ne sont pas égaux. Réalité: La fonction de régression pas à pas dans R,...
Pour analyser les dénombrements d'oiseaux gonflés à zéro, j'aimerais appliquer des modèles de comptage à gonflage nul en utilisant le package R pscl . Cependant, après avoir regardé l'exemple fourni dans la documentation pour l'une des principales fonctions ( ? Zeroinfl ), je commence à douter du...
L'implémentation randomForest ne permet pas l'échantillonnage au-delà du nombre d'observations, même lors d'un échantillonnage avec remplacement. Pourquoi est-ce? Fonctionne bien: rf <- randomForest(Species ~ ., iris, sampsize=c(1, 1, 1), replace=TRUE) rf <- randomForest(Species ~ ., iris,...
Concernant la pondération du score de propension (IPTW) lors de la modélisation du risque proportionnel de Cox des données de survie dans le temps: J'ai des données de registre prospectives où nous souhaitons étudier l'effet du traitement d'un médicament que, dans la plupart des cas, les patients...
Dans l'article fondateur "Rao-Blackwellised Particle Filtering for Dynamic Bayesian Networks" par A. Doucet et. Al. un filtre séquentiel de monte carlo (filtre à particules) est proposé, qui utilise une sous-structure linéaire dans un processus de markov . Par marginalisation de cette structure...
Pour les données de comptage que j'ai collectées, j'utilise la régression de Poisson pour construire des modèles. Je le fais en utilisant la glmfonction dans R, où j'utilise family = "poisson". Pour évaluer les modèles possibles (j'ai plusieurs prédicteurs) j'utilise l'AIC. Jusqu'ici tout va bien....
J'ai un ensemble de données contenant le nombre d'actions effectuées par des individus au cours de 7 jours. L'action spécifique ne devrait pas être pertinente pour cette question. Voici quelques statistiques descriptives pour l'ensemble de données: GammeSignifierVarianceNombre d'observations0 -...
Quelle est la signification de t valueet Pr(>|t|)lors de l'utilisation de la summary()fonction sur le modèle de régression linéaire dans R? Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 10.1595 1.3603 7.469 1.11e-13 *** log(var) 0.3422 0.1597 2.143 0.0322 *...