Here is the code I am writing, however I was unable to double the bar for each feature:
[Xr ,Yr ,feat_names] = load_dataset_sc(\'boavista_covid_dados_abertos_