Salut,
Je regarde une étude qui a été réalisée en R, j'ai réussi à reproduire les résultats du VCA en R en utilisant la syntaxe suivante, VisDil1 <- lmer (l2t ~ (1 | Op / Day / MDCK), Dil1), pour les résultats suivants
Effets aléatoires:
Nom des groupes Variance Std.Dev.
MDCK: (Jour: Op) (Interception) 0,02309 0,15195
Jour: Op (Interception) 0,00861 0,09279
Op (Interception) 0,00000 0,00000
quand je lance ce que je crois être le même modèle dans JMP, je reçois
Effet aléatoire Var Ratio Var composante Pct du total
MDCK [Jour, Op] 0,1904762 0,031746 16 000
Jour [Op] 0 0 0,000
Op 0 0 0,000
Résiduel 0.1666667 84.000
Total 0,1984127 100 000
J'ai utilisé le script suivant
Fit Model(
Y( :l2t ),
Effects,
Random Effects( :MDCK[:Day, :Op], :Day[:Op], :Op ),
NoBounds( 0 ),
Estimate Only Variance Components( 1 ),
Personality( "Standard Least Squares" ),
Method( "REML" ),
Run Model
);
J'ai regardé cela à travers 8 dilutions, et c'est le même résultat dans chaque dilution, R a une composante de variance pour Day: Op, tandis que jmp combine cet effet dans l'effet MDCK [Day, Op], donc je me sens i Je ne spécifie pas le modèle dans JMP de la même manière que dans R pour obtenir le même résultat.
Quelqu'un a-t-il des conseils ???
Merci
Veuillez clarifier la structure de la conception et la randomisation.
Op est au sommet de la hiérarchie.
Les 4 jours sont différents pour les deux op.
Les deux MDCK sont différents pour chaque jour et Op.
Le titre n'est-il pas un effet fixe potentiellement important?
Salut Mark,
Merci pour votre réponse, veuillez trouver ci-joint une copie de la conception imbriquée, je suis d'accord que la syntaxe dosent reflète la conception de l'expérience, je pense qu'un modèle plus approprié serait
RBC
Jour [RBC]
Op [RBC, jour]
Cependant, c'est une autre histoire, ce que j'essaie de comprendre, c'est pourquoi lorsque je lance ce que je crois être le même modèle à la fois dans R et JMP, j'obtiens un composant de variance pour Day dans R, mais JMP semble combiner ce composant dans le MDCK variable.
La variable l2t est le log2 du titre, la variable indépendante.
Merci
