Forum de Statistiques
Vous souhaitez réagir à ce message ? Créez un compte en quelques clics ou connectez-vous pour continuer.
Le Deal du moment : -38%
Ecran PC gaming 23,8″ – ACER KG241Y P3bip ...
Voir le deal
99.99 €

contrastes sur modele lineaire lme

2 participants

Aller en bas

contrastes sur modele lineaire lme Empty contrastes sur modele lineaire lme

Message par fabien4944 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

Revenir en haut Aller en bas

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

Message par Nik 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 : 1606
Date d'inscription : 23/05/2008

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