Description Usage Arguments Details Value Author(s) References See Also Examples
The function performs analysis of polynomial regression in simple designs with quantitative treatments.
1 |
data |
data is a data.frame data frame with two columns, treatments and response (completely randomized design) data frame with three columns, treatments, blocks and response (randomized block design) data frame with four columns, treatments, rows, cols and response (latin square design) data frame with five columns, treatments, square, rows, cols and response (several latin squares) |
design |
1 = completely randomized design 2 = randomized block design 3 = latin square design 4 = several latin squares |
list |
FALSE = a single response variable TRUE = multivariable response |
type |
type is form of obtain sum of squares 1 = a sequential sum of squares 2 = a partial sum of squares |
The response and the treatments must be numeric. Other variables can be numeric or factors.
Returns analysis of variance, models, t test for coefficients and R squared and adjusted R squared.
Emmanuel Arnhold <emmanuelarnhold@yahoo.com.br>
KAPS, M. and LAMBERSON, W. R. Biostatistics for Animal Science: an introductory text. 2nd Edition. CABI Publishing, Wallingford, Oxfordshire, UK, 2009. 504p.
SAMPAIO, I. B. M. Estatistica aplicada a experimentacao animal. 3nd Edition. Belo Horizonte: Editora FEPMVZ, Fundacao de Ensino e Pesquisa em Medicina Veterinaria e Zootecnia, 2010. 264p.
lm, lme(package nlme), ea1(package easyanova), er1
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 | # analysis in completely randomized design
data(data1)
r1=er2(data1)
names(r1)
r1
r1[1]
# analysis in randomized block design
data(data2)
r2=er2(data2, design=2)
r2
# analysis in latin square design
data(data3)
r3=er2(data3, design=3)
r3
# analysis in several latin squares
data(data4)
r4=er2(data4, design=4)
r4
# data
treatments=rep(c(0.5,1,1.5,2,2.5,3), c(3,3,3,3,3,3))
r1=rnorm(18,60,3)
r2=r1*1:18
r3=r1*18:1
r4=r1*c(c(1:10),10,10,10,10,10,10,10,10)
data6=data.frame(treatments,r1,r2,r3, r4)
# use the argument list = TRUE
er2(data6, design=1, list=TRUE)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.