J'ai le modèle suivant: > model1<-lmer(aph.remain~sMFS1+sAG1+sSHDI1+sbare+season+crop +(1|landscape),family=poisson) ... et voici la sortie récapitulative. > summary(model1) Generalized linear mixed model fit by the Laplace approximation Formula: aph.remain ~ sMFS1 + sAG1 + sSHDI1 + sbare...
J'essaie de comprendre comment utiliser l'apprentissage automatique pour prédire la série financière 1 étape ou plus dans le futur. J'ai une série temporelle financière avec des données descriptives et je voudrais former un modèle et ensuite utiliser le modèle pour prédire n étapes à venir. Ce que...
Je voudrais comprendre ce que fait le code suivant. La personne qui a écrit le code ne travaille plus ici et il est presque complètement sans papiers. Une personne qui pense que " c'est un modèle de régression logistique bayésienne " m'a demandé d'enquêter bglm <- function(Y,X) { # Y is a vector...
Le package MARSS en R offre une fonction pour l'analyse factorielle dynamique. Dans ce package, le modèle de facteur dynamique est écrit comme une forme spéciale de modèle d'espace d'état et ils supposent que les tendances communes suivent le processus AR (1). Comme je ne connais pas très bien ces...
Quels sont les paramètres start, etastart, mustartdans la fonction GLM () ? J'ai cherché dans les documents et sur Internet, mais je n'ai trouvé aucune explication claire de ce que cela signifie. Cela ressemble aux "valeurs initiales" bayésiennes pour les chaînes, mais je doute que cela soit lié,...
Comment la matrice d'erreur var / cov est-elle calculée par les progiciels d'analyse statistique dans la pratique? Cette idée m'est claire en théorie. Mais pas en pratique. Je veux dire, si j'ai un vecteur de variables aléatoires , je comprends que la matrice de variance / covariance recevra le...
J'ai effectué une régression logistique multivariée avec la variable dépendante Yétant le décès dans une maison de soins infirmiers dans une certaine période d'entrée et j'ai obtenu les résultats suivants (notez que si les variables commencent dans Ac'est une valeur continue tandis que celles qui...
J'ai un ensemble de données de série chronologique auquel j'essaie d'adapter un modèle de Markov caché (HMM) afin d'estimer le nombre d'états latents dans les données. Mon pseudo-code pour ce faire est le suivant: for( i in 2 : max_number_of_states ){ ... calculate HMM with i states ......
J'adapte un modèle d'effets aléatoires glmerà certaines données d'entreprise. L'objectif est d'analyser la performance commerciale par distributeur en tenant compte des variations régionales. J'ai les variables suivantes: distcode: ID distributeur, avec environ 800 niveaux region: ID géographique...
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 6 ans . J'utilise des fichiers texte pour stocker mes données pour R sans aucun...
Wikipédia est-il faux ... ou je ne le comprends pas? Wikipédia: Les carrés blancs et noirs ("motif d'échecs") sont parfaitement dispersés, donc je serais de -1 pour Moran. Si les carrés blancs étaient empilés sur la moitié du plateau et les carrés noirs sur l'autre, le Moran I serait proche de +1....
Édition majeure: Je voudrais dire un grand merci à Dave et Nick jusqu'à présent pour leurs réponses. La bonne nouvelle est que j'ai réussi à faire fonctionner la boucle (principe emprunté à la publication du professeur Hydnman sur la prévision par lots). Pour consolider les requêtes en attente: a)...
J'ai du code et des sorties, et je voudrais construire un modèle. Je ne sais pas comment construire un modèle en utilisant cette sortie: require("splines") x <- c(0.2, 0.23, 0.26, 0.29, 0.33, 0.46, 0.53 ) y <- c(0.211, 0.2026, 0.2034, 0.2167, 0.2177, 0.19225, 0.182) fit <- lm(y ~ ns(x,3))...
J'avais l'impression que la fonction lmer()dans le lme4package ne produisait pas de valeurs p (voir lmer, valeurs p et tout ça ). J'ai utilisé des valeurs p générées par MCMC à la place selon cette question: effet significatif dans le lme4modèle mixte et cette question: impossible de trouver des...
Je ne suis pas mathématicien. J'ai recherché sur Internet KL Divergence. Ce que j'ai appris, c'est que la divergence KL mesure les informations perdues lorsque nous approchons la distribution d'un modèle par rapport à la distribution d'entrée. Je les ai vues entre deux distributions continues ou...
J'ai du mal à trouver un moyen de calculer la valeur de p pour la zone sous une caractéristique d'opérateur de récepteur (ROC). J'ai une variable continue et un résultat de test de diagnostic. Je veux voir si AUROC est statistiquement significatif. J'ai trouvé de nombreux packages traitant des...
J'ai un ensemble de données qui représente 1000 documents et tous les mots qui y apparaissent. Les lignes représentent donc les documents et les colonnes les mots. Ainsi, par exemple, la valeur dans la cellule représente la fois où le mot apparaît dans le document . Maintenant, je dois trouver des...
J'utilise le package caret pour entraîner un objet randomForest avec 10x10CV. library(caret) tc <- trainControl("repeatedcv", number=10, repeats=10, classProbs=TRUE, savePred=T) RFFit <- train(Defect ~., data=trainingSet, method="rf", trControl=tc, preProc=c("center", "scale")) Après cela, je...
Question Les résultats des tests de trois groupes de personnes sont enregistrés en tant que vecteurs distincts dans R. set.seed(1) group1 <- rnorm(100, mean = 75, sd = 10) group2 <- rnorm(100, mean = 85, sd = 10) group3 <- rnorm(100, mean = 95, sd = 10) Je veux savoir s'il y a une...
Pour un modèle linéaire gaussien où est supposé se trouver dans un espace vectoriel et a la distribution normale standard sur , la statistique du test pour où est un espace vectoriel, est une fonction un à un croissante de la statistique de déviance : Comment savoir que cette statistique fournit le...