Multi-Objective Optimisation using the Bees Algorithm

authors: D.T. Pham, Afshin Ghanbarzadeh

This paper describes the first application of the Bees Algorithm to multi-objective optimisation problems. The Bees Algorithm is a search procedure inspired by the way honey bees forage for food. A standard mechanical design problem, the design of a welded beam structure, was used to benchmark the Bees Algorithm. The paper presents the results obtained showing the robust performance of the Bees Algorithm.

AttachmentSize
Multi-Objective.wmv5.39 MB

a pdf file
Eldukhri's picture
Submitted by Eldukhri on Wed, 04/07/2007 - 9:14am.

Dear Author,

Thank you for submitting your interesting paper on Bees Algorithm. It would be even more stimulating to the discussion if you could upload your presentation. We look forward to hearing you explaining your work and its main outcome.


Vlado's picture
Submitted by Vlado on Wed, 04/07/2007 - 10:19am.

How do you determine the number of patches required to discover a useful Pareto optimal set?

Do you discover the Pareto front concurrently or require a restart of the algorithm to find another solution?


Afshin's picture
Submitted by Afshin on Thu, 05/07/2007 - 9:05am.

Hi

Thanks for your interest, the presentation is there now. Looking forward to see your comments.

Afshin


Afshin's picture
Submitted by Afshin on Thu, 05/07/2007 - 9:17am.

Dear Vlado,

Hi and thanks for your nice question.

The number of selected sites is chosen empirically, but it does not have major influence on the results as we repeat the algorithm enough number of iterations. In the multi-objective Bees Algorithm, the Pareto optimal set will be constructed and updated in each iteration gradually, so the number of selected sites or patches is not a critical parameter.

The Pareto optimal set will be created during the running of algorithm and in the end of each run of the program, the Pareto optimal set is constructed (it is not needed to run the program many times to find non-dominated solutions one by one).

Hope this answer clear more the way of working of the multi-objective Bees Algorithm. Please do not hesitate to ask any other possible question.

Many thanks,

Afshin


Vlado's picture
Submitted by Vlado on Thu, 05/07/2007 - 12:32pm.

I'm not sure I understand your reply.

Can you be more precise in what you mean by iterations in each case. Do you mean iterations/generations of one run of the algorithm, or do you mean multiple iterative runs.

It makes a difference, since the former means that the number of non-dominated points is lower or equal to the number of patches, which in terms means that the number of patches is actually important.

Or alternatively do you mean that you run the program a number of times and accumulate somewhere the results.


Sameh's picture
Submitted by Sameh on Thu, 05/07/2007 - 1:26pm.

Thank you Afshin for your powerpint presentation

As I can notice from your paper, you managed to decrease the number of parameters need to be tuned. The question is how have you done this?
and is it possible to apply it to all other applications or it is specific for this problem?

Thanks

Sameh Otri
PhD Research Student


Afshin's picture
Submitted by Afshin on Thu, 05/07/2007 - 1:31pm.

Hi

Let me explain it in this way.

In each iteration/generation, the non-dominated solutions are candidate to be added to the Pareto optimal solution. Each of these new non-dominated solution which is not dominated by old member of the Pareto optimal set, will be added to the Pareto optimal set. Also, the old members of the Pareto optimal will be checked, if any of them is dominated be the new member, the dominated ones will be removed from the set.

As program is running, the Pareto optima set will grow ( add and remove from it) and in the end of the run the constructed set will be the answer and this way the Pareto optimal set will be created during the running of program.

So, we do not need to run program many times to create the Pareto optimal set, just one run and get the Pareto optimal set.

Regards,

Afshin


Afshin's picture
Submitted by Afshin on Thu, 05/07/2007 - 1:37pm.

Hi and Thanks Sameh,

In the multi-objective version of the Bees Algorithm, all the selected sites for neighbourhood search are in the same level of attention, so we do not need to devide them to elits and other selected sites and hence the recruited bees for all the selected sites will be equal. This theory helped us to decrease the parameters in the Bees Algorithm by 2.

At the moment, this theory is suitable for multi-objective verion but it is a good point to see how can we implement to the single function optimisation version too.

Regards,

Afshin


Sameh's picture
Submitted by Sameh on Thu, 05/07/2007 - 2:02pm.

are they e and nep?

Sameh Otri
PhD Research Student


