Considérez le code et la sortie suivants: par(mfrow=c(3,2)) # generate random data from weibull distribution x = rweibull(20, 8, 2) # Quantile-Quantile Plot for different distributions qqPlot(x, "log-normal") qqPlot(x, "normal") qqPlot(x, "exponential", DB = TRUE) qqPlot(x, "cauchy") qqPlot(x,...
J'ai un problème avec la 2l.normméthode d'imputation multiniveau dans mice. Malheureusement, je ne peux pas publier d'exemple reproductible en raison de la taille de mes données - lorsque je réduis la taille, le problème disparaît. Pour une variable particulière, micegénère les erreurs et...
J'ai un ensemble de données sur 20 ans d'un décompte annuel de l'abondance des espèces pour un ensemble de polygones (~ 200 polygones de forme irrégulière et continue). J'ai utilisé une analyse de régression pour inférer les tendances (changement du nombre par an) pour chaque polygone, ainsi que...
Je suis tombé sur cette image dans un article de blog ici . J'ai été déçu que la lecture de la déclaration n'ait pas suscité la même expression faciale pour moi que pour ce type. Alors, que veut dire l'affirmation selon laquelle l'hypothèse nulle est de savoir comment les fréquentistes expriment un...
J'essaie d'obtenir le prior de Jeffreys pour une distribution binomiale négative. Je ne vois pas où je me trompe, donc si quelqu'un pouvait aider à le souligner, ce serait apprécié. D'accord, la situation est la suivante: je dois comparer les distributions antérieures obtenues à l'aide d'un binôme...
J'ai deux séries chronologiques de données quotidiennes. L'un est sign-upset l'autre terminationsdes abonnements. J'aimerais prédire ce dernier en utilisant les informations contenues dans les deux variables. En regardant le graphique de ces séries, il est évident que les terminaisons sont...
Version courte: existe-t-il une introduction aux écrits de Ronald Fisher (articles et livres) sur les statistiques qui s'adresse à ceux qui ont peu ou pas de connaissances en statistiques? Je pense à quelque chose comme un "lecteur Fisher annoté" destiné aux non-statisticiens. J'explique la...
J'ai lu à propos de l'estimation du maximum de vraisemblance et de l'estimation maximum a posteriori et jusqu'à présent, je n'ai rencontré d'exemples concrets qu'avec l'estimation du maximum de vraisemblance. J'ai trouvé quelques exemples abstraits d'estimation maximale a posteriori, mais rien de...
J'essaie d'utiliser le neuralnetpackage de R (documentation ici ) pour la prédiction. Voici ce que j'essaie de faire: library(neuralnet) x <- cbind(runif(50, min=1, max=500), runif(50, min=1, max=500)) y <- x[, 1] * x[, 2] train <- data.frame(x, y) n <- names(train) f <-...
Je ne sais pas comment ces complots sont appelés et j'ai donc donné à cette question un titre stupide. Disons que j'ai un ensemble de données ordonné comme suit 4253 4262 4270 4383 4394 4476 4635 ... Chaque numéro correspond au nombre de publications qu'un certain utilisateur a contribué à un site...
Existe-t-il une stratégie pour choisir le nombre d'arbres dans un GBM? Plus précisément, l' ntreesargument dans Rla gbmfonction de. Je ne vois pas pourquoi vous ne devriez pas définir ntreesla valeur raisonnable la plus élevée. J'ai remarqué qu'un plus grand nombre d'arbres réduit clairement la...
En utilisant Amelia dans R, j'ai obtenu plusieurs jeux de données imputés. Après cela, j'ai effectué un test de mesures répétées dans SPSS. Maintenant, je veux regrouper les résultats des tests. Je sais que je peux utiliser les règles de Rubin (implémentées via n'importe quel package d'imputation...
J'ai utilisé le code R suivant pour s'adapter à un modèle probit: p1 <- glm(natijeh ~ ., family=binomial(probit), data=data1) stepwise(p1, direction='backward/forward', criterion='BIC') Je veux savoir ce que fait stepwiseet backward/forwardfait exactement et comment sélectionner les...
J'ai une question concernant l'interprétation de l'appel tsboot dans R. J'ai vérifié la documentation du Kendall et du paquet de démarrage, mais je ne suis pas plus intelligent qu'avant. Quand je lance un bootstrap en utilisant par exemple l'exemple dans le paquet Kendall, où la statistique de test...
Dans un modèle hiérarchique de données où il semble typique en pratique de choisir des valeurs ( sorte que la moyenne et la variance de la distribution gamma correspondent à peu près à la moyenne et à la variance des données (par exemple, Clayton et Kaldor, 1987 "Empirical Bayes Estimates of...
J'exécute une simulation rapide pour comparer différentes méthodes de clustering et je suis actuellement confronté à un problème en essayant d'évaluer les solutions de cluster. Je connais différentes métriques de validation (beaucoup se trouvent dans cluster.stats () dans R), mais je suppose que...
J'ai appliqué quelques données pour trouver la meilleure solution de variables du modèle de régression en utilisant la régression de crête dans R. J'ai utilisé lm.ridgeet glmnet(quand alpha=0), mais les résultats sont très différents, surtout quand lambda=0. Il suppose que les deux estimateurs de...
Pour les non-statisticiens comme moi, il est très difficile de saisir l'idée de VImétrique (variation des informations) même après avoir lu l'article pertinent de Marina Melia " Comparing clusterings - An information based distance " (Journal of Multivariate Analysis, 2007). En fait, je ne connais...
J'aimerais vérifier Rsi mes données correspondent aux distributions log-normales ou Pareto. Comment pourrais-je faire ça? Peut ks.test- être pourrais - je m'aider à le faire, mais comment puis-je obtenir les paramètres et pour la distribution de Pareto pour mes
J'ai écrit un programme pour simuler un pronation aléatoire de la carte. Chaque carte est numérotée, avec une couleur allant de CLUBS, DIAMONDS, HEARTS, SPADESet allant de deux à dix, puis Jack, Queen, King et Ace. Ainsi, les deux clubs ont un nombre de 1, les trois clubs un 2 .... L'as des clubs...