I'd like to estimate propensity scores for multiple treatments and try to use the mnps command in Stata (https://www.rand.org/statistics/twang/downloads.html). I have followed the codes suggested by online tutorial. However, I am facing the following error: Array

my Code:
use C:\Program Files (x86)\Stata14\ado\base\6.20.dta,clear
global ado "C:\Program Files (x86)\Stata14\ado\base\adofiles"
adopath + "$ado"
mnps i.insurancedstatus i.race i.sex i.diagnisticperiod i.subtype i.summarystage i.age i.marriedstatus i.ses, ///
ntrees(3000) intdepth(3) shrinkage(0.01) ///
stopmethod(es.mean ks.mean) estimand(ATE) ///
rcmd(C:\Program Files\R\R-4.0.2\bin\Rscript.exe) ///
objpath(C:\Program Files (x86)\Stata14\ado\base\output4) ///
plotname(C:\Program Files (x86)\Stata14\ado\base\mnps_example_plot.pdf)
save C:\Program Files(x86)\Stata14\ado\base\aod_ate_wgts,replace

FYI, the following is the codes suggested by online tutorial (https://www.rand.org/statistics/twan...-tutorial.html).
Code:

mnps treat illact crimjust subprob subdep white, ///
ntrees(3000) intdepth(3) shrinkage(0.01) ///
stopmethod(es.mean ks.mean) estimand(ATE) ///
rcmd(C:\Program Files\R\R-3.3.1\bin\Rscript.exe) ///
objpath(C:\Users\username\twang\output) ///
plotname(C:\users\username\twang\output\mnps_examp le_plot.pdf)
save C:\Users\username\twang\output\aod_ate_wgts, replace

Thank you in advance!