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».
Disons que j'ai des données: x1 <- rnorm(100,2,10) x2 <- rnorm(100,2,10) y <- x1+x2+x1*x2+rnorm(100,1,2) dat <- data.frame(y=y,x1=x1,x2=x2) res <- lm(y~x1*x2,data=dat) summary(res) Je veux tracer l'interaction continue par continue de telle sorte que x1 soit sur l'axe X et x2 soit représenté par 3 lignes, une qui représente x2 à un …
Voici ce que je veux faire, mais il ne semble pas y avoir de predictméthode pour le mlogit. Des idées? library(mlogit) data("Fishing", package = "mlogit") Fish <- mlogit.data(Fishing, varying = c(2:9), shape = "wide", choice = "mode") Fish_fit<-Fish[-1,] Fish_test<-Fish[1,] m <- mlogit(mode ~price+ catch | income, data = Fish_fit) predict(m,newdata=Fish_test)
Je cherche à expliquer (visuellement) la corrélation linéaire simple aux étudiants de première année. La manière classique de visualiser serait de donner un nuage de points Y ~ X avec une droite de régression. Récemment, j'ai eu l'idée d'étendre ce type de graphisme en ajoutant à l'intrigue 3 images supplémentaires, …
En supposant donc qu'il est utile de tester l'hypothèse de normalité pour anova (voir 1 et 2 ) Comment peut-il être testé en R? Je m'attendrais à faire quelque chose comme: ## From Venables and Ripley (2002) p.165. utils::data(npk, package="MASS") npk.aovE <- aov(yield ~ N*P*K + Error(block), npk) residuals(npk.aovE) qqnorm(residuals(npk.aov)) …
Quiconque suit le baseball a probablement entendu parler de la performance de type MVP de nulle part de Jose Bautista de Toronto. Au cours des quatre années précédentes, il a réalisé environ 15 circuits par saison. L'année dernière, il a atteint 54, un nombre dépassé par seulement 12 joueurs dans …
Je veux implémenter (en R) le modèle linéaire dynamique très simple suivant pour lequel j'ai 2 paramètres inconnus variant dans le temps (la variance de l'erreur d'observation et la variance de l'erreur d'état ϵ 2 t ).ϵ1tϵt1\epsilon^1_tϵ2tϵt2\epsilon^2_t Ouitθt + 1==θt+ ϵ1tθt+ ϵ2tYt=θt+ϵt1θt+1=θt+ϵt2 \begin{matrix} Y_t & = & \theta_t + \epsilon^1_t\\ …
Exemple de code: (pc.cr <- princomp(USArrests)) summary(pc.cr) loadings(pc.cr) ## note that blank entries are small but not zero J'obtiens des sorties différentes de chacun, et je ne suis pas sûr de comprendre quelle est la différence. Voici la sortie: > summary(pc.cr) Importance of components: Comp.1 Comp.2 Comp.3 Comp.4 Standard deviation …
Je recherche des corrélations entre les réponses à différentes questions dans une enquête ("euh, voyons si les réponses à la question 11 sont en corrélation avec celles de la question 78"). Toutes les réponses sont catégoriques (la plupart vont de «très malheureux» à «très heureux»), mais quelques-unes ont un ensemble …
Fermé. Cette question est hors sujet . Il n'accepte pas actuellement les réponses. Voulez-vous améliorer cette question? Mettez à jour la question afin qu'elle soit sur le sujet pour la validation croisée. Fermé il y a 2 ans . Je veux calculer un résumé d'une variable dans un data.frame pour …
Au risque de rendre la question spécifique au logiciel, et avec l'excuse de son ubiquité et de ses idiosyncrasies, je veux poser des questions sur la fonction biplot()dans R, et, plus précisément, sur le calcul et le tracé de sa valeur par défaut, des flèches rouges superposées, correspondant aux variables …
Le mgcvpackage pour Ra deux fonctions pour ajuster les interactions des produits tensoriels: te()et ti(). Je comprends la division de base du travail entre les deux (ajustement d'une interaction non linéaire vs décomposition de cette interaction en effets principaux et interaction). Ce que je ne comprends pas, c'est pourquoi te(x1, …
Je me demande comment puis-je simuler un échantillon de n durées de vie de distribution Weibull qui incluent des observations de type I censurées à droite. Par exemple, ayons le n = 3, la forme = 3, l'échelle = 1 et le taux de censure = .15, et le temps …
J'ai un cadre de données qui contient des valeurs sur 4 colonnes: Par exemple: ID, price, click count,rating Ce que je voudrais faire, c'est "diviser" cette trame de données en N groupes différents où chaque groupe aura un nombre égal de lignes avec la même distribution de prix, le nombre …
We use cookies and other tracking technologies to improve your browsing experience on our website,
to show you personalized content and targeted ads, to analyze our website traffic,
and to understand where our visitors are coming from.
By continuing, you consent to our use of cookies and other tracking technologies and
affirm you're at least 16 years old or have consent from a parent or guardian.