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

How to preserve electronics (computers, tablets and phones) for hundreds of years Announcing the arrival of Valued Associate #679: Cesar Manara Planned maintenance scheduled April 23, 2019 at 23:30 UTC (7:30pm US/Eastern) The network's official Twitter account is up and running again. What content…How might modern humans leave a message for 50,000 years?How would an aquatic race develop computers?How do we make computers care for human life?Keeping supplies for 1000 yearsWhat will cell phones become in 100 years?What body armor protects against a laser?Could a society attain 1800s technology if limited to using copper and bronze?Electricity power (magic) and electronicsScavenging metal resources in post-apocalypse after 100 yearsWhy would precursors create devices that can survive and still work after hundreds of years?

155 수학 과학 기타 둘러보기 메뉴eh추가해eh문서를 완성해

해남서초등학교 학교 연혁 학교 동문 외부 링크 둘러보기 메뉴북위 34° 34′ 21″ 동경 126° 35′ 50″ / 북위 34.5726385° 동경 126.5972670°  / 34.5726385; 126.5972670해남서초등학교eh