<div dir="ltr">Alexandre,<div><br></div><div>Salvo se eu me perdi na leitura do seu texto e na sua função, basta você decidir o quê você vai colocar nas células da matriz quando os valores não estiverem como você deseja e usar um teste if com cláusula else nas atribuições, ou uma construção similar à que você já faz para determinar se vai escrever na diagonal superior ou inferior com seus vetores above1 e above2.</div><div><br></div><div>HTH </div></div><div class="gmail_extra"><br><div class="gmail_quote">2015-08-16 23:50 GMT-03:00 ASANTOS <span dir="ltr"><<a href="mailto:alexandresantosbr@yahoo.com.br" target="_blank">alexandresantosbr@yahoo.com.br</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Prezados listeiros,<br>
<br>
No meu CRM abaixo tenho uma matriz de correlação com o valor de p na diagonal superior o e valor r de Pearson da diagonal inferior, no entanto eu gostaria de exibir uma matriz de correlação apenas com os p valores maiores ou igual 0.05 e r maior ou igual a 0.08 isso é possível?<br>
<br>
Segue CRM:<br>
<br>
require(Hmisc)## Pacote<br>
<br>
#Variáveis artificiais<br>
var1<-sort(rnorm(500,1), decreasing = TRUE)<br>
var2<-rnorm(500,10)<br>
var3<-sort(rnorm(500,22))<br>
var4<-sort(rnorm(500,10))<br>
var5<-sort(rnorm(500,0.5), decreasing = TRUE)<br>
var6<-sort(rnorm(500,25), decreasing = TRUE)<br>
var7<-rnorm(500,10)<br>
dados<-cbind(var1,var2,var3,var4,var5,var6,var7)<br>
<br>
#Função que dispõem p valor na diagonal superior e r valor na inferior<br>
<br>
pn <- function(X){crossprod(!<a href="http://is.na" rel="noreferrer" target="_blank">is.na</a>(X))}<br>
<br>
cor.prob <- function(X){<br>
pair.SampSize <- pn(X)<br>
above1 <- row(pair.SampSize) < col(pair.SampSize)<br>
pair.df <- pair.SampSize[above1] - 2<br>
R <- cor(X, use="pair")<br>
above2 <- row(R) < col(R)<br>
r2 <- R[above2]^2<br>
Fstat <- (r2 * pair.df)/(1 - r2)<br>
R[above2] <- 1 - pf(Fstat, 1, pair.df)<br>
R<br>
}<br>
<br>
#Resultado da matriz<br>
cor.prob(dados)<br>
#<br>
<br>
-- <br>
======================================================================<br>
Alexandre dos Santos<br>
Proteção Florestal<br>
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso<br>
Campus Cáceres<br>
Caixa Postal 244<br>
Avenida dos Ramires, s/n<br>
Bairro: Distrito Industrial<br>
Cáceres - MT CEP: 78.200-000<br>
Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)<br>
<a href="mailto:e-mails%3Aalexandresantosbr@yahoo.com.br" target="_blank">e-mails:alexandresantosbr@yahoo.com.br</a><br>
<a href="mailto:alexandre.santos@cas.ifmt.edu.br" target="_blank">alexandre.santos@cas.ifmt.edu.br</a><br>
Lattes: <a href="http://lattes.cnpq.br/1360403201088680" rel="noreferrer" target="_blank">http://lattes.cnpq.br/1360403201088680</a><br>
======================================================================<br>
<br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br>
</blockquote></div><br></div>