116 résultats trouvés

Aller à la recherche avancée

par Matthieu Lesnoff
11 Mai 2016, 06:57
Forum : Questions en cours
Sujet : tableau moyenne par modalité d'une variable qualitative
Réponses : 6
Vues : 4694

Re: tableau moyenne par modalité d'une variable qualitative

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...
par Matthieu Lesnoff
09 Mai 2016, 13:40
Forum : Questions en cours
Sujet : Pré traitement pour analyse spectrale et PLS
Réponses : 1
Vues : 1223

Re: Pré traitement pour analyse spectrale et PLS

Bonjour

- importation de fichiers asd et pré-traitement de spectres : package prospectr

- pls, pcr, etc. : package pls (par exemple)
par Matthieu Lesnoff
22 Avr 2015, 13:20
Forum : Questions en cours
Sujet : test unilateral sur une proportion
Réponses : 3
Vues : 2015

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...
par Matthieu Lesnoff
03 Fév 2012, 10:56
Forum : Questions en cours
Sujet : Résultats ANOVA hiérarchique
Réponses : 9
Vues : 5342

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)...
par Matthieu Lesnoff
02 Nov 2011, 11:21
Forum : Scripts et fonctions R
Sujet : Exemples glmm avec lme4
Réponses : 0
Vues : 9991

Exemples glmm avec lme4

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...
par Matthieu Lesnoff
19 Oct 2011, 11:46
Forum : Questions en cours
Sujet : taille des labels avec xplot()
Réponses : 1
Vues : 2030

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...
par Matthieu Lesnoff
14 Oct 2011, 18:48
Forum : Questions en cours
Sujet : gestion de chaines de caractères
Réponses : 6
Vues : 3644

Re: gestion de chaines de caractères

peut être regarder du côté de l'aide de la fonction grep() (?)
par Matthieu Lesnoff
10 Oct 2011, 10:42
Forum : Questions en cours
Sujet : tirage binomial
Réponses : 1
Vues : 2156

Re: tirage binomial

pour les tirages d'éléments avec ou sans remise, voir la fonction sample() de R
par Matthieu Lesnoff
27 Sep 2011, 17:07
Forum : Questions en cours
Sujet : GLMM
Réponses : 2
Vues : 2692

Re: GLMM

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...
par Matthieu Lesnoff
23 Sep 2011, 18:58
Forum : Questions en cours
Sujet : lme4 : facteur emboîté dans une interaction
Réponses : 7
Vues : 4407

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...
par Matthieu Lesnoff
20 Sep 2011, 23:21
Forum : Questions en cours
Sujet : lme4 : facteur emboîté dans une interaction
Réponses : 7
Vues : 4407

Re: lme4 : facteur emboîté dans une interaction

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' (...
par Matthieu Lesnoff
20 Sep 2011, 22:54
Forum : Questions en cours
Sujet : QQ-PLOT
Réponses : 6
Vues : 7813

Re: QQ-PLOT

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 :...
par Matthieu Lesnoff
22 Juin 2011, 08:58
Forum : Questions en cours
Sujet : Possibilité d'avoir plusieurs graphes à l'écran ?
Réponses : 6
Vues : 3592

Re: Possibilité d'avoir plusieurs graphes à l'écran ?

ou utiliser l'argument "split" de la fonction print.trellis(), mais il faut alors remplacer boxplot() par bwplot() (fonction du package lattice).

voir les exemples dans ?print.trellis
par Matthieu Lesnoff
17 Juin 2011, 17:52
Forum : Questions en cours
Sujet : Calcul intervalle de confiance de paramètres estimés
Réponses : 5
Vues : 4672

Re: Calcul intervalle de confiance de paramètres estimés

A partir de régressions non linéaires, j'essaye de calculer les intervalles de confiance de mes paramètres estimés. Pour calculer les intervalles de confiance je fais: confint(model_pocmat2) La fonction confint calcule des IC de type "profile likelihood" ; si besoin et en alternat...
par Matthieu Lesnoff
04 Avr 2011, 18:33
Forum : Questions en cours
Sujet : objet formula et modification
Réponses : 8
Vues : 3132

les méthodes terms, all.vars, etc fournissent la totalités des termes, y en a t'il une spécifique pour récupérer ces noms de variables? je ne pense pas (?) car "1|id" est aussi une variable, mais vous pouvez y arriver facilement par : > f <- formula(y~x+z+e+(1|id)) > z <- ...

Aller à la recherche avancée