How to run two-way repeated within subject comparisons in R?Post Hoc tests for ANOVARepeated-Measures ANOVA: ezANOVA vs. aov vs. lme syntaxRepeated-measures ANOVA in R: +Error(subject) or + Error(subject/ VI1 * VI2)?How to exclude unwanted comparisons in two-way ANOVA in RHow to set up a 2 way repeated measures ANCOVA using R?post-hoc test for one-way ANOVA with random effectHow to perform Tukey HSD unconfounded comparisons with adjusted p-value in R or PythonRepeated Measures Two-Way Anova using Two Different ListsSpecify within-subjects and between-subjects ANOVA model using lme or lmer, as fixed-effectsR - post-hoc t-tests in mixed ANOVA design and triple interaction

Applications of pure mathematics in operations research

Why are we moving in circles with a tandem kayak?

Does Ubuntu reduce battery life?

What are the cons of stateless password generators?

Would people understand me speaking German all over Europe?

"Valet parking " or "parking valet"

Security measures that could plausibly last 150+ years?

How should I quote American English speakers in a British English essay?

Why don't short runways use ramps for takeoff?

How did astronauts using rovers tell direction without compasses on the Moon?

On the sensitivity conjecture?

What force enables us to walk? Friction or normal reaction?

Should I put my name first, or last in the team members list

Is it okay for me to decline a project on ethical grounds?

How close to the Sun would you have to be to hear it?

How can flights operated by the same company have such different prices when marketed by another?

Why does the Rust compiler not optimize code assuming that two mutable references cannot alias?

Complaints from (junior) developers against solution architects: how can we show the benefits of our work and improve relationships?

Rampant sharing of authorship among colleagues in the name of "collaboration". Is not taking part in it a death knell for a future in academia?

Is it possible for a particle to decay via gravity?

How can Paypal know my card is being used in another account?

How to prevent a single-element caster from being useless against immune foes?

How would a lunar colony attack Earth?

Do the books ever say oliphaunts aren’t elephants?



How to run two-way repeated within subject comparisons in R?


Post Hoc tests for ANOVARepeated-Measures ANOVA: ezANOVA vs. aov vs. lme syntaxRepeated-measures ANOVA in R: +Error(subject) or + Error(subject/ VI1 * VI2)?How to exclude unwanted comparisons in two-way ANOVA in RHow to set up a 2 way repeated measures ANCOVA using R?post-hoc test for one-way ANOVA with random effectHow to perform Tukey HSD unconfounded comparisons with adjusted p-value in R or PythonRepeated Measures Two-Way Anova using Two Different ListsSpecify within-subjects and between-subjects ANOVA model using lme or lmer, as fixed-effectsR - post-hoc t-tests in mixed ANOVA design and triple interaction






.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty margin-bottom:0;








1















I have this following data (3x3 within-subject design):



subj phase cope mean
1 A a ...
1 A b ...
1 A c ...
1 B a ...
1 B b ...
1 B c ...
1 C a ...
1 C b ...
1 C c ...
2 A a ...
2 A b ...
2 A c ...
2 B a ...
2 B b ...
2 B c ...
2 C a ...
2 C b ...
2 C c ...
...


Both phase and cope are categorical and the DV (mean) is continuous.



How do I run a 3-by-3 (ABCxabc) within subject ANOVA/lme with this data and what should I do for the post-hoc test?



After searching on the web, I have found the following 3 potential options:



1. baselinemodel <- lme(m ~ 1, random = ~1 | subj/phase/cope, data = csv_temp, method = "ML")
copemodel <- lme(m ~ cope, random = ~1 | subj/phase/cope, data = csv_temp, method = "ML")
phasemodel <- lme(m ~ phase, random = ~1 | subj/phase/cope, data = csv_temp, method = "ML")
fullmodel <- lme(m ~ phase + cope, random = ~1 | subj/phase/cope, data = csv_temp, method = "ML")
interactionmodel <- lme(m ~ phase * cope, random = ~1 | subj/phase/cope, data = csv_temp, method = "ML")
anova(baselinemodel, copemodel, phasemodel, fullmodel, interactionmodel)

2. m.lme <- lme(m ~ cope * phase, data = csv_temp, random = ~1 | subj)

3. m.lme <- lme(m ~ cope * phase, random = list(subj = pdBlocked(list(~1, pdIdent(~ cope - 1), pdIdent(~ phase - 1)))), data = csv_temp)
anova(m.lme)


They produced different results. Which way is the correct way? And how do I followup with post-hoc tests (at this time, neither method yield significant interactions)?










