Atomic modeling of high-Z partially ionized plasma is essential for simulating radiation hydrodynamics of laser-produced plasma. A collisional-radiative model based upon an average atom model is used to calculate plasma opacity and emissivity. Because line radiations are most dominant in such plasma, the detail configuration accounting (DCA) for electronic state is required. We propose a statistical method to carry out the DCA with the use of the average population of bound electrons. Further modeling of line group made of the same transition from ions in different change states is discussed by considering the detail structure (hierarchy) of the line group.