Pessoal, Desculpem, a anterior fornece uma matriz nula da dimensão de 'm', acabei de observar. A correta é essa: # R help m <- matrix(runif(24), ncol=3) res <- matrix(0, nrow(m), ncol(m)) mg <- mean(m) # média geral mc <- colMeans(m) # média das colunas mr <- rowMeans (m) # médias das linhas for (i in 1:nrow(m)) for (j in 1:ncol(m)) res[i,j] <- m[i,j] + mg - mc[j] - mr[i] res -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Em 14 de abril de 2012 08:11, Jose Claudio Faria <joseclaudio.faria@gmail.com> escreveu:
Pessoal,
Uma pequena correção (apenas no comentário) para dar mais clareza
m <- matrix(1:24, ncol=3) res <- matrix(0, nrow(m), ncol(m)) mg <- mean(m) # média geral mc <- colMeans(m) # média das colunas mr <- rowMeans (m) # médias das linhas for (i in 1:nrow(m)) for (j in 1:ncol(m)) res[i,j] <- m[i,j] + mg - mc[j] - mr[i] res -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
Em 14 de abril de 2012 08:06, Jose Claudio Faria <joseclaudio.faria@gmail.com> escreveu:
Pessoal,
Gostaria de evitar um loop no procedimento abaixo da forma mais elegante possível:
# R help m <- matrix(runif(24), ncol=3) res <- matrix(0, nrow(m), ncol(m)) mg <- mean(m) # média geral mc <- colMeans(m) # média das colunas mr <- rowMeans (m) # médias das variáveis for (i in 1:nrow(m)) for (j in 1:ncol(m)) res[i,j] <- m[i,j] + mg - mc[j] - mr[i] res
Toda ajuda é bem vinda... -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\