Как сгруппировать стандартные ошибки plm на другом уровне, а не по идентификатору или времени?

Я запускаю панельную регрессию по plm. plm может группировать стандартные ошибки только на уровне "группы" или "времени". Но я хочу группировать на уровне страны, а не группы или времени. (Для моих данных группа — это уровень фирмы.)

Я использую следующий код для кластеризации: summary(r, vcov = function(x) vcovHC(x, cluster = "group", method = "arellano")) Но параметры для cluster — это только группа и время. Я хочу кластеризоваться на уровне страны.