# Only one independent series
data(seaLevelTempSOI)
model.dlm = dlm(x = seaLevelTempSOI$LandOcean,
y = seaLevelTempSOI$GMSL , q = 5)
summary(model.dlm)
residuals(model.dlm)
coef(model.dlm)
fitted(model.dlm)
removed = list(x = c(1))
model.dlm = dlm(x = seaLevelTempSOI$LandOcean,
y = seaLevelTempSOI$GMSL , q = 5,
remove = removed)
summary(model.dlm)
# Remove intercept as well
removed = list(x = c(1,-1))
model.dlm = dlm(x = seaLevelTempSOI$LandOcean,
y = seaLevelTempSOI$GMSL , q = 5,
remove = removed)
summary(model.dlm)
removed = list(x = c(0,1))
model.dlm = dlm(x = seaLevelTempSOI$LandOcean,
y = seaLevelTempSOI$GMSL , q = 5,
remove = removed)
summary(model.dlm)
model.dlm = dlm(formula = GMSL ~ LandOcean ,
data = seaLevelTempSOI , q = 5)
summary(model.dlm)
removed = list(x = c(0,1))
model.dlm = dlm(formula = GMSL ~ LandOcean ,
data = seaLevelTempSOI , q = 5,
remove = removed)
summary(model.dlm)
# Remove intercept as well
removed = list(x = c(0,-1))
model.dlm = dlm(formula = GMSL ~ LandOcean ,
data = seaLevelTempSOI , q = 2,
remove = removed)
summary(model.dlm)
# Multiple independent series
data(M1Germany)
data = M1Germany[1:144,]
model.dlm = dlm(formula = logprice ~ interest + logm1,
data = data.frame(data) , q = 4)
summary(model.dlm)
removed = list(interest = c(1,3), logm1 = c(2))
removed
model.dlm = dlm(formula = logprice ~ interest + logm1,
data = data.frame(data) , q = 4 , remove = removed)
summary(model.dlm)
removed = list(interest = c(0,1,3), logm1 = c(0,2))
removed
model.dlm = dlm(formula = logprice ~ interest + logm1,
data = data.frame(data) , q = 4 , remove = removed)
summary(model.dlm)
removed = list( logm1 = c(1,2))
removed
model.dlm = dlm(formula = logprice ~ interest + logm1,
data = data.frame(data) , q = 4 , remove = removed)
summary(model.dlm)
Run the code above in your browser using DataLab