Problem mit <=
Verfasst: Samstag 23. November 2013, 22:42
Hallo ich habe mir folgenden Code überlegt. Ein gesuchtes Ergebnis ist taum1. Das ist das Maximum der Funktion lnl. Dann soll bestimmt werden bei welchem Wert für tau lnl jeweils um 0.5 vom Maximum abweicht. Eigentlich würde ich erwarten, dass in der Variablen df=lnl-lmax>=-0.5 alle Werte der Funktion lnl drin stecken, die grösser als -0.5 sind. Wenn ich mir angucke was in der Variablen df drinnsteckt sind da aber nur noch nullen eingetragen, anstatt zahlen. ich brauche aber alle Zahlen die die Bedingung erfüllen. Für tauf sollte dann ein Wert ungleich null herauskommen. Es wäre schön wenn mir jemand sagen könnte, wo mein Fehler liegt.
Code: Alles auswählen
from pylab import *
import scipy.stats as stats
hold(True)
tautrue=5
t=stats.expon(0.0,tautrue).rvs(1000)
tau = frange(3.5,6.5,0.03)
n=len(t)
lnl =-sum(t)/tau-n* log(tau)
lmax=max(lnl)
taum1=tau[argmax(lnl)]
print taum1
df=lnl-lmax>=-0.5
taut=min(df)
tauh=max(df)
tauf=(taut - tauh)/2