Python stats.circstd方法代碼示例

本文整理匯總了Python中scipy.stats.circstd方法的典型用法代碼示例。如果您正苦於以下問題：Python stats.circstd方法的具體用法？Python stats.circstd怎麽用？Python stats.circstd使用的例子？那麽, 這裏精選的方法代碼示例或許可以為您提供幫助。您也可以進一步了解該方法所在類scipy.stats的用法示例。

在下文中一共展示了stats.circstd方法的15個代碼示例，這些例子默認根據受歡迎程度排序。您可以為喜歡或者感覺有用的代碼點讚，您的評價將有助於係統推薦出更棒的Python代碼示例。

示例1: test_circstd_axis

# 需要導入模塊: from scipy import stats [as 別名]
# 或者: from scipy.stats import circstd [as 別名]
def test_circstd_axis(self):
        x = np.array([[355,5,2,359,10,350],
                      [351,7,4,352,9,349],
                      [357,9,8,358,4,356]])

        S1 = stats.circstd(x, high=360)
        S2 = stats.circstd(x.ravel(), high=360)
        assert_allclose(S1, S2, rtol=1e-11)

        S1 = stats.circstd(x, high=360, axis=1)
        S2 = [stats.circstd(x[i], high=360) for i in range(x.shape[0])]
        assert_allclose(S1, S2, rtol=1e-11)

        S1 = stats.circstd(x, high=360, axis=0)
        S2 = [stats.circstd(x[:,i], high=360) for i in range(x.shape[1])]
        assert_allclose(S1, S2, rtol=1e-11)

開發者ID:ktraunmueller，項目名稱:Computable，代碼行數:18，代碼來源:test_morestats.py

示例2: test_circstd_axis

# 需要導入模塊: from scipy import stats [as 別名]
# 或者: from scipy.stats import circstd [as 別名]
def test_circstd_axis(self):
        x = np.array([[355, 5, 2, 359, 10, 350],
                      [351, 7, 4, 352, 9, 349],
                      [357, 9, 8, 358, 4, 356]])

        S1 = stats.circstd(x, high=360)
        S2 = stats.circstd(x.ravel(), high=360)
        assert_allclose(S1, S2, rtol=1e-11)

        S1 = stats.circstd(x, high=360, axis=1)
        S2 = [stats.circstd(x[i], high=360) for i in range(x.shape[0])]
        assert_allclose(S1, S2, rtol=1e-11)

        S1 = stats.circstd(x, high=360, axis=0)
        S2 = [stats.circstd(x[:, i], high=360) for i in range(x.shape[1])]
        assert_allclose(S1, S2, rtol=1e-11)

開發者ID:Relph1119，項目名稱:GraphicDesignPatternByPython，代碼行數:18，代碼來源:test_morestats.py

示例3: test_circstd

# 需要導入模塊: from scipy import stats [as 別名]
# 或者: from scipy.stats import circstd [as 別名]
def test_circstd(self):
        """ Test custom circular std."""

        ref_std = scistats.circstd(self.test_angles, **self.circ_kwargs)
        test_std = pystats.nan_circstd(self.test_angles, **self.circ_kwargs)

        assert ref_std == test_std

開發者ID:pysat，項目名稱:pysat，代碼行數:9，代碼來源:test_utils_stats.py

示例4: test_circstd_nan

# 需要導入模塊: from scipy import stats [as 別名]
# 或者: from scipy.stats import circstd [as 別名]
def test_circstd_nan(self):
        """ Test custom circular std with NaN."""

        ref_std = scistats.circstd(self.test_angles, **self.circ_kwargs)
        ref_nan = scistats.circstd(self.test_nan, **self.circ_kwargs)
        test_nan = pystats.nan_circstd(self.test_nan, **self.circ_kwargs)

        assert np.isnan(ref_nan)
        assert ref_std == test_nan

開發者ID:pysat，項目名稱:pysat，代碼行數:11，代碼來源:test_utils_stats.py

示例5: test_deprecation_warning_circstd

# 需要導入模塊: from scipy import stats [as 別名]
# 或者: from scipy.stats import circstd [as 別名]
def test_deprecation_warning_circstd(self):
        """Test if circstd in stats is deprecated"""

        with warnings.catch_warnings(record=True) as war:
            try:
                pystats.nan_circstd(None)
            except TypeError:
                # Setting input to None should produce a TypeError after
                # warning is generated
                pass

        assert len(war) >= 1
        assert war[0].category == DeprecationWarning

