[R-br] plotar valores específicos com função image( )

Paulo Justiniano paulojus em leg.ufpr.br
Quarta Junho 6 18:06:44 BRT 2012


Coloque os valores que nao quer plotar como NA e mude a cor de fundo apra 
cor desejada

valores NA sao plotados com a for de fundo (background)


On Wed, 6 Jun 2012, Vitor Aguiar wrote:

> Olá,
> 
> desculpe por estar postando a pergunta novamente, mas não tive resposta da primeira vez. Se alguém puder ajudar, agradeço muito.
> 
> Tenho uma matrix com que estou fazendo um heatmap usando a função image( ). Gostaria de plotar não apenas os valores observados na
> matrix, mas também valores que seriam possíveis mas não são observados. 
> 
> Isso significa o seguinte: 
> Tenho uma matrix 16x16. Valores possíveis, mas não observados, correspondem a "0" (zeros) onde a soma (índice da linha + índice da
> coluna) são iguais ou menores que 15.
> 
> Exemplo: Se eu tenho um valor "0" na posição [14, 3], eu NÃO quero plotar, pois 14 + 3 = 17.
> Se eu tenho um valor "0" na posição [11, 3], eu quero plotar, pois 11 + 3 = 14.
> 
> Eu queria que esses valores fossem plotados em cor diferente do meu heatmap, por exemplo, em azul.
> 
> É possível fazer isso?
> 
> Matrix:
> x = c(3045, 893, 692, 830, 617, 155, 246, 657, 105, 60, 18, 7, 7, 4, 2, 11234, 2985, 2242, 2471, 1575, 366, 503, 1283, 170, 79, 32, 6,
> 4, 1, 3, 19475, 4756, 3233, 3251, 1810, 409, 575, 1210, 139, 41, 11, 4, 2, 0, 0, 20830, 4739, 2990, 2531, 1346, 298, 325, 612, 60, 17,
> 1, 0, 1, 0, 0, 15304, 3196, 1885, 1440, 610, 117, 115, 185, 14, 2, 0, 0, 0, 0, 0, 8026, 1535, 806, 539, 223, 33, 37, 39, 0, 0, 0, 0, 0,
> 0, 0, 3300, 562, 286, 141, 45, 14, 5, 12, 0, 0, 0, 0, 0, 0, 0, 1067, 160, 65, 40, 14, 1, 2, 0, 0, 0, 0, 0, 0, 0, 0, 277, 47, 6, 2,
> 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 72, 6, 1, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 11, 5, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
> 
> xmat = matrix(x, ncol = 12)
> 
> xmat = cbind(xmat, rep(0, 15), rep(0, 15), rep(0, 15), rep(0, 15))
> xmat = rbind(xmat, rep(0, 16))
> dimnames(xmat) = list(0:15, 0:15)
> image(log(xmat))
> 
> 
> Muito Obrigado,
> Vitor
> 
> 
>


Mais detalhes sobre a lista de discussão R-br