student-t-Konfidenzintervall in python

Ich bin daran interessiert, mit python zu berechnen ein Konfidenzintervall von einem student t.

Ich bin mit dem StudentTCI () - Funktion in Mathematica und müssen nun code, der die gleiche Funktion in python http://reference.wolfram.com/mathematica/HypothesisTesting/ref/StudentTCI.html

Ich bin nicht ganz sicher, wie man diese Funktion selber, aber bevor ich begeben sich auf, die, ist diese Funktion in python irgendwo? Wie numpy? (Ich habe nicht verwendet numpy und mein Berater empfohlen, nicht mit numpy, wenn möglich).

Was wäre der einfachste Weg, dieses problem zu lösen? Kann ich kopieren Sie den Quellcode aus dem StudentTCI() in numpy (falls vorhanden) in meinem code, wie die definition einer Funktion?

edit: ich werde bauen müssen die Schüler TCI mit python-code (wenn möglich). Die Installation von scipy hat sich in eine Sackgasse. Ich habe das gleiche problem, alle anderen hat, und es gibt keinen Weg, ich kann verlangen, Scipy für den code, den ich verteilen, wenn es so lange dauern, einzurichten.

Jemand wissen, wie der Blick auf den source-code für den Algorithmus in der scipy-version? Ich denke, ich werde die Umgestaltung in ein python-definition.

InformationsquelleAutor SwimBikeRun | 2013-06-20
  1. 13

    Ich denke, Sie könnte verwenden scipy.stats.t und es ist interval Methode:

    In [1]: from scipy.stats import t
In [2]: t.interval(0.95, 10, loc=1, scale=2)  # 95% confidence interval
Out[2]: (-3.4562777039298762, 5.4562777039298762)
In [3]: t.interval(0.99, 10, loc=1, scale=2)  # 99% confidence interval
Out[3]: (-5.338545334351676, 7.338545334351676)

    Sicher, Sie können Ihre eigene Funktion, wenn Sie möchten. Lasst es so Aussehen, wie in Mathematica:

    from scipy.stats import t


def StudentTCI(loc, scale, df, alpha=0.95):
    return t.interval(alpha, df, loc, scale)

print StudentTCI(1, 2, 10)
print StudentTCI(1, 2, 10, 0.99)

    Ergebnis:

    (-3.4562777039298762, 5.4562777039298762)
(-5.338545334351676, 7.338545334351676)
    • Danke. Dies ist, was ich brauche. Ich bin versucht zu installieren, scipy (win7 64x), und es ist nicht so einfach, wie es sein sollte?? Ich bin mit 32-bit-python also ich glaube, ich möchte die 32 bit-scipy als gut, aber ich bin mir nicht sicher, wie Sie es bekommen. Ich installiere ein Ding namens anaconda jetzt was ich denke, sollte den job zu erledigen? Da gibt es keine Art, einfach zu installieren, scipy ich es nehmen?
    • Anaconda ist getan, und es hat nicht scipy installieren. line 29, in <module> von scipy.stats-import t ImportError: No module named scipy.Statistiken
    • Ich schlage vor, Sie installieren Enthought Baldachin. Es bietet eine einfache und unkomplizierte Installation und Verwaltung von python-packages. Scipy ist in der Express-version, so dass Sie nicht haben, um es zu bezahlen. Es gibt auch pythonxy kommt mit scipy als gut, aber ich habe nie wirklich genutzt.
    • Nur für die künftigen Nutzer, ich verbrachte viel zu viel Zeit bei der Installation von virtual box + linux, und mehrere dual-boot-Linux-Versionen (ubuntu, mint) und jetzt bin ich wieder in windows 7 64 bit aber mit der läuft alles wie 32-bit-tools. Alles funktioniert nun Super.
    • Wo bekommt man die Dokumentation für die t.Intervall-Funktion? Ich kann nicht herausfinden, was die loc-und scale-diesem - docs.scipy.org/doc/scipy/reference/generated/...
    • loc ist ein Lage-parameter oder den Mittelwert μ. Maßstab ist in diesem Fall die Varianz σ. Wenn Sie nicht möchten, zu verwenden SciPy finden Sie den Quellcode für die distribution und es Methoden hier. SciPy jedoch setzt auf NumPy daher die Benutzer haben, um es zu installieren.
    • Ah ich sehe. Ich bin versucht zu kopieren über den Quellcode, so dass meine Benutzer nicht haben, um es zu installieren. So auf Linie 488,489, sehe ich def Abstand(self, alpha): return self.dist.Intervall(alpha -, *selbst.args, **self.kwds). Dann kann ich nicht finden, wo dieses dist.inverval definiert ist. Wo befindet sich
    • github.com/scipy/scipy/blob/v0.13.3/scipy/stats/...
    • Der Rückgabewert ist die zwei-tailed-test. Ist es möglich, den ein tailed-test? Zum Beispiel hier, um die one-tailed test von 0,95 und 10 Grad?

    InformationsquelleAutor abudis
Schreibe einen Kommentar