def METHOD_NAME(y_true_reference: pd.Series, y_pred_reference: pd.Series) -> Tuple:
    """
    Calculate sampling error components for Root Mean Squared Logarithmic Error (RMSLE) using reference data.
    Parameters
    ----------
    y_true_reference: pd.Series
        Target values for the reference dataset.
    y_pred_reference: pd.Series
        Predictions for the reference dataset.
    Returns
    -------
    (std,): Tuple[np.ndarray]
    """
    squared_log_error = (np.log(1 + y_true_reference) - np.log(1 + y_pred_reference)) ** 2
    squared_log_error_std = np.std(squared_log_error)
    squared_log_error_mean = np.mean(squared_log_error)
    return squared_log_error_std, squared_log_error_mean