Nichtlineare Gleichung mit Summierungen in R

Question

Ich habe die schwierigste Zeit versucht, diese Gleichung in einen nichtlinearen Löser in R zu implementieren. Ich probiere sowohl die nleqslv und BB Pakete, aber bisher nichts als Fehler. Ich habe die Dokumentation gesucht und gelesen, bis meine Augen geblutet haben, aber ich kann mein Gehirn nicht darum wickeln. Die Gleichung selbst funktioniert wie folgt:

The Equation

s2 * sum(price^(2*x+2)) - s2.bar * sum(price^(2*x)) = 0

Wo s2, s2.bar und price Vektoren gleicher Länge bekannt sind.

Der letzte Versuch, den ich in BB versuchte, war dies:

gamma = function(x){ 
     n = len(x) 
     f = numeric(n) 
     f[n] = s2*sum(price^(2*x[n]+2)) - s2.bar*sum(price^(2*x[n])) 
     f 
     } 

g0 = rnorm(length(price)) 
results = BBsolve(par=g0, fn=gamma) 

Jede Hilfe sehr

geschätzt

Answer 1

Von Ihnen Beschreibung der verschiedenen Teile in der Funktion können Sie die Formel zu haben scheinen verwendet durcheinander.

Ihre Funktion gamma sollte höchstwahrscheinlich als

gamma <- function(x){ 
    f <- s2*sum(price^(2*x+2)) - s2.bar*sum(price^(2*x)) 
    f 
} 

s2 sind price und s2.bar Vektoren, die aus Ihrer Beschreibung geschrieben werden, so dass die Formel, die Sie gab einen Vektor zurück.

Da Sie keine Daten angegeben haben, können wir nicht testen. Ich habe versucht, mit zufällig generierten Werten für s2, price, s2.bar zu testen. Manchmal bekommt man eine Lösung mit nleqslv und BB aber nicht immer.

Im Fall des Pakets nleqslv funktioniert die Standardmethode nicht immer. Da das Paket verschiedene Methoden hat, sollten Sie die Funktion testnslv aus dem Paket verwenden, um festzustellen, ob eine der bereitgestellten Methoden eine Lösung findet.