fragility_curve.py 文件源码

def sample_damage_state(self, Pr):
        """
        Sample the damage state using a uniform random variable

         Parameters
        -----------
        Pr : pd.Dataframe
            Probability of exceeding a damage state

        Returns
        -------
        damage_state : pd.Series
            The damage state of each element
        """
        p = pd.Series(data = np.random.uniform(size=Pr.shape[0]), index=Pr.index)

        damage_state = pd.Series(data=[None]* Pr.shape[0], index=Pr.index)

        for DS_names in Pr.columns:
            damage_state[p < Pr[DS_names]] = DS_names

        return damage_state