Hola a todos,
Pego un trozo de una práctica de R sobre el tema (como veis está en
formato spin: https://yihui.name/knitr/demo/stitch) por si resulta de
interés:
#' Consideramos datos recogidos en un estudio de mejora de calidad en
una fábrica de semiconductores.
#' Se obtuvo una muestra de obleas que se clasificaron dependiendo de si
se encontraron partículas
#' en la matriz que producía la oblea y de si la calidad de oblea era buena
#' (Para más detalles Hall, 1994. Analysis of defectivity of
semiconductor wafers by contigency table.
#' Proceedings of the Institute of Environmental Sciences 1, 177-183).
n <- c(320, 14, 80, 36)
particulas <- gl(2, 1, 4, labels = c("no","si"))
calidad <- gl(2, 2, labels = c("buena", "mala"))
df <- data.frame(n, particulas, calidad)
df
#' En lugar de estar en el formato de un conjunto de datos (`data.frame`)
#' puede que los datos estén en formato de tabla (`table`, `matrix`):
tabla <- xtabs(n ~ calidad + particulas)
tabla
#' Por ejemplo, si se quiere simular bajo independencia,
#' estimando las probabilidades a partir de la tabla,
#' consideraríamos:
pind <- (rowSums(tabla) %o% colSums(tabla))/(sum(tabla)^2)
matrix(pind, nrow = nrow(tabla))
rtablas <- rmultinom(ntsim, sum(n), pind) # Se fija el total de la tabla...
rtablas[ , 1:5] # Las cinco primeras simulaciones
#' Para realizar el contraste de independencia:
res <- chisq.test(tabla)
res
El 17/04/2017 a las 17:40, Carlos J. Gil Bellosta escribió:
a <- 1:4
b <- 3:6
a %*% t(b)
El 17 de abril de 2017, 17:26, Horacio <horacio9...@gmail.com> escribió:
Buenas tengo una matriz de contingencia de 5x4 donde en la última fila
y última columna tengo las frecuencias marginales, en función de estas
quiero sacar las esperadas, pero cuando hago por ejemplo el producto,,
Esperadas<-Estrategia["Suma",]*Estrategia[,"Suma_e"]
E1 E2 E3 E4 Suma
5550 14160 52650 89600 20720
Warning message:
In Estrategia["Suma", ] * Estrategia[, "Suma_e"] :
longitud de objeto mayor no es múltiplo de la longitud de uno menor
me da un vector y no una matriz,,, como puedo hacer esto? no sé si se
entiende lo que busco
Gracias,,,
_______________________________________________
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es
[[alternative HTML version deleted]]
_______________________________________________
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es
--
Ruben Fernandez Casal
Department of Mathematics
Faculty of Computer Science
Universidade da Coruña
Corporate email: ruben.fcasal <at> udc <dot> es
--
_______________________________________________
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es
