Les posteurs les plus actifs de la semaine
Aytan
contrastes sur modele lineaire lme Vote_lcapcontrastes sur modele lineaire lme Voting_barcontrastes sur modele lineaire lme Vote_rcap 
c@ssoulet
contrastes sur modele lineaire lme Vote_lcapcontrastes sur modele lineaire lme Voting_barcontrastes sur modele lineaire lme Vote_rcap 
Eric Wajnberg
contrastes sur modele lineaire lme Vote_lcapcontrastes sur modele lineaire lme Voting_barcontrastes sur modele lineaire lme Vote_rcap 
Adddd
contrastes sur modele lineaire lme Vote_lcapcontrastes sur modele lineaire lme Voting_barcontrastes sur modele lineaire lme Vote_rcap 
gg
contrastes sur modele lineaire lme Vote_lcapcontrastes sur modele lineaire lme Voting_barcontrastes sur modele lineaire lme Vote_rcap 


contrastes sur modele lineaire lme

Aller en bas

contrastes sur modele lineaire lme Empty contrastes sur modele lineaire lme

Message par fabien4944 le Mar 18 Oct 2011 - 15:02

Bonjour

Je voudrais tester un contraste entre les interations (CS:A) et (CS:B) du modele lineaire aleatoire de la forme
Indice~A + CS+CP+CS:A+CP:B

CS et Cp sont seulement des variables binaires

j'utilise pour le modele lineaire le package nlme et je vois pas comment faire pour tester mon contrastes.

quelqu'un peut il me renseigner

merci d'avance


fabien4944

Nombre de messages : 8
Date d'inscription : 18/08/2011

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

contrastes sur modele lineaire lme Empty Re: contrastes sur modele lineaire lme

Message par Nik le Mar 18 Oct 2011 - 20:15

Salut,

Si j'ai bien compris : les contrastes sous R se gèrent sur les facteurs en eux-même avec la fonction contrasts().
Donc tu dois imposer ton contrasts à tes facteurs en construisant la matrice des contrastes qui t'intéresse :
un exemple
Code:
> test <- rep(c("mod1","mod2"),each=10)
> test <- as.factor(test)
> contrasts(test)
    mod2
mod1    0
mod2    1
> moncontrast <- matrix(c(1,0),nc=1)
> contrasts(test) <- moncontrast
> contrasts(test)
    [,1]
mod1    1
mod2    0
> model.matrix(~test)
  (Intercept) test1
1            1    1
2            1    1
3            1    1
4            1    1
5            1    1
6            1    1
7            1    1
8            1    1
9            1    1
10          1    1
11          1    0
12          1    0
13          1    0
14          1    0
15          1    0
16          1    0
17          1    0
18          1    0
19          1    0
20          1    0
attr(,"assign")
[1] 0 1
attr(,"contrasts")
attr(,"contrasts")$test
    [,1]
mod1    1
mod2    0

tu peux donc mettre ce que tu veux dans tes contrastes pour chaque vecteur.

Nik

Nombre de messages : 1605
Date d'inscription : 23/05/2008

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

Revenir en haut

- Sujets similaires

 
Permission de ce forum:
Vous ne pouvez pas répondre aux sujets dans ce forum