如何总结所有G值？

Question

在我的数据子集的所有行上使用G.test后

apply(datamixG +1 , 1, G.test)

我为每一行看到这样的输出

[[1]]

        G-test for given probabilities

data:  [(newX,,i)
G = 3.9624, df = 1, p-value = 0.04653

我有46行。我需要总结df和G值。有没有办法让R报告不同的G值和/或总结所有的G值和df？

Answer 1

我假设您正在使用G.test包中的RVAideMemoire函数：

# Sample data (always a good idea to post!)
dat <- matrix(1:4, nrow=2)
library(RVAideMemoire)
tests <- apply(dat, 1, G.test)

您可以使用unlist和lapply从列表中的每个元素中提取单个值，并返回结果的向量：

dfs <- unlist(lapply(tests, "[[", "parameter"))
dfs
# df df 
#  1  1 
sum(dfs)
# [1] 2

Gs <- unlist(lapply(tests, "[[", "statistic"))
Gs
#         G         G 
# 1.0464963 0.6795961
sum(Gs)
# [1] 1.726092