predict_example_lrml.R

predict_from_theta <- function(x) {
  
  x <- cbind(x, rep(1, length(x)))
  dot(x, theta)
  
}

predict_from_theta(rnorm(10))
           [,1]
 [1,] -1.530065
 [2,]  8.036859
 [3,]  6.895680
 [4,]  3.170026
 [5,]  4.905467
 [6,]  2.644702
 [7,] 12.555390
 [8,]  1.172425
 [9,]  3.776509
[10,]  4.849211