从ns对象中提取结

时间:2012-02-16 20:42:13

标签: r regression splines

library(splines)
x <- runif(100000)
spline <- ns(x, df=5)

如何判断此脊椎对象的结点?

> spline ['knots']
[1] NA

str让我接近我需要的东西,但我希望能够提取样条曲线的矢量,格式化它,并将它打印在Sweave文档中。

> str(sp)
 ns [1:117542, 1:5] 0.527 0.474 0.455 0.472 0.498 ...
 - attr(*, "dimnames")=List of 2
  ..$ : NULL
  ..$ : chr [1:5] "1" "2" "3" "4" ...
 - attr(*, "degree")= num 3
 - attr(*, "knots")= Named num [1:4] 1.03 1.55 1.99 2.7
  ..- attr(*, "names")= chr [1:4] "20%" "40%" "60%" "80%"
 - attr(*, "Boundary.knots")= num [1:2] 0.0214 4.9999
 - attr(*, "intercept")= logi FALSE
 - attr(*, "class")= chr [1:3] "ns" "basis" "matrix"

1 个答案:

答案 0 :(得分:6)

嗯,这样的东西肯定可以提取这些东西:

> attr(spline, "knots")
      20%       40%       60%       80% 
0.1966951 0.3973261 0.5982386 0.7993613