回到这个thread ,我正在计算两个椭圆之间的交点,例如versicolor 和 virginca 交集的体积: 使用以下 mwe:
data(iris)
log.ir <- log(iris[, 1:4])
ir.species <- iris[, 5]
ir.pca <- prcomp(log.ir, center = TRUE, scale. = TRUE)
library(ggbiplot)
g <- ggbiplot(ir.pca, obs.scale = 1, var.scale = 1,
groups = ir.species, ellipse = TRUE,
circle = TRUE)
g <- g + scale_color_discrete(name = '')
g <- g + theme(legend.direction = 'horizontal',
legend.position = 'top')
print(g)
然后我使用重叠{siar}计算椭球体,如下所示:
library(siar)
setosa <- ir.pca$x[ir.species=="setosa",]
versicolor <- ir.pca$x[ir.species=="versicolor",]
virginica <- ir.pca$x[ir.species=="virginica",]
tmp <- overlap(versicolor[,1], versicolor[,2], virginica[,1], virginica[,2], steps = 5)
virginica.percentage <- round(x=(tmp$overlap/tmp$area2*100), digits = 2)
versicolor.percentage <- round(x=(tmp$overlap/tmp$area1*100), digits = 2)
> virginica.percentage [1] 1.54
> versicolor.percentage[1] 1.56
这显然比图中所示的要少得多 2如上所述。
作为ggbiplot code on line 50表示,ggbiplot 默认计算正态概率椭球,概率为“ellipse.prob = 0.68”,即上图中使用的值。
由于documentation重叠 {siar} 将椭圆计算为“小样本量校正标准椭圆”。
因此,我问自己如何修改重叠{siar}来计算类似于 ggbiplot 中绘制/计算的正态概率椭球体。
最佳答案
最新版本的 SIBER 可从 https://github.com/AndrewLJackson/SIBER 获取包括两个新函数来计算用户定义大小的椭圆重叠(即不同的分位数)。
maxLikOverlap()
和
bayesianOverlap()
我建议您使用
从 master 分支安装 # install.packages("devtools") # install if necessary
devtools::install_github("andrewljackson/SIBER@master",
build_vingettes = TRUE)
library(SIBER)
并遵循小插图“椭圆重叠”。如果没有为您安装小插图,您可以直接从 github 存储库下载它们 https://github.com/AndrewLJackson/SIBER/tree/master/vignettes
关于r - 如何调整重叠{siar}中的椭圆重叠计算以匹配ggbiplot中的正态概率椭圆体?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/37349878/