【问题标题】:VennDiagram for data gives error?数据的 VennDiagram 给出错误?
【发布时间】:2012-05-02 18:12:53
【问题描述】:

我有以下数据:

Set1    : 82
Set2    : 44
Set3    : 56
Set4    : 53
1,2     : 27
1,3     : 37
1,4     : 30
2,3     : 22
2,4     : 14
3,4     : 19
1,2,3   : 18
1,2,4   : 13
1,3,4   : 20
2,3,4   : 11
1,2,3,4 : 11

1,2 表示Set1 . intersection . Set2,以此类推...

当我尝试使用 draw.quad.venn 为该数据集制作 VennDiagram 时,出现以下错误,

Error in draw.quad.venn(82, 44, 56, 53, 27, 37, 30, 22, 14, 19, 18, 13,  : 
Impossible: partial areas negative

我不明白我做错了什么?

更新:

以下是我使用的命令行:

v.all <- draw.quad.venn( 82, 44, 56, 53, 27, 37, 30, 22, 14, 19, 18, 13, 20, 11, 10, category = c( "Set1", "Set2", "Set3", "Set4" ), fill=c( "red", "blue", "orange", "green" ), cex = 0.75, cat.cex=0.85 )

有趣的是,当我使用时

v <- venneuler( c( A=82, B=44, C=56, D=53, "A&B"=27, "A&C"=37, "A&D"=30, "B&C"=22, "B&D"=14, "C&D"=19, "A&B&C"=18, "A&B&D"=13, "A&C&D"=20, "B&C&D"=11, "A&B&C&D"=11 ))
plot(v)

有效!!

所以现在我真的很困惑。

【问题讨论】:

  • 你应该 a) 显示你用来调用 draw.quad.venn 的代码(我们怎么知道你是否调用不正确?)和 b) 使用 dput 显示此代码的可重现版本数据(我们可以放入自己的代码中进行测试)
  • set1 83 是第一组的总空间吗?还是 set1 的总非重叠空间?如果是总空间,则它小于与集合 2、3 和 4 重叠的区域。
  • @user1317221 总空间为 82 并没有错,因为 2,3,4 是非分离的。我敢打赌,没有查看 ?draw.quad.venn 是因为他的参数顺序错误,或者(根据您的提示)该函数需要非重叠空间,而不是总空间。
  • 很高兴它有效。 venneuler 将A= 82 绘制为非重叠空间。您是否检查了您的draw.quad.venn set1 82 是指@Carl 和我指出的不重叠或总空间?这可能是最初的问题。
  • @user1317221 据我所知,draw.quad.venn Set1 82 是指总空间。但我可能是错的。我想我现在明白@Carl 在说什么了。

标签: r venn-diagram


【解决方案1】:

代码中的约束是直截了当的。正是这里;

     if (any(a1 < 0, a2 < 0, a3 < 0, a4 < 0, a5 < 0, a6 < 0, a7 < 
            0, a8 < 0, a9 < 0, a10 < 0, a11 < 0, a12 < 0, a13 < 0, 
            a14 < 0, a15 < 0)) {
            stop("Impossible: partial areas negative")
        }

【讨论】:

    【解决方案2】:

    没有您的代码就很难提供帮助。由于数据应该重叠的方式,可能无法绘制数据。也许您的数据有错误?

    但是,如果是代码问题,使用您的一些数据和不同的包,我可以提供替代路线:

    require(venneuler)
    m <- as.matrix(c(
    0, 27, 37, 30,
    27, 0, 22, 14,
    37 ,22, 0, 19,
    30,14,19,0)byrow=T,nrow=3)
    v = venneuler(m)
    plot(v)
    

    这给了

    【讨论】:

      猜你喜欢
      • 1970-01-01
      • 2015-07-31
      • 1970-01-01
      • 1970-01-01
      • 1970-01-01
      • 1970-01-01
      • 2015-10-31
      • 1970-01-01
      • 1970-01-01
      相关资源
      最近更新 更多