def calc_PearsonR(signalA, signalB, decimals=3, norm=True, fisher=True):
    """ Calculate Pearson R correlation coefficient between
    the time-courses extracted from the given coordinate and given signal

    :param signalA: 1D time-series data array
    :param signalB: 1D time-series data array
    :param norm: perform standard normalization
    :param fisher: perform fisher Z transdormation
    :type signalA: numpy.ndarray
    :type signalB: numpy.ndarray
    :type norm: bool
    :type fisher: bool
    :return: Pearson's correlation coefficient and it's p value
    :rtype return: r, p
    """
    from .signal import standard_norm
    from scipy.stats import pearsonr
    from numpy import arctanh
    if norm:
        signalA = standard_norm(signalA, decimals=decimals)
        signalB = standard_norm(signalB, decimals=decimals)
    r, p = pearsonr(signalA, signalB)
    if fisher:
        r = arctanh(r)
    return r, p
