Displaying 1 result from an estimated 1 matches for "bgls_vib".
2007 Oct 28
0
Request for help with nls error
...le=outfile,"\n",
append=F)
cat("OLS estimate - VIB ",round(bols_vib,6),file=outfile,"\n","\n",append=T)
cat("Estimate of assumed constant SD of y
",round(sigols_vib,6),file=outfile,
"\n","\n","\n",append=T)
bgls_vib <- bols_vib
diss_vib_update2<-diss_vib-30.57897
vib_dat<-data.frame(time_vib,diss_vib_update2)
for (k in 1:cmax){
mu_vib <- unweightfunc_vib(time_vib,bgls_vib[1],bgls_vib[2],bgls_vib[3])
mudot_vib <- prod(mu_vib)**(1/n_vib)
mudot_vib <-rep(mudot_vib,n_vib)
resid_vib <- diss_...