我正在使用镶嵌套件等学习统计数据和R。对于两条尾巴的测试,如果我可以从分布函数中看到(并得到结果)而不必乘以2,那将很有用。
在当前版本的镶嵌中,有什么方法可以做到这一点?
例如。而不是xpnorm(1.96, lower.tail = FALSE) * 2
,我希望得到相同的结果而不用乘以2并同时遮住两条尾巴。
答案 0 :(得分:0)
xpnorm()有意跟在pnorm()之后,因此每个仅处理一个(类型)尾部。但是,您可以指定多个剪切点(例如c(-1.96,1.96))来获得带有两个阴影的图片。但是您仍然需要一些额外的算法来获得两个尾部概率之和。
library(mosaic)
xpnorm(c(-1.96, 1.96))
#>
#> If X ~ N(0, 1), then
#> P(X <= -1.96) = P(Z <= -1.96) = 0.025 P(X <= 1.96) = P(Z <= 1.96) = 0.975
#> P(X > -1.96) = P(Z > -1.96) = 0.975 P(X > 1.96) = P(Z > 1.96) = 0.025
#>
#> [1] 0.0249979 0.9750021
由reprex package(v0.2.0)于2018-08-10创建。
朝另一个方向发展(如果您知道概率并想找到关键值),我们引入了xcnorm()来查找以指定中心概率为边界的端点。这是一个示例:
library(mosaic)
xcnorm(0.90)
#>
#> If X ~ N(0, 1), then
#> P(X <= -1.644854) = 0.05 P(X <= 1.644854) = 0.95
#> P(X > -1.644854) = 0.95 P(X > 1.644854) = 0.05
#>
#> [1] -1.644854 1.644854
由reprex package(v0.2.0)于2018-08-10创建。