Hi there, ============================x=c(145,40,40,120,180, 140,155,90,160,95, 195,150,205,110,160, 45,40,195,65,145, 195,230,115,235,225, 120,55,50,80,45 ) y2=c( rep(as.character(1),5), rep(as.character(2),5), rep(as.character(3),5), rep(as.character(4),5), rep(as.character(5),5), rep(as.character(6),5) ) crd2=data.frame(x,y2) model1=aov(x~y2,data=crd2) TukeyHSD(model1) ============================ I can do the HSD like that, Q1-but how do I do the Fisher's LSD? Q2-and how can I do that table in the attachement (page 8/9)? ( i dont know how to describe that thing, orginally produced in SAS) Q3-in the tukey output, can i ask to to show whether it's significant or not by indicating ***? (attachment page 6) Thanks! casperyc SAS output http://n4.nabble.com/file/n1583134/ws1sols.pdf ws1sols.pdf -- View this message in context: http://n4.nabble.com/Fisher-s-LSD-and-tukey-output-thing-tp1583134p1583134.html Sent from the R help mailing list archive at Nabble.com.