Я использую пакет rlm R и экспериментирую с надежной регрессией, используя функцию Хубера. Вот мой код: myfit= rlm(formula = depvar ~ indep1+indep2, init="ls",data = my_input_data,psi =psi.huber, k=0.99,method = "M", maxit=200) k — это...

Я хотел бы построить прогнозы моей пробит-модели и включить доверительный интервал 95% на основе двухсторонних кластеризованных стандартных ошибок, которые я получил с помощью функции cluster.vcov() из пакета multiwaycov. Однако функции...