루프의 prop.test 함수

Question

루프를 사용하여 R에서 2 샘플 비율 테스트를 수행하고 Health Center와 Measure로 나누기를 원합니다. 다음은 내 데이터가 어떻게 설정되었는지 보여주는 예입니다 (웹 사이트에서 내 데이터 세트의 이미지를 업로드 할 수 없음).

기본적으로 Health Center A와 Measure A의 소품 .test 기능을 사용하고 모든 내 건강 센터 (29 개)와 조치 (14 개)에 대해이 작업을 반복하십시오. 이 코드를 반복하는 코드가 무엇인지 확신 할 수 없기 때문에 원하는 모든 비례 테스트를 수행하고 원하는대로 분할 할 수 있습니다.

도움이 될 것입니다.

나는 모든 0을 삭제했습니다. 그러나 코드가 제대로 작동했지만 원하는대로 데이터를 분할하지 않았습니다. 보건 센터와 측정 단위로 나뉘어지기를 바랬지 만 보건 센터 옆에서 나누었 기 때문에 measure.t가있는 보건 센터 A 대신에 propest.test가 분석하는 28 가지 조치가 있습니다. 아래 출력 예제를 참조하십시오.

Answer 1

이

lapply(split(df, df$Health_Center), function(x) prop.test(as.matrix(cbind(x[,3], x[,4]-x[,3])))) 

출력

$A1 

     2-sample test for equality of proportions with continuity correction 

data: as.matrix(cbind(x[, 3], x[, 4] - x[, 3])) 
X-squared = 1.713, df = 1, p-value = 0.1906 
alternative hypothesis: two.sided 
95 percent confidence interval: 
-0.28136975 0.04846377 
sample estimates: 
    prop 1 prop 2 
0.2307692 0.3472222 


$A2 

     2-sample test for equality of proportions with continuity correction 

data: as.matrix(cbind(x[, 3], x[, 4] - x[, 3])) 
X-squared = 0.12192, df = 1, p-value = 0.727 
alternative hypothesis: two.sided 
95 percent confidence interval: 
-0.1789668 0.1114439 
sample estimates: 
    prop 1 prop 2 
0.4800000 0.5137615 

입력 데이터

,691를 사용하여 분할 및 prop.test을 수행 할 수 있습니다

df <- data.frame(Health_Center=c("A1","A2","A1","A2"), 
       Measure=c("A","B","A","B"), 
       Numerator=c(15,48,25,56), 
       Denominator=c(65,100,72,109), stringsAsFactors=F)