Skip to content

Instantly share code, notes, and snippets.

@s-n-1-0
Created December 22, 2023 10:20
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save s-n-1-0/dac385c3e858e3d5cb4920a24cf6bade to your computer and use it in GitHub Desktop.
Save s-n-1-0/dac385c3e858e3d5cb4920a24cf6bade to your computer and use it in GitHub Desktop.
2要因混合計画ANOVA結果
try({anovakun(df, "AsB", 2, 3)})
[ AsB-Type Design ]
This output was generated by anovakun 4.8.9 under R version 4.3.2.
It was executed on Fri Dec 22 19:19:28 2023.
<< DESCRIPTIVE STATISTICS >>
------------------------------
A B n Mean S.D.
------------------------------
a1 b1 30 5.5082 1.0298
a1 b2 30 5.6037 0.9004
a1 b3 30 5.5916 1.1800
a2 b1 30 5.4325 0.8130
a2 b2 30 5.9707 1.0746
a2 b3 30 6.3453 0.9330
------------------------------
<< SPHERICITY INDICES >>
== Mendoza's Multisample Sphericity Test and Epsilons ==
-------------------------------------------------------------------------
Effect Lambda approx.Chi df p LB GG HF CM
-------------------------------------------------------------------------
B 0.1527 3.6478 5 0.6012 ns 0.5000 0.9981 1.0337 1.0300
-------------------------------------------------------------------------
LB = lower.bound, GG = Greenhouse-Geisser
HF = Huynh-Feldt-Lecoutre, CM = Chi-Muller
<< ANOVA TABLE >>
---------------------------------------------------------
Source SS df MS F-ratio p-value
---------------------------------------------------------
A 5.4600 1 5.4600 5.0517 0.0284 *
s x A 62.6873 58 1.0808
---------------------------------------------------------
B 7.6284 2 3.8142 4.0274 0.0204 *
A x B 5.1672 2 2.5836 2.7280 0.0695 +
s x A x B 109.8598 116 0.9471
---------------------------------------------------------
Total 190.8027 179 1.0659
+p < .10, *p < .05, **p < .01, ***p < .001
<< POST ANALYSES >>
< MULTIPLE COMPARISON for "B" >
== Shaffer's Modified Sequentially Rejective Bonferroni Procedure ==
== The factor < B > is analysed as dependent means. ==
== Alpha level is 0.05. ==
--------------------------
B n Mean S.D.
--------------------------
b1 60 5.4703 0.9207
b2 60 5.7872 1.0002
b3 60 5.9685 1.1210
--------------------------
----------------------------------------------------------
Pair Diff t-value df p adj.p
----------------------------------------------------------
b1-b3 -0.4981 2.8487 58 0.0061 0.0182 b1 < b3 *
b1-b2 -0.3169 1.7476 58 0.0858 0.0858 b1 = b2
b2-b3 -0.1812 1.0253 58 0.3095 0.3095 b2 = b3
----------------------------------------------------------
< SIMPLE EFFECTS for "A x B" INTERACTION >
--------------------------------------------------------------------------
Effect Lambda approx.Chi df p LB GG HF CM
--------------------------------------------------------------------------
B at a1 0.5408 1.1869 2 0.5524 ns 0.5000 0.9602 1.0268 1.0126
B at a2 0.2949 2.3581 2 0.3076 ns 0.5000 0.9253 0.9856 0.9720
--------------------------------------------------------------------------
LB = lower.bound, GG = Greenhouse-Geisser
HF = Huynh-Feldt-Lecoutre, CM = Chi-Muller
----------------------------------------------------------
Source SS df MS F-ratio p-value
----------------------------------------------------------
A at b1 0.0860 1 0.0860 0.0999 0.7531 ns
Er at b1 49.9230 58 0.8607
----------------------------------------------------------
A at b2 2.0201 1 2.0201 2.0556 0.1570 ns
Er at b2 57.0004 58 0.9828
----------------------------------------------------------
A at b3 8.5211 1 8.5211 7.5311 0.0081 **
Er at b3 65.6238 58 1.1314
----------------------------------------------------------
B at a1 0.1624 2 0.0812 0.0826 0.9209 ns
s x B at a1 57.0478 58 0.9836
----------------------------------------------------------
B at a2 12.6332 2 6.3166 6.9371 0.0020 **
s x B at a2 52.8120 58 0.9106
----------------------------------------------------------
+p < .10, *p < .05, **p < .01, ***p < .001
< MULTIPLE COMPARISON for "B at a2" >
== Shaffer's Modified Sequentially Rejective Bonferroni Procedure ==
== The factor < B at a2 > is analysed as dependent means. ==
== Alpha level is 0.05. ==
----------------------------------------------------------
Pair Diff t-value df p adj.p
----------------------------------------------------------
b1-b3 -0.9128 4.3748 29 0.0001 0.0004 b1 < b3 *
b1-b2 -0.5382 2.0230 29 0.0524 0.0524 b1 = b2
b2-b3 -0.3746 1.4388 29 0.1609 0.1609 b2 = b3
----------------------------------------------------------
output is over --------------------///
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment