For a nested model you want to use the coxme function, which is the much superior successor to frailty(). It is currently found in the kinship library. Terry Therneau