0

我有一个 n>1000 的数据框，其中每一行都包含 Year（数字年份）和 Gender（“男性”或“女性”）列的数据。我想计算性别 == “男性”在年份之间成对比例的 t 检验。我已经成功地创建了这个比例的情节，为此我附上了代码和情节。我无法将此扩展到 prop_test() 函数。我无法附加我的数据，但包含示例数据集的代码。

sample_data <- as_tibble(data.frame(Gender = sample(c("Male", "Female"), 1000, replace = TRUE), 
                   Year = sample(c(2016, 2017, 2018, 2019), 1000, replace = TRUE)))

sample_data %>% 
  group_by(Year) %>% 
  summarise(prop = sum(Gender == "Male", na.rm = TRUE) / n()) %>% 
  ggplot(mapping = aes(x = Year, y = prop)) +
  geom_bar(stat = "identity") +
  geom_text(aes(label = round(prop, 2), vjust = -0.25)) +
  labs(y = "prop Male")

由此产生的男性比例图（按年份）

图像

请告知我如何修改用于生成绘图的代码以计算比例的成对 t 检验。我尝试过类似的方法

sample_data %>%    prop_test(Gender ~ Year)

但这会产生错误“rowSums(x) 中的错误：'x' 必须是数字”

4

0 回答 0