R-beta: spline problems(?)

OK, thanks, it looks like I can use the splines library instead (I had no
idea how, but now I do!). Should we call it a bug then?

Question No. 2:

From the same S-plus program, I was forced to change a line that was giving

me Segmentation Error problems (a crash on Solaris 2.5):

 error _ t(y.fit) %*% ( DBD + Lambda ) %*% y.fit +
 rbind(rep(1,n)) %*% F.vec.sq - 2 * t(y.fit) %*% F.vec

Apparently, leaving out the middle term did the trick:
 error _ t(y.fit) %*% ( DBD + Lambda ) %*% y.fit - 2 * t(y.fit) %*% F.vec
, while the program still runs more or less as expected. There's nothing
nasty about this middle term, is there?

Finally, a certain READ.MAT function, that reads in a single-column file,
taking the first two entries as the dimension of a matrix and storing the
remaining values in such a matrix, wouldn't work; anybody can see why?

READ.MAT _function(file, nrow = 0) {
if(nrow > 0) {
    n <- scan(file, n = 2)[2] * nrow
    tmp <- scan(file, n = n + 2)
  } 

## first point of action: simply read in the file in vector 'tmp'
else {
    tmp <- scan(file)
  }

if(length(tmp)<=2) {
    return(NULL) 
  }

## the following clause should do the vector>>matrix conversion

 else {
    x <- matrix(tmp[ - c(1, 2)], ncol = tmp[2], byrow = T)
    return(x)
  }
}

Finally, if anybody is curious to see what is giving me all this trouble,
it is some code for what is termed Mixture Spline models:
http://cm.bell-labs.com/cm/ms/departments/sia/dxsun/speech/mixSplineSpec