Afshin's picture
Submitted by Afshin on Thu, 05/07/2007 - 2:07pm.

They are 'e' and 'nsp', as all selected sites are like elits so we do not have other selected sites and hense 'nsp'.

Regards

Afshin


smekwn's picture
Submitted by smekwn on Fri, 06/07/2007 - 12:09pm.

Hi. A very good paper .. just to let you know at the last paragraph in page 1, the word should be "beam" not "bean".

"To show the performance of the algorithm, two
test functions plus a well known engineering problem,
welded bean design are used and the results are
presented."

Thank you.

K.W. Ng


Afshin's picture
Submitted by Afshin on Fri, 06/07/2007 - 12:33pm.

Hi Ng,

Thanks for your hint. The mistake has been corrected and the corrected paper is re-uploaded.

Many thanks,

Afshin


smekwn's picture
Submitted by smekwn on Fri, 06/07/2007 - 1:07pm.

Hi Afshin,

You're welcome. Now your pdf file cannot be opened.
Error message "File does not begin with %pdf." Is it because
of transition process?

Thank you.

K.W.Ng


soroka's picture
Submitted by soroka on Fri, 06/07/2007 - 1:25pm.

Your paper isn't recognised as a pdf file, you didn't by any chance upload the word version of the file?

Anthony


Afshin's picture
Submitted by Afshin on Fri, 06/07/2007 - 1:32pm.

Hi Anthony,

As far as I remember, I re-uploaded the pdf one, but will upload it again.

Regards,

Afshin


Afshin's picture
Submitted by Afshin on Fri, 06/07/2007 - 1:34pm.

Seems everything is back to normal !!!!

Afshin


zaidimuhamad's picture
Submitted by zaidimuhamad on Fri, 06/07/2007 - 2:38pm.

I think this is a very clear and nice presentation. Well done Afshin.


Sameh's picture
Submitted by Sameh on Fri, 06/07/2007 - 3:02pm.

It is time to invite (recruit) new non-bees to our hive.

please call every one in your emailing list to participate..

Sameh Otri
PhD Research Student


Vlado's picture
Submitted by Vlado on Wed, 11/07/2007 - 10:39am.

So, to clarify if I understand well

You are using the bees to generate solutions with the tendency towards improvement and use a meta-heuristic to seed out dominated solutions. It makes sense.

Do you have any guidance about the choice of parameters to improve the quality of the Pareto optimal set found?

Does the number of Pareto optimal points found correlate to the number of sites maintained?


Afshin's picture
Submitted by Afshin on Wed, 11/07/2007 - 6:19pm.

Hi Vlado,

The parameters are chosen imperically, but it is not like that the certain parameter values will give the acceptable answers, if they are chosen in a normal range we will get acceptable answers but maybe in more iterations.

The number of Pareto optimal points is corelated to the number of selected sites and number of iterations, for example if we decrease the number of selected sites but increase the number of iterations, we will find roughly the same number of Pareto optimal points.

Hope I am giving the answer.

Afshin


mazhari's picture
Submitted by mazhari on Wed, 11/07/2007 - 9:49pm.

Dear Afshin, I read your very interesting paper. I found it really interesting and would like to ask you whether it can be extended for optimisation of the cargo loading in containers in order to reduce the broken stowage of cargo. If so, it will cost a lot to the shipping industry, specially for those planning the container stowage.
Wish you all the best in your carrier and life.

Capt. Shahriar


mazhari's picture
Submitted by mazhari on Wed, 11/07/2007 - 9:55pm.

Dear Afshin,

Would you plaese let me know if the bee algorithem can be used in calculation of the centre of gravity of cargo especially when an specific position is to be set for the centre of gravity of loads and conainers.
If so, then it might be adopted for pre calculation of the centre of gravity of ships prior to dry docking, in order to avoid the ship capsize or unwanted list or trim.

kind regards

Captain Shahriar


Afshin's picture
Submitted by Afshin on Wed, 11/07/2007 - 9:55pm.

Dear Capt. Shahriar

Thanks for your comments and interest in the Bees Algorithm.

This paper is about multi objective function in continious domain, but the bee team are working on the combinatorial problems too and the specific problem you mention is one of the problems they have in their hand. As you have valuable experiences in this field, I am sure my colleagues looking forward to consult about these kind of problems with you.

Regrards,

Ashin


