Два хвоста в распределительных участках

Question

Я изучаю статистику и R, используя пакет мозаики, среди прочего. Для двухстороннего теста было бы полезно, если бы я мог видеть (и получать результаты) из функций распределения без умножения на 2.

Есть ли способ сделать это в текущей версии мозаики?

Например. вместо xpnorm(1.96, lower.tail = FALSE) * 2 Я хотел бы получить тот же результат без умножения на 2 и затенения обоих хвостов.

Answer 1

xpnorm() намеренно следует pnorm(), поэтому каждый обрабатывает только один (тип) tail. Однако вы можете указать несколько точек разреза (скажем, 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

Создано 2018-08-10 пакетом представлением (v0.2.0).

Двигаясь в другом направлении (если вы знаете вероятность и хотите найти критические значения), мы ввели 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

Создано 2018-08-10 пакетом представлением (v0.2.0).