From ac0d1833af3029961a23707fd647562af22e52d4 Mon Sep 17 00:00:00 2001 From: npai9 Date: Sun, 5 Jan 2025 10:29:54 -0500 Subject: [PATCH] fix face bug with missing data and 1 PC --- R/fpca.face.R | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/R/fpca.face.R b/R/fpca.face.R index 70e94f31..770dd829 100644 --- a/R/fpca.face.R +++ b/R/fpca.face.R @@ -318,11 +318,11 @@ function(Y=NULL,ydata=NULL,Y.pred = NULL,argvals=NULL,pve = 0.99, npc = NULL, ######################################### if(imputation) { - A.N <- A[,1:N] + A.N <- A[,1:N, drop = F] d <- Sigma[1:N] sigmahat2 <- max(mean(Y[!Index.miss]^2) -sum(Sigma),0) Xi <- t(A.N)%*%Ytilde - Xi <- t(as.matrix(B%*%(A0%*%((A.N%*%diag(d/(d+sigmahat2/J)))%*%Xi)))) + Xi <- t(as.matrix(B%*%(A0%*%((A.N%*%diag(d/(d+sigmahat2/J), nrow = N))%*%Xi)))) Y <- Y*(1-Index.miss) + Xi*Index.miss if(sum(is.na(Y))>0) print("error")