Boa noite,

estou com um código que necessita o uso de um for dentro de um while. Entretanto, o uso do for aumentou o tempo de execução.

Gostaria de saber se alguém tem alguma sugestão para contornar este problema???? Acredito que a questão seja organizacional, portanto, não estou enviando um CMR mas caso seja necessário será encaminhado!!!!!!!  

   

                          dif <- 1

                          tol <- 1e-05                                                                  

                          i <- 0

                          neff <- c(2)

                      while(dif > tol){

                                    i <- i + 1

                                    LHS <- rBind( cBind(XpX, XpZ), cBind(ZpX, (ZpZ + Ginv)) )

                                    invLHS <- solve(LHS)

                                    est <- solve(LHS,RHS)

                                    b <- est[1:ncol(X.f),]

                                    u <- est[(ncol(X.f) + 1):nrow(est),]

                                    sig2E <- (ypy - t(b)%*%Xpy - t(u)%*%Zpy)/(N - rankX)

                                    sig2U <- numeric(0)                                  

                                            for(j in 1:length(neff)){

                                                 sig2U <- c(sig2U, ( t(est[z[j]:lz[j],])%*%est[z[j]:lz[j],] + sig2E*sum(diag( invLHS[z[j]:lz[j], z[j]:lz[j]] )) )/neff[j])                                                                                                

                                            }                                    

                                    sig2U <- do.call(rBind,lapply(sig2U,as.matrix))                                     

                                    rel <- as.numeric(sig2E)/sig2U

                                         if(!is.null(var)){

                                                lamb <- rel

                                                ratio <- rep(lamb,neff)

                                                newvar <- rBind(sig2U,sig2E)

                                                cds <- sig2U/sum(newvar) 

                                                dif <- as.numeric(max(abs(newvar - ovar)))                                           

                                         }

                                         if(is.null(var)){

                                                lamb <- rel

                                                ratio <- rep(lamb,neff)

                                                newvar <- rBind(sig2U,sig2E)

                                                cds <- sig2U/sum(newvar) 

                                                dif <- as.numeric(max(abs(cds - ocds)))                                                                                           

                                         }                                     

                                     ovar <- newvar

                                     ocds <- cds

                                     oratio <- ratio

                                     olamb <- lamb 

                                     Ginv <- Diagonal(sum(neff),oratio)                                                                                      

                              }

Desde já agradeço por qualquer colaboração.

Atenciosamente,

Tiago.

#################################################################

 

Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES

 

Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas

 

#################################################################