File: contr.fda.R

package info (click to toggle)
r-cran-mda 0.4-10-2
  • links: PTS, VCS
  • area: main
  • in suites: buster
  • size: 780 kB
  • sloc: fortran: 2,597; f90: 15; makefile: 2
file content (13 lines) | stat: -rw-r--r-- 332 bytes parent folder | download | duplicates (4) 
1
2
3
4
5
6
7
8
9
10
11
12
13
 
contr.fda <-
function (p = rep(1, d[1]), contrast.default = contr.helmert(length(p))) 
{
    d <- dim(contrast.default)
    sqp <- sqrt(p/sum(p))
    x <- cbind(1, contrast.default) * outer(sqp, rep(1, d[2] + 
        1))
    qx <- qr(x)
    J <- qx$rank
    qr.qy(qx, diag(d[1])[, seq(2, J)])/outer(sqp, rep(1, J - 
        1))
}