[R-br] Ordenamento de matriz by colnames e by rownames
Luiz Roberto Martins Pinto
luizroberto.uesc em gmail.com
Segunda Agosto 17 17:20:59 BRT 2015
Walmes,
Agradeço as suas considerações.
Depois que enviei a mensagem, consegui resolver o problema da seguinte
forma:
n=14
pv<-matrix(nrow=(n),ncol=(n))
# estimo valores para pv[i,j], que, obrigatoriamente, seguem a seguinte
ordem de entrada na matriz pv
# esta matriz de p-valor é simétrica e tem 1.00 na diagonal principal
rn=c("SSRY_027","SSRY_021","SSRY_141","SSRY_008","SSRY_185","SSRY_324","SSRY_040","SSRY_035","SSRY_235","SSRY_183","SSRY_043","GA_005","GA_136","GA_012")
colnames(pv)=rn
rownames(pv)=rn
### passo-a-passo
pv0=(pv[,order(rownames(pv))])
tpv0=t(pv0)
pv2=tpv0[,order(colnames(tpv0))]
pv2 tem o formato que desejo.
Abraços,
Luiz Roberto.
Luiz Roberto Martins Pinto
Prof. Pleno/DCET/UESC
Laboratório de Estatística Computacional
Universidade Estadual de Santa Cruz
Ilhéus-Bahia-Brasil
luizroberto.uesc em gmail.com
skype: lrmpinto
http://lattes.cnpq.br/2732314327604831
Em 17 de agosto de 2015 16:40, Walmes Zeviani <walmeszeviani em gmail.com>
escreveu:
> O problema é sort() e ordem() ordenam um vetor, olham para uma dimensão.
> Você não pode sisplesmente usar order() nas linhas e depois nas colunas
> porque se fizer isso vai bagungar com as ligações fora da diagonal. Então a
> tarefa é mais complicada porque você tem que ordenar as linhas e colunas
> sem estragar a estrutura de p-valores. Uma vez eu, uns 5 anos atrás, eu fiz
> uma função para ordenar uma matriz de contrastes. Tinha na diagonal a média
> de cada tratamento, na diagonal de cima o contraste entre duas médias e na
> debaixo o p-valor do teste de hipótese sobre o contraste. Levei dias para
> construir e não consegui. Meu código se perdeu. Lembro que eu resolvi
> mudando a ordem dos níveis do fator antes de começar a análise, ao invés de
> ser alfabética, eu ordenei pela média amostral (ou média ajustada de uma
> análise).
>
> Walmes.
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150817/fa870664/attachment.html>
Mais detalhes sobre a lista de discussão R-br