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

Swift 4 - func physicsWorld not invoked on collision? The Next CEO of Stack OverflowHow to call Objective-C code from Swift#ifdef replacement in the Swift language@selector() in Swift?#pragma mark in Swift?Swift for loop: for index, element in array?dispatch_after - GCD in Swift?Swift Beta performance: sorting arraysSplit a String into an array in Swift?The use of Swift 3 @objc inference in Swift 4 mode is deprecated?How to optimize UITableViewCell, because my UITableView lags

Access current req object everywhere in Node.js ExpressWhy are global variables considered bad practice? (node.js)Using req & res across functionsHow do I get the path to the current script with Node.js?What is Node.js' Connect, Express and “middleware”?Node.js w/ express error handling in callbackHow to access the GET parameters after “?” in Express?Modify Node.js req object parametersAccess “app” variable inside of ExpressJS/ConnectJS middleware?Node.js Express app - request objectAngular Http Module considered middleware?Session variables in ExpressJSAdd properties to the req object in expressjs with Typescript