R/fullKF.R

"fullKF" <-
  function(theta, kinscal, kmat, jvec,fixedkmat=FALSE,
           kinscalspecial = list(), kinscalspecialspec = list(),
           nocolsums = FALSE)
  {
    #fullKF
    
    dimk <- nrow(kmat)
    A <- matrix(nrow = dimk, ncol = dimk)
    K <- fillK(theta, kinscal, kmat, fixedkmat, kinscalspecial,
               kinscalspecialspec, nocolsums)
    eigenlijk <- eigen(K, only.values = F)
    V <- eigenlijk$vectors
    gamma <- solve(V) %*% jvec
    for(j in 1:dimk) 
      A[j, ] <- V[ ,j] * gamma[j]
    list(A=A, values = eigenlijk$values)
  }

Try the TIMP package in your browser

Any scripts or data that you put into this service are public.

TIMP documentation built on May 2, 2019, 5:55 p.m.