420 résultats trouvés

Aller à la recherche avancée

par Maxime Hervé
05 Juil 2018, 10:33
Forum : Questions en cours
Sujet : Package emmeans et fonction cld extraction
Réponses : 2
Vues : 90

Re: Package emmeans et fonction cld extraction

Bonjour, si j'ai bien compris et que vous souhaitez récupérer la dernière colonne du tableau renvoyé par cld() : cld(E)$".group" La colonne emmeans correspond bien aux moyennes (ajustées). Attention dans le cas d'un modèle avec fonction de lien il faut "back-transformer" ces moye...
par Maxime Hervé
18 Mai 2018, 06:40
Forum : Questions en cours
Sujet : comment créer un subset à plusieurs conditions ?
Réponses : 4
Vues : 186

Re: comment créer un subset à plusieurs conditions ?

Bonjour,
alternativement, pour raccourcir la commande surtout s'il y a plus de deux possibilités :

Code : Tout sélectionner

soucheBC <- subset(data,souche %in% c("B","C"))

Maxime
par Maxime Hervé
19 Mar 2018, 07:41
Forum : Questions en cours
Sujet : équivalent fonction TukeyHSD pour modèle linéaire mixte
Réponses : 3
Vues : 257

Re: équivalent fonction TukeyHSD pour modèle linéaire mixte

Florent Aubry a écrit :Il faut utiliser les packages multcomp (fonction glht puis summary) ou lsmeans (fonction lsmeans)

Plutôt que lsmeans, emmeans (la nouvelle implémentation de lsmeans, qui disparaîtra à terme).
Maxime
par Maxime Hervé
06 Aoû 2017, 06:31
Forum : Questions en cours
Sujet : comparaison post hoc après modèle linéaire mixte
Réponses : 2
Vues : 454

Re: comparaison post hoc après modèle linéaire mixte

Bonjour,
l'erreur est que vous avez donné votre tableau de données en argument à lsmeans(), alors que la fonction attend un modèle. C'est donc plutôt :

Code : Tout sélectionner

lsmeans(m1,pairwise~Groupe)

Maxime
par Maxime Hervé
23 Juin 2017, 15:50
Forum : Questions en cours
Sujet : Comparaison moyennes suivant un facteur
Réponses : 1
Vues : 339

Re: Comparaison moyennes suivant un facteur

Bonjour,
essayez donc ça :

Code : Tout sélectionner

apply(typo[,2:12],MARGIN=2,function(x) pairwise.perm.t.test(x,typo[,1]))

Maxime
par Maxime Hervé
02 Mar 2017, 14:36
Forum : Questions en cours
Sujet : fonction set.seed()
Réponses : 17
Vues : 16343

Re: fonction set.seed()

Ola j'aurais dû préciser, je ne parlais que de situations simples du type générer une série de données pour un exemple comme plus haut dans le post. Pas d'une véritable étude basée sur une simulation, où l'on génère un grand nombre de séries de données. Dans un cas pareil bien sûr que je ne prendrai...
par Maxime Hervé
02 Mar 2017, 13:27
Forum : Questions en cours
Sujet : fonction set.seed()
Réponses : 17
Vues : 16343

Re: fonction set.seed()

Bonjour,

l'idée est que le résultat soit "aléatoire", donc le choix du nombre doit l'être lui-même. Personnellement je prends toujours l'heure qu'il est, en l'occurrence 14h27 donc 1427.

Maxime
par Maxime Hervé
30 Jan 2017, 15:00
Forum : Questions en cours
Sujet : Modele effets imbriqués
Réponses : 2
Vues : 359

Re: Modele effets imbriqués

Bonjour,
en utilisant le package lsmeans ça donnerait :

Code : Tout sélectionner

lsm <- lsmeans(model_AC,~Operator|Lab)
contrast(lsm,"pairwise")

La syntaxe Operator|Lab précisant que vous souhaitez faire des comparaisons multiples entre niveaux d'Operator au sein de chaque niveau de Lab.
Maxime
par Maxime Hervé
13 Déc 2016, 14:58
Forum : Questions en cours
Sujet : Obtenir Rcmdr en Syrie
Réponses : 4
Vues : 576

Re: Obtenir Rcmdr en Syrie

Une question est : Comment peut-elle récupérer pbkrtest et car ? Serait-il possible par exemple, de lui envoyer des fichiers zippés qu'elle pourrait utiliser ? Oui c'est possible en théorie, vous pouvez télécharger les fichiers zip sur le site du CRAN. Par contre pour pbkrtest il va vous falloir al...
par Maxime Hervé
13 Déc 2016, 14:29
Forum : Questions en cours
Sujet : Obtenir Rcmdr en Syrie
Réponses : 4
Vues : 576

Re: Obtenir Rcmdr en Syrie

Bonjour, la dernière version du package pbkrtest (celle que vous essayez de télécharger) exige une version de R >= 3.2.3. Votre collègue a peut-être une version plus ancienne ? Cela expliquerait que le package n'est pas trouvé. Dans ce cas là la solution est d'installer une version plus récente de R...
par Maxime Hervé
03 Nov 2016, 14:40
Forum : Questions en cours
Sujet : Erreur cercle de corrélation RDA
Réponses : 3
Vues : 553

Re: Erreur cercle de corrélation RDA

La situation s'éclaircit. Il semble qu'au moins une de vos variables à expliquer soit constante. Pourriez-vous le vérifier ?

Si c'est le cas, supprimez la ou les variables en question de l'analyse. De toute façon elle ne sert à rien puisque ne variant pas, elle n'apporte aucune information.

Maxime
par Maxime Hervé
03 Nov 2016, 08:45
Forum : Questions en cours
Sujet : Erreur cercle de corrélation RDA
Réponses : 3
Vues : 553

Re: Erreur cercle de corrélation RDA

Bonjour,

il doit y avoir au moins un coefficient de corrélation qui ne peut pas se calculer. Pouvez-vous afficher le résultat de MVA.cor(RDA,set=2) ?

Maxime
par Maxime Hervé
17 Oct 2016, 08:48
Forum : Questions en cours
Sujet : comment ajouter les facteurs à une RDA
Réponses : 5
Vues : 565

Re: comment ajouter les facteurs à une RDA

Guillaume Foursov a écrit :R studio ne le reconnait pas quand je l'appelle avec library(RVAideMemoire)


Vous voulez dire qu'il ne trouve pas la fonction ? Dans ce cas c'est que vous devez mettre à jour le package.

Maxime
par Maxime Hervé
14 Oct 2016, 16:02
Forum : Questions en cours
Sujet : comment ajouter les facteurs à une RDA
Réponses : 5
Vues : 565

Re: comment ajouter les facteurs à une RDA

Bonjour,
je ne suis pas sûr de totalement saisir la question, mais si vous souhaitez représenter un facteur à la fois (ou une interaction entre plusieurs facteurs) il y a la fonction MVA.plot() du package RVAideMemoire qui gère un assez grand nombre d'analyses multivariées.
Maxime
par Maxime Hervé
28 Sep 2016, 12:38
Forum : Questions en cours
Sujet : FactoMineR: variables centrées-réduites
Réponses : 2
Vues : 2211

Re: FactoMineR: variables centrées-réduites

Bonjour,

une réduction qui utilise n-1 et l'autre n ?

Maxime

Aller à la recherche avancée