In this paper, we study inference methods for regression analysis of clustered current status data with informative cluster sizes. When the correlated failure times of interest arise from a general class of semiparametric transformation frailty models, we develop a nonparametric maximum likelihood estimation based method for regression analysis and conduct an expectation-maximization algorithm to implement it. The asymptotic properties including consistency and asymptotic normality of the proposed estimators are established. Extensive simulation studies are conducted and indicate that the proposed method works well. The developed approach is applied to analyze a real-life data set from a tumorigenicity study.