youYou may want to use "faisalconjoint"faisalconjoint package in R, it is tested with many published and research data, it works perfectly, one on important thing its works without design restriction and rank procedure. It works in all condition and provide accurate estimates.