Le livre de statistiques que je lis recommande l'oméga carré pour mesurer les effets de mes expériences. J'ai déjà prouvé en utilisant un plan de parcelle divisé (mélange de plans intra-sujets et inter-sujets) que mes facteurs intra-sujets sont statistiquement significatifs avec p <0,001 et F =...
J'examine quelques données de couverture génomique qui sont essentiellement une longue liste (quelques millions de valeurs) d'entiers, chacun indiquant dans quelle mesure (ou "profondément") cette position dans le génome est couverte. Je voudrais rechercher dans ces données des "vallées",...
Question: Avec une chaîne MCMC à 10 dimensions, disons que je suis prêt à vous remettre une matrice des tirages: 100 000 itérations (lignes) par 10 paramètres (colonnes), comment identifier au mieux les modes postérieurs? Je suis particulièrement préoccupé par plusieurs modes. Contexte:Je me...
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 2 ans . Je voudrais apprendre à utiliser Rcpp . J'ai parcouru les documents sur...
J'ajuste un modèle linéaire dans R en utilisant glmnet. Le modèle original (non régularisé) a été ajusté en utilisant lmet n'avait pas de terme constant (c'est-à-dire qu'il était sous la forme lm(y~0+x1+x2,data)). glmnetprend une matrice de prédicteurs et un vecteur de réponses. J'ai lu la...
Vous vous demandez si quelqu'un a traversé un package / une fonction dans R qui combinera les niveaux d'un facteur dont la proportion de tous les niveaux d'un facteur est inférieure à un certain seuil? Plus précisément, l'une des premières étapes de la préparation des données que je mène est de...
Existe-t-il un moyen de tracer la ligne de régression d'un modèle par morceaux comme celui-ci, autre que d'utiliser linespour tracer chaque segment séparément ou en utilisant geom_smooth(aes(group=Ind), method="lm", fill=FALSE)? m.sqft <- mean(sqft) model <-
Je comprends la prémisse de l'algorithme kNN pour les données spatiales. Et je sais que je peux étendre cet algorithme à utiliser sur n'importe quelle variable de données continue (ou données nominales avec distance de Hamming). Cependant, quelles stratégies sont utilisées pour traiter des données...
J'ai un ensemble de données avec beaucoup de zéros qui ressemble à ceci: set.seed(1) x <- c(rlnorm(100),rep(0,50)) hist(x,probability=TRUE,breaks = 25) Je voudrais tracer une ligne pour sa densité, mais la density()fonction utilise une fenêtre mobile qui calcule les valeurs négatives de x....
J'ai lu la description de la régression des crêtes dans Applied Linear Statistical Models , 5e Ed chapitre 11. La régression des crêtes est effectuée sur les données de graisse corporelle disponibles ici . Le manuel correspond à la sortie dans SAS, où les coefficients transformés en arrière sont...
J'envisage la construction d' interfaces et R MATLAB® à Ross Quinlan de C5.0 (pour ceux qui ne le connaissent pas, C5.0 est un algorithme d'arbre de décision et progiciel, une extension de C4.5 ), et je suis en train de avoir une idée des composants que j'aurais besoin d'écrire. La seule...
Cette question peut sembler très large, mais voici ce que je recherche. Je sais qu'il existe de nombreux excellents livres sur les méthodes économétriques et de nombreux excellents articles de présentation sur les techniques économétriques. Il existe même d'excellents exemples reproductibles...
J'ai du mal à comprendre un ou deux aspects du package de cluster. Je suis de près l'exemple de Quick-R , mais je ne comprends pas un ou deux aspects de l'analyse. J'ai inclus le code que j'utilise pour cet exemple particulier. ## Libraries library(stats) library(fpc) ## Data mydata =...
J'exécute une simulation sur R et un cluster d'ordinateurs et j'ai le problème suivant. Sur chacun des ordinateurs X, je lance: fxT2 <- function(i) runif(10) nessay <- 100 c(mclapply(1:nessay, fxT2), recursive=TRUE) Il y a 32 ordinateurs, chacun avec 16 cœurs. Cependant, environ 2% des...
Avant de soumettre ma méta-analyse, je veux faire un graphique en entonnoir pour tester l'hétérogénéité et le biais de publication. J'ai la taille d'effet groupée et les tailles d'effet de chaque étude, qui prennent des valeurs de -1 à +1. J'ai les tailles d'échantillon n1, n2 pour les patients et...
Comment dois-je syntaxer la rmafonction du package metafor afin d'obtenir des résultats dans l'exemple réel d'une petite méta-analyse suivante? (statistique aléatoire SMD à effet aléatoire) study, mean1, sd1, n1, mean2, sd2, n2 Foo2000, 0.78, 0.05, 20, 0.82, 0.07, 25 Sun2003, 0.74, 0.08, 30, 0.72,...
Je voudrais connaître les valeurs (x, y)utilisées dans le traçage plot(b, seWithMean=TRUE)dans le package mgcv . Est-ce que quelqu'un sait comment extraire ou calculer ces valeurs? Voici un exemple: library(mgcv) set.seed(0) dat <- gamSim(1, n=400, dist="normal", scale=2) b <- gam(y~s(x0),...
Une question de débutant ici. J'effectue actuellement une régression non paramétrique en utilisant le package np dans R. J'ai 7 fonctionnalités et en utilisant une approche par force brute, j'ai identifié les meilleures 3. Mais, bientôt j'aurai bien plus de 7 fonctionnalités! Ma question est de...
J'ai besoin de dessiner un graphique complexe pour l'analyse visuelle des données. J'ai 2 variables et un grand nombre de cas (> 1000). Par exemple (le nombre est 100 si la dispersion est moins "normale"): x <- rnorm(100,mean=95,sd=50) y <- rnorm(100,mean=35,sd=20) d <-...
Disons que je régresse Y sur X1 et X2, où X1 est une variable numérique et X2 est un facteur à quatre niveaux (A: D). Existe-t-il un moyen d'écrire la fonction de régression linéaire lm(Y ~ X1 + as.factor(X2))afin que je puisse choisir un niveau particulier de X2 - disons, B - comme ligne de...