Hi,
To extract the p values alone:
lapply(split(dat1,dat1$date),function(x)
TukeyHSD(aov(value~name,data=x),"name",ordered=TRUE)[[1]][4])
#$`2008-10-05`
#[1] 0.1402197
#$`2008-11-05`
#[1] 0.004488184
#$`2008-12-05`
#[1] 0.5441151
A.K.
----- Original Message -----
From: arun <smartpink111 at yahoo.com>
To: raz <barvazduck at gmail.com>
Cc: R help <r-help at r-project.org>
Sent: Tuesday, November 13, 2012 12:57 PM
Subject: Re: [R] Tukey test for subgroups in a data frame
HI,
Try this:
dat1<-read.table(text="
? date name? value
? 10/5/2008????? L???? 26
? 10/5/2008????? L???? 30
? 10/5/2008????? L???? 54
? 10/5/2008????? L???? 25
? 10/5/2008???? M???? 18
? 10/5/2008???? M???? 21
? 10/5/2008???? M???? 29
? 10/5/2008???? M???? 17
? 11/5/2008????? L???? 32
? 11/5/2008????? L???? 35
? 11/5/2008????? L???? 44
? 11/5/2008????? L???? 28
? 11/5/2008???? M???? 19
? 11/5/2008???? M???? 22
? 11/5/2008???? M???? 18
? 11/5/2008???? M???? 16
? 12/5/2008????? L???? 24
? 12/5/2008????? L???? 33
? 12/5/2008????? L???? 35
? 12/5/2008????? L???? 28
? 12/5/2008???? M???? 29
? 12/5/2008???? M???? 23
? 12/5/2008???? M???? 39
? 12/5/2008???? M???? 42
? ",sep="",header=TRUE,stringsAsFactors=FALSE)
?dat1$date<-as.Date(dat1$date,format="%m/%d/%Y")
dat1<-within(dat1,{name<-factor(name)})
res<- lapply(split(dat1,dat1$date),function(x)
TukeyHSD(aov(value~name,data=x),"name",ordered=TRUE))
?res
#$`2008-10-05`
?# Tukey multiple comparisons of means
? #? 95% family-wise confidence level
?? # factor levels have been ordered
#Fit: aov(formula = value ~ name, data = x)
#$name
? #? diff?????? lwr????? upr???? p adj
#L-M 12.5 -5.501854 30.50185 0.1402197
#$`2008-11-05`
?# Tukey multiple comparisons of means
? #? 95% family-wise confidence level
?? # factor levels have been ordered
#Fit: aov(formula = value ~ name, data = x)
#$name
?# ? diff???? lwr????? upr???? p adj
#L-M?? 16 7.13522 24.86478 0.0044882
----------------------------------------
A.K.
----- Original Message -----
From: raz <barvazduck at gmail.com>
To: r-help at r-project.org
Cc:
Sent: Tuesday, November 13, 2012 11:27 AM
Subject: [R] Tukey test for subgroups in a data frame
Hello,
I have a data frame with the following columns:
"date","name","value"
the name is the same for each date
I would like to get TukeyHSD p-value for the differences of "value"
between
"name"s in each "date" separately I tried different ANOVA
(aov()) but can
only get either tukey by "name" or by "data" but not
"name by date"
any suggestions?
Thanks
--
\m/
??? [[alternative HTML version deleted]]
______________________________________________
R-help at r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
______________________________________________
R-help at r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.