Isoklinen und/oder Richtungsfeld zum Plotten hinzufügen

Question

Ich löste Differentialgleichungen mit R, als ich auf ein Problem stieß: Ich musste meinem Plot Isoclines und ein Richtungsfeld hinzufügen, wusste aber nicht wie. Welches Paket soll ich installieren/welche Funktion soll ich aufrufen/sollte ich das alles manuell machen? Hier

Answer 1

ist keine Antwort, sondern nur eine Zusammenfassung der in den Kommentaren gegebenen Antworten:

• die LVCompGames in primer Paket gibt dieses Grundstück:

1. Diese ist die Lösung von @Ben Bolker. Die Lösung basiert auf dem graphics Paket und der filled.contour Funktion.

• Dann Vorschlag @ Josh O'Brien das rasterVis Paket mit vectorplot Funktion

Ich denke, alle Antworten angepasst werden können, wenn Sie uns weitere Informationen geben von dem, was Sie als Ausgabe wünschen.

Answer 2

Mit der Funktion my.symbols im TeachingDemos-Paket können Sie eigene Symbole definieren, die zu vorhandenen Plots hinzugefügt werden. Die Unterstützungsfunktion ms.arrows codiert bereits grundlegende Pfeile, die zu einem Plot hinzugefügt werden. Sie können sie einfach (oder die anderen ms. Funktionen) modifizieren, um ein anderes Symbol zum Hinzufügen zu erstellen.

Answer 3

Es gibt ein großes Paket namens phaseR. Es verwendet den ODE-Solver aus dem Paket deSolve und lässt Sie problemlos Nulllinien, Trajektorien und so weiter hinzufügen. Geben Sie nur die Parameter Ihrer ODE-Funktion "t", "y" und "parameter" an. Andernfalls werden die Phaser-Funktionen einen Fehler geben:

Error in deriv(t = 0, y = c(x[i], y[j]), parameters = parameters) :
unused arguments (y = c(x[i], y[j]), parameters = parameters)

Hier ist ein funktionierendes Codebeispiel für die Lotka-Volterra-Räuber-Beute-System:

require(deSolve) 
require(phaseR) 

model.LV <- function(t, y, parameters){ 
    with(as.list(parameters),{ 

    N<-y[1] 
    P<-y[2] 
    dN <- a*N - b*N*P 
    dP <- c*N*P -d*P 

    list(c(dN,dP)) 
    }) 
} 

params.LV<-c(a=0.4, b=0.3, c=0.1, d=0.2) 

data.LV<-as.data.frame(lsoda(c(N=1,P=1),seq(1,250,by=0.5), model.LV, params.LV)) 
# plot the time series of both populations 
plot(data.LV$time,data.LV$N, main="Time series of L-V equations", xlab="time", 
    ylab="Population densities N, P", 
    type="l", col="green", 
    ylim=c(0,max(data.LV$N,data.LV$C))) 
lines(data.LV$time,data.LV$P,col="red") 

# plot the trajectories of the system 
plot(data.LV$N, data.LV$P, type="l", col="blue", main="Trajectory of L-V equations", 
    xlab="Prey density N", ylab="Predator density P", xlim=c(0,5), ylim=c(0,3)) 
#add Nullclines 
nullclines(model.LV, x.lim=c(0.1,5),y.lim=c(0.1,3), parameters=params.LV, system="two.dim", colour=c("green","red"), add=TRUE)