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]]
