确保样条回归中的节点位于正确的位置

Question

我有时间序列数据：

library(xts)
library(splines)
set.seed(123)
time <- seq(as.POSIXct("2015-09-01"),as.POSIXct("2015-09-01 23:59:59"),by="hour")
ob <- xts(rnorm(length(time),150,5),time))

对象ob是每小时的时间序列对象。现在，我想对它进行样条回归。我想在上午7点和下午4点打结。 R中的以下语句是否确保此

ns(ob,knots = c(7,16)) # 7 corresponds to 7 AM and 16 corresponds to 4 PM

另外，我应该如何检查在上述时间放置结？

Answer 1

你有点走错了路。您似乎希望按时回归观察，因此您应该真正提供时间索引而不是观察ob到ns。

y <- as.vector(ob)    ## observations
x <- 1:24    ## 24 hourse

然后考虑一个模型：

y ~ ns(x, knots = c(7, 16))

正如您所看到的，确实没有必要使用＆＃34; xts＆＃34;对象在这里。

ns生成设计矩阵。检查

X <- ns(x, knots = c(7, 16))

您将看到属性：

#attr(,"degree")
#[1] 3
#attr(,"knots")
#[1]  7 16
#attr(,"Boundary.knots")
#[1]  1 24
#attr(,"intercept")
#[1] FALSE
#attr(,"class")
#[1] "ns"     "basis"  "matrix"

＆＃34;结＆＃34;字段为您提供内部结的位置信息。