[R] problem with optim: (list) object cannot be coerced to 'double'

Uwe Ligges ligges at statistik.uni-dortmund.de
Tue Mar 14 16:17:58 CET 2006


Zebouni, Stephane (Exchange) wrote:

> Hi,
> 
>  
> 
> I am trying to use optim to solve a heavy calibration problem. I supply
> the parameters in vector form. But before entering my target 
> 
>  
> 
> The call is simply:
> 
>  
> 
> optim(par = parameters, fn = SumLSQ, method = "Nelder-Mead")
> 
>  
> 
> the function SumLSQ is simply:
> 
>  
> 
> SumLSQ<-function(parameters, data = timeseries){
> 
>  
> 
> print("sumLSQ")
> 
>       
> 
>  
> 
>       nbseries = dim(timeseries)[2]/3
> 
>  
> 
>       SumLSQ = 0
> 
>       
> 
>       for (i in (1:nbseries)){
> 
>             SumLSQ = SumLSQ +
> LSQ(parameters,timeseries[,((i-1)*3+1):(i*3)])
> 
>       }
> 
>  
> 
> }
> 
>  
> 
> I actually never enter the objective function SumLSQ. I always receive
> the error:
> 
>  
> 
> (list) object cannot be coerced to 'double'
> 
>  
> 
> And I really don't know where it comes from. II thought it was the
> format of "par" argument but I think it is correct to supply a numerical
> vector  ...

Yes, it is. Please type
   str(parameters)
an tell us the result.

Uwe Ligges



> 
>  
> 
> Thnaks a lot for any help!
> 
>  
> 
> Stephane
> 
>  
> 
>  
> 
> 
> 
> ------------------------------------------------------------------------
> 
> 
> 
> ***********************************************************************
> Bear Stearns is not responsible for any recommendation, solicitation, 
> offer or agreement or any information about any transaction, customer 
> account or account activity contained in this communication.
> ***********************************************************************
> 
> 
> 
> ------------------------------------------------------------------------
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html




More information about the R-help mailing list