[R-br] Dúvida com quando o array tem NA

[Victor Eduardo]

Consegui resolver o problema, pessoal.


:)

Em 14 de junho de 2012 10:38, Victor Eduardo <victorduca08 em gmail.com>escreveu:

> Galera, estava com um array de tamanho 150 onde  tinham matrizes 1x1.
> Sendo que em alguns tinham NA. Dentro de minha função criei uma variável
> auxliar para identificar esses NA. Fiz do seguinte modo:
>
> aux<-NULL
> for(t in 1:n){
> if(is.na(Y[,,t])==TRUE){
>  aux[t]<-0
> }else{
> aux[t]<-1
>  }
> }
>
>
> Ao identificar, o aux terá tamanho de 150, de 0 e 1.
>
> Feito essa verificação, estou pulando para a parte debaixo que é assim:
>
> 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'
>
>
>
>
> Abraços!
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120614/72c148c2/attachment.html>