#Data
v<-c(1,2,2,5,6,5,2,1)
TR<-c(1,1,1,2,2,2,3,3)
d<-data.frame(v,TR)
d$TR<-as.factor(d$TR)
str(d)
## 'data.frame': 8 obs. of 2 variables:
## $ v : num 1 2 2 5 6 5 2 1
## $ TR: Factor w/ 3 levels "1","2","3": 1 1 1 2 2 2 3 3
#Do the Anova Table
d.fit<-aov(v~TR,data=d)
summary(d.fit)
## Df Sum Sq Mean Sq F value Pr(>F)
## TR 2 26.167 13.083 35.68 0.0011 **
## Residuals 5 1.833 0.367
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
You can also embed plots, for example:
Note that the echo = FALSE parameter was added to the code chunk to prevent printing of the R code that generated the plot.