Afshin's picture
Submitted by Afshin on Wed, 11/07/2007 - 10:12pm.

Dear Captain Shahriar,

Thanks again for your nice comments and question. I beleive the Bees Algorithm could be used to optimise the location of the containers in the cargo ship in order to have minimum difference between real C.G. and desired C.G. It would be a very intersting area and I think the next paper will be in this area !!!!

It might be interesting as well to optimise the location of C.G. and broken stowage, using the multi objective Bees Algorithm.

Regards,

Afshin


Vlado's picture
Submitted by Vlado on Thu, 12/07/2007 - 8:59am.

Does this mean that the sites don't get stuck at a local optimum? Or you refer to the effect what the scouts have on the whole process


Sameh's picture
Submitted by Sameh on Thu, 12/07/2007 - 11:06am.

Dear Shahriar,

Optimisations of the cargo loading in containers and calculating the centre of gravity of cargo seem an interesting application and a promising area of research for the bees Algorithm.
As my bee-college Afshin mentioned earlier it is a similar to some of the applications we are dealing with at the moment.
If you have any matrail to read and data sets to run our algorithm on, please email them to me at otris1@cf.ac.uk

Thanks again for your interest and
good luck, captain, in your journey among the bees world!

King-Bee
Sameh
Sameh Otri
PhD Research Student


Ahmed Haj Darwish's picture
Submitted by Ahmed Haj Darwish on Thu, 12/07/2007 - 11:37am.

How many objective functions can the Bees Algorithm work with?


Afshin's picture
Submitted by Afshin on Thu, 12/07/2007 - 11:46am.

Hi Vlado,

You are absolutely right, the algorithm will not stuck in local minimum. In the basic version, this problem was existing and when the algorithm reaches a local optimum, it was stucking there till finding better solution from random search. But, in the new version, when there is no new information to discover in a patch (means bee is located in the local or global optimum), the algorithm records that position and move the bee to other positions as it does not help to keep searching in the neighbourhood of that site.

Regards,

Afshin


Afshin's picture
Submitted by Afshin on Thu, 12/07/2007 - 11:54am.

Hi Ahmed,

Theoritically, there is no resteriction about the number of objective functions but the more the number of objective function, the more number of population and selected sites and iteration needed. This means, increasing the number of objective functions, increases the computation to get Pareto optimal set.

Regards,

Afshin


MASOUD's picture
Submitted by MASOUD on Thu, 12/07/2007 - 11:55am.

Hi there
Sounds very interesting, I was just curious to ask if there is any ready-to use statistical package or syntax commands for the optimization processing.
In addition, I was thinking to ask if the complete control over independent factors is an absolute prerequisite of your model. A “Yes” on “No” answer to this question could play a key role on your model replicability, for instance, on replication of bee algorithm on data from human sources and mental constructs like creativity, intelligence or addiction and …).

tnx,
Masoud


Afshin's picture
Submitted by Afshin on Thu, 12/07/2007 - 12:11pm.

Dear Masoud,

Thanks for your comment. There are some packages which include optimisatio toolbox like Matlab. In these packages you can find toolboxes for optimisation with algorithms like Genetic Algorithm, Simulated Anealing, and ... . But not the Bees Algorithm at the moment maybe they have not discovered how good is the Bees Algorithm !!!. I am sure they will include the Bees Algorithm in their packages very soon, meantime you are so welcome the software provided in Cardiff Bay Bees to optimise single and multi objective functions.

Regarding the second question, the Bees Algorithm does not need any specific information about the problem unless the objective function. But in the end the algorithm will come up with some values for the parameters which withoes parameter value we will recieve the optimum objective function. If some there is a resteriction about the parameters values, these limitations can be stated as constraints in the algorithm and algorithm avoids to choose a value in unfeasiable regions.

Hope could give you some explanations and please do not hesitate to contact if there is any other issues.

Regards,

Afshin


Sameh's picture
Submitted by Sameh on Thu, 12/07/2007 - 12:21pm.

How do you judge if there is new information or not?

What criteria do u use?

Sameh Otri
PhD Research Student


Afshin's picture
Submitted by Afshin on Thu, 12/07/2007 - 12:28pm.

Dear Sameh,

