请求计算R语言数据框各列频率的乘积

咱们可以分步骤来实现这个需求，用R的基础函数就能轻松搞定，具体方法如下：

方法一：分步实现（逻辑更清晰） #

1. 先把每列的频数分布存下来——你已经用apply(df, 2, table)拿到了这个结果，我们直接存起来方便后续处理：

# 存储每列的频数表
freq_tables <- apply(df, 2, table)

2. 计算总样本量（也就是数据框的总行数）：

total_obs <- nrow(df)

3. 遍历每个频数表，先把频数转成频率（每个类别出现次数除以总样本量），再计算该列所有频率的乘积：

# 用lapply得到列表格式的结果
col_freq_products <- lapply(freq_tables, function(counts) {
  frequencies <- counts / total_obs
  prod(frequencies)
})

# 如果想要向量格式的结果，换成sapply就行
col_freq_products <- sapply(freq_tables, function(counts) prod(counts / total_obs))

方法二：一步到位（代码更简洁） #

如果不想存中间的频数表，也可以直接对每列做处理，一行逻辑搞定：

col_freq_products <- sapply(df, function(column) {
  prod(table(column) / nrow(df))
})

结果说明 #

运行完代码后，col_freq_products会是一个带列名的向量，每个元素对应原数据框中一列的频率乘积。比如针对你给出的示例数据，列C的频率乘积就是(62/364)*(65/364)*(65/364)*(56/364)*(60/364)*(56/364)的计算结果。

内容的提问来源于stack exchange，提问作者Omry Atia