{"id":329852,"date":"2022-02-19T15:00:38","date_gmt":"2022-02-19T15:00:38","guid":{"rendered":"http:\/\/savepearlharbor.com\/?p=329852"},"modified":"-0001-11-30T00:00:00","modified_gmt":"-0001-11-29T21:00:00","slug":"","status":"publish","type":"post","link":"https:\/\/savepearlharbor.com\/?p=329852","title":{"rendered":"\u0420\u0435\u0433\u0440\u0435\u0441\u0441\u0438\u044f \u0438 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 \u0441 \u043d\u0435\u0443\u0441\u0442\u0440\u0430\u043d\u0438\u043c\u044b\u043c\u0438 \u0440\u0430\u0437\u0440\u044b\u0432\u0430\u043c\u0438 \u043f\u0435\u0440\u0432\u043e\u0433\u043e \u0440\u043e\u0434\u0430<\/span>"},"content":{"rendered":"
<\/div>\n
\n
\n
\n
\n

\u041e \u043f\u0430\u043a\u0435\u0442\u0435 BinSeqBstrap<\/p>\n

\u041f\u043e\u0441\u0442\u0430\u043d\u043e\u0432\u043a\u0430 \u0437\u0430\u0434\u0430\u0447\u0438<\/h2>\n

\u0414\u043e\u043f\u0443\u0441\u0442\u0438\u043c, \u0443 \u043d\u0430\u0441 \u0435\u0441\u0442\u044c \u043a\u0430\u043a\u0430\u044f \u043a\u0443\u0441\u043e\u0447\u043d\u043e-\u0433\u043b\u0430\u0434\u043a\u0430\u044f \u0444\u0443\u043d\u043a\u0446\u0438\u044f f(x), \u043a \u043a\u043e\u0442\u043e\u0440\u044b\u0439 \u043f\u0440\u0438\u0431\u0430\u0432\u043b\u0435\u043d \u043d\u0435\u043a\u0438\u0439 \u0441\u043b\u0443\u0447\u0430\u0439\u043d\u044b\u0439 \u0448\u0443\u043c, \u0441\u043e\u043e\u0442\u0432\u0435\u0442\u0441\u0442\u0432\u0443\u044e\u0449\u0438\u0439 \u0443\u0441\u043b\u043e\u0432\u0438\u044f\u043c \u0413\u0430\u0443\u0441\u0441\u0430-\u041c\u0430\u0440\u043a\u043e\u0432\u0430. \u0418 \u0432\u0441\u0435 \u0445\u043e\u0440\u043e\u0448\u043e, \u0442\u043e\u043b\u044c\u043a\u043e \u044d\u0442\u0430 \u0441\u0430\u043c\u0430\u044f \u0444\u0443\u043d\u043a\u0446\u0438\u044f f(x) \u2013 \u0444\u0443\u043d\u043a\u0446\u0438\u044f \u0441 \u043d\u0435\u0443\u0441\u0442\u0440\u0430\u043d\u0438\u043c\u044b\u043c(-\u0438)\u00a0 \u0440\u0430\u0437\u0440\u044b\u0432\u043e\u043c(-\u0430\u043c\u0438) \u043f\u0435\u0440\u0432\u043e\u0433\u043e \u0440\u043e\u0434\u0430, \u0442\u043e \u0435\u0441\u0442\u044c \u0432 \u043a\u0430\u043a\u043e\u0439-\u0442\u043e \u0442\u043e\u0447\u043a\u0435 \u043b\u0435\u0432\u044b\u0439 \u0438 \u043f\u0440\u0430\u0432\u044b\u0439 \u043f\u0440\u0435\u0434\u0435\u043b \u044d\u0442\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 \u0440\u0430\u0432\u043d\u044b \u0440\u0430\u0437\u043d\u044b\u043c \u0447\u0438\u0441\u043b\u0430\u043c, \u0430 \u0443 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 \u0435\u0441\u0442\u044c \u0441\u043a\u0430\u0447\u043e\u043a. \u0417\u0430\u0434\u0430\u0447\u0430 \u2013 \u043a\u0430\u043a-\u0442\u043e \u043d\u0443\u0436\u043d\u043e \u043d\u0430\u0443\u0447\u0438\u0442\u044c \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c \u0440\u0430\u0441\u043f\u043e\u0437\u043d\u0430\u0432\u0430\u0442\u044c \u044d\u0442\u043e\u0442 \u0441\u043a\u0430\u0447\u043e\u043a.<\/p>\n

\u041c\u0438\u043d\u0443\u0442\u043a\u0430 \u0442\u0435\u043e\u0440\u0438\u0438<\/h2>\n

\u0422\u0435\u043e\u0440\u0435\u0442\u0438\u0447\u0435\u0441\u043a\u0438\u0435 \u043e\u0441\u043d\u043e\u0432\u044b \u0438\u0437\u043b\u043e\u0436\u0435\u043d\u044b \u0432 \u0432\u0438\u043d\u044c\u0435\u0442\u043a\u0435, \u043d\u0430\u043f\u0438\u0441\u0430\u043d\u043d\u043e\u0439 \u041a\u044d\u0442\u0438 \u041c\u0430\u043a\u0414\u044d\u0439\u0434 \u0438 \u0424\u043b\u043e\u0440\u0438\u0430\u043d\u0430 \u041f\u044d\u0439\u043d\u0430 \u0438\u0437 \u041a\u044d\u043c\u0431\u0440\u0438\u0434\u0436\u0430, \u043e\u043f\u0443\u0431\u043b\u0438\u043a\u043e\u0432\u0430\u043d\u043d\u043e\u0439 \u0432\u043e\u0442 \u0442\u0443\u0442 https:\/\/cran.r-project.org\/web\/packages\/BinSegBstrap\/vignettes\/BinSegBstrap.pdf<\/a><\/p>\n

\u0415\u0441\u043b\u0438 \u0432\u043a\u0440\u0430\u0442\u0446\u0435, \u0442\u043e \u0444\u0443\u043d\u043a\u0446\u0438\u043e\u043d\u0430\u043b \u043f\u0430\u043a\u0435\u0442\u0430 \u043f\u043e\u0437\u0432\u043e\u043b\u044f\u0435\u0442 \u0440\u0435\u0448\u0438\u0442\u044c 4 \u0440\u0430\u0437\u043d\u044b\u0435 \u043f\u043e\u0434\u0437\u0430\u0434\u0430\u0447\u0438, \u0431\u0430\u0437\u043e\u0432\u0430\u044f \u0438\u0437 \u043a\u043e\u0442\u043e\u0440\u044b\u0445 — \u043e\u043f\u0440\u0435\u0434\u0435\u043b\u0435\u043d\u0438\u0435 \u043c\u0435\u0441\u0442\u0430 \u0438 \u0440\u0430\u0437\u043c\u0435\u0440\u0430 \u0435\u0434\u0438\u043d\u0441\u0442\u0432\u0435\u043d\u043d\u043e\u0433\u043e \u0441\u043a\u0430\u0447\u043a\u0430 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 \u043f\u0440\u0438 \u0438\u0437\u0432\u0435\u0441\u0442\u043d\u043e\u0439 \u0448\u0438\u0440\u0438\u043d\u0435 \u043e\u043a\u043d\u0430 h<\/p>\n

\u041e\u0441\u043d\u043e\u0432\u0430: \u0432\u043e\u0442 \u044d\u0442\u0430 \u0432\u043e\u0442 \u0434\u043b\u0438\u043d\u043d\u0430\u044f \u0444\u043e\u0440\u043c\u0443\u043b\u0430<\/p>\n

<\/figcaption><\/figure>\n

\u0420\u0435\u0437\u0443\u043b\u044c\u0442\u0430\u0442 \u044d\u0442\u043e\u0439 \u0444\u043e\u0440\u043c\u0443\u043b\u044b \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u0435\u0442\u0441\u044f \u0434\u043b\u044f \u0432\u044b\u0447\u0438\u0441\u043b\u0435\u043d\u0438\u044f \u043c\u0435\u0441\u0442\u0430 \u0441\u043a\u0430\u0447\u043a\u0430 <\/p>\n

\"n
n \u2013 \u043e\u0431\u0449\u0435\u0435 \u043a\u043e\u043b\u0438\u0447\u0435\u0441\u0442\u0432\u043e \u0442\u043e\u0447\u0435\u043a<\/figcaption><\/figure>\n

\u0421\u0430\u043c \u043c\u0435\u0442\u043e\u0434 \u043e\u0441\u043d\u043e\u0432\u0430\u043d \u043d\u0430 \u0438\u0434\u0435\u044f\u0445 \u043d\u0435\u043f\u0430\u0440\u0430\u043c\u0435\u0442\u0440\u0438\u0447\u0435\u0441\u043a\u043e\u0439 \u0440\u0435\u0433\u0440\u0435\u0441\u0441\u0438\u0438. \u0414\u043b\u044f \u043e\u043f\u0440\u0435\u0434\u0435\u043b\u0435\u043d\u0438\u044f \u043c\u0435\u0441\u0442\u043e\u043f\u043e\u043b\u043e\u0436\u0435\u043d\u0438\u044f \u0442\u043e\u0447\u0435\u043a \u0440\u0430\u0437\u0440\u044b\u0432\u0430 \u0432\u044b\u0431\u0438\u0440\u0430\u0435\u0442\u0441\u044f \u043e\u043a\u043d\u043e \u0438\u0437 \u043d\u0435\u043a\u043e\u0442\u043e\u0440\u043e\u0433\u043e \u043a\u043e\u043b\u0438\u0447\u0435\u0441\u0442\u0432\u0430 \u0442\u043e\u0447\u0435\u043a, \u0438 \u0441\u0447\u0438\u0442\u0430\u0435\u0442\u0441\u044f \u0440\u0430\u0437\u043d\u0438\u0446\u0430 \u043c\u0435\u0436\u0434\u0443 \u0441\u0443\u043c\u043c\u0430\u043c\u0438 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0439 \u0441\u043b\u0435\u0432\u0430 \u0438 \u0441\u043f\u0440\u0430\u0432\u0430. \u0414\u043b\u044f \u043e\u043f\u0440\u0435\u0434\u0435\u043b\u0435\u043d\u0438\u044f \u0440\u0430\u0437\u043c\u0435\u0440\u0430 \u0441\u043a\u0430\u0447\u043a\u0430 \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u0435\u0442\u0441\u044f \u0440\u0430\u0437\u043d\u0438\u0446\u0430 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0439 \u0441\u0433\u043b\u0430\u0436\u0435\u043d\u043d\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 \u0432 \u043e\u043f\u0440\u0435\u0434\u0435\u043b\u0435\u043d\u043d\u043e\u0439 \u0442\u043e\u0447\u0435\u043a \u0440\u0430\u0437\u0440\u044b\u0432\u0430.<\/p>\n

\u041f\u0440\u0430\u043a\u0442\u0438\u043a\u0430<\/h2>\n

\u0412 \u043a\u043e\u0434\u0435 \u044d\u0442\u043e \u0432\u044b\u0433\u043b\u044f\u0434\u0438\u0442 \u0442\u0430\u043a:<\/p>\n

library(BinSegBstrap) ## \u0427\u0430\u0441\u0442\u044c 1 set.seed(1) n <- 1:100 signal <- 0.1*n-5 signal[51:100] <- signal[51:100] + 5  y <- rnorm(n) + signal # \u041f\u0440\u0438\u0434\u0443\u043c\u044b\u0432\u0430\u0435\u043c \u0438\u0441\u043a\u0443\u0441\u0441\u0442\u0432\u0435\u043d\u043d\u044b\u0435 \u0434\u0430\u043d\u043d\u044b\u0435 est <- estimateSingleCp(y = y, bandwidth = 0.1) est$cp # \u043e\u043f\u0440\u0435\u0434\u0435\u043b\u044f\u0435\u043c \u0442\u043e\u0447\u043a\u0443 \u0440\u0430\u0437\u0440\u044b\u0432\u0430 est$size # \u043e\u043f\u0440\u0435\u0434\u0435\u043b\u044f\u0435\u043c \u0441\u043a\u0430\u0447\u043e\u043a \u0444\u0443\u043d\u043a\u0446\u0438\u0438 plot(y, pch = 16, col = \"grey30\") lines(signal) lines(est$est, col = \"red\") # est$est - \u043f\u043e\u0434\u043e\u0433\u043d\u0430\u043d\u043d\u044b\u0435 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u044f<\/code><\/pre>\n
\"\u041c\u0435\u0441\u0442\u043e\u043f\u043e\u043b\u043e\u0436\u0435\u043d\u0438\u0435
\u041c\u0435\u0441\u0442\u043e\u043f\u043e\u043b\u043e\u0436\u0435\u043d\u0438\u0435 \u0441\u043a\u0430\u0447\u043a\u0430 \u0442\u043e\u0447\u043d\u043e \u0443\u0441\u0442\u0430\u043d\u043e\u0432\u043b\u0435\u043d\u043e, \u0435\u0433\u043e \u0440\u0430\u0437\u043c\u0435\u0440 \u2013 \u043f\u0440\u0438\u043c\u0435\u0440\u043d\u043e.<\/figcaption><\/figure>\n
\"\u0410
\u0410 \u0432\u043e\u0442 \u0433\u0440\u0430\u0444\u0438\u043a\u0430. \u041a\u0440\u0430\u0441\u043d\u0430\u044f \u043b\u0438\u043d\u0438\u044f — \u044d\u0442\u043e \u0441\u0433\u043b\u0430\u0436\u0435\u043d\u043d\u0430\u044f \u0440\u0435\u0433\u0440\u0435\u0441\u0441\u0438\u044f, \u0447\u0435\u0440\u043d\u0430\u044f — \u043f\u043e\u0434\u043e\u0431\u0440\u0430\u043d\u043d\u044b\u0435 \u043b\u0438\u043d\u0438\u0438 \u0443\u0440\u0430\u0432\u043d\u0435\u043d\u0438\u0439<\/figcaption><\/figure>\n
\u0427\u0430\u0441\u0442\u044c 2. \u041e\u043f\u0440\u0435\u0434\u0435\u043b\u0435\u043d\u0438\u0435 \u043f\u0430\u0440\u0430\u043c\u0435\u0442\u0440\u043e\u0432 \u0441\u043a\u0430\u0447\u043a\u0430 \u043f\u0440\u0438 \u043d\u0435\u0438\u0437\u0432\u0435\u0441\u0442\u043d\u043e\u0439 \u0432\u0435\u043b\u0438\u0447\u0438\u043d\u0435 <\/p>\n
\u0421\u043e\u0431\u0441\u0442\u0432\u0435\u043d\u043d\u043e, \u043a \u0438\u0441\u0445\u043e\u0434\u043d\u043e\u0439 \u0444\u043e\u0440\u043c\u0443\u043b\u0435 \u0434\u043e\u0441\u0442\u0440\u0430\u0438\u0432\u0430\u0435\u0442\u0441\u044f \u0442\u043e\u043b\u044c\u043a\u043e \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c \u043a\u0440\u043e\u0441\u0441-\u0432\u0430\u043b\u0438\u0434\u0430\u0446\u0438\u0438, \u043a\u043e\u0442\u043e\u0440\u044b\u0439 \u043f\u043e\u0437\u0432\u043e\u043b\u044f\u0435\u0442 \u043d\u0430\u0439\u0442\u0438 \u043b\u0443\u0447\u0448\u0438\u0439 \u0440\u0430\u0437\u043c\u0435\u0440 \u043e\u043a\u043d\u0430, \u0432 \u043a\u043e\u0442\u043e\u0440\u043e\u043c \u0438\u0441\u043a\u0430\u0442\u044c \u0441\u043a\u0430\u0447\u043e\u043a.<\/p>\n
\u0417\u0430\u0434\u0430\u0434\u0438\u043c \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c\u0443 \u0437\u0430\u0434\u0430\u0447\u043a\u0443 \u043f\u043e\u0441\u043b\u043e\u0436\u043d\u0435\u0435:<\/p>\n
set.seed(1) n <- 1:100 signal1 <- 5-0.1*n signal2 <- 0.01*n^2-20 signal<-c(signal1[1:50],signal2[51:100]) y <- rnorm(n) + signal est <- estimateSingleCp(y = y) est$bandwidth # \u043f\u043e\u0434\u043e\u0431\u0440\u0430\u043d\u043d\u0430\u044f \u0448\u0438\u0440\u0438\u043d\u0430 \u043e\u043a\u043d\u0430 est$cp # \u043e\u043f\u0440\u0435\u0434\u0435\u043b\u044f\u0435\u043c \u0442\u043e\u0447\u043a\u0443 \u0440\u0430\u0437\u0440\u044b\u0432\u0430 est$size # \u043e\u043f\u0440\u0435\u0434\u0435\u043b\u044f\u0435\u043c \u0441\u043a\u0430\u0447\u043e\u043a \u0444\u0443\u043d\u043a\u0446\u0438\u0438 plot(y, pch = 16, col = \"grey30\") # plot of observations, true and estimated signal lines(signal) lines(est$est, col = \"red\")<\/code><\/pre>\n
\"\u0421\u043d\u043e\u0432\u0430
\u0421\u043d\u043e\u0432\u0430 \u0432\u0441\u0435 \u043e\u043f\u0440\u0435\u0434\u0435\u043b\u0435\u043d\u043e \u0432\u0435\u0440\u043d\u043e<\/figcaption><\/figure>\n
\u0427\u0430\u0441\u0442\u044c 3. \u041f\u0440\u043e\u0432\u0435\u0440\u043a\u0430 \u0441\u0442\u0430\u0442\u0438\u0441\u0442\u0438\u0447\u0435\u0441\u043a\u043e\u0439 \u0433\u0438\u043f\u043e\u0442\u0435\u0437\u044b \u043e \u043d\u0430\u043b\u0438\u0447\u0438\u0438 \u0441\u043a\u0430\u0447\u043a\u0430<\/p>\n
\u041e\u0431\u044a\u044f\u0432\u043b\u044f\u0435\u043c \u043d\u0443\u043b\u0435\u0432\u0443\u044e \u0433\u0438\u043f\u043e\u0442\u0435\u0437\u0443 \u043e \u0442\u043e\u043c, \u0447\u0442\u043e \u0440\u0430\u0437\u043c\u0435\u0440 \u0441\u043a\u0430\u0447\u043a\u0430 \u041a \u0440\u0430\u0432\u0435\u043d 0 (\u0430\u043b\u044c\u0442\u0435\u0440\u043d\u0430\u0442\u0438\u0432\u043d\u0430\u044f \u2013 \u043d\u0435 \u0440\u0430\u0432\u0435\u043d).<\/p>\n
\u0414\u0430\u043b\u044c\u0448\u0435 \u0441 \u043f\u043e\u043c\u043e\u0449\u044c\u044e \u0431\u0443\u0442\u0441\u0442\u0440\u044d\u043f\u0430 \u043f\u043e\u043b\u0443\u0447\u0430\u0435\u043c \u043c\u043d\u043e\u0436\u0435\u0441\u0442\u0432\u043e \u043e\u0446\u0435\u043d\u043e\u043a \u0432\u0435\u043b\u0438\u0447\u0438\u043d\u044b \u044d\u0442\u043e\u0433\u043e \u0441\u043a\u0430\u0447\u043a\u0430 (\u043c\u043d\u043e\u0433\u043e-\u043c\u043d\u043e\u0433\u043e \u0440\u0430\u0437 \u0440\u0435\u0448\u0430\u0435\u043c \u0447\u0443\u0442\u044c \u0440\u0430\u0437\u043d\u044b\u0435 \u0437\u0430\u0434\u0430\u0447\u0438 \u0438\u0437 \u043f.2) \u0438 \u0441\u043c\u043e\u0442\u0440\u0438\u043c, \u0432\u044b\u043f\u043e\u043b\u043d\u044f\u0435\u0442\u0441\u044f \u0433\u0438\u043f\u043e\u0442\u0435\u0437\u0430 \u0438\u043b\u0438 \u043d\u0435\u0442.<\/p>\n
\u0412 \u043a\u043e\u0434\u0435 \u044d\u0442\u043e \u0432\u044b\u0433\u043b\u044f\u0434\u0438\u0442 \u0442\u0430\u043a:<\/p>\n
test <- BstrapTest(y = y) test$outcome # \u0420\u0435\u0437\u0443\u043b\u044c\u0442\u0430\u0442 \u043f\u0440\u043e\u0432\u0435\u0440\u043a\u0438 \u0433\u0438\u043f\u043e\u0442\u0435\u0437\u044b \u043e \u043d\u0443\u043b\u0435\u0432\u043e\u043c \u0441\u043a\u0430\u0447\u043a\u0435 test$pValue<\/code><\/pre>\n
\"\u0418\u0442\u043e\u0433\u043e
\u0418\u0442\u043e\u0433\u043e \u2013 \u043d\u0443\u043b\u0435\u0432\u0430\u044f \u0433\u0438\u043f\u043e\u0442\u0435\u0437\u0430 \u043d\u0430 \u0443\u0440\u043e\u0432\u043d\u0435 \u0437\u043d\u0430\u0447\u0438\u043c\u043e\u0441\u0442\u0438 0.05 \u043e\u0442\u0432\u0435\u0440\u0433\u0430\u0435\u0442\u0441\u044f<\/figcaption><\/figure>\n
\u0427\u0430\u0441\u0442\u044c 4. \u041f\u043e\u0438\u0441\u043a \u0432\u0441\u0435\u0445 \u0432\u043e\u0437\u043c\u043e\u0436\u043d\u044b\u0445 \u0442\u043e\u0447\u0435\u043a \u0440\u0430\u0437\u0440\u044b\u0432\u0430<\/p>\n
\u0421\u043e\u0441\u0442\u0430\u0432\u0438\u043c \u0444\u0443\u043d\u043a\u0446\u0438\u044e, \u0432 \u043a\u043e\u0442\u043e\u0440\u043e\u0439 \u0435\u0441\u0442\u044c \u0434\u0432\u0435 \u0442\u043e\u0447\u043a\u0438 \u0440\u0430\u0437\u0440\u044b\u0432\u0430, \u0438 \u043f\u043e\u0441\u043c\u043e\u0442\u0440\u0438\u043c, \u0441\u043f\u0440\u0430\u0432\u0438\u0442\u0441\u044f \u043b\u0438 \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c \u0441 \u043d\u0435\u0439<\/p>\n
set.seed(1) signal1<-abs(-8-seq(from=-10, to=-5, length.out=50)) signal2<-3*exp(seq(from=-5, to=0, length.out=50)) signal3<-sqrt(seq(from=0, to=10, length.out=50))-4 signal <- c(signal1,signal2,signal3) y <- rnorm(150) + signal est <- BinSegBstrap(y = y) est$cps # \u043d\u0430\u0445\u043e\u0434\u0438\u043c \u0442\u043e\u0447\u043a\u0438 \u0440\u0430\u0437\u0440\u044b\u0432\u0430 plot(y, pch = 16, col = \"grey30\") # plot of observations, true and estimated signal lines(signal) lines(est$est, col = \"red\")<\/code><\/pre>\n
\"\u041d\u0430
\u041d\u0430 \u0443\u0434\u0438\u0432\u043b\u0435\u043d\u0438\u0435, \u0432\u0441\u0435 \u043d\u0430\u0448\u043b\u043e\u0441\u044c<\/figcaption><\/figure>\n
\u0412\u043c\u0435\u0441\u0442\u043e \u0437\u0430\u043a\u043b\u044e\u0447\u0435\u043d\u0438\u044f<\/h2>\n
\u0412\u0435\u0441\u044c\u043c\u0430 \u0438\u043d\u0442\u0435\u0440\u0435\u0441\u043d\u0430\u044f \u0437\u0430\u0434\u0443\u043c\u043a\u0430. \u0414\u043b\u044f \u0440\u0435\u0448\u0435\u043d\u0438\u044f \u043f\u0440\u0430\u043a\u0442\u0438\u0447\u0435\u0441\u043a\u0438\u0445 \u0437\u0430\u0434\u0430\u0447 \u0432 \u0442\u0435\u043a\u0443\u0449\u0435\u043c \u0441\u043e\u0441\u0442\u043e\u044f\u043d\u0438\u0438 \u043f\u0430\u043a\u0435\u0442 \u043f\u0440\u0438\u043c\u0435\u043d\u0438\u043c \u043c\u0430\u043b\u043e, \u043a \u0442\u043e\u043c\u0443 \u0436\u0435 \u0440\u0435\u0434\u043a\u043e \u043f\u0440\u0438\u0445\u043e\u0434\u0438\u0442\u0441\u044f \u0434\u0430\u043d\u043d\u044b\u0435 \u0430\u043f\u043f\u0440\u043e\u043a\u0441\u0438\u043c\u0438\u0440\u043e\u0432\u0430\u0442\u044c \u0444\u0443\u043d\u043a\u0446\u0438\u044f\u043c\u0438 \u0441 \u0440\u0430\u0437\u0440\u044b\u0432\u0430\u043c\u0438, \u043d\u043e \u043a\u0430\u043a \u0438\u043d\u0441\u0442\u0440\u0443\u043c\u0435\u043d\u0442 \u0432\u0441\u043f\u043e\u043c\u043e\u0433\u0430\u0442\u0435\u043b\u044c\u043d\u043e\u0433\u043e \u0430\u043d\u0430\u043b\u0438\u0437\u0430 \u043c\u043e\u0436\u0435\u0442 \u043f\u0440\u0438\u0433\u043e\u0434\u0438\u0442\u044c\u0441\u044f. \u041d\u0443 \u0438 \u0438\u043d\u0442\u0435\u0440\u0435\u0441\u043d\u043e \u043f\u043e\u0442\u0435\u0441\u0442\u0438\u0442\u044c, \u0434\u043e \u043a\u0430\u043a\u0438\u0445 \u043f\u0440\u0435\u0434\u0435\u043b\u043e\u0432 \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c \u0440\u0430\u0431\u043e\u0442\u0430\u0435\u0442 \u0445\u043e\u0440\u043e\u0448\u043e, \u0430 \u043a\u043e\u0433\u0434\u0430 \u043d\u0430\u0447\u043d\u0435\u0442 \u043f\u0440\u043e\u043f\u0443\u0441\u043a\u0430\u0442\u044c \u0440\u0430\u0437\u0440\u044b\u0432\u044b \/ \u043d\u0430\u0445\u043e\u0434\u0438\u0442\u044c \u043b\u0438\u0448\u043d\u0438\u0435.<\/p>\n<\/div>\n<\/div>\n<\/div>\n
<\/div>\n<\/div>\n
  
 \u0441\u0441\u044b\u043b\u043a\u0430 \u043d\u0430 \u043e\u0440\u0438\u0433\u0438\u043d\u0430\u043b \u0441\u0442\u0430\u0442\u044c\u0438  https:\/\/habr.com\/ru\/post\/652491\/<\/a><\/p>\n","protected":false},"excerpt":{"rendered":"
<\/div>\n
\n
\n
\n
\n
\u041e \u043f\u0430\u043a\u0435\u0442\u0435 BinSeqBstrap<\/p>\n
\u041f\u043e\u0441\u0442\u0430\u043d\u043e\u0432\u043a\u0430 \u0437\u0430\u0434\u0430\u0447\u0438<\/h2>\n
\u0414\u043e\u043f\u0443\u0441\u0442\u0438\u043c, \u0443 \u043d\u0430\u0441 \u0435\u0441\u0442\u044c \u043a\u0430\u043a\u0430\u044f \u043a\u0443\u0441\u043e\u0447\u043d\u043e-\u0433\u043b\u0430\u0434\u043a\u0430\u044f \u0444\u0443\u043d\u043a\u0446\u0438\u044f f(x), \u043a \u043a\u043e\u0442\u043e\u0440\u044b\u0439 \u043f\u0440\u0438\u0431\u0430\u0432\u043b\u0435\u043d \u043d\u0435\u043a\u0438\u0439 \u0441\u043b\u0443\u0447\u0430\u0439\u043d\u044b\u0439 \u0448\u0443\u043c, \u0441\u043e\u043e\u0442\u0432\u0435\u0442\u0441\u0442\u0432\u0443\u044e\u0449\u0438\u0439 \u0443\u0441\u043b\u043e\u0432\u0438\u044f\u043c \u0413\u0430\u0443\u0441\u0441\u0430-\u041c\u0430\u0440\u043a\u043e\u0432\u0430. \u0418 \u0432\u0441\u0435 \u0445\u043e\u0440\u043e\u0448\u043e, \u0442\u043e\u043b\u044c\u043a\u043e \u044d\u0442\u0430 \u0441\u0430\u043c\u0430\u044f \u0444\u0443\u043d\u043a\u0446\u0438\u044f f(x) \u2013 \u0444\u0443\u043d\u043a\u0446\u0438\u044f \u0441 \u043d\u0435\u0443\u0441\u0442\u0440\u0430\u043d\u0438\u043c\u044b\u043c(-\u0438)\u00a0 \u0440\u0430\u0437\u0440\u044b\u0432\u043e\u043c(-\u0430\u043c\u0438) \u043f\u0435\u0440\u0432\u043e\u0433\u043e \u0440\u043e\u0434\u0430, \u0442\u043e \u0435\u0441\u0442\u044c \u0432 \u043a\u0430\u043a\u043e\u0439-\u0442\u043e \u0442\u043e\u0447\u043a\u0435 \u043b\u0435\u0432\u044b\u0439 \u0438 \u043f\u0440\u0430\u0432\u044b\u0439 \u043f\u0440\u0435\u0434\u0435\u043b \u044d\u0442\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 \u0440\u0430\u0432\u043d\u044b \u0440\u0430\u0437\u043d\u044b\u043c \u0447\u0438\u0441\u043b\u0430\u043c, \u0430 \u0443 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 \u0435\u0441\u0442\u044c \u0441\u043a\u0430\u0447\u043e\u043a. \u0417\u0430\u0434\u0430\u0447\u0430 \u2013 \u043a\u0430\u043a-\u0442\u043e \u043d\u0443\u0436\u043d\u043e \u043d\u0430\u0443\u0447\u0438\u0442\u044c \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c \u0440\u0430\u0441\u043f\u043e\u0437\u043d\u0430\u0432\u0430\u0442\u044c \u044d\u0442\u043e\u0442 \u0441\u043a\u0430\u0447\u043e\u043a.<\/p>\n
\u041c\u0438\u043d\u0443\u0442\u043a\u0430 \u0442\u0435\u043e\u0440\u0438\u0438<\/h2>\n
\u0422\u0435\u043e\u0440\u0435\u0442\u0438\u0447\u0435\u0441\u043a\u0438\u0435 \u043e\u0441\u043d\u043e\u0432\u044b \u0438\u0437\u043b\u043e\u0436\u0435\u043d\u044b \u0432 \u0432\u0438\u043d\u044c\u0435\u0442\u043a\u0435, \u043d\u0430\u043f\u0438\u0441\u0430\u043d\u043d\u043e\u0439 \u041a\u044d\u0442\u0438 \u041c\u0430\u043a\u0414\u044d\u0439\u0434 \u0438 \u0424\u043b\u043e\u0440\u0438\u0430\u043d\u0430 \u041f\u044d\u0439\u043d\u0430 \u0438\u0437 \u041a\u044d\u043c\u0431\u0440\u0438\u0434\u0436\u0430, \u043e\u043f\u0443\u0431\u043b\u0438\u043a\u043e\u0432\u0430\u043d\u043d\u043e\u0439 \u0432\u043e\u0442 \u0442\u0443\u0442 https:\/\/cran.r-project.org\/web\/packages\/BinSegBstrap\/vignettes\/BinSegBstrap.pdf<\/a><\/p>\n
\u0415\u0441\u043b\u0438 \u0432\u043a\u0440\u0430\u0442\u0446\u0435, \u0442\u043e \u0444\u0443\u043d\u043a\u0446\u0438\u043e\u043d\u0430\u043b \u043f\u0430\u043a\u0435\u0442\u0430 \u043f\u043e\u0437\u0432\u043e\u043b\u044f\u0435\u0442 \u0440\u0435\u0448\u0438\u0442\u044c 4 \u0440\u0430\u0437\u043d\u044b\u0435 \u043f\u043e\u0434\u0437\u0430\u0434\u0430\u0447\u0438, \u0431\u0430\u0437\u043e\u0432\u0430\u044f \u0438\u0437 \u043a\u043e\u0442\u043e\u0440\u044b\u0445 — \u043e\u043f\u0440\u0435\u0434\u0435\u043b\u0435\u043d\u0438\u0435 \u043c\u0435\u0441\u0442\u0430 \u0438 \u0440\u0430\u0437\u043c\u0435\u0440\u0430 \u0435\u0434\u0438\u043d\u0441\u0442\u0432\u0435\u043d\u043d\u043e\u0433\u043e \u0441\u043a\u0430\u0447\u043a\u0430 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 \u043f\u0440\u0438 \u0438\u0437\u0432\u0435\u0441\u0442\u043d\u043e\u0439 \u0448\u0438\u0440\u0438\u043d\u0435 \u043e\u043a\u043d\u0430 h<\/p>\n
\u041e\u0441\u043d\u043e\u0432\u0430: \u0432\u043e\u0442 \u044d\u0442\u0430 \u0432\u043e\u0442 \u0434\u043b\u0438\u043d\u043d\u0430\u044f \u0444\u043e\u0440\u043c\u0443\u043b\u0430<\/p>\n
<\/figcaption><\/figure>\n
\u0420\u0435\u0437\u0443\u043b\u044c\u0442\u0430\u0442 \u044d\u0442\u043e\u0439 \u0444\u043e\u0440\u043c\u0443\u043b\u044b \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u0435\u0442\u0441\u044f \u0434\u043b\u044f \u0432\u044b\u0447\u0438\u0441\u043b\u0435\u043d\u0438\u044f \u043c\u0435\u0441\u0442\u0430 \u0441\u043a\u0430\u0447\u043a\u0430 <\/p>\n
n \u2013 \u043e\u0431\u0449\u0435\u0435 \u043a\u043e\u043b\u0438\u0447\u0435\u0441\u0442\u0432\u043e \u0442\u043e\u0447\u0435\u043a<\/figcaption><\/figure>\n
\u0421\u0430\u043c \u043c\u0435\u0442\u043e\u0434 \u043e\u0441\u043d\u043e\u0432\u0430\u043d \u043d\u0430 \u0438\u0434\u0435\u044f\u0445 \u043d\u0435\u043f\u0430\u0440\u0430\u043c\u0435\u0442\u0440\u0438\u0447\u0435\u0441\u043a\u043e\u0439 \u0440\u0435\u0433\u0440\u0435\u0441\u0441\u0438\u0438. \u0414\u043b\u044f \u043e\u043f\u0440\u0435\u0434\u0435\u043b\u0435\u043d\u0438\u044f \u043c\u0435\u0441\u0442\u043e\u043f\u043e\u043b\u043e\u0436\u0435\u043d\u0438\u044f \u0442\u043e\u0447\u0435\u043a \u0440\u0430\u0437\u0440\u044b\u0432\u0430 \u0432\u044b\u0431\u0438\u0440\u0430\u0435\u0442\u0441\u044f \u043e\u043a\u043d\u043e \u0438\u0437 \u043d\u0435\u043a\u043e\u0442\u043e\u0440\u043e\u0433\u043e \u043a\u043e\u043b\u0438\u0447\u0435\u0441\u0442\u0432\u0430 \u0442\u043e\u0447\u0435\u043a, \u0438 \u0441\u0447\u0438\u0442\u0430\u0435\u0442\u0441\u044f \u0440\u0430\u0437\u043d\u0438\u0446\u0430 \u043c\u0435\u0436\u0434\u0443 \u0441\u0443\u043c\u043c\u0430\u043c\u0438 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0439 \u0441\u043b\u0435\u0432\u0430 \u0438 \u0441\u043f\u0440\u0430\u0432\u0430. \u0414\u043b\u044f \u043e\u043f\u0440\u0435\u0434\u0435\u043b\u0435\u043d\u0438\u044f \u0440\u0430\u0437\u043c\u0435\u0440\u0430 \u0441\u043a\u0430\u0447\u043a\u0430 \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u0435\u0442\u0441\u044f \u0440\u0430\u0437\u043d\u0438\u0446\u0430 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0439 \u0441\u0433\u043b\u0430\u0436\u0435\u043d\u043d\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 \u0432 \u043e\u043f\u0440\u0435\u0434\u0435\u043b\u0435\u043d\u043d\u043e\u0439 \u0442\u043e\u0447\u0435\u043a \u0440\u0430\u0437\u0440\u044b\u0432\u0430.<\/p>\n
\u041f\u0440\u0430\u043a\u0442\u0438\u043a\u0430<\/h2>\n
\u0412 \u043a\u043e\u0434\u0435 \u044d\u0442\u043e \u0432\u044b\u0433\u043b\u044f\u0434\u0438\u0442 \u0442\u0430\u043a:<\/p>\n
library(BinSegBstrap) ## \u0427\u0430\u0441\u0442\u044c 1 set.seed(1) n <- 1:100 signal <- 0.1*n-5 signal[51:100] <- signal[51:100] + 5  y <- rnorm(n) + signal # \u041f\u0440\u0438\u0434\u0443\u043c\u044b\u0432\u0430\u0435\u043c \u0438\u0441\u043a\u0443\u0441\u0441\u0442\u0432\u0435\u043d\u043d\u044b\u0435 \u0434\u0430\u043d\u043d\u044b\u0435 est <- estimateSingleCp(y = y, bandwidth = 0.1) est$cp # \u043e\u043f\u0440\u0435\u0434\u0435\u043b\u044f\u0435\u043c \u0442\u043e\u0447\u043a\u0443 \u0440\u0430\u0437\u0440\u044b\u0432\u0430 est$size # \u043e\u043f\u0440\u0435\u0434\u0435\u043b\u044f\u0435\u043c \u0441\u043a\u0430\u0447\u043e\u043a \u0444\u0443\u043d\u043a\u0446\u0438\u0438 plot(y, pch = 16, col = \"grey30\") lines(signal) lines(est$est, col = \"red\") # est$est - \u043f\u043e\u0434\u043e\u0433\u043d\u0430\u043d\u043d\u044b\u0435 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u044f<\/code><\/pre>\n
\u041c\u0435\u0441\u0442\u043e\u043f\u043e\u043b\u043e\u0436\u0435\u043d\u0438\u0435 \u0441\u043a\u0430\u0447\u043a\u0430 \u0442\u043e\u0447\u043d\u043e \u0443\u0441\u0442\u0430\u043d\u043e\u0432\u043b\u0435\u043d\u043e, \u0435\u0433\u043e \u0440\u0430\u0437\u043c\u0435\u0440 \u2013 \u043f\u0440\u0438\u043c\u0435\u0440\u043d\u043e.<\/figcaption><\/figure>\n
\u0410 \u0432\u043e\u0442 \u0433\u0440\u0430\u0444\u0438\u043a\u0430. \u041a\u0440\u0430\u0441\u043d\u0430\u044f \u043b\u0438\u043d\u0438\u044f — \u044d\u0442\u043e \u0441\u0433\u043b\u0430\u0436\u0435\u043d\u043d\u0430\u044f \u0440\u0435\u0433\u0440\u0435\u0441\u0441\u0438\u044f, \u0447\u0435\u0440\u043d\u0430\u044f — \u043f\u043e\u0434\u043e\u0431\u0440\u0430\u043d\u043d\u044b\u0435 \u043b\u0438\u043d\u0438\u0438 \u0443\u0440\u0430\u0432\u043d\u0435\u043d\u0438\u0439<\/figcaption><\/figure>\n
\u0427\u0430\u0441\u0442\u044c 2. \u041e\u043f\u0440\u0435\u0434\u0435\u043b\u0435\u043d\u0438\u0435 \u043f\u0430\u0440\u0430\u043c\u0435\u0442\u0440\u043e\u0432 \u0441\u043a\u0430\u0447\u043a\u0430 \u043f\u0440\u0438 \u043d\u0435\u0438\u0437\u0432\u0435\u0441\u0442\u043d\u043e\u0439 \u0432\u0435\u043b\u0438\u0447\u0438\u043d\u0435 <\/p>\n
\u0421\u043e\u0431\u0441\u0442\u0432\u0435\u043d\u043d\u043e, \u043a \u0438\u0441\u0445\u043e\u0434\u043d\u043e\u0439 \u0444\u043e\u0440\u043c\u0443\u043b\u0435 \u0434\u043e\u0441\u0442\u0440\u0430\u0438\u0432\u0430\u0435\u0442\u0441\u044f \u0442\u043e\u043b\u044c\u043a\u043e \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c \u043a\u0440\u043e\u0441\u0441-\u0432\u0430\u043b\u0438\u0434\u0430\u0446\u0438\u0438, \u043a\u043e\u0442\u043e\u0440\u044b\u0439 \u043f\u043e\u0437\u0432\u043e\u043b\u044f\u0435\u0442 \u043d\u0430\u0439\u0442\u0438 \u043b\u0443\u0447\u0448\u0438\u0439 \u0440\u0430\u0437\u043c\u0435\u0440 \u043e\u043a\u043d\u0430, \u0432 \u043a\u043e\u0442\u043e\u0440\u043e\u043c \u0438\u0441\u043a\u0430\u0442\u044c \u0441\u043a\u0430\u0447\u043e\u043a.<\/p>\n
\u0417\u0430\u0434\u0430\u0434\u0438\u043c \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c\u0443 \u0437\u0430\u0434\u0430\u0447\u043a\u0443 \u043f\u043e\u0441\u043b\u043e\u0436\u043d\u0435\u0435:<\/p>\n
set.seed(1) n <- 1:100 signal1 <- 5-0.1*n signal2 <- 0.01*n^2-20 signal<-c(signal1[1:50],signal2[51:100]) y <- rnorm(n) + signal est <- estimateSingleCp(y = y) est$bandwidth # \u043f\u043e\u0434\u043e\u0431\u0440\u0430\u043d\u043d\u0430\u044f \u0448\u0438\u0440\u0438\u043d\u0430 \u043e\u043a\u043d\u0430 est$cp # \u043e\u043f\u0440\u0435\u0434\u0435\u043b\u044f\u0435\u043c \u0442\u043e\u0447\u043a\u0443 \u0440\u0430\u0437\u0440\u044b\u0432\u0430 est$size # \u043e\u043f\u0440\u0435\u0434\u0435\u043b\u044f\u0435\u043c \u0441\u043a\u0430\u0447\u043e\u043a \u0444\u0443\u043d\u043a\u0446\u0438\u0438 plot(y, pch = 16, col = \"grey30\") # plot of observations, true and estimated signal lines(signal) lines(est$est, col = \"red\")<\/code><\/pre>\n
\u0421\u043d\u043e\u0432\u0430 \u0432\u0441\u0435 \u043e\u043f\u0440\u0435\u0434\u0435\u043b\u0435\u043d\u043e \u0432\u0435\u0440\u043d\u043e<\/figcaption><\/figure>\n
\u0427\u0430\u0441\u0442\u044c 3. \u041f\u0440\u043e\u0432\u0435\u0440\u043a\u0430 \u0441\u0442\u0430\u0442\u0438\u0441\u0442\u0438\u0447\u0435\u0441\u043a\u043e\u0439 \u0433\u0438\u043f\u043e\u0442\u0435\u0437\u044b \u043e \u043d\u0430\u043b\u0438\u0447\u0438\u0438 \u0441\u043a\u0430\u0447\u043a\u0430<\/p>\n
\u041e\u0431\u044a\u044f\u0432\u043b\u044f\u0435\u043c \u043d\u0443\u043b\u0435\u0432\u0443\u044e \u0433\u0438\u043f\u043e\u0442\u0435\u0437\u0443 \u043e \u0442\u043e\u043c, \u0447\u0442\u043e \u0440\u0430\u0437\u043c\u0435\u0440 \u0441\u043a\u0430\u0447\u043a\u0430 \u041a \u0440\u0430\u0432\u0435\u043d 0 (\u0430\u043b\u044c\u0442\u0435\u0440\u043d\u0430\u0442\u0438\u0432\u043d\u0430\u044f \u2013 \u043d\u0435 \u0440\u0430\u0432\u0435\u043d).<\/p>\n
\u0414\u0430\u043b\u044c\u0448\u0435 \u0441 \u043f\u043e\u043c\u043e\u0449\u044c\u044e \u0431\u0443\u0442\u0441\u0442\u0440\u044d\u043f\u0430 \u043f\u043e\u043b\u0443\u0447\u0430\u0435\u043c \u043c\u043d\u043e\u0436\u0435\u0441\u0442\u0432\u043e \u043e\u0446\u0435\u043d\u043e\u043a \u0432\u0435\u043b\u0438\u0447\u0438\u043d\u044b \u044d\u0442\u043e\u0433\u043e \u0441\u043a\u0430\u0447\u043a\u0430 (\u043c\u043d\u043e\u0433\u043e-\u043c\u043d\u043e\u0433\u043e \u0440\u0430\u0437 \u0440\u0435\u0448\u0430\u0435\u043c \u0447\u0443\u0442\u044c \u0440\u0430\u0437\u043d\u044b\u0435 \u0437\u0430\u0434\u0430\u0447\u0438 \u0438\u0437 \u043f.2) \u0438 \u0441\u043c\u043e\u0442\u0440\u0438\u043c, \u0432\u044b\u043f\u043e\u043b\u043d\u044f\u0435\u0442\u0441\u044f \u0433\u0438\u043f\u043e\u0442\u0435\u0437\u0430 \u0438\u043b\u0438 \u043d\u0435\u0442.<\/p>\n
\u0412 \u043a\u043e\u0434\u0435 \u044d\u0442\u043e \u0432\u044b\u0433\u043b\u044f\u0434\u0438\u0442 \u0442\u0430\u043a:<\/p>\n
test <- BstrapTest(y = y) test$outcome # \u0420\u0435\u0437\u0443\u043b\u044c\u0442\u0430\u0442 \u043f\u0440\u043e\u0432\u0435\u0440\u043a\u0438 \u0433\u0438\u043f\u043e\u0442\u0435\u0437\u044b \u043e \u043d\u0443\u043b\u0435\u0432\u043e\u043c \u0441\u043a\u0430\u0447\u043a\u0435 test$pValue<\/code><\/pre>\n
\u0418\u0442\u043e\u0433\u043e \u2013 \u043d\u0443\u043b\u0435\u0432\u0430\u044f \u0433\u0438\u043f\u043e\u0442\u0435\u0437\u0430 \u043d\u0430 \u0443\u0440\u043e\u0432\u043d\u0435 \u0437\u043d\u0430\u0447\u0438\u043c\u043e\u0441\u0442\u0438 0.05 \u043e\u0442\u0432\u0435\u0440\u0433\u0430\u0435\u0442\u0441\u044f<\/figcaption><\/figure>\n
\u0427\u0430\u0441\u0442\u044c 4. \u041f\u043e\u0438\u0441\u043a \u0432\u0441\u0435\u0445 \u0432\u043e\u0437\u043c\u043e\u0436\u043d\u044b\u0445 \u0442\u043e\u0447\u0435\u043a \u0440\u0430\u0437\u0440\u044b\u0432\u0430<\/p>\n
\u0421\u043e\u0441\u0442\u0430\u0432\u0438\u043c \u0444\u0443\u043d\u043a\u0446\u0438\u044e, \u0432 \u043a\u043e\u0442\u043e\u0440\u043e\u0439 \u0435\u0441\u0442\u044c \u0434\u0432\u0435 \u0442\u043e\u0447\u043a\u0438 \u0440\u0430\u0437\u0440\u044b\u0432\u0430, \u0438 \u043f\u043e\u0441\u043c\u043e\u0442\u0440\u0438\u043c, \u0441\u043f\u0440\u0430\u0432\u0438\u0442\u0441\u044f \u043b\u0438 \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c \u0441 \u043d\u0435\u0439<\/p>\n
set.seed(1) signal1<-abs(-8-seq(from=-10, to=-5, length.out=50)) signal2<-3*exp(seq(from=-5, to=0, length.out=50)) signal3<-sqrt(seq(from=0, to=10, length.out=50))-4 signal <- c(signal1,signal2,signal3) y <- rnorm(150) + signal est <- BinSegBstrap(y = y) est$cps # \u043d\u0430\u0445\u043e\u0434\u0438\u043c \u0442\u043e\u0447\u043a\u0438 \u0440\u0430\u0437\u0440\u044b\u0432\u0430 plot(y, pch = 16, col = \"grey30\") # plot of observations, true and estimated signal lines(signal) lines(est$est, col = \"red\")<\/code><\/pre>\n
\u041d\u0430 \u0443\u0434\u0438\u0432\u043b\u0435\u043d\u0438\u0435, \u0432\u0441\u0435 \u043d\u0430\u0448\u043b\u043e\u0441\u044c<\/figcaption><\/figure>\n
\u0412\u043c\u0435\u0441\u0442\u043e \u0437\u0430\u043a\u043b\u044e\u0447\u0435\u043d\u0438\u044f<\/h2>\n
\u0412\u0435\u0441\u044c\u043c\u0430 \u0438\u043d\u0442\u0435\u0440\u0435\u0441\u043d\u0430\u044f \u0437\u0430\u0434\u0443\u043c\u043a\u0430. \u0414\u043b\u044f \u0440\u0435\u0448\u0435\u043d\u0438\u044f \u043f\u0440\u0430\u043a\u0442\u0438\u0447\u0435\u0441\u043a\u0438\u0445 \u0437\u0430\u0434\u0430\u0447 \u0432 \u0442\u0435\u043a\u0443\u0449\u0435\u043c \u0441\u043e\u0441\u0442\u043e\u044f\u043d\u0438\u0438 \u043f\u0430\u043a\u0435\u0442 \u043f\u0440\u0438\u043c\u0435\u043d\u0438\u043c \u043c\u0430\u043b\u043e, \u043a \u0442\u043e\u043c\u0443 \u0436\u0435 \u0440\u0435\u0434\u043a\u043e \u043f\u0440\u0438\u0445\u043e\u0434\u0438\u0442\u0441\u044f \u0434\u0430\u043d\u043d\u044b\u0435 \u0430\u043f\u043f\u0440\u043e\u043a\u0441\u0438\u043c\u0438\u0440\u043e\u0432\u0430\u0442\u044c \u0444\u0443\u043d\u043a\u0446\u0438\u044f\u043c\u0438 \u0441 \u0440\u0430\u0437\u0440\u044b\u0432\u0430\u043c\u0438, \u043d\u043e \u043a\u0430\u043a \u0438\u043d\u0441\u0442\u0440\u0443\u043c\u0435\u043d\u0442 \u0432\u0441\u043f\u043e\u043c\u043e\u0433\u0430\u0442\u0435\u043b\u044c\u043d\u043e\u0433\u043e \u0430\u043d\u0430\u043b\u0438\u0437\u0430 \u043c\u043e\u0436\u0435\u0442 \u043f\u0440\u0438\u0433\u043e\u0434\u0438\u0442\u044c\u0441\u044f. \u041d\u0443 \u0438 \u0438\u043d\u0442\u0435\u0440\u0435\u0441\u043d\u043e \u043f\u043e\u0442\u0435\u0441\u0442\u0438\u0442\u044c, \u0434\u043e \u043a\u0430\u043a\u0438\u0445 \u043f\u0440\u0435\u0434\u0435\u043b\u043e\u0432 \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c \u0440\u0430\u0431\u043e\u0442\u0430\u0435\u0442 \u0445\u043e\u0440\u043e\u0448\u043e, \u0430 \u043a\u043e\u0433\u0434\u0430 \u043d\u0430\u0447\u043d\u0435\u0442 \u043f\u0440\u043e\u043f\u0443\u0441\u043a\u0430\u0442\u044c \u0440\u0430\u0437\u0440\u044b\u0432\u044b \/ \u043d\u0430\u0445\u043e\u0434\u0438\u0442\u044c \u043b\u0438\u0448\u043d\u0438\u0435.<\/p>\n<\/div>\n<\/div>\n<\/div>\n
<\/div>\n<\/div>\n
  
 \u0441\u0441\u044b\u043b\u043a\u0430 \u043d\u0430 \u043e\u0440\u0438\u0433\u0438\u043d\u0430\u043b \u0441\u0442\u0430\u0442\u044c\u0438  https:\/\/habr.com\/ru\/post\/652491\/<\/a>
<\/br><\/br><\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[],"tags":[],"class_list":["post-329852","post","type-post","status-publish","format-standard","hentry"],"_links":{"self":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/posts\/329852","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=329852"}],"version-history":[{"count":0,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/posts\/329852\/revisions"}],"wp:attachment":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=329852"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=329852"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=329852"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}