####ANOVA model > mc<-aov(rating~line+rep+plant %in% rep, data=c) > summary(mc) Df Sum Sq Mean Sq F value Pr(>F) line 27 93931 3479 8.709 <2e-16 *** rep 1 1230 1230 3.080 0.0809 . rep:plant 6 2573 429 1.074 0.3798 Residuals 189 75497 399 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 ####LSD results > print(LSD.test(mc,'line')) $statistics Mean CV MSerror LSD 48.84375 40.91895 399.4542 19.71249 $parameters Df ntr t.value alpha test name.t 189 28 1.972595 0.05 Fisher-LSD line $means rating std r LCL UCL Min Max PI223 0.375 0.5175492 8 -13.56383 14.31383 0 1 PI39 0.375 0.5175492 8 -13.56383 14.31383 0 1 top_bunch 51.875 21.0335616 8 37.93617 65.81383 10 80 topper 35.625 22.7466638 8 21.68617 49.56383 10 70 v045 38.125 27.8948920 8 24.18617 52.06383 5 80 v046 56.250 19.2260983 8 42.31117 70.18883 20 80 v047 54.375 25.6956778 8 40.43617 68.31383 15 85 v048 57.500 20.0000000 8 43.56117 71.43883 20 80 v049 44.375 23.8203844 8 30.43617 58.31383 15 80 v050 46.875 25.2045206 8 32.93617 60.81383 15 90 v051 56.250 15.2947236 8 42.31117 70.18883 40 80 v052 51.875 18.1142209 8 37.93617 65.81383 20 75 v053 27.500 10.6904497 8 13.56117 41.43883 10 40 v054 48.125 11.9335960 8 34.18617 62.06383 35 65 v055 44.375 19.8993001 8 30.43617 58.31383 15 75 v056 96.250 3.5355339 8 82.31117 110.18883 90 100 v057 24.375 13.7418807 8 10.43617 38.31383 5 45 v058 46.250 24.6040066 8 32.31117 60.18883 15 85 v059 69.375 14.9851117 8 55.43617 83.31383 45 85 v060 49.375 24.2660109 8 35.43617 63.31383 15 85 v061 41.875 25.7650788 8 27.93617 55.81383 20 90 v062 44.375 27.9588345 8 30.43617 58.31383 15 85 v063 86.250 14.3302876 8 72.31117 100.18883 65 100 v064 70.625 13.9993622 8 56.68617 84.56383 50 90 v065 75.000 12.5356634 8 61.06117 88.93883 55 90 v066 48.750 28.3788150 8 34.81117 62.68883 15 95 v067 58.125 32.9434038 8 44.18617 72.06383 10 100 v068 43.125 15.1038075 8 29.18617 57.06383 30 70 $comparison NULL $groups trt means M 1 v056 96.250 a 2 v063 86.250 ab 3 v065 75.000 bc 4 v064 70.625 bcd 5 v059 69.375 bcd 6 v067 58.125 cde 7 v048 57.500 cdef 8 v046 56.250 cdef 9 v051 56.250 cdef 10 v047 54.375 defg 11 top_bunch 51.875 defg 12 v052 51.875 defg 13 v060 49.375 efg 14 v066 48.750 efg 15 v054 48.125 efg 16 v050 46.875 efgh 17 v058 46.250 efgh 18 v049 44.375 efgh 19 v055 44.375 efgh 20 v062 44.375 efgh 21 v068 43.125 efghi 22 v061 41.875 efghi 23 v045 38.125 fghi 24 topper 35.625 ghi 25 v053 27.500 hi 26 v057 24.375 i 27 PI223 0.375 j 28 PI39 0.375 j