Plotten Einfluss von 2 Variablen auf Ausgang

Question

Hallo Art und weise Internet-Freunde,

Ich habe mit R aus verschiedenen Winkeln liebäugelt, aber ich scheine zu wenig Fortschritte gemacht zu haben. Dies ist wahrscheinlich einfach, aber meine Meinung und Erfahrung ist so ...

Ich mag schließlich den Einfluss unterschiedlichen m plotten und n auf dem Ausgang (b) in der Gleichung:

b = (0,15 m + 0,15 n)/n. Wo m und n haben Bereiche von -1 bis 1.

Ich stelle mir ein Konturdiagramm, um dies zu visualisieren, aber ich bin auf der Stufe des Erhaltens von Werten von b mit entsprechenden m und n Eingaben stecken.

Die neueste Ansatz beteiligt Bootstrapping Schätzungen für m und n b-Werte zurück, aber soweit ich weiß, gibt es keine Möglichkeit, den entsprechenden Eingang m und n Werte zu erhalten:

m<-seq(from=-1,1,length.out=100) 
n<-seq(from=-1,1,length.out=100) 
z<-rnorm(100) 
b<-((0.15*(sample(m,1,replace=T))+0.15*(sample(n,1,replace=T)))/(sample(n,1,replace=T))) 

library("boot") 
bfunc<-function(m,n){ 
    (0.15*(sample(m,1,replace=T))+0.15*(sample(n,1,replace=T)))/(sample(n,1,replace=T)) 
} 
bootb<-boot(data=z,statistic=(bfunc),R=1000) 
bootb$t 

Meine Frage: Wie bekomme ich meine Ausgabe (b) und entsprechende Eingaben (m und n), damit ich die Daten plotten kann? Modifikationen der oben genannten oder ganz andere Wege sind alle willkommen ... Ich muss lernen!

Jede Hilfe wird sehr geschätzt, danke.

Answer 1

Ich sehe nicht, warum Sie das Bootstrap

dataset <- expand.grid(
    m = seq(-1, 1, length.out = 101), 
    n = seq(-1, 1, length.out = 101) 
) 
dataset$b <- (0.15 * dataset$m + 0.15 * dataset$n)/dataset$n 

library(ggplot2) 
ggplot(dataset, aes(x = m, y = n, z = b)) + 
    geom_contour(aes(colour = ..level..)) + 
    scale_colour_gradient2() 
ggplot(dataset, aes(x = m, y = n, fill = b)) + 
    geom_tile() + 
    scale_fill_gradient2()