Le coefficient de corrélation de Pearson est calculé à l'aide de la formule r =c o v ( X, Y)v a r ( X)√v a r ( Y)√r=cov(X,Y)var(X)var(Y)r = \frac{cov(X,Y)}{\sqrt{var(X)} \sqrt{var(Y)}}. Comment cette formule contient-elle les informations selon lesquelles les deux variables et sont corrélées ou...
Je voudrais obtenir les simulations postérieures des composantes de variance d'un modèle lmer () avec la fonction mcmcsamp (). Comment faire ? Par exemple, ci-dessous est le résultat d'un ajustement lmer (): > fit Linear mixed model fit by REML Formula: y ~ 1 + (1 | Part) + (1 | Operator) + (1 |...
J'essaie d'adapter une distribution gamma inverse à trois paramètres à mes données en R ou en Python. Je voudrais le faire en utilisant l'estimation du maximum de vraisemblance (MLE). Le pdf du gamma inverse à trois paramètres est donné par: Où Γ est la fonction gamma, ρ est la forme, α est...
Contexte : Israël (et le Moyen-Orient en général) est rempli de protestations. Je suis curieux, lorsqu'on me donne une photo, d'estimer le nombre de personnes sur cette photo (souvent une photo d'une grande foule). Quelle modélisation peut offrir une solution à ce problème? (et, bien sûr, cela...
En tant qu'institution financière, nous nous heurtons souvent à l'analyse de données chronologiques. Souvent, nous finissons par faire une régression en utilisant des variables de séries chronologiques. Dans ce cas, nous rencontrons souvent des résidus avec une structure de séries chronologiques...
J'ai le plmpaquet et je voudrais exécuter des tests de racine unitaire sur certaines variables. J'obtiens l'erreur suivante: > purtest(data$tot.emp) Error in data.frame(baldwin = c(59870, 61259, 60397, 58919, 57856, 57227, : arguments imply differing number of rows: 14, 19, 11, 12, 1, 20, 18,...
Comment puis-je compléter le carré à partir du point où je me suis arrêté, et est-ce correct jusqu'à présent? J'ai un avant normal pour ββ\beta de la forme p ( β|σ2) ∼ N( 0 ,σ2V)p(β|σ2)∼N(0,σ2V)p(\beta|\sigma^2)\sim \mathcal{N}(0,\sigma^2V), obtenir: p ( β|σ2) = ( 2 πσ2V)p2exp[...
J'ai énormément de difficulté à évaluer avec le package en R. Dans mon cas, les valeurs de , , sont toujours des nombres réels positifs. Même ainsi, la fonction hypergéométrique est incroyablement sensible à leurs valeurs. Je ne recherche pas une précision extrême; Je peux utiliser Excel pour...
Lors du calcul de modèles de régression avec R, j'utilise régulièrement la fonction relevel pour que mon modèle me donne également des résultats pour l'autre niveau. J'ai remarqué que parfois, mais pas souvent, cela changeait le modèle dans le sens où les niveaux d'autres facteurs qui étaient...
J'analyse un certain ensemble de données et j'ai besoin de comprendre comment choisir le meilleur modèle qui correspond à mes données. J'utilise R. Un exemple de données dont je dispose est le suivant: corr <- c(0, 0, 10, 50, 70, 100, 100, 100, 90, 100, 100) Ces chiffres correspondent au...
Étant donné que l'estimation postérieure de σ′2σ′2\sigma'^{2} d'une vraisemblance normale et d'un gamma inverse antérieur sur σ2σ2\sigma^2 est: σ′2∼IG(α+n2,β+∑ni=1(yi−μ)22)σ′2∼IG(α+n2,β+∑i=1n(yi−μ)22)\sigma'^{2}\sim\textrm{IG}\left(\alpha + \frac{n}{2}, \beta...
Une variable de réponse y est une fonction non linéaire d'un certain nombre de variables prédictives X (dans mes données réelles, la réponse est distribuée de façon binomiale, mais ici j'utilise une valeur normalement distribuée pour plus de simplicité). Je peux modéliser les relations entre les...
Je veux que R affiche les données qu'il me donne de la summary()fonction dans un tableau afin que je puisse facilement les partager. Je suis actuellement en train de faire summary()dans la console, puis de prendre une capture d'écran, mais je préfère que cela soit généré comme un joli tableau, tout...
Je dois faire quelques comparaisons moyennes simples entre les groupes (tests F ANOVA de base) sur des données avec des valeurs manquantes. J'utilise le package de souris dans R pour l'imputation multiple, mais je ne peux regrouper les résultats que pour les coefficients du modèle linéaire, ou le...
Je calcule des probabilités conditionnelles et des intervalles de confiance à 95% associés. Pour bon nombre de mes cas, j'ai un décompte simple des xsuccès des nessais (à partir d'un tableau de contingence), donc je peux utiliser un intervalle de confiance binomial, tel que celui fourni par...
Supposons que je construis un classificateur de régression logistique qui prédit si quelqu'un est marié ou célibataire. (1 = marié, 0 = célibataire) Je veux choisir un point sur la courbe précision-rappel qui me donne au moins 75% de précision, donc je veux choisir les seuils et , de sorte...
J'ai un problème de régression multiple, que j'ai essayé de résoudre en utilisant une régression multiple simple: model1 <- lm(Y ~ X1 + X2 + X3 + X4 + X5, data=data) Cela semble expliquer les 85% de variance (selon R au carré), ce qui semble assez bon. Cependant, ce qui m'inquiète, c'est...
Mise à jour : Puisque je sais maintenant que mon problème est appelé séparation quasi-complète, j'ai mis à jour la question pour refléter cela (merci à Aaron). J'ai un ensemble de données provenant d'une expérience dans laquelle 29 participants humains (facteur code) ont travaillé sur un ensemble...
J'ai 2 facteurs Aet B(5 × 3) et une covariable Xdans une conception intra-sujet. Voici comment je spécifie mon modèle global: lme.out = lme(y~ A*B*X, random=~1|Subject, data=mydata) Mon interprétation est que je regarde un graphique y~x, où la pente change en raison de la covariable, et la ligne se...
Je ne sais pas vraiment ce qui est possible et je voudrais un pointeur dans la bonne direction. J'ai des mesures de temps et de position qui pourraient être n'importe quoi d'une personne marchant, un véhicule sur une route, ou dans un parking, ou une imprimante dans un bureau. Je dois calculer les...