所有可能的结果

Question

我的数据框是这样的

gen<-c("A","B","C")
prob<-c("0.95","0.82","0.78")
mw<-c("10","20","50")
df<-data.frame(gen,prob,mw)

 gen prob mw
1   A 0.95 10
2   B 0.82 20
3   C 0.78 50

现在我想要 (A,B,C), (A,B),(A,C),(B,C),(A),(B),(C) 的所有可能结果,(NONE)，例如 (A,B,C)=0.95*0.82*0.78= 0.60762 的概率。

trials <- data.frame(matrix(nrow=0,ncol=length(gen)))
for(i in 1:length(gen)){
  trials.tmp <- t(combn(gen,i))
  trials <- rbind(trials,cbind(trials.tmp, matrix(nrow=nrow(trials.tmp),      
  ncol=length(gen)-i ) ))
  }
 trials

V1   V2   V3
1  A <NA> <NA>
2  B <NA> <NA>
3  C <NA> <NA>
4  A    B <NA>
5  A    C <NA>
6  B    C <NA>
7  A    B   C 

但我仍然缺少组合 (NA,NA,NA)。如何制作具有所有结果和概率的新数据框。

Answer 1

你可以试试zx8754提到的combn，例如

a=0.95
b=0.82
c=0.78
x <- c(a,b,c)
df <- rbind(t(combn(x, 3)), cbind(t(combn(x, 2)), NA), cbind(t(combn(x, 1)), NA, NA))
apply(df, 1, function(x) prod(x[!is.na(x)]))
[1] 0.60762 0.77900 0.74100 0.63960 0.95000 0.82000 0.78000