Difference between revisions of "Bsplines.R"
From Organic Design wiki
(R code for plotting B-splines) |
(version using 'bs' function) |
||
| Line 15: | Line 15: | ||
X | X | ||
# p41 | # p41 | ||
| − | knots <- c( | + | par(mfrow=c(2,1)) |
| + | knots <- c(1:3,5,7,8,10) | ||
atx <- seq(0:11, by=0.1) | atx <- seq(0:11, by=0.1) | ||
for(ord in 1:4) { | for(ord in 1:4) { | ||
| − | B <- | + | B <- bs(x=atx, degree=ord, knots=knots, intercept=TRUE, Boundary.knots=c(0,11)) |
matplot(atx, B[,1], type="l", ylim=0:1, lty=2) | matplot(atx, B[,1], type="l", ylim=0:1, lty=2) | ||
matlines(atx, B[,-1], col=1, lty=1) | matlines(atx, B[,-1], col=1, lty=1) | ||
Revision as of 02:57, 5 September 2007
Code snipits and programs written in R, S or S-PLUS
- p38
pos <- function(x) ifelse(x>0, x, 0) x <- 1:7 y <- c(8,3,8,5,9,14,11) knot <- 4 plot(x,y) X <- cbind(1,x,x^2, x^3, pos(x-knot)^3) fit <- lm(y~X-1) xx <- seq(1,7, length=200) XX <- cbind(1,xx, xx^2, xx^3, pos(xx-knot)^3) lines(xx, predict(fit, as.data.frame(XX))) # Bug here abline(v=knot, lty=2) X
- p41
par(mfrow=c(2,1)) knots <- c(1:3,5,7,8,10) atx <- seq(0:11, by=0.1) for(ord in 1:4) { B <- bs(x=atx, degree=ord, knots=knots, intercept=TRUE, Boundary.knots=c(0,11)) matplot(atx, B[,1], type="l", ylim=0:1, lty=2) matlines(atx, B[,-1], col=1, lty=1) title(paste("B splines of order", ord)) abline(v=knots, lty=2) }
