objF {SPPr} | R Documentation |
Objective function (Least squares)
objF(p)
p |
p : fitting parameters: epsilon real, imaginary, thickness, scaling factor. |
Sum of squares (to be minimized by optim)
baptiste
## Not run: parameters <- list(layer2fit=1,epsr=c(-12),epsi=c(1.2),thick=c(40), cel=3e8,lambda=632.8e-9,nPrism=1.46,anglePrism=pi/3,nWater=1, zeroShift=-0.4,normalisation=23,pIni=c(-12,1.2,40,1)) attach(parameters,2,warn.conflicts = FALSE) data(scanSPP) head(scanSPP) dataNorm<-data2norm(scanSPP,parameters) #save(file="dataNorm.rda",dataNorm) data(dataNorm) plotData(dataNorm=dataNorm,int=TRUE,ylim=c(0,max(1,max(dataNorm$N)))) plotModel(p=pIni,dataNorm=dataNorm,parameters= parameters,int=TRUE,col=2) result <-optim(pIni, objF) yFin<-reflModel(result$par,dataNorm,parameters) yErr<-abs(yFin-dataNorm$N) plotModel(p=result$par,dataNorm=dataNorm,parameters= parameters,int=TRUE,col=3) ## End(Not run)