R/gen.Data.lmFormula.R

`gen.Data.lmFormula` <-
function(fmla, pars, var.Frame, data.Sigma=1, data.Size, train.Frame, test.Size, ...){
  if(dim(train.Frame)[1] > dim(var.Frame)[1]) stop ("Screwed up training frame / var.Frame")
  Terms <- delete.response(terms(fmla))
  m <- model.frame(Terms, var.Frame)
  X <- model.matrix(Terms, m)
  temp <- X %*% pars
  temp <- temp + rnorm(length(temp), sd = data.Sigma)
  train.Size <- dim(train.Frame)[1]
  train.Inds = 1:train.Size
  test.Inds = (train.Size + 1):dim(var.Frame)[1]
  if(dim(train.Frame)[1] == dim(var.Frame)[1]) test.Inds = NULL
  list(data.Vector = temp[train.Inds], test.Vector = temp[test.Inds])
}

Try the CombMSC package in your browser

Any scripts or data that you put into this service are public.

CombMSC documentation built on May 2, 2019, 2:32 p.m.