cor.prob <- function(X, dfr=nrow(X)-2){
R <- cor(X,use='na.or.complete')
above <- row(R) < col(R)
r2 <- R[above]^2
Fstat <- r2*dfr/(1-r2)
R[above] <- 1-pf(Fstat, 1, dfr)
R
}
\begin{signature}
<<>>=
Prof. Dr. Ivan Bezerra Allaman
Universidade Estadual de Santa Cruz
Departamento de Ciências Exatas e Tecnológicas
Ilhéus/BA - Brasil
Fone: +55 73 3680-5596
E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com
@
\end{signature}