[R-br] Matriz de correlação entre vetores de diferentes tamanhos [RESOLVIDO]

[Éder Comunello]

Alexandre, boa tarde!

Só uma curiosidade... O resultado final que você deseja é realmente uma
matriz mista (correlações e probalidades)?

​
================================================
Éder Comunello
Agronomist (UEM), MSc in Environ. Sciences (UEM)
DSc in Agricultural Systems Engineering (USP/Esalq)
Brazilian Agricultural Research Corporation (Embrapa)
Dourados, MS, Brazil |<O>|
================================================
GEO, -22.2752, -54.8182, 408m
UTC-04:00 / DST: UTC-03:00




Em 21 de março de 2016 11:02, ASANTOS <alexandresantosbr em yahoo.com.br>
escreveu:

> Muito obrigado Éder,
>
>       Foi muito útil sim, já estava largando mão,
>
> Redobrados agradecimentos,
>
> --
> ======================================================================
> Alexandre dos Santos
> Proteção Florestal
> IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
> Campus Cáceres
> Caixa Postal 244
> Avenida dos Ramires, s/n
> Bairro: Distrito Industrial
> Cáceres - MT                      CEP: 78.200-000
> Fone: (+55) 65 8132-8112 (TIM)   (+55) 65 9686-6970 (VIVO)e-mails:alexandresantosbr em yahoo.com.br
>         alexandre.santos em cas.ifmt.edu.br
> Lattes: http://lattes.cnpq.br/1360403201088680
> OrcID: orcid.org/0000-0001-8232-6722
> Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10
> LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635
> ======================================================================
>
> Em 21/03/2016 09:21, Éder Comunello escreveu:
>
> Alexandre, bom dia!
>
> Se ainda lhe for útil, segue uma sugestão...
>
> ### <code r>
> var1 <- rnorm(200,5,0.25)
> var2 <- 1:500
> var3 <- rnorm(100,5,0.25)
> var4 <- 500:1
>
> cor.prob <- function(X){
>       pair.SampSize <- crossprod(!is.na(X))
>       above1 <- row(pair.SampSize) < col(pair.SampSize)
>       pair.df <- pair.SampSize[above1] - 2
>
>       R <- cor(X, use="pair")
>       above2 <- row(R) < col(R)
>       r2 <- R[above2]^2
>       Fstat <- (r2 * pair.df)/(1 - r2)
>       R[above2] <- 1 - pf(Fstat, 1, pair.df)
>       R
>       }
>
> tam   <- 100 # amostra
> n_sim <-  99 # simulações
>
> RES <- NULL
> for(i in 1:n_sim){
>      dados   <- sapply(c("var1", "var2", "var3", "var4"), function(x)
> sample(get(x), 100))
>      correla <- round(cor.prob(dados),4)
>      res     <- cbind(i, correla)
>      RES     <- rbind(RES, res)
> }
> #
>
> RES[c(1:5,392:396),]
> #       i    var1    var2    var3   var4
> # var1  1  1.0000  0.8025  0.8451 0.9178
> # var2  1 -0.0253  1.0000  0.6609 0.0106
> # var3  1 -0.0198  0.0444  1.0000 0.4850
> # var4  1 -0.0104 -0.2545 -0.0706 1.0000
> # var1  2  1.0000  0.5200  0.1704 0.8461
> # var4 98  0.1309  0.0332 -0.1361 1.0000
> # var1 99  1.0000  0.1433  0.2047 0.3260
> # var2 99 -0.1474  1.0000  0.5391 0.9524
> # var3 99 -0.1279 -0.0621  1.0000 0.2885
> # var4 99 -0.0992 -0.0060  0.1072 1.0000
> ### </code>
>
>
>
> ​
> ================================================
> Éder Comunello
> Agronomist (UEM), MSc in Environ. Sciences (UEM)
> DSc in Agricultural Systems Engineering (USP/Esalq)
> Brazilian Agricultural Research Corporation (Embrapa)
> Dourados, MS, Brazil |<O>|
> ================================================
> GEO, -22.2752, -54.8182, 408m
> UTC-04:00 / DST: UTC-03:00
>
>
>
>
>
> _______________________________________________
> R-br mailing listR-br em listas.c3sl.ufpr.brhttps://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/20160321/fc683e36/attachment.html>