I would use aov( ) mod<-aov(Score~Group, data=x) summary(mod) JJ -- View this message in context: http://n4.nabble.com/t-test-tp1470539p1470558.html Sent from the R help mailing list archive at Nabble.com.