share|improve this question


























  • This might get a better response on stats.stackexchange.com

    – Stedy
    Mar 26 at 22:17

















1















I have this following data (3x3 within-subject design):



subj phase cope mean
1 A a ...
1 A b ...
1 A c ...
1 B a ...
1 B b ...
1 B c ...
1 C a ...
1 C b ...
1 C c ...
2 A a ...
2 A b ...
2 A c ...
2 B a ...
2 B b ...
2 B c ...
2 C a ...
2 C b ...
2 C c ...
...


Both phase and cope are categorical and the DV (mean) is continuous.



How do I run a 3-by-3 (ABCxabc) within subject ANOVA/lme with this data and what should I do for the post-hoc test?



After searching on the web, I have found the following 3 potential options:



1. baselinemodel <- lme(m ~ 1, random = ~1 | subj/phase/cope, data = csv_temp, method = "ML")
copemodel <- lme(m ~ cope, random = ~1 | subj/phase/cope, data = csv_temp, method = "ML")
phasemodel <- lme(m ~ phase, random = ~1 | subj/phase/cope, data = csv_temp, method = "ML")
fullmodel <- lme(m ~ phase + cope, random = ~1 | subj/phase/cope, data = csv_temp, method = "ML")
interactionmodel <- lme(m ~ phase * cope, random = ~1 | subj/phase/cope, data = csv_temp, method = "ML")
anova(baselinemodel, copemodel, phasemodel, fullmodel, interactionmodel)

2. m.lme <- lme(m ~ cope * phase, data = csv_temp, random = ~1 | subj)

3. m.lme <- lme(m ~ cope * phase, random = list(subj = pdBlocked(list(~1, pdIdent(~ cope - 1), pdIdent(~ phase - 1)))), data = csv_temp)
anova(m.lme)


They produced different results. Which way is the correct way? And how do I followup with post-hoc tests (at this time, neither method yield significant interactions)?










share|improve this question


























  • This might get a better response on stats.stackexchange.com

    – Stedy
    Mar 26 at 22:17













1












1








1


1






I have this following data (3x3 within-subject design):



subj phase cope mean
1 A a ...
1 A b ...
1 A c ...
1 B a ...
1 B b ...
1 B c ...
1 C a ...
1 C b ...
1 C c ...
2 A a ...
2 A b ...
2 A c ...
2 B a ...
2 B b ...
2 B c ...
2 C a ...
2 C b ...
2 C c ...
...


Both phase and cope are categorical and the DV (mean) is continuous.



How do I run a 3-by-3 (ABCxabc) within subject ANOVA/lme with this data and what should I do for the post-hoc test?



After searching on the web, I have found the following 3 potential options:



1. baselinemodel <- lme(m ~ 1, random = ~1 | subj/phase/cope, data = csv_temp, method = "ML")
copemodel <- lme(m ~ cope, random = ~1 | subj/phase/cope, data = csv_temp, method = "ML")
phasemodel <- lme(m ~ phase, random = ~1 | subj/phase/cope, data = csv_temp, method = "ML")
fullmodel <- lme(m ~ phase + cope, random = ~1 | subj/phase/cope, data = csv_temp, method = "ML")
interactionmodel <- lme(m ~ phase * cope, random = ~1 | subj/phase/cope, data = csv_temp, method = "ML")
anova(baselinemodel, copemodel, phasemodel, fullmodel, interactionmodel)

2. m.lme <- lme(m ~ cope * phase, data = csv_temp, random = ~1 | subj)

3. m.lme <- lme(m ~ cope * phase, random = list(subj = pdBlocked(list(~1, pdIdent(~ cope - 1), pdIdent(~ phase - 1)))), data = csv_temp)
anova(m.lme)


They produced different results. Which way is the correct way? And how do I followup with post-hoc tests (at this time, neither method yield significant interactions)?










share|improve this question
















I have this following data (3x3 within-subject design):



subj phase cope mean
1 A a ...
1 A b ...
1 A c ...
1 B a ...
1 B b ...
1 B c ...
1 C a ...
1 C b ...
1 C c ...
2 A a ...
2 A b ...
2 A c ...
2 B a ...
2 B b ...
2 B c ...
2 C a ...
2 C b ...
2 C c ...
...


Both phase and cope are categorical and the DV (mean) is continuous.



How do I run a 3-by-3 (ABCxabc) within subject ANOVA/lme with this data and what should I do for the post-hoc test?



After searching on the web, I have found the following 3 potential options:



