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...