search for: f_vib

Displaying 1 result from an estimated 1 matches for "f_vib".

2007 Oct 28
0
Request for help with nls error
...tfile,"\n",append=TRUE) n_vib <- length(time_vib) p_vib <- 3 cmax <- 30 WLSiter <- 500 beta0_vib <- list(Wmax_star= 30, tau_star=15, gamma_star=.8) t0 <- 20 tol <- 10^-8 meanfunc_vib <- function(time_vib,Wmax_star,tau_star,gamma_star){ diff<-time_vib-t0 f_vib<-Wmax_star*(1-exp(-log(2)*((diff/tau_star)^gamma_star))) f_vib # meanfunc_vib_deriv <- function(time_vib,Wmax_star,tau_star,gamma_star) meangrad_vib <- array(0,c(n_vib,p_vib)) meangrad_vib[,1] <- 1-exp(-log(2)*(diff/tau_star)^gamma_star) meangrad_vib[,2] <- -(Wmax_star*log...