Добавление нескольких ковариат в «rdplot» в R

Question

Я пытаюсь добавить несколько ковариат в функцию 'rdplot', которая строит модель прерывности регрессии в R. Добавление только одной ковариаты работает нормально следующим образом:

rdplot(dataset$CONTINUANCE, 
       dataset$T,
       p=scenario.polynomial_order,
       nbins=c(scenario.bandwidth, scenario.bandwidth),
       kernel=scenario.kernel,
       weights=dataset$USERS,
       h=h,
       covs=cbind(dataset$VIEW_COUNT))

... тогда как добавление нескольких ковариат делает не:

rdplot(dataset$CONTINUANCE, 
       dataset$T,
       p=scenario.polynomial_order,
       nbins=c(scenario.bandwidth, scenario.bandwidth),
       kernel=scenario.kernel,
       weights=dataset$USERS,
       h=h,
       covs=cbind(dataset$VIEW_COUNT, dataset$COMMENT_COUNT))

Ошибка:

Error in crossprod(U_p_l[, colsZ], invG_p_l %*% U_p_l) : 
  non-conformable arguments
Calls: rdplot -> crossprod

Я что-то не так делаю или есть ошибка в функции 'rdplot'?