r语言 - 确保样条回归中的结被放置在正确的位置

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为小时时间序列对象。现在，我想对它进行样条回归。我想在7a处打结。上午和下午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))

如您所见，这里确实没有必要使用"xts"对象。

---

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"

"knots"字段给出了内部结的位置信息。