MSQC :: mult.chart игнорирует аргумент xlim

Question

Я вычислил Hotelling T2 по нормализованному набору данных и получил следующий график в R:

На графике я хочу изучить интервал50-100 по оси X тесно.Есть ли в R функция или метод, с помощью которых это можно сделать?Благодарю.

Answer 1

К сожалению, MSQC::mult.chart не поддерживает аргумент xlim.Но вы можете извлечь значения для построения и воспроизведения вручную.Смотрите пример ниже.

library(MSQC)
data(dowel1)
# default
mult.chart(dowel1, type = "chi", alpha = 0.05)
#> [[1]]
#> [1] "Chi-squared Control Chart"
#> 
#> $ucl
#> [1] 5.99
#> 
#> $t2
#>       [,1]
#>  [1,] 1.62
#> ...
#> 
#> $Xmv
#> [1] 0.5 1.0
#> 
#> $covariance
#>         [,1]    [,2]
#> [1,] 4.9e-05 8.6e-05
#> [2,] 8.6e-05 4.2e-04

# manual
mc <- mult.chart(dowel1, type = "chi", alpha = 0.05)
plot(seq_along(mc$t2), mc$t2, ylim = c(0, mc$ucl), type = "l")
points(seq_along(mc$t2), mc$t2)
abline(h = mc$ucl, col = 2)

# restricted
plot(seq_along(mc$t2), mc$t2, ylim = c(0, mc$ucl), type = "l", xlim = c(5, 20))
points(seq_along(mc$t2), mc$t2)
abline(h = mc$ucl, col = 2)

^{Создано в 2019-02-10 при представлении пакета (v0.2.1)}