Nous avons un ensemble de données avec deux covariables et une variable de regroupement catégorique et nous voulons savoir s'il existe des différences significatives entre la pente ou l'ordonnée à l'origine parmi les covariables associées aux différentes variables de regroupement. Nous avons...
@whuber a montré comment simuler des résultats multivariés ( , y 2 ety1y1y_1y2y2y_2 ) pour un point dans le temps.y3y3y_3 Comme nous le savons, les données longitudinales se produisent souvent dans les études médicales. Ma question est de savoir comment simuler des résultats multivariés de mesures...
J'utilise le Root Mean Squared Error(RMSE) pour mesurer la précision des valeurs prédites à l'aide d'un modèle. Je comprends que la valeur retournée utilise les unités de mes mesures (plutôt qu'un pourcentage). Cependant, je voudrais citer mes valeurs en pourcentage. L'approche que j'ai adoptée...
Le modèle Bradley – Terry – Luce (BTL) indique que , où est la probabilité que l'objet soit jugé "meilleur", plus lourd, etc., que l'objet , et et sont des paramètres.pj i= l o gje t- 1( δj- δje)pjje=logjet-1(δj-δje)p_{ji} = logit^{-1}(\delta_j - \delta_i)pje
Mes données sont une série chronologique de la population occupée, L, et la période de temps, l'année. n.auto=auto.arima(log(L),xreg=year) summary(n.auto) Series: log(L) ARIMA(2,0,2) with non-zero mean Coefficients: ar1 ar2 ma1 ma2 intercept year 1.9122 -0.9567 -0.3082 0.0254 -3.5904 0.0074 s.e....
Je lis A. Agresti (2007), An Introduction to Categorical Data Analysis , 2nd. édition, et je ne sais pas si je comprends bien ce paragraphe (p.106, 4.2.1) (bien que cela devrait être facile): Dans le tableau 3.1 sur le ronflement et les maladies cardiaques du chapitre précédent, 254 sujets ont...
Je me demandais si quelqu'un avait de l'expérience avec la fonction souris, comme décrit dans Souris: Imputation multivariée par équations chaînées dans R (JSS 2011 45 (3))? J'ai un ensemble de données avec un certain nombre de variables, chacune avec différents degrés de données manquantes. Ma...
J'utilise le package randomForest dans R (R version 2.13.1, randomForest version 4.6-2) pour la régression et j'ai remarqué un biais significatif dans mes résultats: l'erreur de prédiction dépend de la valeur de la variable de réponse. Les valeurs élevées sont sous-prédites et les valeurs faibles...
Les greffes suivantes sont extraites de cet article . Je suis novice dans le bootstrap et j'essaie d'implémenter le bootstrap paramétrique, semi-paramétrique et non paramétrique pour le modèle mixte linéaire avec le R bootpackage. Code R Voici mon Rcode: library(SASmixed) library(lme4)...
Comment puis-je estimer les intervalles de confiance à 95% en utilisant le profilage des paramètres estimés en maximisant une fonction log-vraisemblance en utilisant optim dans R? Je sais que je peux estimer asymptotiquement la matrice de covariance en inversant la toile de jute , mais je crains...
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 2 ans . Le code suivant PredictNew <- predict (glm.fit, newdata =...
J'ai un ensemble de données médicales avec environ 200 variables. L'une des variables est un bio-marqueur (concentration d'une enzyme particulière). Sa distribution est asymétrique, et le problème est que les valeurs au-dessus d'un certain niveau sont censurées / coupées à ce niveau. Ainsi, alors...
J'ai un problème qui, je pense, devrait être simple, mais je n'arrive pas à le comprendre. Je regarde la pollinisation des graines, j'ai des plantes (n = 36) qui fleurissent en grappes, j'échantillon 3 grappes de fleurs de chaque plante et 6 gousses de graines de chaque grappe (18 gousses de...
J'essaie d'écrire un script R pour simuler l'interprétation des expériences répétées d'un intervalle de confiance à 95%. J'ai trouvé qu'il surestime la proportion de fois où la valeur réelle de la population d'une proportion est contenue dans l'IC à 95% de l'échantillon. Pas une grande différence -...
J'essaie d'analyser l'effet de l'année sur logInd variable pour un groupe particulier d'individus (j'ai 3 groupes). Le modèle le plus simple: > fix1 = lm(logInd ~ 0 + Group + Year:Group, data = mydata) > summary(fix1) Call: lm(formula = logInd ~ 0 + Group + Year:Group, data = mydata)...
Voici le résumé-sortie du modèle de Coxph que j'ai utilisé (j'ai utilisé R et la sortie est basée sur le meilleur modèle final, c'est-à-dire que toutes les variables explicatives significatives et leurs interactions sont incluses): coxph(formula = Y ~ LT + Food + Temp2 + LT:Food + LT:Temp2 +...
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 2 ans . R semble être capable de produire de jolis tracés récapitulatifs à...
J'ai ajusté le modèle ARIMA (5,1,2) en utilisant la auto.arima()fonction dans R et en regardant l'ordre, nous pouvons dire que ce n'est pas le meilleur modèle à prévoir. Si des valeurs aberrantes existent dans la série de données, quelle est la méthode pour adapter un modèle à ces...
J'exécute des statistiques sur 5 IVs (5 traits de personnalité, extraversion, agrément, conscience, névrosisme, ouverture) contre 3 DVs Attitude envers PCT, Attitude envers CBT, Attitude envers PCT vs CBT. J'ai également ajouté l'âge et le sexe pour voir quels sont les autres effets. Je teste pour...
J'ai une probabilité assez faible conduisant l'échantillonneur Metropolis-Hastings à se déplacer dans l'espace des paramètres de manière très irrégulière, c'est-à-dire qu'aucune convergence ne peut être atteinte quels que soient les paramètres de distribution de la proposition (dans mon cas, il est...