R/generate_gaussian_quad_points_surv.R

In JiyuanHu/JointMM: A joint modeling tool for analyzing zero-inflated longitudinal proportions and time-to-event data

generate_gaussian_quad_points_surv <-
function(quad.n){
  gherm <- gauss.quad(quad.n, kind = "hermite")
  nodes = gherm$nodes
  weights = gherm$weights
  weights = as.numeric(outer(weights,weights))
  nodes = cbind(rep(nodes,quad.n),rep(nodes,each = quad.n))
  bivar.gherm = list(nodes = nodes, weights = weights)
  return(bivar.gherm)
}