search for: mudot_vib

...star))/tau_star) weightgrad_vib[,"gamma_star"] <- w12_vib*(Wmax_star*log(2)*((diff/tau_star)^gamma_star)*log(diff/tau_star)*exp(-log(2)*(diff/tau_star)^gamma_star)) attr(weightf_vib,"gradient") <- weightgrad_vib weightf_vib } trkfunc_vib <- function(resid_vib,mudot_vib,mu_vib,theta_vib){ trk_vib <- resid_vib*((mudot_vib/mu_vib)**theta_vib) trkgrad_vib <- array(0,c(length(mu_vib),1),list(,NULL,c("theta_vib"))) trkgrad_vib[,"theta_vib"] <- trk_vib*log(mudot_vib/mu_vib) attr(trk_vib,"gradient") <- trkgrad_vib #...