prior.b {emulator} | R Documentation |
Gives the fitted regression coefficients corresponding to the specified regression model.
prior.b(H, Ainv, d, b0 = NULL, B0 = NULL) prior.B(H , Ainv , B0=NULL)
H |
Regression basis function (eg that returned by regressor.multi() ) |
Ainv |
inv(A) where A is a correlation matrix (eg that
returned by corr.matrix() ) |
d |
Vector of data points |
b0 |
prior constant |
B0 |
prior coefficients |
Robin K. S. Hankin
# example has 10 observations on 6 dimensions. # function is just sum( (1:6)*x) where x=c(x_1, ... , x_2) data(toy) val <- toy d <- apply(val,1,function(x){sum((1:6)*x)}) #add some noise: d <- jitter(d) A <- corr.matrix(val,scales=rep(1,ncol(val)), power=1.5) Ainv <- solve(A) H <- regressor.multi(val) prior.b(H,Ainv,d) prior.B(H,Ainv)