兩個比例的多重比較prop.test

Question

我有大量的治療和對照組，我需要提供人口比例的比較。我正在尋找一種方法來循環訪問每個類別提供測試的data.frame。

樣本數據：

test_data <- data.frame(
Category = c("A","A","B","B"), 
Churn = c(56,46,83,58), 
Other = c(180,555,144,86)) 

例如，比較類別A（56/180至五百五十五分之四十六）等等。

我的初始溶液：

by(test_data, test_data$Category, 
function(x) prop.test(test_data$Churn, test_data$Other)) 

的問題：該溶液由輸出類別但提供了一個4樣本測試，而不是一個雙樣品檢驗。我發現很多解決方案可以很好地遍歷行，但不是很多的類。現在輸出列表很好。

真的很感謝這個幫助！

Answer 1

您的by()功能不正確。您沒有使用傳入的x值。通過使用原始變量名稱（test_data），每個by()調用都沒有數據是子集。嘗試

by(test_data, test_data$Category, 
function(x) prop.test(x$Churn, x$Other))