estimate_shifts: Estimate shifts

In abichat/zazou: Z-Scores As Ornstein-Uhlenbeck

Description

Estimate shifts

Usage

estimate_shifts(
  zscores,
  tree,
  alphaOU,
  lambda = NULL,
  method = c("lasso", "scaled lasso", "L-BFGS-B"),
  criterion = c("bic", "pbic"),
  beta0 = rep(0, length(tree$edge.length)),
  constraint_type = c("beta", "yhat", "none"),
  ...
)

Arguments

zscores	z-scores (size m)
tree	the tree used to compute the incidence and covariance matrices (if not specified)
alphaOU	a grid of positive alpha parameters
lambda	a grid of positive regularization parameters
method	method to use for the optimization. One of lasso, scaled lasso, or L-BFGS-B.
criterion	criterion on which the selection is done.
beta0	initial position (size n+m)
constraint_type	Constrains on shifts.
...	further arguments passed to or from other methods.

Value

an object of class shiftpunct

abichat/zazou documentation built on Sept. 8, 2021, 6:53 a.m.