This is an R HTML document. When you click the Knit HTML button a web page will be generated that includes both content as well as the output of any embedded R code chunks within the document. You can embed an R code chunk like this:

way <- read.table("3way.txt",h=T)
a <- with(way,(aov(way$x~way$color*way$length*way$order
                   +Error(way$sub / (way$color))
                   )
               )
          )
a
## 
## Call:
## aov(formula = way$x ~ way$color * way$length * way$order + Error(way$sub/(way$color)))
## 
## Grand Mean: 0.5473647
## 
## Stratum 1: way$sub
## 
## Terms:
##                  way$order
## Sum of Squares  0.01966049
## Deg. of Freedom          1
## 
## 3 out of 4 effects not estimable
## Estimated effects are balanced
## 
## Stratum 2: way$sub:way$color
## 
## Terms:
##                  way$color
## Sum of Squares  0.01374355
## Deg. of Freedom          1
## 
## 3 out of 4 effects not estimable
## Estimated effects are balanced
## 
## Stratum 3: Within
## 
## Terms:
##                 way$color way$length way$order way$color:way$length
## Sum of Squares    0.00073   46.00434   0.12368              0.00675
## Deg. of Freedom         1          1         1                    1
##                 way$color:way$order way$length:way$order
## Sum of Squares              0.08389              0.02541
## Deg. of Freedom                   1                    1
##                 way$color:way$length:way$order Residuals
## Sum of Squares                         0.00172  19.14238
## Deg. of Freedom                              1       458
## 
## Residual standard error: 0.2044397
## Estimated effects may be unbalanced
summary(a)
## 
## Error: way$sub
##           Df  Sum Sq Mean Sq
## way$order  1 0.01966 0.01966
## 
## Error: way$sub:way$color
##           Df  Sum Sq Mean Sq
## way$color  1 0.01374 0.01374
## 
## Error: Within
##                                 Df Sum Sq Mean Sq  F value Pr(>F)    
## way$color                        1   0.00    0.00    0.017 0.8952    
## way$length                       1  46.00   46.00 1100.698 <2e-16 ***
## way$order                        1   0.12    0.12    2.959 0.0861 .  
## way$color:way$length             1   0.01    0.01    0.162 0.6879    
## way$color:way$order              1   0.08    0.08    2.007 0.1572    
## way$length:way$order             1   0.03    0.03    0.608 0.4359    
## way$color:way$length:way$order   1   0.00    0.00    0.041 0.8394    
## Residuals                      458  19.14    0.04                    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

You can also embed plots, for example: