The R Project SVN R

Rev

Rev 59039 | Rev 66444 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 59039 Rev 61160
Line 84... Line 84...
84 
  \code{\link{glm}},
 84 
  \code{\link{glm}},
85 
  \code{\link{lm}}.
 85 
  \code{\link{lm}}.
86 
}
 86 
}
87 
\examples{
 87 
\examples{
88 
(cH <- contr.helmert(4))
 88 
(cH <- contr.helmert(4))
89 
apply(cH, 2,sum) # column sums are 0
 89 
apply(cH, 2, sum) # column sums are 0
90 
crossprod(cH) # diagonal -- columns are orthogonal
 90 
crossprod(cH) # diagonal -- columns are orthogonal
91 
contr.helmert(4, contrasts = FALSE) # just the 4 x 4 identity matrix
 91 
contr.helmert(4, contrasts = FALSE) # just the 4 x 4 identity matrix
92 
 
 92 
 
93 
(cT <- contr.treatment(5))
 93 
(cT <- contr.treatment(5))
94 
all(crossprod(cT) == diag(4)) # TRUE: even orthonormal
 94 
all(crossprod(cT) == diag(4)) # TRUE: even orthonormal
95 
 
 95 
 
96 
(cT. <- contr.SAS(5))
 96 
(cT. <- contr.SAS(5))
97 
all(crossprod(cT.) == diag(4)) # TRUE
 97 
all(crossprod(cT.) == diag(4)) # TRUE
98 
 
 98 
 
99 
zapsmall(cP <- contr.poly(3)) # Linear and Quadratic
 99 
zapsmall(cP <- contr.poly(3)) # Linear and Quadratic
100 
zapsmall(crossprod(cP), digits=15) # orthonormal up to fuzz
 100 
zapsmall(crossprod(cP), digits = 15) # orthonormal up to fuzz
101 
}
 101 
}
102 
\keyword{design}
 102 
\keyword{design}
103 
\keyword{regression}
 103 
\keyword{regression}
104 
\keyword{array}
 104 
\keyword{array}