In the new version of the Bees Algorithm two parts have been added (detals are in the Lapmamap 2007 paper titled Preliminary Design using the Bees Algorithm) shrinking and abondonment. If after some iterations, the representative bee in a site is not changing, it means that the neighbourhood search can not find any other better solution in that specific patch size, so we shrink the neighbourhood size to do more detal searching and if after some shrinkings still the representative bee is unchanged, it is assumed that the bee is located on a peak ( which can bee local or global optimum ).

Regards,

Afshin


Sameh's picture
Submitted by Sameh on Thu, 12/07/2007 - 1:06pm.

What if these criteria (number of shrinking and the shrinking step) were not right? [i.e. you stopped shrinking after 20 times but if u made the shrink number 21st you may find better results or when you select the shrinking step it was slightly greater the position of the optimum]

In these cases, the Bees Algorithm may miss the global optimum as well (if we consider the patch as a search space and zoom in it) as you will stop search around that patch?

So is there a real guarantee that the Bees Algorithm will not stuck in local optimum?

Thanks again

Sameh Otri
PhD Research Student


urazy's picture
Submitted by urazy on Thu, 12/07/2007 - 1:23pm.

Uraz Yavanoglu

This paper is really good value. However i dont know other optimization results on other algorithms. Do you provide a comparison chart between bee algorithms results vs. other algorithms

Thank You


Afshin's picture
Submitted by Afshin on Thu, 12/07/2007 - 1:30pm.

Thanks Sameh for your detailed questions.

If we do not choose the parameters, defenitaly we will not get good results but the range of parameters which are acceptable or give good results are very wide. In the case of number iteration before abondonment, it should have a good chance to find a better solution in neighbourhood if there is any by trying enough number of iterations. If the optimum is outside the patch size, when algorithm finds a better solution and choose that one as representative, in the next attempt, the patch will be shifted as the new representative be in the centre so the patch will move without changing the size of course. In this case, as algorithm runs, the representative will be more close to the optimum (local or global) and after some iterations, the optimum will be in the patch.

Please let me know if you need more explanation or it raised another question.

Regards,

Afshin


Afshin's picture
Submitted by Afshin on Thu, 12/07/2007 - 1:33pm.

Thanks for your comment Mr Uraz.

There is a comparison with GA (the only result could be found in the litrature).The GA and the Bees Algorithm results are illustrated and can be seen that the Bees Algorithm will give longer Pareto front and more Pareto optimal points.

Regards,

Afshin


MASOUD's picture
Submitted by MASOUD on Thu, 12/07/2007 - 1:44pm.

Dear Afshin
Thanks very much for prompt answer and the compliments for software use.
Just with few more Qs:

Is the accuracy of algorithm's output a function of sample size? What about the measurement erros?

thanks again,
Masoud


Afshin's picture
Submitted by Afshin on Thu, 12/07/2007 - 1:52pm.

Thanks again for your comments.

About the accuracy, I would like to refer you to the Bees paper in IPROMS 2006, which you can find a very good comparison about the accuracy and speed of the algorithm with other optimisation algorithms. But here the comparison can only made by visuality as unfortunately the results we from GA are only shown in the graph. Accordings the graphs illustrated in the paper for results obtained from the GA and the Bees Algorithm, it can be seen that the Bees Algorithm can give slightly more accurate answers (the answers are closer to 0,0).

Regards,

Afshin


mass845's picture
Submitted by mass845 on Fri, 13/07/2007 - 7:48am.

Hi Afshin,

Thanks for the paper. I have some curiosity. In you paper, you used Bees Algorithm in minimising problem. I wonder if we can do the same method for minimising and maximising problem in the same time. If so how?

Are there any related techniques (conventional for example) for the same problem? And why you use optimisation rather than these techniques?

I'm agree with one of the comments that there should have a comparison with GA.

Ady


ang's picture
Submitted by ang on Fri, 13/07/2007 - 10:49am.

Dear authors,

From the article, I understand that you were comparing your work with the work done by Deb[12], but from the figure 4, without quoting the reference, I assumed that you have also produced your own version of NSGA and NSGAII tests and compared with the MO bees algorithm results.

Please clarify.

Thanks

Ang


ang's picture
Submitted by ang on Fri, 13/07/2007 - 11:14am.

Dear authors,

I have a design question here. Your case study is very interesting and basically is to find the "best" parameters, i.e. the setting of weld thickness, weld length, beam thickness and beam width. When the results are generated, from designer perspective, what's the "best" setting (only one is needed) and what is the acceptable beam deflection and the lowest cost based on your graph? I believed designer needs a final answer to carry on with his work. I can't find the best result in the paper. Is it possible to know the best result for MO bees algorithm and the best result for NSGA and NSGAII?

