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

