I know ivprobit function can be implemented well in cross-sectional data, and I see we can use cluster() command if we have panel data in Stata. However, what can I do in R in this case?