1. baselinemodel <- lme(m ~ 1, random = ~1 | subj/phase/cope, data = csv_temp, method = "ML")
copemodel <- lme(m ~ cope, random = ~1 | subj/phase/cope, data = csv_temp, method = "ML")
phasemodel <- lme(m ~ phase, random = ~1 | subj/phase/cope, data = csv_temp, method = "ML")
fullmodel <- lme(m ~ phase + cope, random = ~1 | subj/phase/cope, data = csv_temp, method = "ML")
interactionmodel <- lme(m ~ phase * cope, random = ~1 | subj/phase/cope, data = csv_temp, method = "ML")
anova(baselinemodel, copemodel, phasemodel, fullmodel, interactionmodel)

2. m.lme <- lme(m ~ cope * phase, data = csv_temp, random = ~1 | subj)

3. m.lme <- lme(m ~ cope * phase, random = list(subj = pdBlocked(list(~1, pdIdent(~ cope - 1), pdIdent(~ phase - 1)))), data = csv_temp)
anova(m.lme)


They produced different results. Which way is the correct way? And how do I followup with post-hoc tests (at this time, neither method yield significant interactions)?







r






share|improve this question















share|improve this question













share|improve this question




share|improve this question








edited Mar 26 at 22:01







Jiang

















asked Mar 26 at 21:20









JiangJiang

62 bronze badges




62 bronze badges















  • This might get a better response on stats.stackexchange.com

    – Stedy
    Mar 26 at 22:17

















  • This might get a better response on stats.stackexchange.com

    – Stedy
    Mar 26 at 22:17
















This might get a better response on stats.stackexchange.com

– Stedy
Mar 26 at 22:17





This might get a better response on stats.stackexchange.com

– Stedy
Mar 26 at 22:17












0






active

oldest

votes










Your Answer






StackExchange.ifUsing("editor", function ()
StackExchange.using("externalEditor", function ()
StackExchange.using("snippets", function ()
StackExchange.snippets.init();
);
);
, "code-snippets");

StackExchange.ready(function()
var channelOptions =
tags: "".split(" "),
id: "1"
;
initTagRenderer("".split(" "), "".split(" "), channelOptions);

StackExchange.using("externalEditor", function()
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled)
StackExchange.using("snippets", function()
createEditor();
);

else
createEditor();

);

function createEditor()
StackExchange.prepareEditor(
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader:
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
,
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
);



);













draft saved

draft discarded


















StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f55366359%2fhow-to-run-two-way-repeated-within-subject-comparisons-in-r%23new-answer', 'question_page');

);

Post as a guest















Required, but never shown

























0






active

oldest

votes








0






active

oldest

votes









active

oldest

votes






active

oldest

votes




Is this question similar to what you get asked at work? Learn more about asking and sharing private information with your coworkers using Stack Overflow for Teams.







Is this question similar to what you get asked at work? Learn more about asking and sharing private information with your coworkers using Stack Overflow for Teams.



















draft saved

draft discarded
















































Thanks for contributing an answer to Stack Overflow!


  • Please be sure to answer the question. Provide details and share your research!

But avoid


  • Asking for help, clarification, or responding to other answers.

  • Making statements based on opinion; back them up with references or personal experience.

To learn more, see our tips on writing great answers.




draft saved


draft discarded














StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f55366359%2fhow-to-run-two-way-repeated-within-subject-comparisons-in-r%23new-answer', 'question_page');

);

Post as a guest















Required, but never shown





















































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown

































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown







Popular posts from this blog

Kamusi Yaliyomo Aina za kamusi | Muundo wa kamusi | Faida za kamusi | Dhima ya picha katika kamusi | Marejeo | Tazama pia | Viungo vya nje | UrambazajiKuhusu kamusiGo-SwahiliWiki-KamusiKamusi ya Kiswahili na Kiingerezakuihariri na kuongeza habari

SQL error code 1064 with creating Laravel foreign keysForeign key constraints: When to use ON UPDATE and ON DELETEDropping column with foreign key Laravel error: General error: 1025 Error on renameLaravel SQL Can't create tableLaravel Migration foreign key errorLaravel php artisan migrate:refresh giving a syntax errorSQLSTATE[42S01]: Base table or view already exists or Base table or view already exists: 1050 Tableerror in migrating laravel file to xampp serverSyntax error or access violation: 1064:syntax to use near 'unsigned not null, modelName varchar(191) not null, title varchar(191) not nLaravel cannot create new table field in mysqlLaravel 5.7:Last migration creates table but is not registered in the migration table

은진 송씨 목차 역사 본관 분파 인물 조선 왕실과의 인척 관계 집성촌 항렬자 인구 같이 보기 각주 둘러보기 메뉴은진 송씨세종실록 149권, 지리지 충청도 공주목 은진현