{"id":188244,"date":"2013-07-30T01:31:03","date_gmt":"2013-07-29T21:31:03","guid":{"rendered":"http:\/\/savepearlharbor.com\/?p=188244"},"modified":"-0001-11-30T00:00:00","modified_gmt":"-0001-11-29T21:00:00","slug":"","status":"publish","type":"post","link":"https:\/\/savepearlharbor.com\/?p=188244","title":{"rendered":"<span class=\"post_title\">\u0421\u043a\u0440\u044b\u0442\u044b\u0435 \u0446\u0435\u043f\u0438 \u041c\u0430\u0440\u043a\u043e\u0432\u0430, \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c \u0411\u0430\u0443\u043c\u0430-\u0412\u0435\u043b\u0448\u0430<\/span>"},"content":{"rendered":"<div class=\"content html_format\">   \t\u0421\u043a\u0440\u044b\u0442\u044b\u0435 \u043c\u043e\u0434\u0435\u043b\u0438\/\u0446\u0435\u043f\u0438 \u041c\u0430\u0440\u043a\u043e\u0432\u0430 \u043e\u0434\u043d\u0438 \u0438\u0437 \u043f\u043e\u0434\u0445\u043e\u0434\u043e\u0432 \u043a \u043f\u0440\u0435\u0434\u0441\u0442\u0430\u0432\u043b\u0435\u043d\u0438\u044e \u0434\u0430\u043d\u043d\u044b\u0445. \u041c\u043d\u0435 \u043e\u0447\u0435\u043d\u044c \u043f\u043e\u043d\u0440\u0430\u0432\u0438\u043b\u043e\u0441\u044c \u043a\u0430\u043a \u043e\u0431\u043e\u0431\u0449\u0430\u0435\u0442\u0441\u044f \u043c\u043d\u043e\u0436\u0435\u0441\u0442\u0432\u043e \u0442\u0430\u043a\u0438\u0445 \u043f\u043e\u0434\u0445\u043e\u0434\u043e\u0432 \u0432 <a href=\"http:\/\/habrahabr.ru\/company\/surfingbird\/blog\/177889\/\">\u044d\u0442\u043e\u0439 \u0441\u0442\u0430\u0442\u044c\u0435<\/a>.<\/p>\n<p>  \u0412 \u043f\u0440\u043e\u0434\u043e\u043b\u0436\u0435\u043d\u0438\u0435 \u0436\u0435 \u043c\u043e\u0435\u0439 <a href=\"http:\/\/habrahabr.ru\/post\/180109\/\">\u043f\u0440\u0435\u0434\u044b\u0434\u0443\u0449\u0435\u0439 \u0441\u0442\u0430\u0442\u044c\u0438<\/a> \u043e\u043f\u0438\u0441\u0430\u043d\u0438\u044f \u0441\u043a\u0440\u044b\u0442\u044b\u0445 \u043c\u043e\u0434\u0435\u043b\u0435\u0439 \u041c\u0430\u0440\u043a\u043e\u0432\u0430, \u0437\u0430\u0434\u0430\u043c\u0441\u044f \u0432\u043e\u043f\u0440\u043e\u0441\u043e\u043c: \u043e\u0442\u043a\u0443\u0434\u0430 \u0432\u0437\u044f\u0442\u044c \u0445\u043e\u0440\u043e\u0448\u0443\u044e \u043c\u043e\u0434\u0435\u043b\u044c? \u041e\u0442\u0432\u0435\u0442 \u0434\u043e\u0441\u0442\u0430\u0442\u043e\u0447\u043d\u043e \u0441\u0442\u0430\u043d\u0434\u0430\u0440\u0442\u0435\u043d, \u0432\u0437\u044f\u0442\u044c \u043d\u0435\u043f\u043b\u043e\u0445\u0443\u044e \u043c\u043e\u0434\u0435\u043b\u044c \u0438 \u0441\u0434\u0435\u043b\u0430\u0442\u044c \u0438\u0437 \u043d\u0435\u0435 \u0445\u043e\u0440\u043e\u0448\u0443\u044e.<\/p>\n<p>  \u041d\u0430\u043f\u043e\u043c\u043d\u044e \u043f\u0440\u0438\u043c\u0435\u0440: \u043d\u0430\u043c \u043d\u0443\u0436\u043d\u043e \u0440\u0435\u0430\u043b\u0438\u0437\u043e\u0432\u0430\u0442\u044c \u0434\u0435\u0442\u0435\u043a\u0442\u043e\u0440 \u043b\u0436\u0438, \u043a\u043e\u0442\u043e\u0440\u044b\u0439 \u043f\u043e \u043f\u043e\u0434\u0440\u0430\u0433\u0438\u0432\u0430\u043d\u0438\u044e \u0440\u0443\u043a \u0447\u0435\u043b\u043e\u0432\u0435\u043a\u0430, \u043e\u043f\u0440\u0435\u0434\u0435\u043b\u044f\u0435\u0442, \u0433\u043e\u0432\u043e\u0440\u0438\u0442 \u043e\u043d \u043f\u0440\u0430\u0432\u0434\u0443 \u0438\u043b\u0438 \u043d\u0435\u0442. \u0414\u043e\u043f\u0443\u0441\u0442\u0438\u043c, \u043a\u043e\u0433\u0434\u0430 \u0447\u0435\u043b\u043e\u0432\u0435\u043a \u043b\u0436\u0435\u0442, \u0440\u0443\u043a\u0438 \u0442\u0440\u044f\u0441\u0443\u0442\u0441\u044f \u0447\u0443\u0442\u044c \u0431\u043e\u043b\u044c\u0448\u0435, \u043d\u043e \u043d\u0430\u043c \u043d\u0435 \u0438\u0437\u0432\u0435\u0441\u0442\u043d\u043e \u043d\u0430 \u0441\u043a\u043e\u043b\u044c\u043a\u043e \u0438\u043c\u0435\u043d\u043d\u043e. \u0412\u043e\u0437\u044c\u043c\u0435\u043c \u043c\u043e\u0434\u0435\u043b\u044c \u043d\u0430\u043e\u0431\u0443\u043c, \u043f\u0440\u043e\u0433\u043e\u043d\u0438\u043c \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c \u0412\u0438\u0442\u0435\u0440\u0431\u0438 \u0438\u0437 \u043f\u0440\u0435\u0434\u044b\u0434\u0443\u0449\u0435\u0439 \u0441\u0442\u0430\u0442\u044c\u0438 \u0438 \u043f\u043e\u043b\u0443\u0447\u0438\u043c \u0434\u043e\u0432\u043e\u043b\u044c\u043d\u043e \u0441\u0442\u0440\u0430\u043d\u043d\u044b\u0435 \u0440\u0435\u0437\u0443\u043b\u044c\u0442\u0430\u0442\u044b:<br \/>  <img decoding=\"async\" src=\"http:\/\/habrastorage.org\/storage2\/e4c\/622\/5f1\/e4c6225f162e1efd9387f0191c1e2b86.png\"\/><br \/>  <a name=\"habracut\"><\/a><br \/>  \u041d\u0438\u043a\u0430\u043a\u043e\u0433\u043e \u0441\u043c\u044b\u0441\u043b\u0430 \u0432 \u043f\u043e\u043b\u0443\u0447\u0438\u0432\u0448\u0435\u043c\u0441\u044f \u0440\u0430\u0437\u0431\u0438\u0435\u043d\u0438\u0438 \u043d\u0435\u0442; \u0442\u043e\u0447\u043d\u0435\u0435, \u0435\u0441\u0442\u044c, \u043d\u043e \u043d\u0435 \u0442\u043e \u043a\u043e\u0442\u043e\u0440\u044b\u0439 \u0445\u043e\u0442\u0435\u043b\u043e\u0441\u044c \u0431\u044b \u0432\u0438\u0434\u0435\u0442\u044c. \u041f\u0435\u0440\u0432\u044b\u0439 \u0432\u043e\u043f\u0440\u043e\u0441: \u043a\u0430\u043a \u043e\u0446\u0435\u043d\u0438\u0442\u044c \u0432\u0435\u0440\u043e\u044f\u0442\u043d\u043e\u0441\u0442\u044c \u043f\u043e\u044f\u0432\u043b\u0435\u043d\u0438\u044f \u0438\u043c\u0435\u044e\u0449\u0435\u0439\u0441\u044f \u0432\u044b\u0431\u043e\u0440\u043a\u0438 \u0432 \u0432\u044b\u0431\u0440\u0430\u043d\u043d\u043e\u0439 \u043c\u043e\u0434\u0435\u043b\u0438. \u0418 \u0432\u0442\u043e\u0440\u043e\u0439 \u0432\u043e\u043f\u0440\u043e\u0441, \u043a\u0430\u043a \u043d\u0430\u0439\u0442\u0438 \u0442\u0430\u043a\u0443\u044e \u043c\u043e\u0434\u0435\u043b\u044c, \u0432 \u043a\u043e\u0442\u043e\u0440\u043e\u0439 \u0438\u043c\u0435\u044e\u0449\u0430\u044f\u0441\u044f \u0432\u044b\u0431\u043e\u0440\u043a\u0430 \u0431\u044b\u043b\u0430 \u0431\u044b \u043d\u0430\u0438\u0431\u043e\u043b\u0435\u0435 \u043f\u0440\u0430\u0432\u0434\u043e\u043f\u043e\u0434\u043e\u0431\u043d\u0430, \u0442\u043e \u0435\u0441\u0442\u044c \u043a\u0430\u043a \u043e\u0446\u0435\u043d\u0438\u0442\u044c \u043a\u0430\u0447\u0435\u0441\u0442\u0432\u043e \u0438 \u043a\u0430\u043a \u043e\u0440\u0433\u0430\u043d\u0438\u0437\u043e\u0432\u0430\u0442\u044c \u0430\u0434\u0430\u043f\u0442\u0438\u0432\u043d\u043e\u0441\u0442\u044c \u043a \u043d\u043e\u0432\u044b\u043c \u0434\u0430\u043d\u043d\u044b\u043c. \u041a\u0430\u043a \u0440\u0430\u0437 \u043d\u0430 \u043d\u0438\u0445 \u0438 \u043e\u0442\u0432\u0435\u0447\u0430\u044e\u0442 the forward-backward procedure \u0438 the Baum-Welch reestimation procedure.<\/p>\n<p>  The forward-backward procedure \u043f\u043e\u0437\u0432\u043e\u043b\u044f\u0435\u0442 \u043e\u0446\u0435\u043d\u0438\u0442\u044c \u0432\u0435\u0440\u043e\u044f\u0442\u043d\u043e\u0441\u0442\u044c \u0442\u043e\u0433\u043e, \u0447\u0442\u043e \u0432 \u0434\u0430\u043d\u043d\u043e\u0439 \u043c\u043e\u0434\u0435\u043b\u0438 HMM \u043f\u043e\u044f\u0432\u0438\u0442\u0441\u044f \u0432\u044b\u0431\u043e\u0440\u043a\u0430 \u043d\u0430\u0431\u043b\u044e\u0434\u0430\u0435\u043c\u044b\u0445 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0439 <img decoding=\"async\" src=\"http:\/\/habr.habrastorage.org\/post_images\/36b\/656\/2a8\/36b6562a88c03a8fac4b242f1b324f9e.png\" alt=\"O=(o_1,\\ldots,o_{T_m})\"\/> (\u041e\u0431\u043e\u0437\u043d\u0430\u0447\u0435\u043d\u0438\u044f \u0432\u0437\u044f\u0442\u044b \u0438\u0437 \u043f\u0440\u0435\u0434\u044b\u0434\u0443\u0449\u0435\u0439 \u0441\u0442\u0430\u0442\u044c\u0438). \u041e\u0431\u043e\u0437\u043d\u0430\u0447\u0438\u043c \u043e\u0446\u0435\u043d\u043e\u0447\u043d\u0443\u044e \u0444\u0443\u043d\u043a\u0446\u0438\u044e \u043a\u0430\u043a <img decoding=\"async\" src=\"http:\/\/habr.habrastorage.org\/post_images\/2b0\/bd8\/821\/2b0bd8821cfba0ca0973a23f06ee3882.png\" alt=\"P(O|HMM)\"\/>. \u0418 \u0432\u0432\u0435\u0434\u0435\u043c \u0434\u0432\u0435 \u0434\u043e\u043f\u043e\u043b\u043d\u0438\u0442\u0435\u043b\u044c\u043d\u044b\u0435 \u0444\u0443\u043d\u043a\u0446\u0438\u0438:<\/p>\n<pre><code>AlphaFunc&lt;-function(S,I,T,P,O,Tm) {   la&lt;-I*Es(P,O[1])   alpha&lt;-cbind(la)   for(t in 2:Tm) {     la&lt;-as.vector((la*Es(P,O[t]))%*%T)     alpha&lt;-cbind(alpha,la)   }   return (alpha) }  BetaFunc&lt;-function(S,I,T,P,O,Tm) {   lb&lt;-rep(1,length(S))   beta&lt;-cbind(lb)   for(t in Tm:2) {     lb&lt;-as.vector(T%*%(Es(P,O[t])*lb))     beta&lt;-cbind(lb, beta)   }   return (beta) } <\/code><\/pre>\n<p>  \u0424\u0443\u043d\u043a\u0446\u0438\u0438 <img decoding=\"async\" src=\"http:\/\/habr.habrastorage.org\/post_images\/120\/241\/fe4\/120241fe47b639bdfc2972c5f20b35bf.png\" alt=\"\\alpha_t (s)\"\/>, <img decoding=\"async\" src=\"http:\/\/habr.habrastorage.org\/post_images\/417\/07d\/007\/41707d0075722e7faebd372015105e94.png\" alt=\"\\beta_t (s)\"\/> \u043e\u0431\u043b\u0430\u0434\u0430\u044e\u0442 \u0441\u043b\u0435\u0434\u0443\u044e\u0449\u0438\u043c\u0438 \u0438\u043d\u0442\u0435\u0440\u0435\u0441\u043d\u044b\u043c\u0438 \u0441\u0432\u043e\u0439\u0441\u0442\u0432\u0430\u043c\u0438:<\/p>\n<p>  <img decoding=\"async\" src=\"http:\/\/habr.habrastorage.org\/post_images\/5ad\/210\/975\/5ad210975f48774a6e7cf05905bec343.png\" alt=\"\\sum_{s\\in S} \\beta_1(s)\\pi_s f_s(o_1) = \\sum_{s\\in S} \\alpha_{T_m}(s) = P(O|HMM)=\\sum_{s\\in S} \\alpha_t(s)\\beta_t(s),\\forall t \\in \\overline{1,T_m}.\"\/><br \/>  \u0418\u0442\u0430\u043a, \u043c\u044b \u0443\u0436\u0435 \u0443\u043c\u0435\u0435\u043c \u043e\u0446\u0435\u043d\u0438\u0432\u0430\u0442\u044c \u043f\u043e\u044f\u0432\u043b\u0435\u043d\u0438\u0435 \u0432\u044b\u0431\u043e\u0440\u043a\u0438 O \u0432 \u043c\u043e\u0434\u0435\u043b\u0438 HMM \u0438 \u043d\u0430\u0445\u043e\u0434\u0438\u0442\u044c \u043f\u043e\u0441\u043b\u0435\u0434\u043e\u0432\u0430\u0442\u0435\u043b\u044c\u043d\u043e\u0441\u0442\u044c \u0441\u043a\u0440\u044b\u0442\u044b\u0445 \u0441\u043e\u0441\u0442\u043e\u044f\u043d\u0438\u0439 <img decoding=\"async\" src=\"http:\/\/habr.habrastorage.org\/post_images\/8a2\/b60\/ebf\/8a2b60ebfacb96a0f5706e9203496976.png\" alt=\"Q^*\"\/>. \u041e\u0442\u043a\u0440\u044b\u0442\u044b\u043c \u043e\u0441\u0442\u0430\u0435\u0442\u0441\u044f \u0432\u043e\u043f\u0440\u043e\u0441, \u043a\u0430\u043a \u043d\u0430\u0439\u0442\u0438 \u043c\u043e\u0434\u0435\u043b\u044c, \u0432 \u043a\u043e\u0442\u043e\u0440\u043e\u0439 \u0432\u0435\u0440\u043e\u044f\u0442\u043d\u043e\u0441\u0442\u044c \u043f\u043e\u044f\u0432\u043b\u0435\u043d\u0438\u044f \u0438\u043c\u0435\u044e\u0449\u0435\u0439\u0441\u044f \u0432\u044b\u0431\u043e\u0440\u043a\u0438 \u0431\u0443\u0434\u0435\u0442 \u043c\u0430\u043a\u0441\u0438\u043c\u0430\u043b\u044c\u043d\u043e\u0439, \u0442\u043e \u0435\u0441\u0442\u044c <img decoding=\"async\" src=\"http:\/\/habr.habrastorage.org\/post_images\/87c\/29f\/339\/87c29f339efbd8e4dbcdeea69c6e7a69.png\" alt=\"P(O|HMM^*)=\\max\\limits_{HMM\\in\\Omega}P(O|HMM)\"\/>, \u0433\u0434\u0435 <img decoding=\"async\" src=\"http:\/\/habr.habrastorage.org\/post_images\/26e\/1ba\/d5a\/26e1bad5aa12d05a634e00702403776e.png\" alt=\"\\Omega\"\/> \u043d\u0435\u043a\u043e\u0442\u043e\u0440\u043e\u0435 \u043c\u043d\u043e\u0436\u0435\u0441\u0442\u0432\u043e \u0434\u043e\u043f\u0443\u0441\u0442\u0438\u043c\u044b\u0445 \u043c\u043e\u0434\u0435\u043b\u0435\u0439. \u041a\u0430\u043a \u0443\u0442\u0432\u0435\u0440\u0436\u0434\u0430\u0435\u0442\u0441\u044f, \u0444\u043e\u0440\u043c\u0430\u043b\u044c\u043d\u043e\u0433\u043e \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c\u0430 \u043d\u0430\u0445\u043e\u0436\u0434\u0435\u043d\u0438\u044f \u043e\u043f\u0442\u0438\u043c\u0443\u043c\u0430 \u043d\u0435\u0442, \u043d\u043e \u0437\u0430\u0447\u0430\u0441\u0442\u0443\u044e \u043c\u043e\u0436\u043d\u043e \u0445\u043e\u0442\u044f \u0431\u044b \u043d\u0430 \u043a\u0430\u043f\u0435\u043b\u044c\u043a\u0443 \u043f\u0440\u0438\u0431\u043b\u0438\u0437\u0438\u0442\u044c\u0441\u044f \u043a \u00ab\u0441\u043e\u0432\u0435\u0440\u0448\u0435\u043d\u0441\u0442\u0432\u0443\u00bb, \u0437\u0430 \u0441\u0447\u0451\u0442 \u043d\u0435\u043a\u043e\u0442\u043e\u0440\u043e\u0439 \u043c\u043e\u0434\u0438\u0444\u0438\u043a\u0430\u0446\u0438\u0438 \u043c\u043e\u0434\u0435\u043b\u0438. \u0418\u043c\u0435\u043d\u043d\u043e \u0442\u0430\u043a \u0438 \u0440\u0430\u0431\u043e\u0442\u0430\u044e\u0442 \u0438\u0442\u0435\u0440\u0430\u0446\u0438\u043e\u043d\u043d\u044b\u0435 <a href=\"http:\/\/ru.wikipedia.org\/wiki\/EM-%E0%EB%E3%EE%F0%E8%F2%EC\">EM-\u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c\u044b<\/a>.<\/p>\n<p>  \u0418\u0442\u0430\u043a, \u043f\u0443\u0441\u0442\u044c \u0435\u0441\u0442\u044c \u043d\u0430\u0447\u0430\u043b\u044c\u043d\u0430\u044f \u043c\u043e\u0434\u0435\u043b\u044c <img decoding=\"async\" src=\"http:\/\/habr.habrastorage.org\/post_images\/fd6\/b43\/898\/fd6b4389832a0e5f709efe6c0f7bb46a.png\" alt=\"HMM^k\"\/>\u043d\u0435\u043e\u0431\u0445\u043e\u0434\u0438\u043c\u043e \u043d\u0430\u0439\u0442\u0438 \u0442\u0430\u043a\u0443\u044e \u043c\u043e\u0434\u0435\u043b\u044c <img decoding=\"async\" src=\"http:\/\/habr.habrastorage.org\/post_images\/36b\/745\/9a8\/36b7459a8e508df22d3aa8af1eacc778.png\" alt=\"HMM^{k+1}\"\/>, \u0447\u0442\u043e <img decoding=\"async\" src=\"http:\/\/habr.habrastorage.org\/post_images\/750\/117\/d31\/750117d31106be143d7f87b079b3980b.png\" alt=\"P(O|HMM^{k+1} )\\geq P(O|HMM^{k})\"\/>. \u041f\u0440\u0438\u0441\u0442\u0443\u043f\u0438\u043c.<\/p>\n<p>  \u0412\u0432\u0435\u0434\u0435\u043c \u0435\u0449\u0435 \u0434\u0432\u0435 \u0432\u0441\u043f\u043e\u043c\u043e\u0433\u0430\u0442\u0435\u043b\u044c\u043d\u044b\u0435 \u0444\u0443\u043d\u043a\u0446\u0438\u0438:<\/p>\n<pre><code>xiGamma&lt;-function(t,alpha,beta,S,T,P,O) {   xi&lt;-matrix(0,nrow=length(S),ncol=length(S))   denom&lt;-alpha[,t]%*%beta[,t]   for(sin in S) {     for(sout in S) {       nom&lt;-alpha[sin,t]*T[sin,sout]*Es(P,O[t+1])[sout]*beta[sout,t+1]       xi[sin,sout]&lt;-nom\/denom     }   }   return(list(xi=xi,g=rowSums(xi))) }  XiGammas&lt;-function(Tm,alpha,beta,S,T,P,O) {   xis&lt;-c()   gs&lt;-c()   for(t in 1:(Tm-1)) {     xg&lt;-xiGamma(t,alpha,beta,S,T,P,O)     xis&lt;-cbind(xis,as.vector(xg$xi))     gs&lt;-cbind(gs,as.vector(xg$g))   }   return(list(xi=xis,g=gs)) } <\/code><\/pre>\n<p>  <img decoding=\"async\" src=\"http:\/\/habr.habrastorage.org\/post_images\/67d\/c25\/dd4\/67dc25dd4c0c20055f44bdc1a9cf7ea4.png\" alt=\"\\xi(s,s&#39;)\"\/> \u0438\u043d\u0442\u0435\u0440\u043f\u0440\u0435\u0442\u0438\u0440\u0443\u0435\u0442\u0441\u044f \u043a\u0430\u043a \u0434\u043e\u043b\u044f \u043f\u0435\u0440\u0435\u0445\u043e\u0434\u043e\u0432 \u0438\u0437 \u0441\u043e\u0441\u0442\u043e\u044f\u043d\u0438\u044f s \u0432 \u0441\u043e\u0441\u0442\u043e\u044f\u043d\u0438\u0435 s\u2019 \u0432 \u043c\u043e\u043c\u0435\u043d\u0442 \u0432\u0440\u0435\u043c\u0435\u043d\u0438 t, \u0430 <img decoding=\"async\" src=\"http:\/\/habr.habrastorage.org\/post_images\/cb1\/2e9\/4b3\/cb12e94b3d6f4dd224a89c1505a15aab.png\" alt=\"\\gamma_t (s)\"\/> \u2013 \u043a\u0430\u043a \u043e\u0431\u0449\u0430\u044f \u0434\u043e\u043b\u044f \u043f\u0435\u0440\u0435\u0445\u043e\u0434\u043e\u0432 \u0438\u0437 \u0441\u043e\u0441\u0442\u043e\u044f\u043d\u0438\u044f s \u0432 \u043c\u043e\u043c\u0435\u043d\u0442 \u0432\u0440\u0435\u043c\u0435\u043d\u0438 t.<\/p>\n<p>  \u041f\u043e\u0434\u043f\u0440\u0430\u0432\u043b\u0435\u043d\u043d\u044b\u0435 \u043f\u0430\u0440\u0430\u043c\u0435\u0442\u0440\u044b \u043c\u043e\u0434\u0435\u043b\u0438 <img decoding=\"async\" src=\"http:\/\/habr.habrastorage.org\/post_images\/757\/251\/016\/757251016f7c7a0e009d3d478f9fb789.png\" alt=\"HMM^{k+1}=\\langle M^{k+1},S^{k+1},I^{k+1},T^{k+1},E^{k+1}\\rangle\"\/> \u0431\u0443\u0434\u0443\u0442 \u0438\u043c\u0435\u0442\u044c \u0441\u043b\u0435\u0434\u0443\u044e\u0449\u0438\u0439 \u0432\u0438\u0434:<\/p>\n<pre><code>BaumWelshFunc&lt;-function(S,I,T,P,O,Tm) {   a&lt;-alphaLog(S,I,T,P,O,Tm)   b&lt;-betaLog(S,I,T,P,O,Tm)      xg&lt;-XiGammas(Tm, a,b, S,T,P,O)   xi&lt;-rowSums(xg$xi)   dim(xi)&lt;-c(length(S),length(S))   g&lt;-rowSums(xi)   new_I &lt;- g\/sum(g)    small_t&lt;-rep(g,length(S))   dim(small_t)&lt;-c(length(S),length(S))   new_T &lt;- xi\/small_t   norm_params&lt;-NormalReestimation(xg$g,O[1:(Tm-1)],S)   new_P&lt;-c()   for(s in S) {     new_P&lt;-c(new_P,Curry(&quot;dnorm&quot;,mean=norm_params$mu[s], sd=sqrt(norm_params$sigma2[s])))   }   return(list(S=S,I=new_I,T=new_T,P=new_P,N=norm_params))   } <\/code><\/pre>\n<p>  \u041f\u0435\u0440\u0435\u043d\u0430\u0441\u0442\u0440\u043e\u0439\u043a\u0430 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 <img decoding=\"async\" src=\"http:\/\/habr.habrastorage.org\/post_images\/1e4\/83c\/c5d\/1e483cc5dd367c0055b9d9c8de914deb.png\" alt=\"f_s(o)\"\/> (NormalReestimation \u0432 \u043d\u0430\u0448\u0435\u043c \u0441\u043b\u0443\u0447\u0430\u0435) \u0437\u0430\u0432\u0438\u0441\u0438\u0442 \u043e\u0442 \u043a\u043e\u043d\u043a\u0440\u0435\u0442\u043d\u043e\u0439 \u0441\u0438\u0442\u0443\u0430\u0446\u0438\u0438. \u041a\u043b\u0430\u0441\u0441\u0438\u0447\u0435\u0441\u043a\u0430\u044f \u0440\u0435\u0430\u043b\u0438\u0437\u0430\u0446\u0438\u044f \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c\u0430 \u043f\u0440\u0435\u0434\u043f\u043e\u043b\u0430\u0433\u0430\u0435\u0442, \u0447\u0442\u043e \u043d\u0430\u0431\u043b\u044e\u0434\u0435\u043d\u0438\u044f \u043e\u043f\u0438\u0441\u044b\u0432\u0430\u044e\u0442\u0441\u044f \u043a\u043e\u043d\u0435\u0447\u043d\u044b\u043c\u0438 \u0434\u0438\u0441\u043a\u0440\u0435\u0442\u043d\u044b\u043c\u0438 \u0440\u0430\u0441\u043f\u0440\u0435\u0434\u0435\u043b\u0435\u043d\u0438\u044f\u043c\u0438, \u0447\u0442\u043e \u043d\u0430\u043c \u043d\u0435 \u043f\u043e\u0434\u0445\u043e\u0434\u0438\u0442. \u041f\u0440\u0435\u0434\u043b\u0430\u0433\u0430\u0435\u0442\u0441\u044f \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u043e\u0432\u0430\u0442\u044c \u0441\u043b\u0435\u0434\u0443\u044e\u0449\u0438\u0439 \u043c\u0435\u0442\u043e\u0434 \u043e\u0446\u0435\u043d\u043a\u0438 \u043f\u0430\u0440\u0430\u043c\u0435\u0442\u0440\u043e\u0432 \u0434\u043b\u044f \u043d\u0430\u0431\u043b\u044e\u0434\u0430\u0435\u043c\u044b\u0445 \u0432\u0435\u043b\u0438\u0447\u0438\u043d:<\/p>\n<pre><code>NormalReestimation&lt;-function(g,O,S) {   denom&lt;-rowSums(g)   mu&lt;-g%*%O\/denom   sigma2&lt;-c()   for(s in S) {     s2&lt;-g[s,]%*%(O-mu[s])^2     sigma2&lt;-c(sigma2,s2)   }   return(list(mu=as.vector(mu),sigma2=sigma2\/denom)) } <\/code><\/pre>\n<p>  \u0420\u0430\u0441\u0441\u043c\u043e\u0442\u0440\u0438\u043c \u0437\u0430\u0434\u0430\u0447\u0443 \u043e \u0434\u0435\u0442\u0435\u043a\u0442\u043e\u0440\u0435 \u043b\u0436\u0438 \u0432 \u0431\u043e\u043b\u0435\u0435 \u043e\u0431\u0449\u0435\u043c \u0432\u0430\u0440\u0438\u0430\u043d\u0442\u0435. \u0423 \u043d\u0430\u0441 \u0435\u0441\u0442\u044c \u043e\u0431\u0449\u0430\u044f \u043c\u043e\u0434\u0435\u043b\u044c \u0442\u043e\u0433\u043e, \u043a\u0430\u043a \u043d\u0430\u0448 \u0434\u0435\u0442\u0435\u043a\u0442\u043e\u0440 \u043b\u0436\u0438 \u0440\u0435\u0430\u0433\u0438\u0440\u0443\u0435\u0442 \u043d\u0430 \u043b\u043e\u0436\u044c \u0438 \u0438\u0441\u043a\u0440\u0435\u043d\u043d\u043e\u0441\u0442\u044c \u0441\u0443\u0431\u044a\u0435\u043a\u0442\u043e\u0432 \u0438 \u0435\u0441\u0442\u044c \u043a\u043e\u043d\u043a\u0440\u0435\u0442\u043d\u044b\u0439 \u0438\u043d\u0434\u0438\u0432\u0438\u0434, \u043a\u043e\u0442\u043e\u0440\u044b\u0439 \u0438\u043d\u043e\u0433\u0434\u0430 \u0432\u0440\u0435\u0442 \u0438\u043d\u043e\u0433\u0434\u0430 \u0433\u043e\u0432\u043e\u0440\u0438\u0442 \u043f\u0440\u0430\u0432\u0434\u0443. \u041c\u044b \u0445\u043e\u0442\u0438\u043c \u043d\u0430\u0441\u0442\u0440\u043e\u0438\u0442\u044c \u043d\u0430\u0448\u0443 \u043c\u043e\u0434\u0435\u043b\u044c \u043d\u0430 \u0440\u0430\u0431\u043e\u0442\u0443 \u0441 \u044d\u0442\u0438\u043c, \u043a\u043e\u043d\u043a\u0440\u0435\u0442\u043d\u044b\u043c \u0438\u043d\u0434\u0438\u0432\u0438\u0434\u043e\u043c.<\/p>\n<p>  \u0414\u0430\u043d\u043e:<\/p>\n<ul>\n<li>\u041c\u043e\u0434\u0435\u043b\u044c. \u0418\u043c\u0435\u044e\u0442\u0441\u044f \u0434\u0432\u0430 \u0441\u043a\u0440\u044b\u0442\u044b\u0445 \u0441\u043e\u0441\u0442\u043e\u044f\u043d\u0438\u044f, \u043a\u043e\u0433\u0434\u0430 \u0441\u0444\u0435\u0440\u0438\u0447\u0435\u0441\u043a\u0438\u0439 \u0447\u0435\u043b\u043e\u0432\u0435\u043a \u0433\u043e\u0432\u043e\u0440\u0438\u0442 \u043f\u0440\u0430\u0432\u0434\u0443 \u0438 \u043a\u043e\u0433\u0434\u0430 \u043b\u0436\u0435\u0442. \u0412 \u0442\u043e\u043c \u0438 \u0432 \u0434\u0440\u0443\u0433\u043e\u043c \u0441\u043b\u0443\u0447\u0430\u0435 \u043f\u043e\u0434\u0451\u0440\u0433\u0438\u0432\u0430\u043d\u0438\u0435 \u0440\u0443\u043a \u043c\u043e\u0436\u043d\u043e \u043f\u0440\u0435\u0434\u0441\u0442\u0430\u0432\u0438\u0442\u044c \u043a\u0430\u043a \u0431\u0435\u043b\u044b\u0439 \u0448\u0443\u043c, \u043d\u043e \u0434\u0438\u0441\u043f\u0435\u0440<\/li>\n<li>\u0412\u044b\u0431\u043e\u0440\u043a\u0430 \u0438\u0437 10000 \u043f\u043e\u0441\u043b\u0435\u0434\u043e\u0432\u0430\u0442\u0435\u043b\u044c\u043d\u044b\u0445 \u043d\u0430\u0431\u043b\u044e\u0434\u0435\u043d\u0438\u0439 \u0437\u0430 \u043d\u0430\u0448\u0438\u043c \u0438\u043d\u0434\u0438\u0432\u0438\u0434\u043e\u043c.<\/li>\n<\/ul>\n<p>  \u041d\u0435\u043e\u0431\u0445\u043e\u0434\u0438\u043c\u043e:<\/p>\n<ol>\n<li> \u041e\u0446\u0435\u043d\u0438\u0442\u044c \u0434\u0438\u0441\u043f\u0435\u0440\u0441\u0438\u044e \u043f\u043e\u0434\u0440\u0430\u0433\u0438\u0432\u0430\u043d\u0438\u044f \u0440\u0443\u043a \u0432 \u0441\u043b\u0443\u0447\u0430\u0435 \u0438\u0441\u043a\u0440\u0435\u043d\u043d\u043e\u0441\u0442\u0438\/\u043b\u0436\u0438 \u044d\u0442\u043e\u0433\u043e \u0438\u043d\u0434\u0438\u0432\u0438\u0434\u0430.<\/li>\n<li> \u041e\u043f\u0440\u0435\u0434\u0435\u043b\u0438\u0442\u044c \u043c\u043e\u043c\u0435\u043d\u0442\u044b \u0441\u043c\u0435\u043d\u044b \u044d\u0442\u0438\u0445 \u0441\u043e\u0441\u0442\u043e\u044f\u043d\u0438\u0439 \u0434\u043b\u044f \u043d\u0435\u0433\u043e.<\/li>\n<\/ol>\n<p>  \u0420\u0435\u0448\u0435\u043d\u0438\u0435:<\/p>\n<ol>\n<li>\u0417\u0430\u0434\u0430\u0434\u0438\u043c \u043d\u0430\u0448\u0443 \u043f\u044f\u0442\u0435\u0440\u043a\u0443:\n<ol>\n<li>\u0427\u0438\u0441\u043b\u043e \u0441\u043e\u0441\u0442\u043e\u044f\u043d\u0438\u0439 <img decoding=\"async\" src=\"http:\/\/habr.habrastorage.org\/post_images\/69f\/d82\/e61\/69fd82e61c861fd32742038b16e7db21.png\" alt=\"M\\leftarrow 2\"\/>;<\/li>\n<li>\u0421\u043e\u0441\u0442\u043e\u044f\u043d\u0438\u044f <img decoding=\"async\" src=\"http:\/\/habr.habrastorage.org\/post_images\/ed8\/ea6\/a94\/ed8ea6a944c0b542bad0960ab16881d9.png\" alt=\"S\\leftarrow \\{0,1\\}\"\/>;<\/li>\n<li>\u041d\u0430\u0447\u0430\u043b\u044c\u043d\u043e\u0435 \u0440\u0430\u0441\u043f\u0440\u0435\u0434\u0435\u043b\u0435\u043d\u0438\u0435 <img decoding=\"async\" src=\"http:\/\/habr.habrastorage.org\/post_images\/d9a\/207\/e51\/d9a207e5196877c249f59cac599b6f98.png\" alt=\"I\\leftarrow(1\/2,1\/2)\"\/>;<\/li>\n<li>\u041f\u043e \u0443\u043c\u043e\u043b\u0447\u0430\u043d\u0438\u044e \u043f\u0440\u0435\u0434\u043f\u043e\u043b\u043e\u0436\u0438\u043c, \u0447\u0442\u043e \u0441\u043e\u0441\u0442\u043e\u044f\u043d\u0438\u044f \u0440\u0430\u0432\u043d\u043e\u043f\u0440\u0430\u0432\u043d\u044b \u0438 \u0441\u043a\u0430\u0436\u0435\u043c \u0441\u043c\u0435\u043d\u0430 \u043f\u0440\u043e\u0438\u0441\u0445\u043e\u0434\u0438\u0442 \u0432 \u0441\u0440\u0435\u0434\u043d\u0435\u043c, \u0440\u0430\u0437 \u0432 999 \u0442\u0430\u043a\u0442\u043e\u0432 \u0432\u0440\u0435\u043c\u0435\u043d\u0438:<br \/>  <img decoding=\"async\" src=\"http:\/\/habr.habrastorage.org\/post_images\/97b\/c3f\/c65\/97bc3fc658d404cf4b2a2a2672e9a812.png\" alt=\"T\\leftarrow\\left(\\begin{matrix}0.999&amp;0.001\\\\0.001&amp;0.999\\end{matrix}\\right)\"\/>;<\/li>\n<li><img decoding=\"async\" src=\"http:\/\/habr.habrastorage.org\/post_images\/9ce\/490\/747\/9ce490747f8d1b5d8f091d0d7c1b8419.png\" alt=\"E=(N(0,1\/2\\cdot\\sigma^2 ),N(0,2 \\cdot\\sigma^2 ))\"\/>, \u0433\u0434\u0435 <img decoding=\"async\" src=\"http:\/\/habr.habrastorage.org\/post_images\/e4d\/efe\/2d4\/e4defe2d4b733f962204a7e16174c4d7.png\" alt=\"\\sigma^2=var(X)\"\/>.<br \/>  \u041f\u043e\u044f\u0441\u043d\u044e, \u0447\u0442\u043e <img decoding=\"async\" src=\"http:\/\/habr.habrastorage.org\/post_images\/c60\/398\/7ea\/c603987eafb6ba3a939c3642f2887bba.png\" alt=\"var(X)\"\/> \u2013 \u044d\u0442\u043e \u0432\u044b\u0431\u043e\u0440\u043e\u0447\u043d\u0430\u044f \u0434\u0438\u0441\u043f\u0435\u0440\u0441\u0438\u044f \u0434\u043b\u044f \u0432\u0441\u0435\u0439 \u0432\u044b\u0431\u043e\u0440\u043a\u0438 \u0446\u0435\u043b\u0438\u043a\u043e\u043c. \u041d\u0430\u0447\u0430\u043b\u044c\u043d\u044b\u0435 \u043e\u0446\u0435\u043d\u043e\u0447\u043d\u044b\u0435 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 \u0434\u043b\u044f <img decoding=\"async\" src=\"http:\/\/habr.habrastorage.org\/post_images\/475\/be5\/4e6\/475be54e6435fd900c9766fbf59b4990.png\" alt=\"E_1, E_2\"\/> \u0432\u044b\u0431\u0438\u0440\u0430\u044e\u0442\u0441\u044f \u0442\u0430\u043a\u0438\u043c \u043e\u0431\u0440\u0430\u0437\u043e\u043c, \u0447\u0442\u043e\u0431\u044b \u043f\u0435\u0440\u0432\u043e\u0435 \u0441\u043e\u0441\u0442\u043e\u044f\u043d\u0438\u0435 \u043e\u0442\u0432\u0435\u0447\u0430\u043b\u043e \u043c\u0435\u043d\u044c\u0448\u0438\u043c \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u044f\u043c \u0434\u0438\u0441\u043f\u0435\u0440\u0441\u0438\u0438, \u0430 \u0432\u0442\u043e\u0440\u043e\u0435 \u2013 \u0431\u043e\u043b\u044c\u0448\u0438\u043c. \u0421\u0443\u0449\u0435\u0441\u0442\u0432\u0435\u043d\u043d\u043e, \u0447\u0442\u043e \u043e\u043d\u0438 \u0440\u0430\u0437\u043d\u044b\u0435.<\/li>\n<\/ol>\n<\/li>\n<li>\u0411\u0443\u0434\u0435\u043c \u043f\u0435\u0440\u0435\u043d\u0430\u0441\u0442\u0440\u0430\u0438\u0432\u0430\u0442\u044c \u043f\u0430\u0440\u0430\u043c\u0435\u0442\u0440\u044b \u0441 \u043f\u043e\u043c\u043e\u0449\u044c\u044e \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c\u0430 \u0411\u0430\u0443\u043c\u0430-\u0412\u0435\u043b\u0448\u0430 \u0434\u043e \u0442\u0435\u0445 \u043f\u043e\u0440, \u043f\u043e\u043a\u0430 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u043e\u0446\u0435\u043d\u043e\u0447\u043d\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 \u043d\u0435 \u0441\u0442\u0430\u0431\u0438\u043b\u0438\u0437\u0438\u0440\u0443\u0435\u0442\u0441\u044f.<\/li>\n<li>\u041d\u0430\u0439\u0434\u0435\u043c \u043d\u0430\u0438\u0431\u043e\u043b\u0435\u0435 \u043f\u0440\u0430\u0432\u0434\u043e\u043f\u043e\u0434\u043e\u0431\u043d\u044b\u0435 \u0441\u043a\u0440\u044b\u0442\u044b\u0435 \u0441\u043e\u0441\u0442\u043e\u044f\u043d\u0438\u044f \u0434\u043b\u044f \u043a\u0430\u0436\u0434\u043e\u0433\u043e \u043c\u043e\u043c\u0435\u043d\u0442\u0430 \u0432\u0440\u0435\u043c\u0435\u043d\u0438 \u0441 \u043f\u043e\u043c\u043e\u0449\u044c\u044e \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c\u0430 \u0412\u0438\u0442\u0435\u0440\u0431\u0438.<\/li>\n<\/ol>\n<p>  \u041e\u0442\u0432\u0435\u0442 \u043f\u0440\u043e\u0438\u043b\u043b\u044e\u0441\u0442\u0440\u0438\u0440\u0443\u0435\u043c \u043d\u0430 \u0434\u0430\u043d\u043d\u044b\u0445 \u0441 \u043f\u0435\u0440\u0432\u043e\u0439 \u043a\u0430\u0440\u0442\u0438\u043d\u043a\u0438: \u0441\u0432\u0435\u0440\u0445\u0443 \u0441\u043e\u0441\u0442\u043e\u044f\u043d\u0438\u044f, \u043a\u0430\u043a \u043e\u043d\u0438 \u0431\u044b\u043b\u0438 \u0437\u0430\u0434\u0430\u043d\u044b \u043f\u0440\u0438 \u043c\u043e\u0434\u0435\u043b\u0438\u0440\u043e\u0432\u0430\u043d\u0438\u0438, \u0441\u043d\u0438\u0437\u0443 \u0441\u043e\u0441\u0442\u043e\u044f\u043d\u0438\u044f, \u043a\u043e\u0442\u043e\u0440\u044b\u0435 \u0431\u044b\u043b\u0438 \u0432\u044b\u0434\u0435\u043b\u0435\u043d\u044b \u0441 \u043f\u043e\u043c\u043e\u0449\u044c\u044e \u043c\u0435\u0442\u043e\u0434\u0430 \u0441\u043a\u0440\u044b\u0442\u044b\u0445 \u0446\u0435\u043f\u0435\u0439 \u041c\u0430\u0440\u043a\u043e\u0432\u0430.<\/p>\n<p>  <img decoding=\"async\" src=\"http:\/\/habrastorage.org\/storage2\/f67\/7e3\/a88\/f677e3a88a7e702bfade2e921691e402.png\"\/><br \/>  \u041d\u0430 \u0433\u0440\u0430\u0444\u0438\u043a\u0435 \u043d\u0438\u0436\u0435 \u043f\u0440\u0438\u0432\u0435\u0434\u0435\u043c \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u044f \u043b\u043e\u0433\u0430\u0440\u0438\u0444\u043c\u0430 \u043e\u0446\u0435\u043d\u043e\u0447\u043d\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 \u0432 \u0437\u0430\u0432\u0438\u0441\u0438\u043c\u043e\u0441\u0442\u0438 \u043e\u0442 \u0447\u0438\u0441\u043b\u0430 \u0438\u0442\u0435\u0440\u0430\u0446\u0438\u0439 \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c\u0430 \u0411\u0430\u0443\u043c\u0430-\u0412\u0435\u043b\u0448\u0430. \u0417\u0435\u043b\u0435\u043d\u044b\u043c \u0446\u0432\u0435\u0442\u043e\u043c \u043e\u0442\u043c\u0435\u0447\u0435\u043d\u0430 \u043e\u0446\u0435\u043d\u043a\u0430, \u0441\u043e\u043e\u0442\u0432\u0435\u0442\u0441\u0442\u0432\u0443\u044e\u0449\u0430\u044f \u043f\u0430\u0440\u0430\u043c\u0435\u0442\u0440\u0430\u043c \u043c\u043e\u0434\u0435\u043b\u0438\u0440\u043e\u0432\u0430\u043d\u0438\u044f.<\/p>\n<p>  <img decoding=\"async\" src=\"http:\/\/habrastorage.org\/storage2\/f99\/654\/96d\/f9965496dd86cb496f94929d4b8bf506.png\"\/><br \/>  \u041e\u0448\u0438\u0431\u043a\u0430 \u0432 \u043e\u0446\u0435\u043d\u043a\u0435 \u0441\u043a\u0440\u044b\u0442\u044b\u0445 \u0441\u043e\u0441\u0442\u043e\u044f\u043d\u0438\u0439 \u043e\u0442\u043d\u043e\u0441\u0438\u0442\u0435\u043b\u044c\u043d\u043e \u043c\u0430\u043b\u0435\u043d\u044c\u043a\u0430\u044f, \u043d\u043e \u0442\u043e\u043b\u044c\u043a\u043e \u043a\u043e\u0433\u0434\u0430 \u043d\u0430\u0447\u0430\u043b\u044c\u043d\u043e\u0435 \u043f\u0440\u0438\u0431\u043b\u0438\u0436\u0435\u043d\u0438\u0435 \u043d\u0435 \u0441\u0438\u043b\u044c\u043d\u043e \u043e\u0442\u043b\u0438\u0447\u0430\u0435\u0442\u0441\u044f \u043e\u0442 \u0442\u043e\u0433\u043e, \u0447\u0442\u043e \u043c\u044b \u0445\u043e\u0442\u0438\u043c \u043d\u0430\u0439\u0442\u0438. \u041d\u0443\u0436\u043d\u043e \u043f\u043e\u043d\u0438\u043c\u0430\u0442\u044c, \u0447\u0442\u043e \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c \u0411\u0430\u0443\u043c\u0430-\u0412\u0435\u043b\u0448\u0430 \u0432\u0441\u0435\u0433\u043e \u043b\u0438\u0448\u044c \u0448\u0430\u0433 \u043f\u043e\u0438\u0441\u043a\u0430 \u043b\u043e\u043a\u0430\u043b\u044c\u043d\u043e\u0433\u043e \u043e\u043f\u0442\u0438\u043c\u0443\u043c\u0430, \u0430 \u0437\u043d\u0430\u0447\u0438\u0442, \u043c\u043e\u0436\u0435\u0442 \u043d\u0438\u043a\u043e\u0433\u0434\u0430 \u043d\u0435 \u0434\u043e\u0441\u0442\u0438\u0433\u043d\u0443\u0442\u044c \u0433\u043b\u043e\u0431\u0430\u043b\u044c\u043d\u043e\u0433\u043e \u044d\u043a\u0441\u0442\u0440\u0435\u043c\u0443\u043c\u0430.<\/p>\n<p>  \u0424\u043e\u0440\u043c\u0443\u043b\u044b \u043f\u0440\u0435\u0432\u0440\u0430\u0449\u0435\u043d\u044b \u0432 \u043a\u0430\u0440\u0442\u0438\u043d\u043a\u0438 \u0441 \u043f\u043e\u043c\u043e\u0449\u044c\u044e <a href=\"http:\/\/mathurl.com\">mathurl.com<\/a>. \u041d\u0430\u0434\u0435\u044e\u0441\u044c, \u043a\u043e\u0433\u0434\u0430-\u043d\u0438\u0431\u0443\u0434\u044c, \u0437\u0434\u0435\u0441\u044c \u043c\u043e\u0436\u043d\u043e \u0431\u0443\u0434\u0435\u0442 \u0432\u0441\u0442\u0430\u0432\u043b\u044f\u0442\u044c \u0444\u043e\u0440\u043c\u0443\u043b\u044b \u0431\u043e\u043b\u0435\u0435 \u043f\u0440\u043e\u0441\u0442\u044b\u043c \u0438 \u043f\u0440\u0438\u044f\u0442\u043d\u044b\u043c \u0433\u043b\u0430\u0437\u0443 \u0441\u043f\u043e\u0441\u043e\u0431\u043e\u043c \u0438, \u043a\u0441\u0442\u0430\u0442\u0438, \u043e\u0447\u0435\u043d\u044c \u043d\u0435 \u0445\u0432\u0430\u0442\u0430\u0435\u0442 \u043f\u043e\u0434\u0441\u0432\u0435\u0442\u043a\u0438 \u0441\u0438\u043d\u0442\u0430\u043a\u0441\u0438\u0441\u0430 \u0434\u043b\u044f <a href=\"http:\/\/r-project.org\">R<\/a>.    \t<\/p>\n<div class=\"clear\"><\/div>\n<\/p><\/div>\n<p> \u0441\u0441\u044b\u043b\u043a\u0430 \u043d\u0430 \u043e\u0440\u0438\u0433\u0438\u043d\u0430\u043b \u0441\u0442\u0430\u0442\u044c\u0438 <a href=\"http:\/\/habrahabr.ru\/post\/188244\/\"> http:\/\/habrahabr.ru\/post\/188244\/<\/a><\/p>\n","protected":false},"excerpt":{"rendered":"<div class=\"content html_format\">   \t\u0421\u043a\u0440\u044b\u0442\u044b\u0435 \u043c\u043e\u0434\u0435\u043b\u0438\/\u0446\u0435\u043f\u0438 \u041c\u0430\u0440\u043a\u043e\u0432\u0430 \u043e\u0434\u043d\u0438 \u0438\u0437 \u043f\u043e\u0434\u0445\u043e\u0434\u043e\u0432 \u043a \u043f\u0440\u0435\u0434\u0441\u0442\u0430\u0432\u043b\u0435\u043d\u0438\u044e \u0434\u0430\u043d\u043d\u044b\u0445. \u041c\u043d\u0435 \u043e\u0447\u0435\u043d\u044c \u043f\u043e\u043d\u0440\u0430\u0432\u0438\u043b\u043e\u0441\u044c \u043a\u0430\u043a \u043e\u0431\u043e\u0431\u0449\u0430\u0435\u0442\u0441\u044f \u043c\u043d\u043e\u0436\u0435\u0441\u0442\u0432\u043e \u0442\u0430\u043a\u0438\u0445 \u043f\u043e\u0434\u0445\u043e\u0434\u043e\u0432 \u0432 <a href=\"http:\/\/habrahabr.ru\/company\/surfingbird\/blog\/177889\/\">\u044d\u0442\u043e\u0439 \u0441\u0442\u0430\u0442\u044c\u0435<\/a>.<\/p>\n<p>  \u0412 \u043f\u0440\u043e\u0434\u043e\u043b\u0436\u0435\u043d\u0438\u0435 \u0436\u0435 \u043c\u043e\u0435\u0439 <a href=\"http:\/\/habrahabr.ru\/post\/180109\/\">\u043f\u0440\u0435\u0434\u044b\u0434\u0443\u0449\u0435\u0439 \u0441\u0442\u0430\u0442\u044c\u0438<\/a> \u043e\u043f\u0438\u0441\u0430\u043d\u0438\u044f \u0441\u043a\u0440\u044b\u0442\u044b\u0445 \u043c\u043e\u0434\u0435\u043b\u0435\u0439 \u041c\u0430\u0440\u043a\u043e\u0432\u0430, \u0437\u0430\u0434\u0430\u043c\u0441\u044f \u0432\u043e\u043f\u0440\u043e\u0441\u043e\u043c: \u043e\u0442\u043a\u0443\u0434\u0430 \u0432\u0437\u044f\u0442\u044c \u0445\u043e\u0440\u043e\u0448\u0443\u044e \u043c\u043e\u0434\u0435\u043b\u044c? \u041e\u0442\u0432\u0435\u0442 \u0434\u043e\u0441\u0442\u0430\u0442\u043e\u0447\u043d\u043e \u0441\u0442\u0430\u043d\u0434\u0430\u0440\u0442\u0435\u043d, \u0432\u0437\u044f\u0442\u044c \u043d\u0435\u043f\u043b\u043e\u0445\u0443\u044e \u043c\u043e\u0434\u0435\u043b\u044c \u0438 \u0441\u0434\u0435\u043b\u0430\u0442\u044c \u0438\u0437 \u043d\u0435\u0435 \u0445\u043e\u0440\u043e\u0448\u0443\u044e.<\/p>\n<p>  \u041d\u0430\u043f\u043e\u043c\u043d\u044e \u043f\u0440\u0438\u043c\u0435\u0440: \u043d\u0430\u043c \u043d\u0443\u0436\u043d\u043e \u0440\u0435\u0430\u043b\u0438\u0437\u043e\u0432\u0430\u0442\u044c \u0434\u0435\u0442\u0435\u043a\u0442\u043e\u0440 \u043b\u0436\u0438, \u043a\u043e\u0442\u043e\u0440\u044b\u0439 \u043f\u043e \u043f\u043e\u0434\u0440\u0430\u0433\u0438\u0432\u0430\u043d\u0438\u044e \u0440\u0443\u043a \u0447\u0435\u043b\u043e\u0432\u0435\u043a\u0430, \u043e\u043f\u0440\u0435\u0434\u0435\u043b\u044f\u0435\u0442, \u0433\u043e\u0432\u043e\u0440\u0438\u0442 \u043e\u043d \u043f\u0440\u0430\u0432\u0434\u0443 \u0438\u043b\u0438 \u043d\u0435\u0442. \u0414\u043e\u043f\u0443\u0441\u0442\u0438\u043c, \u043a\u043e\u0433\u0434\u0430 \u0447\u0435\u043b\u043e\u0432\u0435\u043a \u043b\u0436\u0435\u0442, \u0440\u0443\u043a\u0438 \u0442\u0440\u044f\u0441\u0443\u0442\u0441\u044f \u0447\u0443\u0442\u044c \u0431\u043e\u043b\u044c\u0448\u0435, \u043d\u043e \u043d\u0430\u043c \u043d\u0435 \u0438\u0437\u0432\u0435\u0441\u0442\u043d\u043e \u043d\u0430 \u0441\u043a\u043e\u043b\u044c\u043a\u043e \u0438\u043c\u0435\u043d\u043d\u043e. \u0412\u043e\u0437\u044c\u043c\u0435\u043c \u043c\u043e\u0434\u0435\u043b\u044c \u043d\u0430\u043e\u0431\u0443\u043c, \u043f\u0440\u043e\u0433\u043e\u043d\u0438\u043c \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c \u0412\u0438\u0442\u0435\u0440\u0431\u0438 \u0438\u0437 \u043f\u0440\u0435\u0434\u044b\u0434\u0443\u0449\u0435\u0439 \u0441\u0442\u0430\u0442\u044c\u0438 \u0438 \u043f\u043e\u043b\u0443\u0447\u0438\u043c \u0434\u043e\u0432\u043e\u043b\u044c\u043d\u043e \u0441\u0442\u0440\u0430\u043d\u043d\u044b\u0435 \u0440\u0435\u0437\u0443\u043b\u044c\u0442\u0430\u0442\u044b:<br \/>  <img decoding=\"async\" src=\"http:\/\/habrastorage.org\/storage2\/e4c\/622\/5f1\/e4c6225f162e1efd9387f0191c1e2b86.png\"\/>  <\/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-188244","post","type-post","status-publish","format-standard","hentry"],"_links":{"self":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/posts\/188244","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=188244"}],"version-history":[{"count":0,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/posts\/188244\/revisions"}],"wp:attachment":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=188244"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=188244"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=188244"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}