Questions marquées «r»

Utilisez cette balise pour toute question * sur le sujet * qui (a) implique «R» en tant que partie critique de la question ou réponse attendue, et (b) n'est pas * seulement * sur la façon d'utiliser «R».

4
Comment obtenir les valeurs utilisées dans plot.gam en mgcv?
Je voudrais connaître les valeurs (x, y)utilisées dans le traçage plot(b, seWithMean=TRUE)dans le package mgcv . Est-ce que quelqu'un sait comment extraire ou calculer ces valeurs? Voici un exemple: library(mgcv) set.seed(0) dat <- gamSim(1, n=400, dist="normal", scale=2) b <- gam(y~s(x0), data=dat) plot(b, seWithMean=TRUE)
4
Étant donné une chaîne 10MC MCMC, comment puis-je déterminer son ou ses modes postérieurs dans R?
Question: Avec une chaîne MCMC à 10 dimensions, disons que je suis prêt à vous remettre une matrice des tirages: 100 000 itérations (lignes) par 10 paramètres (colonnes), comment identifier au mieux les modes postérieurs? Je suis particulièrement préoccupé par plusieurs modes. Contexte:Je me considère comme un statisticien averti en …
5
Omega au carré pour la mesure de l'effet dans R?
Le livre de statistiques que je lis recommande l'oméga carré pour mesurer les effets de mes expériences. J'ai déjà prouvé en utilisant un plan de parcelle divisé (mélange de plans intra-sujets et inter-sujets) que mes facteurs intra-sujets sont statistiquement significatifs avec p <0,001 et F = 17. Maintenant, je cherche …
1
Sortie du modèle logistique en R
J'essaie d'interpréter le type de modèle logistique suivant: mdl <- glm(c(suc,fail) ~ fac1 + fac2, data=df, family=binomial) Est la sortie des predict(mdl)chances de réussite attendues pour chaque point de données? Existe-t-il un moyen simple de tabuler les cotes pour chaque niveau de facteur du modèle, plutôt que pour tous les …
3
Rapport de vraisemblance vs test de Wald
D'après ce que j'ai lu, entre autres sur le site du groupe de consultation en statistiques de l' UCLA, les tests de rapport de vraisemblance et les tests wald sont assez similaires pour tester si deux modèles glm montrent une différence significative dans l'ajustement pour un ensemble de données (excusez-moi …
1
Test de Johansen pour la cointégration
Je teste la cointégration en utilisant le test de Johansen. J'ai vu des questions comme comment interpréter les résultats des tests, mais quand j'interprète les miens, j'ai des doutes. Dans mes résultats r = 3depuis 4.10 < 10.49, je ne peux donc pas former une série stationnaire. Il est le …
2
lmer avec multiplier les données imputées
Comment puis-je obtenir des effets aléatoires groupés pour lmer après une imputation multiple? J'utilise des souris pour imputer plusieurs trames de données. Et lme4 pour un modèle mixte avec interception aléatoire et pente aléatoire. La mise en commun de lmer se passe bien, sauf qu'elle ne regroupe pas les effets …
2
Interprétation de la corrélation CCF dans R
J'utilise ccfpour trouver une corrélation entre 2 séries chronologiques. Je reçois un complot qui ressemble à ça: Notez que je suis principalement intéressé par la corrélation pour le décalage = 0. Des questions: Interprétez-vous correctement qu'il existe une corrélation croisée pour le décalage = 0, car pour ce décalage, la …
En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.