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: