Re: [R-br] Matriz de correlação entre vetores de diferentes tamanhos
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
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@yahoo.com.br alexandre.santos@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 <http://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 list R-br@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.
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@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@yahoo.com.br alexandre.santos@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@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.
Oi Éder, É isso mesmo, estou calculando uma matriz mista média com correlação e probabilidade, agora estou tentando criar um matriz média: #Matriz média através na soma de todas as matrizes e divisão pelo número delas RESmat<-NULL for(i in 1:n_sim){ (RES[RES[,2]==i, 3:6]+RES[RES[,2]==i+1, 3:6])/length(n_sim) } E depois vou comparar essa matriz média com a informação completa para a matriz de covariância por ML, Abraços, -- ====================================================================== 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@yahoo.com.br alexandre.santos@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 15:52, Éder Comunello escreveu:
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@yahoo.com.br <mailto:alexandresantosbr@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 <tel:%28%2B55%29%2065%208132-8112> (TIM)(+55) 65 9686-6970 <tel:%28%2B55%29%2065%209686-6970> (VIVO) e-mails:alexandresantosbr@yahoo.com.br <mailto:e-mails:alexandresantosbr@yahoo.com.br> alexandre.santos@cas.ifmt.edu.br <mailto:alexandre.santos@cas.ifmt.edu.br> Lattes:http://lattes.cnpq.br/1360403201088680 OrcID:orcid.org/0000-0001-8232-6722 <http://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 <http://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 list R-br@listas.c3sl.ufpr.br <mailto:R-br@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.
participantes (2)
-
ASANTOS -
Éder Comunello