Based on the method of determination for relative stability of each phase from the difference among the interaction parameters of the phases consisting the mixed layer, the types of interactions between layers were specified and interaction parameter between layers in ordered domain was analytically derived as a function parameter between layers in ordered domain was analytically derived as a function of not only temperature and mole fraction of layers but also ordering parameter. Interaction parameter between the different layers in ordered phase, L is as follows:{{{{ {L }_{1 } (X,Q,T)= { C} over { Q} -4(1-2Q) { L}^{2 } - { RT} over {2} ln { 1} over {2 } - { 2RT} over { { X}_{ s} } ln { { 4QX}`_{s } ^{2 } } over {(1- { X}_{s }- { QX}_{s })( { X}_{s }- {QX }_{s } ) } }}}}L2 is the interaction parameter between ordered and disordered phase in domain and is the mole fraction of the domain which represent the infinite length of mixed layer mineral and Q and C are the reaction progress parameter and arbitrary constant, respectively. This equation was used for the I/S mixed layer clay minerals to infer the relative stability of R1 type I/S mixed layer in the temperature range from 373K to 450K. The result of calculation suggest that, owing to the decrease in interaction parameter with increasing temperature. The interaction parameter decreases more rapidly with decreasing mole fraction of smectite in domain, which is consistent with the fact that the probability of finding the series smectite layer is lo in the domain with small mole fraction of smectite layers in natural system.