Moi-même, j'utiliserais toujours la moyenne géométrique pour estimer une médiane lognormale. Cependant, dans le monde de l'industrie, l'utilisation de la médiane de l'échantillon donne parfois de meilleurs résultats. La question est donc la suivante: existe-t-il un intervalle / point de coupure à...
Quelqu'un connaît-il une méthode pour comparer deux variables avec un test du chi carré si les variables proviennent d'enquêtes différentes avec des svydesign()énoncés différents ? Je cherche à tester une différence dans une distribution variable sur deux vagues d'une enquête, mais l'...
J'ai effectué un regroupement de points de coordonnées (longitude, latitude) et trouvé des résultats surprenants et négatifs à partir de critères de regroupement pour le nombre optimal de clusters. Les critères sont tirés du clusterCrit()package. Les points que j'essaie de regrouper sur un...
J'ai un ensemble de données météorologiques quotidiennes, qui a, sans surprise, un effet saisonnier très fort. J'ai adapté un modèle ARIMA à cet ensemble de données en utilisant la fonction auto.arima du package de prévision. À ma grande surprise, la fonction n'applique aucune opération saisonnière...
Ma conception est la suivante. yyy est la réponse de Bernoulli x1x1x_1 est une variable continue x2x2x_2 est une variable catégorielle (facteur) à deux niveaux L'expérience est complètement dans les sujets. Autrement dit, chaque sujet reçoit chaque combinaison dex1x1x_1 et x2x2x_2. Il s'agit d'une...
Je veux calculer les erreurs standard d'une distribution hyperbolique ajustée. Dans ma notation, la densité est donnée par J'utilise le package HyperbolicDistr dans R. J'évalue les paramètres via la commande
J'ai quelques questions de base concernant l'ACP (analyse en composantes principales) et LDA (analyse discriminante linéaire): Dans l'ACP, il existe un moyen de calculer la proportion de la variance expliquée. Est-il également possible pour LDA? Si c'est le cas, comment? La sortie «Proportion de...
http://cran.r-project.org/web/packages/quadprog/quadprog.pdf Le package R ne quadprogsemble pouvoir résoudre le problème de programmation quadratique que lorsque la matrice est définie positive.DDD Cependant, il existe un cas où la matrice n'est pas définie positive. tel queDDD...
Est-il possible d'effectuer une sélection approximative entièrement bayésienne (1) d'hyper-paramètres (par exemple l'échelle de covariance) avec le code GPML, au lieu de maximiser la vraisemblance marginale (2)? Je pense que l'utilisation de méthodes MCMC pour résoudre les intégrales impliquant des...
Ainsi, il me semble que la fonction de pondération en lm donne aux observations plus de poids plus la valeur de «poids» de l'observation associée est grande, tandis que la fonction lme en lme fait exactement le contraire. Cela peut être vérifié avec une simulation simple. #make 3 vectors- c is used...
J'ai deux questions, Question 1: Comment puis-je montrer que la distribution postérieure est une distribution bêta si la probabilité est binomiale et l'a priori est une bêta Question 2: Comment les choix des paramètres antérieurs affectent-ils le postérieur? Ne devraient-ils pas tous être les...
Je sais qu'il est possible lorsque vous réalisez une régression linéaire multiple de calculer l'intervalle de confiance sur le R au carré et sur le R au carré ajusté. Est-ce que quelqu'un sait comment le faire avec
J'ai un ensemble de données qui contient quelques centaines de transactions provenant de trois fournisseurs opérant dans plus de 100 pays sur une période de trois ans. Nous avons constaté que le pays de vente n'est pas un facteur significatif dans les prix atteints (les produits sont des produits...
Cette question a été soulevée dans un autre fil que j'ai commencé, alors j'ai pensé que j'obtiendrais les opinions de plus de gens à ce sujet. Ma question est Est le résiduel, e, un estimateur de l'erreur, ϵϵ\epsilon? La raison que je demande est la suivante. Dans OLS, la variance des résidus,RSS(...
Ceci est un spin-off de cette question: comment comparer deux groupes avec plusieurs mesures pour chaque individu avec R? Dans les réponses (si j'ai bien compris), j'ai appris que la variance intra-sujet n'affecte pas les inférences faites sur les moyennes de groupe et il est correct de simplement...
J'ai le type de données suivant (codé en R): v.a = c('cat', 'dog', 'dog', 'goat', 'cat', 'goat', 'dog', 'dog') v.b = c(1, 2, 1, 2, 1, 2, 1, 2) v.c = c('blue', 'red', 'blue', 'red', 'red', 'blue', 'yellow', 'yellow') set.seed(12) v.d = rnorm(8) aov(v.a ~ v.b + v.c + v.d) # Error Je voudrais savoir...
J'essaie de faire du SVM à une classe dans R. J'ai essayé d'utiliser le package kernlab e1071 / ksvm. Mais je ne sais pas si je le fais correctement. Existe-t-il un exemple de travail pour SVM à une classe dans R? Aussi, Je donne une grande matrice de prédicteurs sous forme de X. Puisqu'il est...
J'utilise un CFA et j'obtiens de bons indices d'ajustement (CFI = .99, RMSEA = .01) pour une échelle unidimensionnelle. Cependant, lorsque je teste la cohérence interne, j'obtiens les s de Cronbach pauvres ( ). J'ai tout essayé, de la suppression des valeurs aberrantes à la suppression d'éléments,...
J'effectue des régressions quantiles dans R en utilisant le paquet quantreg. Mon ensemble de données comprend 12 328 observations allant de 0,12 à 330. Les points temporels de mes données ne sont pas exactement continus; toutes les données se trouvent dans l'une des quelques dizaines de casiers...
J'ai calculé les valeurs ICC pour deux groupes et je voudrais maintenant comparer les valeurs ICC pour déterminer si les groupes diffèrent dans leur répétabilité. Dans la littérature, les gens ont simplement utilisé des tests t pour comparer la répétabilité, mais je ne sais pas comment procéder....