Hallo,
ich möchte gerne Parameter eines Modells an Messdaten anfitten, stoße
aber auf Probleme, offenbar weil das Modell NIntegrate enthält. Ich
hänge ein Notebook an, hier die wesentlichen Punkte. Das Modell sieht
so aus:
Ecalc[\[Theta]_, gua_] :=
NIntegrate[
P[r]EvonR[r, \[Theta], n[gua]],
{r, 0, lc}
]
Dabei ist gua ein Parameter, der verschiedene Datensätze unterscheidet
(sie wurden also bei unterschiedlichen Konzentrationen der Chemikalie
Guanidiniumchlorid gemessen), und \theta ist der Parameter, der in einem
Datensatz variiert wurde (die Temperatur in Grad Celsius). Die
eigentlichen Fitparameter stecken in der Definition von EvonR.
Ich habe nun so zu fitten versucht:
NonlinearRegress[Data[0], Ecalc[\[Theta], n[0]], {\[Theta]}, FitParams]
Aber das liefert nur "Integrand is not numerical". Statt EvonR die
Funktion explizit hinzuschreiben ändert wie erwartet nichts.
Und eigentlich ist es noch komplizierter, weil ich mehrere Datensätze
bei unterschiedlichen Werten von gua gemessen habe und diese global
fitten möchte. Ich habe dazu bissher immer selber ein ChiSquared
definiert (mit Map[(#[[2]]-Ecalc[#[[1]],gua})^,Data]) und dann NMinimize
verwendet, aber das steigt genauso aus.
Gibt es eine Möglichkeit, dieses Problem zu lösen?
Danke im Voraus,
Frank
--
Frank Küster
Single Molecule Spectroscopy, Protein Folding @ Inst. f. Biochemie, Univ. Zürich
Debian Developer (teTeX/TeXLive)
Minimize-Integral.nb
Description: Mathematica Notebook document
