de manière générale, éviter rbind dès qu'il y a beaucoup d'itérations Préférer "rbindlist" (très rapide) du package data.table (qu'il faut installer) n <- 3 ; p <- 5 ; nrep <- 4 # List of data.frames mylist <- lapply(1:nrep, function(x) as.data.frame(matrix(runif(n * p), nrow = n))) # tran...
Bonjour La fonction 'aggregate' du package de base de R est très pratique car elle renvoie directement un data.frame : res <- aggregate(v1 ~ f1, data = dat, FUN = mean) # on peut spécifier plusieurs variables, plusieurs facteurs # et n'importe quelle autre fonction que la moyenne, p.ex. médiane ci-d...
P=p0 vs p<p0 Mathématiquement oui la statistique du test unilatéral classique est définie sous H0 p = p0, même si intuitivement on poserait plutôt "p >= p0". Mais ce n'est qu'une nuance : en pratique dans ce type de test, on cherche surtout à savoir si H1, par exemple p < p0, est vraisemb...
en complément : dans le cas d'un modèle hiérarchique A/B où B est aléatoire, et lorsque le plan est équlibré, une manière simple de tester A en évitant les package nlme ou lme4 est d'utiliser aov : fm <- aov(y ~ A + Error(B), data = x) summary(fm) # même chose (pour A)...
Quelques scripts pouvant aider pour des analyses glmm avec lme4. Je me suis largement inspiré des informations disponibles sur le site http://glmm.wikidot.com . Certains des scripts seront vraisemblablement bientôt obsolètes, car lme4 évolue et fournira sûrement des fonctions plus efficaces. Par ail...
une approche est de placer "cex" dans la liste "x" (ou "y") de la liste de l'argument "scales", p.ex. : z <- c("J", "A", "S", "O", "N", "D", "J", "F", "M", "A&qu...
Salut, quel est le package meilleur selon vous pour les GLMM? Sara Pour les méthodes ML, deux standards sont lme4 et glmmML. lme4 permet d'ajuster plusieurs facteurs aléatoires emboités ou des facteurs croisés, glmmML non (1 seul facteur). Il ya aussi glmm.admb, + les packages utilisant les approch...
Ma remarque était juste par rapport à : je voudrais emboîter le facteur "dil" dans l'interaction "Opérateur*Equipement"" Sous lme4, "(1 | A/B") est équivalent à "(1 | A) + (1 | A:B)", c'est à dire effet aléatoire "A" + effet aléatoire "B em...
Bonjour, Dans l'exemple ci-dessous je voudrais emboîter le facteur "dil" dans l'interaction "Opérateur*Equipement" : SAS n'a pas forcément la même syntaxe que lme4. Sous lme4, '(1|dil/Operator:Equipement)' signifie : effet 'dil' + effet 'Operator:Equipement' emboité dans 'dil' (...
Je cherche à tracer un qq-plot avec ma distribution empirique et une loi de pareto Un exemple de code avec qqmath de lattice et une fonction construite qchisq01 (mélange de lois de Chi2). Pour votre cas, il suffit de remplacer qchisq01 par la fonction définissant les quantiles de la loi de Pareto :...