Quelle est la valeur donnée dans le résumé d'un modèle coxph dans R? Par exemple,R2R2R^2 Rsquare= 0.186 (max possible= 0.991 ) Je l'ai insensément inclus un manuscrit en tant que valeur et le critique a sauté dessus en disant qu'il n'était pas au courant d'un analogue de la statistique de la...
Pour ma thèse de maîtrise, je travaille sur le développement d'un modèle statistique des transitions entre différents états, défini par le statut sérologique. Pour l'instant, je ne donnerai pas trop de détails dans ce contexte, car ma question est plus générale / théorique. Quoi qu'il en soit, mon...
Je découvre le monde merveilleux de ces "modèles de Markov cachés", également appelés "modèles à changement de régime". Je souhaite adapter un HMM en R pour détecter les tendances et les tournants. Je voudrais construire le modèle le plus générique possible afin de pouvoir le tester sur de nombreux...
J'ai des données de patients traités avec 2 types de traitements différents pendant la chirurgie. J'ai besoin d'analyser son effet sur la fréquence cardiaque. La mesure de la fréquence cardiaque est prise toutes les 15 minutes. Étant donné que la durée de la chirurgie peut être différente pour...
Considérez le modèle linéaire simple: yy=X′ββ+ϵyy=X′ββ+ϵ\pmb{y}=X'\pmb{\beta}+\epsilon où et , et contient une colonne des constantes.ϵi∼i.i.d.N(0,σ2)ϵi∼i.i.d.N(0,σ2)\epsilon_i\sim\mathrm{i.i.d.}\;\mathcal{N}(0,\sigma^2)X∈Rn×pX∈Rn×pX\in\mathbb{R}^{n\times p}p≥2p≥2p\geq2XXX Ma question est, étant...
J'ai des données de densité de poisson que j'essaie de comparer entre plusieurs techniques de collecte différentes, les données ont beaucoup de zéros et l'histogramme semble approprié pour une distribution de poisson, sauf que, comme densités, ce ne sont pas des données entières. Je suis...
Je suis vraiment stupéfait par le fait que le GLM de Poisson accepte les nombres non entiers! Regardez: Données (contenu de data.txt): 1 2001 0.25 1 1 2002 0.5 1 1 2003 1 1 2 2001 0.25 1 2 2002 0.5 1 2 2003 1 1 Script R: t <- read.table("data.txt") names(t) <- c('site', 'year', 'count',...
Après avoir monté un coxmodel, il est possible de faire des prédictions et de récupérer le risque relatif de nouvelles données. Ce que je ne comprends pas, c'est comment le risque relatif est calculé pour un individu et à quoi est-il relatif (c'est-à-dire la moyenne de la population)? Des...
Lorsque vous faites un GLM et que vous obtenez l'erreur "non défini en raison des singularités" dans la sortie anova, comment peut-on contrer cette erreur? Certains ont suggéré que cela est dû à la colinéarité entre les covariables ou que l'un des niveaux n'est pas présent dans l'ensemble de...
La régression bêta (c'est-à-dire GLM avec distribution bêta et généralement la fonction de lien logit) est souvent recommandée pour traiter la réponse aka variable dépendante prenant des valeurs comprises entre 0 et 1, telles que les fractions, les ratios ou les probabilités: régression pour un...
Nous travaillons avec quelques régressions logistiques et nous avons réalisé que la probabilité moyenne estimée est toujours égale à la proportion de celles de l'échantillon; c'est-à-dire que la moyenne des valeurs ajustées est égale à la moyenne de l'échantillon. Quelqu'un peut-il m'expliquer la...
J'ai formé un modèle de régression linéaire, en utilisant un ensemble de variables / fonctionnalités. Et le modèle a de bonnes performances. Cependant, j'ai réalisé qu'il n'y a pas de variable avec une bonne corrélation avec la variable prédite. Comment est-ce
Mise à jour : Désolé pour une autre mise à jour, mais j'ai trouvé des solutions possibles avec des polynômes fractionnaires et le package de risques concurrents avec lesquels j'ai besoin d'aide. Le problème Je ne trouve pas de moyen facile de faire une analyse de coefficient dépendant du temps dans...
J'ai appris le modèle de risque proportionnel de Cox. J'ai beaucoup d'expérience dans l'ajustement de modèles de régression logistique, et donc pour construire l'intuition, j'ai comparé des modèles ajustés à l'aide coxphde la «survie» R avec des modèles de régression logistique ajustés à l'aide...
Dites que vous êtes dans la bibliothèque de votre département de statistique et que vous tombez sur un livre avec l'image suivante en première page. Vous penserez probablement qu'il s'agit d'un livre sur les choses de régression linéaire. Quelle serait l'image qui vous ferait penser aux modèles...
On peut effectuer une régression logit dans R en utilisant un tel code: > library(MASS) > data(menarche) > glm.out = glm(cbind(Menarche, Total-Menarche) ~ Age, + family=binomial(logit), data=menarche) > coefficients(glm.out) (Intercept) Age -21.226395 1.631968 Il semble que l'algorithme...
J'ai six variables dépendantes (données de comptage) et plusieurs variables indépendantes, je vois que dans un MMR le script va comme ceci: my.model <- lm(cbind(DV1,DV2,DV3,DV4,DV5,DV6) ~ IV1 + IV2 + ... + IVn) Mais, comme mes données sont des nombres, je veux utiliser un modèle linéaire...
En pratique, l'utilisation d'un test T standard pour vérifier la signification d'un coefficient de régression linéaire est une pratique courante. La mécanique du calcul a du sens pour moi. Pourquoi la distribution T peut-elle être utilisée pour modéliser la statistique de test standard utilisée...
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 =...
Je suis AW van der Vaart, statistiques asymptotiques (1998). Il parle d'expériences statistiques, affirmant qu'elles sont différentes d'un modèle statistique, mais il ne définit ni l'une ni l'autre. Ma question: Qu'est-ce que (1) une expérience statistique, (2) un modèle statistique et (3) quel est...