Вопрос:
Мне нужно использовать эту формулу для вычисления верхнего предела и нижнего предела для нижнего вектора. Пожалуйста помоги..
> a <- c(2150, 2261, 2018, 2217, 2074, 1986, 1662, 2324, 1591, 1700, 1809, 1824, 2026, 2169, 2072, 1463, 1484, 841, 851, 3742, 889, 905, 1203, 967, 1776, 1075, 1031) > a [1] 2150 2261 2018 2217 2074 1986 1662 2324 1591 1700 1809 1824 2026 2169 2072 [16] 1463 1484 841 851 3742 889 905 1203 967 1776 1075 1031 Лучший ответ:
Это просто простая формула R:
UL = mean(a) + 2.66*sum(abs(diff(a)))/length(a) LL = mean(a) — 2.66*sum(abs(diff(a)))/length(a)
Трудная часть | X_i – X_ (i-1) | и поэтому мы используем несколько встроенных функций R:
- diff дает X_i – X_ (i-1)
- abs дает абсолютное значение
первая часть формулы – это среднее значение, а + / – просто разделяется.
Ответ №1 ctrl.limits <- c( LCL = mean(a) — 2.66*sum(abs(diff(a)))/length(a), UCL = mean(a) + 2.66*sum(abs(diff(a)))/length(a) ) #> ctrl.limits # LCL UCL # 495.1133 2920.4422
Почему 2,66? Он немного запутан, но сводится к тому, что он равен 3 (количество заданных стандартных отклонений в любом направлении), умноженное на отношение ожидаемого диапазона размера выборки 2, взятого от стандартного нормального до его стандартного отклонения (= 1). ‘d2’ является функцией размера выборки. E (Ri) = d2 (n) σ
(См. Стр. 55 (с надписью 47) www.win.tue.nl/~adibucch/2WS10/SPClecturenotes.pdf для вывода.)