Pour une régression linéaire avec plusieurs groupes (groupes naturels définis a priori), est-il acceptable d'exécuter deux modèles différents sur le même ensemble de données pour répondre aux deux questions suivantes? Chaque groupe a-t-il une pente non nulle et une intersection non nulle et quels...
J'essaie de déterminer si de simples probabilités fonctionneront pour mon problème ou s'il vaudra mieux utiliser (et en savoir plus) des méthodes plus sophistiquées comme la régression logistique. La variable de réponse dans ce problème est une réponse binaire (0, 1). J'ai un certain nombre de...
J'essaie de reproduire le travail d'un collègue et je déplace l'analyse de Stata vers R. Les modèles qu'elle utilise invoquent l'option "cluster" dans la fonction nbreg pour regrouper les erreurs standard. Voir http://repec.org/usug2007/crse.pdf pour une description assez complète du quoi et du...
J'essaie d'analyser le décalage entre les séries chronologiques de deux cours boursiers. Dans l'analyse de séries chronologiques régulières, nous pouvons faire Cross Correlaton, VECM (Granger Causality). Cependant, comment peut-on gérer la même chose dans des séries temporelles irrégulièrement...
Je souhaite modéliser deux variables temporelles différentes, dont certaines sont fortement colinéaires dans mes données (âge + cohorte = période). Ce faisant, j'ai rencontré des problèmes avec lmeret et les interactions de poly(), mais ce n'est probablement pas limité à lmer, j'ai obtenu les mêmes...
L'image ci-dessous est tirée de cet article dans Psychological Science . Un collègue a souligné deux choses inhabituelles à ce sujet: Selon la légende, les barres d'erreur indiquent «± 2,04 erreurs standard, l'intervalle de confiance à 95%». Je n'ai vu que ± 1,96 SE utilisé pour l'IC à 95%, et je...
Je fais une ANOVA à sens unique (par espèce) avec des contrastes personnalisés. [,1] [,2] [,3] [,4] 0.5 -1 0 0 0 5 1 -1 0 0 12.5 0 1 -1 0 25 0 0 1 -1 50 0 0 0 1 où je compare l'intensité 0,5 contre 5, 5 contre 12,5 et ainsi de suite. Ce sont les données sur lesquelles je travaille avec les...
J'ai remarqué que lorsque j'essaie de trouver les valeurs critiques pour le Mann-Whitney U en utilisant R, les valeurs sont toujours 1 + valeur critique. Par exemple, pour , la valeur critique (bilatérale) est 8, tandis que pour , la (bilatérale) ) la valeur critique est 22 (consultez les tableaux...
Dans l'exemple suivant, j'ai une trame de données qui se compose d'une série chronologique de mesures de la température de l'eau enregistrées à 5 profondeurs dans l'océan où chaque valeur en Tempcorrespond à la date DateTimeet à la profondeur Depth. set.seed(1) Temp <- rnorm(43800,sd=20) AirT...
Les modèles de pièces biaisées ont généralement un paramètre . Une façon d'estimer la partir d'une série de tirages est d'utiliser une distribution bêta a priori et de calculer la distribution postérieure avec vraisemblance binomiale.θθ=P(Head|θ)θ=P(Head|θ)\theta = P(\text{Head} | \theta)θθ\theta...
J'ai adapté un modèle lognormal utilisant R avec un ensemble de données. Les paramètres résultants étaient: meanlog = 4.2991610 sdlog = 0.5511349 J'aimerais transférer ce modèle sur Scipy, que je n'ai jamais utilisé auparavant. En utilisant Scipy, j'ai pu obtenir une forme et une échelle de 1 et...
Structure de données > str(data) 'data.frame': 6138 obs. of 10 variables: $ RT : int 484 391 422 516 563 531 406 500 516 578 ... $ ASCORE : num 5.1 4 3.8 2.6 2.7 6.5 4.9 2.9 2.6 7.2 ... $ HSCORE : num 6 2.1 7.9 1 6.9 8.9 8.2 3.6 1.7 8.6 ... $ MVMNT : Factor w/ 2 levels "_Withd","Appr": 2 2 1 1 2...
Je n'ai pas la différence entre le rfobject$importanceet importance(rfobject)dans la colonne MeanDecreaseAccuracy. Exemple: > data("iris") > fit <- randomForest(Species~., data=iris, importance=TRUE) > fit$importance setosa versicolor virginica MeanDecreaseAccuracy MeanDecreaseGini...
Mes questions portent sur les GAM dans le package mgcv R. En raison de la petite taille de l'échantillon, je souhaite déterminer l'erreur de prédiction à l'aide de la validation croisée avec laisser-un-out. Est-ce raisonnable? Existe-t-il un package ou un code pour y parvenir? La errorest()fonction...
J'ai une expérience qui est exécutée sur des centaines d'ordinateurs répartis dans le monde entier qui mesure les événements de certains événements. Les événements dépendent chacun les uns des autres pour que je puisse les commander dans l'ordre croissant et ensuite calculer la différence de temps....
J'ai une question concernant la comparaison des modèles à l'aide des facteurs Bayes. Dans de nombreux cas, les statisticiens sont intéressés à utiliser une approche bayésienne avec des a priori impropres (par exemple certains a priori de Jeffreys et a priori de référence). Ma question est, dans les...
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...
Je voudrais faire une analyse de puissance pour un seul échantillon à partir de données binomiales, avec , contre , où est la proportion de succès dans la population. Si , je pourrais utiliser soit l'approximation normale du binôme, soit -test, mais avec , les deux échouent. J'aimerais savoir s'il...
Que se passe-t-il lorsque vous n'avez aucune idée de la distribution des paramètres? Quelle approche devrions-nous utiliser? La plupart du temps, nous visons à sous-estimer si une certaine variable a une influence sur la présence / absence d'une certaine espèce, et si la variable est acceptée ou...
Je me sens dépassé après avoir tenté de creuser dans la littérature sur la façon d'exécuter mon analyse de modèle mixte, en suivant l'utilisation d'AIC pour sélectionner le ou les meilleurs modèles. Je ne pense pas que mes données soient si compliquées, mais je recherche la confirmation que ce que...