写了以下一段程序:
- a <- c(77.7,43.2,42.0,79.1,52.1,43.8,88.2,85.4)
- b <- c(81.9,54.5,70.2,83.9,83.9,57.3,92.6,91.9)
- na <- 1752
- nb <- 1720
- a1 <- round(a*na/100)
- b1 <- round(b*nb/100)
- a0 <- na-a1
- b0 <- nb-b1
- pv <- function(df) {
- pv <- array()
- for (i in 1:length(a)) {
- pv[i] <- chisq.test(matrix(as.vector(df[i,]),2,2))$p.value
- return(pv)
- }
- }
试着改了一点,还是不行:
- > pv <- function(df) {
- + pv <- array(rep(0,length(a)),c(1,length(a)))
- + for (i in 1:length(a)) {
- + pv[i] <- chisq.test(matrix(as.vector(df[i,]),2,2))$p.value
- + return(pv)
- + }
- + }
- > pv(df)
- [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8]
- [1,] 0.002177379 0 0 0 0 0 0 0


雷达卡





京公网安备 11010802022788号







