[R-sig-ME] p-values from lme::anova VS fixed-effects of lme
K Imran M
drki.musa at gmail.com
Sat Dec 10 10:57:05 CET 2016
Dear all,
Firstly, my apology if this question is trivial. This question is related
to this post at
http://stats.stackexchange.com/questions/7185/what-is-the-difference-between-using-aov-and-lme-in-analyzing-a-longitudinal
.
My I ask a question is: If I rerun these codes (lme)
#### LME ######
res <- lme(distance ~ age*Sex, random = ~ 1 | Subject, data = Orthodont)
summary(res)
################
I would obtain this
#####################################################
Fixed effects: distance ~ age * Sex
Value Std.Error DF t-value p-value
(Intercept) 16.340625 0.9813122 79 16.651810 0.0000
age 0.784375 0.0775011 79 10.120823 0.0000
SexFemale 1.032102 1.5374208 25 0.671321 0.5082
age:SexFemale -0.304830 0.1214209 79 -2.510520 0.0141
########################################################
Next, if I run 'anova(res)', I would obtain this
###ANOVA####
numDF denDF F-value p-value
(Intercept) 1 79 4123.156 <.0001
age 1 79 122.450 <.0001
Sex 1 25 9.292 0.0054
age:Sex 1 79 6.303 0.0141
###############
I do not understand why the p-values for age and age:Sex (from anova) are
similar with (lme::anova) BUT for sex: it is different (anova, p-val =
0.0054 vs lme::anova, p-val = 0.508).
Your kind attention is very much appreciated.
Thanks and best wishes,
Kamarul
--
Dr. Kamarul Imran Musa (MD MCommunityMed)
Associate Professor (Epidemiology and Biostatistics) &
Public Health Physician,
Address:
Dept of Community Medicine,
School of Medical Sciences,
Universiti Sains Malaysia,
16150 Kbg Kerian Kelantan
MALAYSIA
ResearcherID: http://www.researcherid.com/rid/N-3198-2015
Google-scholar: 'Kamarul Imran Musa' at https://goo.gl/D3o3y6
ORCID ID: orcid.org/0000-0002-3708-0628
ScopusID: 18634847200
Personal web: www.myanalytics.com.my
Alt emails : drkamarul at usm.my , k.musa at lancaster.ac.uk
[[alternative HTML version deleted]]
More information about the R-sig-mixed-models
mailing list