Existe-t-il de toute façon que je peux exécuter LASSO avec une régression binomiale négative sur R? J'effectue une régression binomiale négative sur mon ensemble de données car les données sont trop dispersées pour imposer une régression de poisson. Pendant ce temps, je suis également confronté à...
J'ai le tableau suivant dans R df <- structure(list(x = structure(c(12458, 12633, 12692, 12830, 13369, 13455, 13458, 13515), class = "Date"), y = c(6080, 6949, 7076, 7818, 0, 0, 10765, 11153)), .Names = c("x", "y"), row.names = c("1", "2", "3", "4", "5", "6", "8", "9"), class = "data.frame")...
J'ai le modèle linéaire généralisé suivant. L'objet glmDVest modélisé comme une proportion de succès sur le total des essais. Les objets x_isont des variables continues. À quoi cela ressemble-t-il en notation mathématique? winp.glm = glm(glmDV ~ x1 + x2 + x3 + x4 + x5 + x6 + x7, data=myData,...
J'ai un ensemble de données sur lequel j'essaie de faire une régression, et échoue. La situation: Des milliers d'opérateurs de robots de combat se battent entre eux à l'aide de robots de combat. Certains robots de combat sont forts et puissants, et d'autres sont faibles; les plus forts gagnent plus...
Mon projet actuel peut m'obliger à construire un modèle pour prédire le comportement d'un certain groupe de personnes. l'ensemble de données de formation ne contient que 6 variables (id est uniquement à des fins d'identification): id, age, income, gender, job category, monthly spend dans laquelle...
J'ai environ 400 individus et> 10 000 points temporels chacun (résultats de simulation). J'aimerais pouvoir suivre leur évolution au fil du temps. Tracer tous les individus est trop compliqué, tracer la moyenne + -sd, min / max ou quantiles est trop peu d'informations à mon goût. Je me demande...
J'utilise la méthode arima du package de statistiques de R avec ma série temporelle de 17376 éléments. Mon objectif est d'obtenir la valeur du critère AIC, j'ai observé lors de mon premier test ceci: ts <- arima(serie[,1], order = c(2,1,1), seasonal = list(order=c(2,0,1),period = 24), method =...
Quelles sont les approches typiques de traitement des données incomplètes dans le filtre de Kalman? Je parle de la situation où certains éléments du vecteur observé sont manquants, distinct du cas où un vecteur observé entier est manquant. Une autre manière de penser serait que la dimension du...
Dans Bayesian Data Analysis , chapitre 13, page 317, deuxième paragraphe complet, dans les approximations modales et distributionnelles, Gelman et al. écrire: Si le plan est de résumer l'inférence par le mode postérieur de [le paramètre de corrélation dans une distribution normale bivariée], nous...
J'ai un problème avec l'interprétation des interactions à 2 et 3 voies dans lmer. Mon DV est la hauteur qui est une variable continue. Tous les IV sont des variables catégorielles. Le premier facteur est animal, rat ou lion. Le deuxième facteur est le sexe, masculin ou féminin. Le troisième facteur...
J'utilise lmerTestpour exécuter des modèles mixtes linéaires (LMM) pour obtenir les valeurs p. Cependant, dans les articles que j'ai écrits, ils montreront à la fois lme4et lmerTest. Ensuite, je ne sais pas si j'ai également citer lme4quand j'utilise uniquement lmerTestpour exécuter des modèles...
Il existe plusieurs méthodes pour faire des prévisions de séries chronologiques équidistantes (par exemple Holt-Winters, ARIMA, ...). Cependant, je travaille actuellement sur l'ensemble de données espacées irrégulières suivant, qui a un nombre variable de points de données par an et aucun...
Après avoir fait quelques recherches sur le sujet, j'ai remarqué un déficit surprenant de packages d'inférence et de bibliothèques qui reposent sur des méthodes de transmission de messages ou d'optimisation pour Python et R. À ma connaissance, ces méthodes sont extrêmement utiles. Par exemple, pour...
Je travaille sur des modèles de données de comptage gonflés à zéro en utilisant le psclpackage. Je me demande simplement pourquoi il n'y a pas de développement de modèles pour les modèles de données de comptage gonflés! Aussi pourquoi il n'y a pas de développement de modèles de données de comptage...
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é il y a 2 ans . ma série chronologique est évidemment périodique, mais la décomposition...
Je suis un peu novice en R et en sélection de fonctionnalités, et j'ai essayé le package Boruta pour sélectionner (diminuer) mon nombre de variables (n = 40). Je pensais que cette méthode prenait également en compte la corrélation possible entre les variables, cependant, deux (parmi les 20...
Supposons que j'ai une fonction telle que: f <- function(x){ exp(x) / (1 + exp(x)) } il est censé fonctionner pour toute valeur réelle de x, mais en fait, il renvoie NaN lorsque x est 710 ou plus. Je me demande quelle est la bonne façon de gérer ce problème. Je me rends compte qu'il est facile...
Étant donné la trame de données suivante: df <- data.frame(x1 = c(26, 28, 19, 27, 23, 31, 22, 1, 2, 1, 1, 1), x2 = c(5, 5, 7, 5, 7, 4, 2, 0, 0, 0, 0, 1), x3 = c(8, 6, 5, 7, 5, 9, 5, 1, 0, 1, 0, 1), x4 = c(8, 5, 3, 8, 1, 3, 4, 0, 0, 1, 0, 0), x5 = c(1, 1, 1, 1, 1, 0, 1, 0, 0, 0, 0, 0), x6 = c(2,...
J'utilise le paquet bnlearn en R pour apprendre la structure de mon réseau bayésien et ses paramètres. Ce que je veux faire, c'est "prédire" la valeur d'un nœud étant donné la valeur des autres nœuds comme preuve (évidemment, à l'exception du nœud dont nous prédisons les valeurs). J'ai des...
J'ai exécuté une courbe de survie de censure d'intervalle avec R, JMP et SAS. Ils m'ont tous deux donné des graphiques identiques, mais les tableaux différaient un peu. Voici le tableau que JMP m'a donné. Start Time End Time Survival Failure SurvStdErr . 14.0000 1.0000 0.0000 0.0000 16.0000 21.0000...