開發者ID:pysat，項目名稱:pysat，代碼行數:15，代碼來源:test_utils_stats.py

示例6: test_circfuncs

# 需要導入模塊: from scipy import stats [as 別名]
# 或者: from scipy.stats import circstd [as 別名]
def test_circfuncs(self):
        x = np.array([355,5,2,359,10,350])
        M = stats.circmean(x, high=360)
        Mval = 0.167690146
        assert_allclose(M, Mval, rtol=1e-7)

        V = stats.circvar(x, high=360)
        Vval = 42.51955609
        assert_allclose(V, Vval, rtol=1e-7)

        S = stats.circstd(x, high=360)
        Sval = 6.520702116
        assert_allclose(S, Sval, rtol=1e-7)

開發者ID:ktraunmueller，項目名稱:Computable，代碼行數:15，代碼來源:test_morestats.py

示例7: test_circfuncs_small

# 需要導入模塊: from scipy import stats [as 別名]
# 或者: from scipy.stats import circstd [as 別名]
def test_circfuncs_small(self):
        x = np.array([20,21,22,18,19,20.5,19.2])
        M1 = x.mean()
        M2 = stats.circmean(x, high=360)
        assert_allclose(M2, M1, rtol=1e-5)

        V1 = x.var()
        V2 = stats.circvar(x, high=360)
        assert_allclose(V2, V1, rtol=1e-4)

        S1 = x.std()
        S2 = stats.circstd(x, high=360)
        assert_allclose(S2, S1, rtol=1e-4)

開發者ID:ktraunmueller，項目名稱:Computable，代碼行數:15，代碼來源:test_morestats.py

示例8: test_circfuncs_array_like

# 需要導入模塊: from scipy import stats [as 別名]
# 或者: from scipy.stats import circstd [as 別名]
def test_circfuncs_array_like(self):
        x = [355,5,2,359,10,350]
        assert_allclose(stats.circmean(x, high=360), 0.167690146, rtol=1e-7)
        assert_allclose(stats.circvar(x, high=360), 42.51955609, rtol=1e-7)
        assert_allclose(stats.circstd(x, high=360), 6.520702116, rtol=1e-7)

開發者ID:ktraunmueller，項目名稱:Computable，代碼行數:7，代碼來源:test_morestats.py

示例9: test_empty

# 需要導入模塊: from scipy import stats [as 別名]
# 或者: from scipy.stats import circstd [as 別名]
def test_empty(self):
        assert_(np.isnan(stats.circmean([])))
        assert_(np.isnan(stats.circstd([])))
        assert_(np.isnan(stats.circvar([])))

開發者ID:ktraunmueller，項目名稱:Computable，代碼行數:6，代碼來源:test_morestats.py

示例10: test_circfuncs

# 需要導入模塊: from scipy import stats [as 別名]
# 或者: from scipy.stats import circstd [as 別名]
def test_circfuncs(self):
        x = np.array([355, 5, 2, 359, 10, 350])
        M = stats.circmean(x, high=360)
        Mval = 0.167690146
        assert_allclose(M, Mval, rtol=1e-7)

        V = stats.circvar(x, high=360)
        Vval = 42.51955609
        assert_allclose(V, Vval, rtol=1e-7)

        S = stats.circstd(x, high=360)
        Sval = 6.520702116
        assert_allclose(S, Sval, rtol=1e-7)

開發者ID:Relph1119，項目名稱:GraphicDesignPatternByPython，代碼行數:15，代碼來源:test_morestats.py

示例11: test_circfuncs_small

# 需要導入模塊: from scipy import stats [as 別名]
# 或者: from scipy.stats import circstd [as 別名]
def test_circfuncs_small(self):
        x = np.array([20, 21, 22, 18, 19, 20.5, 19.2])
        M1 = x.mean()
        M2 = stats.circmean(x, high=360)
        assert_allclose(M2, M1, rtol=1e-5)

        V1 = x.var()
        V2 = stats.circvar(x, high=360)
        assert_allclose(V2, V1, rtol=1e-4)

        S1 = x.std()
        S2 = stats.circstd(x, high=360)
        assert_allclose(S2, S1, rtol=1e-4)

