File: fix.theta.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 (8 lines) | stat: -rw-r--r-- 166 bytes parent folder | download | duplicates (4)
1
2
3
4
5
6
7
8
fix.theta <-
function (theta, Q) 
{
    M <- t(theta) %*% Q %*% theta
    eM <- eigen(M, symmetric = TRUE)
    scale(theta %*% eM$vectors, FALSE, sqrt(eM$values))
}