Displaying 1 result from an estimated 1 matches for "resid_vib".
2007 Oct 28
0
Request for help with nls error
...ar)^gamma_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") <- trkgr...