Hi fellow R users,
I have a question about something that I've spent some time on but couldn't figure out how to solve. Specifically, as the title of the thread suggests, it is about producing clustered standard errors for ordered logit.
While I know how it is done in Stata (the command is ologit y x, vce(cluster variable_to_cluster_by)), when I use the command polr to run the ordered logit model, I get the same point estimates as what Stata produces, but the standard errors are off, and I can't seem to figure out how to do the same thing in R. Any help is greatly appreciated. Thank you so much!