vrstvy = QgsProject.instance().mapLayersByName('kraje_sjtsk')
vrstva = vrstvy[0]
vrstva.selectByExpression('"NM2"=\'Košický\'')

kraj = "D:/PHD/výučba/Python/cvičenia/vrstvy_python/ke_kraj_sjtsk.shp"
zapis = QgsVectorFileWriter.writeAsVectorFormat(vrstva, kraj, "utf-8", 
driverName = "ESRI shapefile", onlySelected=True)

iface.addVectorLayer(kraj, "", "OGR")