플롯이 너무 많이 점프 할 때마다 방향을 반대로 바꿉니다. 효과적인 기준 중 하나는 모든 구성 요소의 총 점프 량을 계산하는 것입니다. 다음 고유 벡터가 부정되면 총 점프 량을 계산합니다. 후자가 적 으면 다음 고유 벡터를 무효화하십시오.

구현은 다음과 같습니다. (나는 익숙하지 않다zoo 더 우아한 솔루션을 허용 할 수있는 .)

require(zoo)
amend <- function(result) {
  result.m <- as.matrix(result)
  n <- dim(result.m)[1]
  delta <- apply(abs(result.m[-1,] - result.m[-n,]), 1, sum)
  delta.1 <- apply(abs(result.m[-1,] + result.m[-n,]), 1, sum)
  signs <- c(1, cumprod(rep(-1, n-1) ^ (delta.1 <= delta)))
  zoo(result * signs)
}

예를 들어, 직교 그룹에서 임의의 보행을 실행하고 관심을 끌기 위해 약간의 지터를 만들어 봅시다.

random.rotation <- function(eps) {
  theta <- rnorm(3, sd=eps)
  matrix(c(1, theta[1:2], -theta[1], 1, theta[3], -theta[2:3], 1), 3)
}
set.seed(17)
n.times <- 1000
x <- matrix(1., nrow=n.times, ncol=3)
for (i in 2:n.times) {
  x[i,] <- random.rotation(.05) %*% x[i-1,]
}

롤링 PCA는 다음과 같습니다.

window <- 31
data <- zoo(x)
result <- rollapply(data, window,
  function(x) summary(princomp(x))$loadings[, 1], by.column = FALSE, align = "right")
plot(result)

이제 고정 버전 :

plot(amend(result))

답변