R bs()函数计算(手算模拟)

原创

于 2019-11-05 22:12:14 发布 · 3.9k 阅读

16 ·

CC 4.0 BY-SA版权

文章标签：

#R B样条

本文探讨R语言中bs()函数的使用，重点在于理解其计算B样条基矩阵的过程。通过手动计算模拟bs()函数，揭示了函数内部节点的设置，并通过实例验证了bs()函数的正确性。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

此处需要注意的是R中的默认degree=3，此处的1：5是并不代表内部节点，而是在最后的输出结果的矩阵里，B样条基函数会返回它们在1，2，3，4，5处的值,所以返回的矩阵才会是5行，下面有一个简单实现bs()计算基矩阵的代码，我忘了是从哪里找到的了，可以帮助理解R中bs()的计算过程

rm(list=ls())
# x <- c(0.0, 11.0, 17.9, 49.3, 77.4)
#这是一个递归的程序，注意一下
basis <- function(x, degree, i, knots) {
  browser()
  if(degree == 0){
    B <- ifelse((x >= knots[i]) & (x < knots[i+1]), 1, 0)
  } else {
    if((knots[degree+i] - knots[i]) == 0) {
      alpha1 <- 0
    } else {
      alpha1 <- (x - knots[i])/(knots[degree+i] - knots[i])
    }
    if((knots[i+degree+1] - knots[i+1]) == 0) {
      alpha2 <- 0
    } else {
      alpha2 <- (knots[i+degree+1] - x)/(knots[i+degree+1] - knots[i+1])
    }
    B <-