One-way random effects analysis of variance vetyou=read.table('vetyou.txt') #Now fit the model you=lm(log(cnt+1)~reg,vetyou) #Obtain the ANOVA table anova.lm(you)

Gives something like this:

Analysis of Variance Table Response: log(cnt + 1) Df Sum Sq Mean Sq F value Pr(>F) reg 4 80.649 20.162 10.353 4.961e-06 *** Residuals 45 87.637 1.947 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1