開發者ID:Relph1119，項目名稱:GraphicDesignPatternByPython，代碼行數:15，代碼來源:test_morestats.py

示例12: test_circfuncs_array_like

# 需要導入模塊: from scipy import stats [as 別名]
# 或者: from scipy.stats import circstd [as 別名]
def test_circfuncs_array_like(self):
        x = [355, 5, 2, 359, 10, 350]
        assert_allclose(stats.circmean(x, high=360), 0.167690146, rtol=1e-7)
        assert_allclose(stats.circvar(x, high=360), 42.51955609, rtol=1e-7)
        assert_allclose(stats.circstd(x, high=360), 6.520702116, rtol=1e-7)

開發者ID:Relph1119，項目名稱:GraphicDesignPatternByPython，代碼行數:7，代碼來源:test_morestats.py

示例13: test_circular_standard_deviation_1d

# 需要導入模塊: from scipy import stats [as 別名]
# 或者: from scipy.stats import circstd [as 別名]
def test_circular_standard_deviation_1d(data):
    high = 8
    low = 4
    assert np.allclose(
        _circular_standard_deviation(data, high=high, low=low), circstd(data, high=high, low=low),
    )

開發者ID:arviz-devs，項目名稱:arviz，代碼行數:8，代碼來源:test_stats_utils.py

示例14: radial_stddev

# 需要導入模塊: from scipy import stats [as 別名]
# 或者: from scipy.stats import circstd [as 別名]
def radial_stddev(L):
    scipy = circstd(L, 360,0)
    return scipy

開發者ID:protwis，項目名稱:protwis，代碼行數:5，代碼來源:models.py

示例15: _mc_error

# 需要導入模塊: from scipy import stats [as 別名]
# 或者: from scipy.stats import circstd [as 別名]
def _mc_error(ary, batches=5, circular=False):
    """Calculate the simulation standard error, accounting for non-independent samples.

    The trace is divided into batches, and the standard deviation of the batch
    means is calculated.

    Parameters
    ----------
    ary : Numpy array
        An array containing MCMC samples
    batches : integer
        Number of batches
    circular : bool
        Whether to compute the error taking into account `ary` is a circular variable
        (in the range [-np.pi, np.pi]) or not. Defaults to False (i.e non-circular variables).

    Returns
    -------
    mc_error : float
        Simulation standard error
    """
    _numba_flag = Numba.numba_flag
    if ary.ndim > 1:

        dims = np.shape(ary)
        trace = np.transpose([t.ravel() for t in ary])

        return np.reshape([_mc_error(t, batches) for t in trace], dims[1:])

    else:
        if _not_valid(ary, check_shape=False):
            return np.nan
        if batches == 1:
            if circular:
                if _numba_flag:
                    std = _circular_standard_deviation(ary, high=np.pi, low=-np.pi)
                else:
                    std = stats.circstd(ary, high=np.pi, low=-np.pi)
            else:
                if _numba_flag:
                    std = np.float(_sqrt(svar(ary), np.zeros(1)))
                else:
                    std = np.std(ary)
            return std / np.sqrt(len(ary))

        batched_traces = np.resize(ary, (batches, int(len(ary) / batches)))

        if circular:
            means = stats.circmean(batched_traces, high=np.pi, low=-np.pi, axis=1)
            if _numba_flag:
                std = _circular_standard_deviation(means, high=np.pi, low=-np.pi)
            else:
                std = stats.circstd(means, high=np.pi, low=-np.pi)
        else:
            means = np.mean(batched_traces, 1)
            if _numba_flag:
                std = _sqrt(svar(means), np.zeros(1))
            else:
                std = np.std(means)

        return std / np.sqrt(batches)

開發者ID:arviz-devs，項目名稱:arviz，代碼行數:63，代碼來源:diagnostics.py

注：本文中的scipy.stats.circstd方法示例由純淨天空整理自Github/MSDocs等開源代碼及文檔管理平台，相關代碼片段篩選自各路編程大神貢獻的開源項目，源碼版權歸原作者所有，傳播和使用請參考對應項目的License；未經允許，請勿轉載。