Modeling Binary Correlated Responses using SAS, SPSS and R