search for: y_est

...5]-(x/number[2])^number[7])^2)+ sum((y- number[6]-(x/number[3])^number[7])^2)+ for minimising this sum, I compute the function "sce": sce<-function(param){ sce_yest<-matrix(nrow=3,ncol=1) for( i in 1:3){ yy<-(y[((long[i]+1):long[i+1])]) xx<-x[(long[i]+1):(long[i+1])] y_est<-(param[i+2]-(xx/param[i])^param[(2*3)+1]) sce_yest[i,]<-sum((yy-y_est)^2) } return(sum(sce_yest)) } Then, I use the fonction optim for obtaining a vector of 7 parameters which will minimise the fonction sce. I use initial parameters at random that I don't have a answer relating to a...

...1):((long[i]+1)+(length(vv)=length(yy))-1))] fl<-optim(p=pinit,sce,method="L-BFGS-B",lower=pinf,upper=psup, control=list(maxit=200000,temp=200)) co<-fl$par coe<-fl$value coeff<-c(co,coe) } coeff #graphical verification plot(x_b,y_b) y_est<-coeff[6]-(x_b/coeff[2])^coeff[7] points(x_b,y_est) -------------- next part -------------- An embedded message was scrubbed... From: Kervahu Anne <anne.kervahu at adria.tm.fr> Subject: syntax problem Date: Fri, 8 Apr 2005 14:31:04 +0200 Size: 2442 Url: https://stat.ethz.ch/piperm...