Ao identificar, o aux terá tamanho de 150, de 0 e 1.
for(t in 1:n){
if(aux[t]==1){
B<-matrix(Z[,,t],ncol=m)
U[,,t]<-Y[,,t]-B%*%At[,,t]-D[,,t]
F[,,t]<-B%*%Pt[,,t]%*%t(B)+H[,,t]
K<-T[,,t]%*%Pt[,,t]%*%t(B)%*%ginv(F[,,t])
L[,,t]<-T[,,t]-K%*%B
At[,,t+1]<-T[,,t]%*%At[,,t]+K%*%U[,,t]+C[,,t]
Pt[,,t+1]<-T[,,t]%*%Pt[,,t]%*%t(L[,,t])+R[,,t]%*%Q[,,t]%*%t(R[,,t])
}else{
At[,,t+1]<-T[,,t]%*%At[,,t]+C[,,t]
Pt[,,t+1]<-T[,,t]%*%Pt[,,t]%*%t(L[,,t])+R[,,t]%*%Q[,,t]%*%t(R[,,t])
}
}
return(list(At,Pt))
Só que o R está dando um erro que não estou conseguindo resolver, o erro é esse: Erro em svd(X) : infinite or missing values in 'x'