Logo Questions Linux Laravel Mysql Ubuntu Git Menu
 

Implementation of Particle Swarm Optimization Algorithm in R

I'm checking a simple moving average crossing strategy in R. Instead of running a huge simulation over the 2 dimenional parameter space (length of short term moving average, length of long term moving average), I'd like to implement the Particle Swarm Optimization algorithm to find the optimal parameter values. I've been browsing through the web and was reading that this algorithm was very effective. Moreover, the way the algorithm works fascinates me...

Does anybody of you guys have experience with implementing this algorithm in R? Are there useful packages that can be used?

Thanks a lot for your comments.

Martin

like image 404
Eva Avatar asked Oct 17 '10 08:10

Eva


People also ask

How does particle swarm optimization algorithm work?

Particle swarm is a population-based algorithm. In this respect it is similar to the genetic algorithm. A collection of individuals called particles move in steps throughout a region. At each step, the algorithm evaluates the objective function at each particle.

What is the initialization step of the particle swarm optimization method?

The first step of the PSO algorithm is to initialize the swarm and control parameters. In the context of the basic PSO, the acceleration constants, c1 and c2, the initial veloc- ities, particle positions and personal best positions need to be specified.


1 Answers

Well, there is a package available on CRAN called pso, and indeed it is a particle swarm optimizer (PSO).

I recommend this package.

It is under actively development (last update 22 Sep 2010) and is consistent with the reference implementation for PSO. In addition, the package includes functions for diagnostics and plotting results.

It certainly appears to be a sophisticated package yet the main function interface (the function psoptim) is straightforward--just pass in a few parameters that describe your problem domain, and a cost function.

More precisely, the key arguments to pass in when you call psoptim:

  • dimensions of the problem, as a vector (par);

  • lower and upper bounds for each variable (lower, upper); and

  • a cost function (fn)

There are other parameters in the psoptim method signature; those are generally related to convergence criteria and the like).

Are there any other PSO implementations in R?

There is an R Package called ppso for (parallel PSO). It is available on R-Forge. I do not know anything about this package; i have downloaded it and skimmed the documentation, but that's it.

Beyond those two, none that i am aware of. About three months ago, I looked for R implementations of the more popular meta-heuristics. This is the only pso implementation i am aware of. The R bindings to the Gnu Scientific Library GSL) has a simulated annealing algorithm, but none of the biologically inspired meta-heuristics.

The other place to look is of course the CRAN Task View for Optimization. I did not find another PSO implementation other than what i've recited here, though there are quite a few packages listed there and most of them i did not check other than looking at the name and one-sentence summary.

like image 123
doug Avatar answered Oct 16 '22 13:10

doug