Je passe en revue les vidéos du cours gratuit d' apprentissage automatique en ligne d' Andrew Ng à Stanford. Il décrit Gradient Descent comme un algorithme pour résoudre la régression linéaire et écrire des fonctions dans Octave pour l'exécuter. Vraisemblablement, je pourrais réécrire ces fonctions...
Sur la base d'un échantillon de temps de survie, je voudrais estimer la probabilité de survivre au temps , pour certains spécifiques , à l'aide de l'estimateur de Kaplan-Meier. Est-il possible de le faire en ? S'il vous plaît, notez que est pas nécessairement un temps
Considérez le biplot PCA suivant: library(mvtnorm) set.seed(1) x <- rmvnorm(2000, rep(0, 6), diag(c(5, rep(1,5)))) x <- scale(x, center=T, scale=F) pc <- princomp(x) biplot(pc) Il y a un tas de flèches rouges tracées, que signifient-elles? Je savais que la première flèche étiquetée "Var1"...
J'ai réalisé une expérience où j'ai élevé différentes familles provenant de deux populations sources différentes. Chaque famille a reçu l'un des deux traitements. Après l'expérience, j'ai mesuré plusieurs traits sur chaque individu. Pour tester un effet du traitement ou de la source ainsi que leur...
J'ai un ensemble de données sur les cas incidents par saison d'une maladie rare. Par exemple, disons qu'il y a eu 180 cas au printemps, 90 en été, 45 en automne et 210 en hiver. Je me demande s'il est approprié de joindre des erreurs standard à ces chiffres. Les objectifs de recherche sont...
J'ai quelques données dichotomiques, seulement des variables binaires, et mon patron m'a demandé d'effectuer une analyse factorielle en utilisant la matrice de corrélations tétrachoriques. J'ai déjà pu m'apprendre à exécuter différentes analyses sur la base des exemples ici et sur le site de...
Considérons un problème d'analyse de données classique où vous avez un résultat YiYiY_{i} et comment elle est liée à un certain nombre de facteurs prédictifs Xi1,...,XipXi1,...,XipX_{i1}, ..., X_{ip} . Le type d'application de base à l'esprit ici est que YiYiY_{i} est un résultat au niveau du...
J'ai les données suivantes et j'aimerais y adapter un modèle de croissance exponentielle négative: Days <- c( 1,5,12,16,22,27,36,43) Emissions <- c( 936.76, 1458.68, 1787.23, 1840.04, 1928.97, 1963.63, 1965.37, 1985.71) plot(Days, Emissions) fit <- nls(Emissions ~ a* (1-exp(-b*Days)),...
J'utilise le package nnet dans R pour tenter de construire un ANN pour prédire les prix immobiliers des condos (projet personnel). Je suis nouveau dans ce domaine et je n'ai pas de formation en mathématiques, alors s'il vous plaît, mettez-moi à nu. J'ai des variables d'entrée qui sont à la fois...
Je lisais sur le calcul de l'estimation non biaisée de l'écart-type et la source que j'ai lue a déclaré (...) sauf dans certaines situations importantes, la tâche a peu de pertinence pour les applications de la statistique car son besoin est évité par des procédures standard, telles que...
Je me demandais quel progiciel statistique de logiciels recommandez-vous pour effectuer l'inférence bayésienne. Par exemple, je sais que vous pouvez exécuter openBUGS ou winBUGS de manière autonome ou vous pouvez également les appeler à partir de R. Mais R a également plusieurs de ses propres...
mgcv::gamPage d'aide de lecture : des intervalles de confiance / crédibles sont facilement disponibles pour toute quantité prévue à l'aide d'un modèle ajusté Cependant, je ne peux pas trouver un moyen d'en obtenir un. Je pensais predict.gamavoir un paramètre type=confidenceet un levelmais ce n'est...
J'ai besoin d'estimer la fonction de densité sur la base d'un ensemble d'observations à l'aide de l'estimateur de densité du noyau. Sur la base du même ensemble d'observations, j'ai également besoin d'estimer les première et deuxième dérivées de la densité en utilisant les dérivées de l'estimateur...
J'essaie de tracer un QQ-plot avec deux ensembles de données d'environ 1,2 million de points, dans R (en utilisant qqplot et en introduisant les données dans ggplot2). Le calcul est assez facile, mais le graphique résultant est douloureusement lent à charger, car il y a tellement de points. J'ai...
J'ai un tas de données provenant de deux échantillons (témoins et traités), contenant chacun plusieurs milliers de valeurs qui doivent subir des tests de signification dans R. Théoriquement, les valeurs devraient être continues, mais en raison de l'arrondissement effectué par le logiciel de mesure,...
J'ai récemment implémenté un filtre de Kalman sur l'exemple simple de mesure d'une position de particules avec une vitesse et une accélération aléatoires. J'ai trouvé que le filtre de Kalman fonctionnait bien, mais je me suis alors demandé quelle était la différence entre cela et simplement faire...
Je viens de commencer à apprendre à utiliser Stan et rstan. À moins que je ne sois toujours confus sur le fonctionnement de JAGS / BUGS, je pensais que vous deviez toujours définir une distribution préalable d'une sorte pour chaque paramètre du modèle à partir duquel tirer. Il semble que vous...
J'ai un vecteur de nombres que j'ai téléchargé ici (... / code / MyData.Rdata) en utilisant dput. Je voudrais obtenir le bca ci donc j'ai écrit ce code: my.mean <- function(dat, idx){ return (mean(dat[idx], na.rm = TRUE)) } boot.out<-boot(data=my.data, statistic = my.mean, R=1000) Mais quand...
J'ai utilisé randomForest pour classer 6 comportements d'animaux (par exemple, debout, marcher, nager, etc.) en fonction de 8 variables (différentes postures corporelles et mouvements). Le MDSplot dans le package randomForest me donne cette sortie et j'ai des problèmes pour interpréter le résultat....
Je fais une validation croisée en utilisant la méthode du congé-un. J'ai une réponse binaire et j'utilise le package de démarrage pour R et la fonction cv.glm . Mon problème est que je ne comprends pas bien la partie "coût" de cette fonction. D'après ce que je peux comprendre, c'est la fonction qui...