you could try "nlsur" (combining a linear and a censored equation) and vce(cluster ..) for your panel variable... however, I would recommend you to have look in the manual and in the code of nlsur with regard to what the cluster option exactly does with this command.