MSDS 596 Homework 10 Due November 28 2017
Notes. The lowest grade among all eleven homework will be dropped, so NO late submission will be
accepted. All homework assignment must be written on standard 8.5 by 11 paper and stapled together.
Computer generated output without detailed explanations and remarks will not receive any credit. You may
type out your answers, but make sure to use different fonts to distinguish your own words with computer
output. Only hard copies are accepted, except under special circumstances. For the simulation and data
analysis problems, keep the code you develop as you may be asked to present your work later.
1 (80 pts). Use scan("lt.txt") to read in the data from the file lt.txt, which is a vector of length 500.
Consider a local trend model for the data
yt = st + et, et ∼ N(0, 0.25);
st+1 = st + ηt, ηt ∼ N(0, 0.01), s0 ∼ N(0.2, 2.25).
Use your own program for part (a), (b), and (c).
(a) Write your own program to implement the Kalman filter. Calculate the exact log likelihood of the
data. [Hint. Use s1|0 = 0.2, and Σ1|0 = 2.26 as the initial values.]
(b) Plot the predicted state variables st|t−1 for 1 ≤ t ≤ T. Also plot the 95% confidence intervals
st|t−1 ± 2pΣt|t−1. Your plot should look similar to Figure 11.4 of the textbook.
(c) Plot the filtered state variables st|
for 1 ≤ t ≤ T. Also plot the 95% confidence intervals st|t±2pΣt|t.
(d) Pretend that you do not know σ2e and σ2η, but you do have the prior knowledge that s0 ∼ N(0.2, 2.25),
use the dlm package to find the MLE of σ2e and σ2η.
(e) Reproduce the plots of part (b) and (c), using the MLE you obtained from part (d). Also, plot the
smoothed state variable st|T , together with the 95% confidence intervals st|T ± 2
p
Σt|T , for every t.
(f) Repeat part (d), using the R function StructTS() with type="level".
2 (20 pts). Rewrite the model (1−.3B)(1−.5B)(1−.6B +.4B2
)(rt −1.3) = at in the original form without
using the back-shift operator B.
1

因为专业,所以值得信赖。如有需要,请加QQ:99515681 或 微信:codehelp

MSDS 596 Homework的更多相关文章

  1. CSCI 1100 — Computer Science 1 Homework

    CSCI 1100 — Computer Science 1 Homework 8CS1 Multiverse: ClassesOverviewThis homework is worth 100 p ...

  2. bzoj 4320: ShangHai2006 Homework

    4320: ShangHai2006 Homework Time Limit: 10 Sec Memory Limit: 128 MB Description 1:在人物集合 S 中加入一个新的程序员 ...

  3. HDU 1789 Doing Homework again(贪心)

    Doing Homework again 这只是一道简单的贪心,但想不到的话,真的好难,我就想不到,最后还是看的题解 [题目链接]Doing Homework again [题目类型]贪心 & ...

  4. NYOJ题目596谁是最好的Coder

    aaarticlea/png;base64,iVBORw0KGgoAAAANSUhEUgAAAscAAAMaCAIAAADlQ3w8AAAgAElEQVR4nO3dO3LbvN4H4G8T7r2Q1F

  5. hdu-1789-Doing Homework again

    /* Doing Homework again Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Oth ...

  6. nyist 596 谁是最好的Coder

    http://acm.nyist.net/JudgeOnline/problem.php?pid=596 谁是最好的Coder 时间限制:1000 ms  |  内存限制:65535 KB 难度:0 ...

  7. HDU 1789 Doing Homework again (贪心)

    Doing Homework again http://acm.hdu.edu.cn/showproblem.php?pid=1789 Problem Description Ignatius has ...

  8. Doing Homework 状态压缩DP

    Doing Homework 题目抽象:给出n个task的name,deadline,need.  每个任务的罚时penalty=finish-deadline;   task不可以同时做.问按怎样的 ...

  9. 机器学习 —— 概率图模型(Homework: Exact Inference)

    在前三周的作业中,我构造了概率图模型并调用第三方的求解器对器进行了求解,最终获得了每个随机变量的分布(有向图),最大后验分布(双向图).本周作业的主要内容就是自行编写概率图模型的求解器.实际上,从根本 ...

随机推荐

  1. 简单使用Moq框架

    Moq框架简单使用   系列目录 Moq库简介及安装 Moq简介 Moq是.net平台下的一个非常流行的模拟库,只要有一个接口它就可以动态生成一个对象,底层使用的是Castle的动态代理功能. 它的流 ...

  2. 使用Docker快速搭建Zookeeper和kafka集群

    使用Docker快速搭建Zookeeper和kafka集群 镜像选择 Zookeeper和Kafka集群分别运行在不同的容器中zookeeper官方镜像,版本3.4kafka采用wurstmeiste ...

  3. layui switch 确定之后才变更状态

    let x = data.elem.checked; data.elem.checked = !x; form.render(); 完整代码 form.on('switch(is_enable)', ...

  4. C++ 智能指针 shared_ptr 分析

    引文: C++对指针的管理提供了两种解决问题的思路: 1.不允许多个对象管理一个指针 2.允许多个对象管理一个指针,但仅当管理这个指针的最后一个对象析构时才调用delete ps:这两种思路的共同点就 ...

  5. 用LabVIEW做声源定位系统

    前一阵子,研发部举办了为期三天的第一届Innovation Day,让大家用3天时间去完成工作之外的一些创意.有人做微信小程序,有人继续研究一些工作中用到的Tool,有人把一直想解决而没时间解决的老b ...

  6. JS Web API 拖拽对话框案例

    <style> .login-header { width: 100%; text-align: right; height: 30px; font-size: 24px; line-he ...

  7. 全栈项目|小书架|微信小程序-点赞功能实现

    微信小程序端的点赞功能其实没什么好介绍的,无非就是调用接口改变点赞状态和点赞数量.需要注意的是取消点赞时的处理,我这里为了减少服务器接口的调用,直接本地存一个变量,修改这里的变量值即可. 由于源码都相 ...

  8. C# vb .net实现黑白效果特效滤镜

    在.net中,如何简单快捷地实现Photoshop滤镜组中的黑白效果效果呢?答案是调用SharpImage!专业图像特效滤镜和合成类库.下面开始演示关键代码,您也可以在文末下载全部源码: 设置授权 第 ...

  9. 关于.Net使用企业库访问MySql数据库

    关于.Net使用企业库访问MySql数据库 在网上看了很多又重写又加WebConfig中的内容,其实不用那么麻烦 企业库5.0访问MySql数据库只需要在Web服务器安装mysql-connector ...

  10. Mybatis中使用association及collection进行自关联示例(含XML版与注解版)

    XML版本: 实体类: @Data @ToString @NoArgsConstructor public class Dept { private Integer id; private Strin ...