By only comparing the number of non-dominated solutions found may not justify fully that the MO bees algorithm is better. The quality of the results need to be emphasized. Do you agree?

Thanks

Ang


Afshin's picture
Submitted by Afshin on Fri, 13/07/2007 - 11:20am.

Thanks for your comments. Maximising f(x) is equal to minimising -f(x) and vice versa. So we can change any maximising problem to another minimising problem and vice versa.

Actually, in the Bees Algorithm maximising problems will be solved. so if we have minimising problem, first we change them to maximising problem then solve them.

There is acomparison with GA and can bee seen that the Bees Algorithm can find longer Pareto front than GA.

This problem is an optimisation problem and there are other techniques and algorithms for solving multi-objective function optimisation problems (Multi-Objective Optimization using Evolutionary Algorithms, by Deb). Here the authors are introducing a new algorithm for solving multi-objective function problems called the Bees Algorithm for multi-objective function.

Regards,

Afshin


Afshin's picture
Submitted by Afshin on Fri, 13/07/2007 - 11:23am.

You are right, the reults are compared with results reported in Deb[12] and we did not test the NSGA and NSGAII.

Regrads,

Afshin


Afshin's picture
Submitted by Afshin on Fri, 13/07/2007 - 11:37am.

Thanks for your comments. In multi-objective problems, there are two stages; first optimisation which will come up with Pareto optimal set or set of non-dominated solutions, and second, decision making or finding the most suitable solution among the Pareto optimal set. In the paper, the first stage was explained and Multi-objective Bees Algorithm is developed to generate the Pareto optimal set. Actually, decision making stage is not optimisation problem any more and in this stage according the specification of the problem and other factors like availabilty of materials, one of the non-dominated solutions will be chosen.

In the specified problem, the goal is to produced as much as possible more non-dominated solutions and of course the Pareto front to be as long as possible.

There is another paper about the same problem but with only onr objective function, and in that case the algorithm came up with one solution rather than a set. This paper is submitted to publish in IMCHE, and has a very good numerical comparison with GA and other techniques.

Regards,

Afshin


ang's picture
Submitted by ang on Fri, 13/07/2007 - 12:36pm.

Dear Afshin,

Currently I have no access to the paper of Deb[12]. But I believe you have the paper, so I have a number of questions here:
Do you know what is the population size being used in their test using NSGA and NSGAII?
Do they use the similar technique like the one you use where you collect pareto-front solutions and update them whenever there's a better one?

Thanks

Ang


Afshin's picture
Submitted by Afshin on Fri, 13/07/2007 - 12:54pm.

Hi Ang,

The refernce number 12 is a Technical Report and you can find it online with this address: http://www.iitk.ac.in/kangal/reports.shtml#2000
Technical Report No. 200002

As you can see in the report, unfortunately there is no information about the population size and other parameters used in thd GAs. The technique for creating the Pareto optimal set is also a bit different as the technique we used in the Multi-Objective Bees Algorithm is based on the carasteristics of the Bees Algorithm.

Regards,

Afshin


Afshin's picture
Submitted by Afshin on Sun, 15/07/2007 - 11:30pm.

I would like to thank IPROMS coference chairman and all the people involved in the conference, to give me this opportunity to present my paper and all the participants who left comments or questions to clarify the paper.

Finally I would like to thank the Co-chair, Sameh Bee for his encouragement and support.

I will keep watching this page for any possible question or comment frequently but feel free to send me email as well on the email addres : sceag1@cf.ac.uk

Regards

Regards,

Afshin


Sameh's picture
Submitted by Sameh on Mon, 16/07/2007 - 6:58am.

Thank you Afshin for your paper, presentation and comments

Sameh Otri
PhD Research Student


bingG's picture
Submitted by bingG on Mon, 16/07/2007 - 8:52am.

Dear author:
Thank you for the great bees application paper. I'm wondering is there any mutation performance like GA in bees during optimization progress? Is that possible to make a small scout of bees (missing their ways) to search in an undesigned route for possibly faster searching speed?


Comment viewing options

Select your preferred way to display the comments and click "Save settings" to activate your changes.

Who's online

There are currently 0 users and 108 guests online.