J'ai le vecteur x <- c(1,2,3,4,5,5,5,6,6,6,6, 7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7, 7,7,7,7,7,7,7,7,8,8,8,8,9,9,9,10) (mon vecteur réel a une longueur> 10 000) et je voudrais trouver les intervalles où se trouvent les 90% de la densité. Est-ce quantile(x, probs=c(0.05,0.95), type=5)le plus...
J'ai le type de données suivant. J'ai évalué 10 individus chacun répété 10 fois. J'ai une matrice de relations 10x10 (relation entre toutes les combinaisons des individus). set.seed(1234) mydata <- data.frame (gen = factor(rep(1:10, each = 10)), repl = factor(rep(1:10, 10)), yld = rnorm(10, 5,...
Y a-t-il une différence entre les termes suivants ou ils sont identiques? Biais Biais systématique Erreurs systématiques S'il existe des différences, veuillez les expliquer. Ces erreurs peuvent-elles être réduites lorsque l'on augmente la taille de l'échantillon? MISE À JOUR: Mon domaine d'intérêt...
Dans R, comment spécifier le modèle lmer sans effet fixe global? Par exemple, si je dis quelque chose comme lmer(y ~ (1 | group) + (0 + x | group), data = my_df) le modèle équipé sera yje j= a +αje+βjeXje jyij=a+αi+βixijy_{ij} = a + \alpha_{i} + \beta_{i} x_{ij} Comment puis-je adapter le modèle...
J'essaie de visualiser certaines données de consommation, qui ont 4 catégories. Les utilisateurs sont libres de basculer entre les différentes catégories. Je voudrais visualiser les trois ou quatre derniers commutateurs pour chaque individu. Nous commencerions donc par un tracé avec une colonne...
J'ai des parcelles de 13 groupes que je montre dans une parcelle. Les groupes ont des populations déséquilibrées et ne sont pas normalement distribués. Je veux montrer quelles paires sont statistiquement similaires (c'est-à-dire, ont une valeur de p kruskal.test <0,05) en mettant a, b, c, etc....
Si j'ai une variable à 4 niveaux, en théorie, je dois utiliser 3 variables fictives. En pratique, comment cela se fait-il réellement? Dois-je utiliser 0-3, dois-je utiliser 1-3 et laisser les 4 vides? Aucune suggestion? REMARQUE: je vais travailler dans R. MISE À JOUR: Que se passerait-il si...
Disons que je génère la probabilité d'un résultat basé sur un certain facteur et trace la courbe de ce résultat. Existe-t-il un moyen d'extraire l'équation de cette courbe de R? > mod = glm(winner~our_bid, data=mydat, family=binomial(link="logit")) > summary(mod) Call: glm(formula = winner ~...
J'ai des données de séries chronologiques où la variable mesurée est des entiers positifs discrets (comptes). Je veux tester s'il y a une tendance à la hausse dans le temps (ou non). La variable indépendante (x) est comprise entre 0 et 500 et la variable dépendante (y) est comprise entre 0 et 8....
Fermé. Cette question est hors sujet . Il n'accepte pas actuellement les réponses. Vous souhaitez améliorer cette question? Mettez à jour la question afin qu'elle soit sur le sujet pour la validation croisée. Fermé il y a 9 mois . Le problème mentionné dans cette question est résolu dans la...
J'ai créé ma propre version légèrement améliorée du termplot que j'utilise dans cet exemple, vous pouvez le trouver ici . J'ai déjà posté sur SO mais plus j'y pense, plus je pense que cela est probablement plus lié à l'interprétation du modèle des risques proportionnels de Cox qu'au codage réel. Le...
Je cherche un exemple de code simple sur la façon d'exécuter un filtre à particules dans R. Le package pomp semble prendre en charge le bit mathématique de l'espace d'état, mais les exemples sont un peu difficiles à suivre par programme pour un développeur OO simple comme moi, en particulier...
Existe-t-il un principe général permettant de calculer la corrélation de Pearson pour deux variables aléatoires X et Y avant de prendre leur transformation logarithmique ou après? Existe-t-il une procédure de test qui est plus appropriée? Ils donnent des valeurs similaires mais différentes, car la...
Fermé. Cette question est hors sujet . Il n'accepte pas actuellement les réponses. Vous souhaitez améliorer cette question? Mettez à jour la question afin qu'elle soit sur le sujet pour la validation croisée. Fermé il y a 3 ans . Je suis très intéressé à apprendre à créer des rapports récurrents à...
En tant que passe-temps secondaire, j'ai exploré les séries chronologiques de prévision (en particulier, en utilisant R). Pour mes données, j'ai le nombre de visites par jour, pour chaque jour remontant à près de 4 ans. Dans ces données, il existe des modèles distincts: Du lundi au vendredi, il y a...
Disons que vous disposez de données de survie comme celle-ci: obs <- data.frame( time = c(floor(runif(100) * 30), floor((runif(100)^2) * 30)), status = c(rbinom(100, 1, 0.2), rbinom(100, 1, 0.7)), group = gl(2,100) ) Pour effectuer un test de classement de journal standard, on peut utiliser...
J'utilise le package kernlab dans R pour construire un SVM pour classer certaines données. Le SVM fonctionne bien en ce qu'il fournit des `` prédictions '' d'une précision décente, mais ma liste de variables d'entrée est plus longue que je ne le souhaiterais et je ne suis pas sûr de l'importance...
J'ai des problèmes pour utiliser les fonctions cor()et cor.test(). J'ai juste deux matrices (seulement des valeurs numériques, et le même nombre de lignes et de colonnes) et je veux avoir le numéro de corrélation et la valeur p correspondante. Lorsque j'utilise, cor(matrix1, matrix2)j'obtiens les...
J'ai besoin d'un package qui peut me donner l'équation d'un modèle SVM linéaire. Actuellement, j'utilise e1071 comme ceci: library(e1071) m = svm(data, labels, type='C', kernel='linear', cost=cost, probability=FALSE, scale=scale) w = t(m$coefs) %*% data[m$index,] #Weight vector b = -model$rho...
Je sais que c'est probablement une question fondamentale ... Mais je ne semble pas trouver la réponse. Je monte un GLM avec une famille Poisson, puis j'ai essayé de jeter un coup d'œil aux prédictions, mais le décalage semble être pris en considération: