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».
Envisager des mesures répétées d'ANOVA (RM-ANOVA) avec un facteur intra-sujet Aet plusieurs mesures par sujet pour chaque niveau de A. Elle est étroitement liée à l'ANOVA bidirectionnelle avec deux facteurs: Aet subject. Ils utilisent la décomposition identique de la somme des carrés en quatre parties: A, subject, A⋅subjectet residual. Cependant, …
J'ai un jeu de données à 1 dimension et j'utilise la boxplotfonction pour faire un diagramme en boîte. Ensuite, je peux voir que j'ai des valeurs aberrantes. Les valeurs aberrantes comptent-elles lorsque les quantiles sont déterminés? Y a-t-il une bonne / mauvaise façon ou les deux sont-elles correctes tant que …
Je voulais estimer l'intervalle de confiance pour l'écart-type de certaines données. Le code R ressemble à ceci: library(boot) sd_boot <- function (x, ind) { res <- sd(x$ReadyChange[ind], na.rm = TRUE) return(res) } data_boot <- boot::boot(data, statistic = sd_boot, R = 10000) plot(data_boot) Et j'ai l'intrigue suivante: Je suis coincé avec …
J'ai effectué un lasso, puis laissé une validation croisée cv<-cv.glmnet(df, df$Price, nfolds = 1500) Lorsque je trace un cv, j'obtiens ce qui suit: J'ai aussi remarqué que j'obtenais 2 lambdas différents: lambda.minetlambda.1se Quelle est la différence entre ces lambdas? Que puis-je comprendre de l'intrigue ci-dessus en général (à quoi servent …
J'ai des données de collecte à long terme et j'aimerais tester si le nombre d'animaux collectés est influencé par les effets météorologiques. Mon modèle ressemble à ci-dessous: glmer(SumOfCatch ~ I(pc.act.1^2) +I(pc.act.2^2) + I(pc.may.1^2) + I(pc.may.2^2) + SampSize + as.factor(samp.prog) + (1|year/month), control=glmerControl(optimizer="bobyqa", optCtrl=list(maxfun=1e9,npt=5)), family="poisson", data=a2) Explication des variables utilisées: SumOfCatch: …
J'essaie de déterminer si les chenilles qui ont une alimentation naturelle (fleur de singe) sont plus résistantes aux prédateurs (fourmis) que les chenilles qui ont une alimentation artificielle (un mélange de germe de blé et de vitamines). J'ai fait une étude d'essai avec un petit échantillon (20 chenilles; 10 par …
Considérez la régression binomiale suivante: # Create some data set.seed(10) n <- 500 x <- runif(n,0,100) y <- x + rnorm(n,sd=100) < 0 # Fit a binomial regression model model <- glm(y ~ x, family="binomial") summary(model) La summaryfonction renvoie une valeur de p de 1.03e-05. Lors de l'utilisation anova.glm, on …
Je n'arrive pas à interpréter ce graphique. Ma variable dépendante est le nombre total de billets de cinéma qui seront vendus pour un spectacle. Les variables indépendantes sont le nombre de jours restants avant le spectacle, les variables factices saisonnières (jour de la semaine, mois de l'année, vacances), le prix, …
Je veux régresser l'économie de carburant sur la cylindrée du moteur, le type de carburant, la transmission 2 vs 4 roues motrices, la puissance, la transmission manuelle vs automatique et le nombre de vitesses. Mon ensemble de données ( lien ) contient des véhicules de 2012 à 2014. fuelEconomy en …
Il semble que je puisse très bien apprendre les paramètres et trouver les probabilités postérieures pour les données d'entraînement, mais je n'ai aucune idée de comment faire de nouvelles prédictions sur de nouvelles données. Le problème vient en particulier des probabilités de transition changeant sur les covariables, il n'est donc …
Fermé. Cette question est hors sujet . Il n'accepte pas actuellement de 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é l'année dernière . J'essaie d'adapter le modèle linéaire généralisé pour la famille weibull, mais quand je l'essaie …
J'ai lu les excellents commentaires sur la façon de traiter les valeurs manquantes avant d'appliquer SVD, mais j'aimerais savoir comment cela fonctionne avec un exemple simple: Movie1 Movie2 Movie3 User1 5 4 User2 2 5 5 User3 3 4 User4 1 5 User5 5 1 5 Étant donné la matrice …
Soit un processus stochastique formé en concaténant les tirages iid d'un processus AR (1), où chaque tirage est un vecteur de longueur 10. En d'autres termes, sont des réalisations d'un processus AR (1); sont tirés du même processus, mais sont indépendants des 10 premières observations; etc.{Xt}{Xt}\left\{X_t\right\}{X1,X2, … ,Xdix}{X1,X2,…,Xdix}\left\{X_1, X_2, \ldots, …
J'ajuste quelques modèles additifs généralisés en utilisant le mgcvpackage en R, et je veux tester entre deux modèles; si je peux supprimer un terme ou non. Cependant, j'obtiens des résultats contradictoires (pour autant que je sache). Un modèle, m1avec un terme lisse pour xajouté, semble donner un meilleur ajustement en …
Quelle est la difference entre (1|DNA.concentration/mouse.id)et (DNA.concentration|mouse.id)? Que signifient les symboles |et la /signification dans la syntaxe de l'effet aléatoire?
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.