marta Garcia
2013-Feb-06 22:41 UTC
[R-es] Construcción de data.frame para 3 y mas factores
Hola estimada comunidad, con el código de abajo hago la data.frame A1 A2 y Ny n 1 1 1 55 105 160 2 1 2 56 104 160 3 1 3 89 71 160 4 2 1 78 82 160 5 2 2 96 64 160 6 2 3 65 95 160 A1 <- factor(c(1,2)) # factor 1, con 2 niveles contrasts(A1) <- contr.treatment(2,1) # fijando el primer nivel 1 de referencia A2 <- factor(c(1,2,3)) # factor 2 con 3 niveles contrasts(A2) <- contr.treatment(3,1)# fijando el nivel 1 de referencia F1 <- merge(A2,A1) names(F1) <- c("A2", "A1") mFor <- formula(cbind(y, Ny) ~ A1 + A2) y <- as.vector(cbind(55,56,89,78,96,65)) # exitos y [1] 55 56 89 78 96 65 n <- as.vector(cbind(160,160,160,160,160,160)) # muestra para la distintas combinaciones de niveles Ny <- n-y # fracasos Dat <- data.frame(A1=F1$A1, A2=F1$A2, y, Ny, n) Dat A1 A2 y Ny n 1 1 1 55 105 160 2 1 2 56 104 160 3 1 3 89 71 160 4 2 1 78 82 160 5 2 2 96 64 160 6 2 3 65 95 160 Lo que quiero es hacer una data.frame para 3 factores y en general para cualquier número pretendo hacer un logit sin usar glm Gracias de antemano [[alternative HTML version deleted]]