-- ====================================================================== 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 ======================================================================
Alexandre, boa tarde!
Não sei se entendi direito, mas fiz uma tentativa. Se a ideia servir dá pra colocar na forma de uma função.
### <code r>treat <- gl(4, 15, labels = paste("t", 1:4, sep="")); treat # Fatores
# Variablesset.sed(124)sp <- cbind(c(rnorm(10, 5, 0.25), rnorm(50, 2.5, 0.25)), rnorm(60, 2.5, 0.25),c(rnorm(10, 12, 0.25), rnorm(50, 2.5, 0.25)), rnorm(60, 2.5, 0.25))colnames(sp) <- c("sp1", "sp2", "sp3", "sp4")head(sp)
#Criação da matriz de contrastesTI <- model.matrix(~ treat-1)head(TI)
f <- length(levels(treat))comb <- t(combn(1:f, 2))n <- nrow(comb)
contr2 <- NULLfor (x in 1:n) {i <- comb[x, 1]j <- comb[x, 2]tmp <- list(TI[,i] - TI[,j]); names(tmp) <- paste0("TI",i, "_", j)contr2 <- c(contr2, tmp)}contr2
head(do.call(cbind, contr2))# TI1_2 TI1_3 TI1_4 TI2_3 TI2_4 TI3_4# 1 1 1 1 0 0 0# 2 1 1 1 0 0 0# 3 1 1 1 0 0 0# 4 1 1 1 0 0 0# 5 1 1 1 0 0 0# 6 1 1 1 0 0 0### </code>
================================================
Éder ComunelloAgronomist (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.