since you basically do a glm, the treatment can be always have one than one group if I understand it correctly. Thus all propensity score matching packages should be able to deal with it.

Look at the MatchIt package for example.

A list of some further packages are listed here:

http://jason.bryer.org/talks/psaworkshop.html

Thank you for your advice! I found this r package too:

https://cran.r-project.org/web/packages/twang/vignettes/mnps.pdf

Hello

You can use the StatMach Package in Rcran site.

Be successfull

Zahra ghayoumi