Propensity Score Matching for more than 2 groups I'm new to propensity score matching (PSM). So, my questions can be bit trivial.
1) Suppose I've 3 treatment levels and want to check the effectiveness of the treatment levels. Treatment levels are taking drug on time, not taking drug on time and not taking drug at regular interval. For this I need to do multinomial logistic regression. 
But in PSM we do case-control study. So, how we are going to define which will be the case and which will be the control? Will it be the case that we will use one group as control for each time and other 2 groups as case?
2) Can anyone please tell me which package to use for multilevel group in R. I checked this link. But this link is old. I also checked this package which seems to do multilevel. But is there any other option for packages?
 A: Check out the WeightIt package. You can simply provide a factor treatment variable and covariates and it will estimate balancing weights for that treatment. It provides an interface to other packages and methods that do this using a unified syntax.
Currently, it provides support for estimating balancing weighting for multinomial treatments using multinomial logistic and probit regression propensity scores, generalized boosted modeling propensity scores (through twang), covariate balancing propensity score (through CBPS), entropy balancing (through ebal), optimization-based weights (through optweight), empirical balancing calibration weights (through ATE), and SuperLearner propensity scores (through SuperLearner). The syntax is the following:
w.out <- weightit(treat ~ cov1 + cov2 + cov3, data = data, method = "ps")

You can change the estimation method with the method parameter and others.
A: Propensity score calculation and subsequent paired analysis is possible in several ways. There are already some overlapping Q&A in CV that you might wish to look at:


*

*Propensity Score Matching in R with Multiple Treatments

*Software that matches 6 groups by propensity score? 

*Comparing
two or more treatments with inverse probablity of treatment
weighting
My advice would be to use the twang R package.
A: I have mostly used PSM for 2 class problems. We predict the probability of the treatment. And then compare effect of treatment vs control in same decile of our probability scores. Customers in same decile typically are similar and so comparable. 
So you can repeat the PSM twice. Control each time being people who have taken medicine on time. Though i am not sure if doing the test will lead to higher type 1 errors like in multiple t tests. 
For PSM in r for 2 variables there is a package matchit but have not used it. 
