R/maxStep.r

Defines functions maxStep

Documented in maxStep

maxStep <- function(beta, beta_new, V, eps = 0){

# function t on page 6
prod <- t(V) %*% beta    
prod_new <- t(V) %*% beta_new    
    
ratio <- - prod / (prod_new - prod)
ratio <- ratio[prod_new > eps]

t <- min(1, ratio)
    
return(t)}

Try the OrdFacReg package in your browser

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

OrdFacReg documentation built on May 1, 2019, 10:06 p.m.