Entspricht der Funktion „expand.grid“ in RVerehrte Gäste,
Ein Kunde bittet darum, wenn möglich die Konturen für seine Modellvorhersagen und/oder Vorhersageintervalle über den gesamten Bereich möglicher Werte für meine unabhängigen Variablen darzustellen.
Er erhält es in R (mithilfe der Funktion expand.grid), wobei die gepunktete Linie das untere 95 %-Vorhersageintervall darstellt (Anhang 1).
Das Beste, was er (und ich auch) tun konnten, ist Anhang 2, und kurz gesagt, er und ich waren nicht in der Lage, die gepunktete Linie zu erreichen.(Ich habe auch einen Dummy-Datensatz mit der Grafik)
Er hat die Frage auch in der Community gepostet.
Hast du irgendwelche Vorschläge?
Danke
Betreff: Äquivalent der Funktion „expand.grid“ in RHallo @MarcoSalfi ,
Die Funktion expand.grid verhält sich wie FacetGrid in Seaborn für Python oder die Funktion facet_grid in Rs ggplot2. Grundsätzlich werden alle möglichen Kombinationen eines oder mehrerer Faktoren berücksichtigt, die als Formel oder Matrix eingegeben werden. Was ist die vollständige R-Codezeile, die mit der Funktion expand.grid verwendet wird? Ich vermute, dass der Rest des Plotcodes das Diagramm für die Dichte mit dem gestrichelten Linienfeld liefert.
Betreff: Äquivalent der Funktion „expand.grid“ in RHallo Marco,
Eine Möglichkeit könnte darin bestehen, ein raumfüllendes Design zu erstellen, um die Bandbreite der Möglichkeiten, die Ihr Benutzer erkunden möchte, angemessen abzudecken, und dann die Formel für alle Kombinationen im raumfüllenden Design zu bewerten.
Ich würde jedoch versuchen zu fragen, welche Informationen er durch die Erstellung der Handlung erhalten möchte. Eine bessere Alternative könnte ein Profiler-Graph der Formel oder des Regressionsmodells sein – ohne und mit Extrapolationskontrolle. Der Profiler ist eine viel bessere Möglichkeit, die Komplexität der hohen Dimensionalität dieses Problems zu bändigen. Möglicherweise möchten wir die übliche Warnung vor der Gefahr einer Extrapolation außerhalb des X-Raums (oder Designraums) aussprechen, wenn empirische Modelle verwendet werden.
Malcolm
Betreff: Äquivalent der Funktion „expand.grid“ in RHi @Chris_Kirchberg , folgend dem vom Kunden verwendeten Code:
## Variablen transformieren
log.AlphaLISA = log(AlphaLISA)
log.rabbit = log(Rabbit.Ab)
log.mouse = log(Mouse.Ab)
log.lysate = log(Lysat)
fit = lm( Log.AlphaLISA ~ Ab*(log.rabbit+log.mouse+log.lysate)^2 +
Ab*I(log.rabbit^2) + Ab*I(log.mouse^2) +
Ab*I(log.lysate^2) )
d.grid = expand.grid(
Ab=c("R1-M1", "R1-M2", "R2-M1", "R2-M2", "R3-M1", "R3-M2"),
log.rabbit=seq(log(0.3), log(3), length=5),
log.mouse=seq(log(0.3), log(3), length=10),
log.lysate=seq(log(0,25), log(2), Länge=10)
)
## 95 % niedrigeres Vorhersageintervall
d.grid = cbind(d.grid,
Vorhersage(fit, newdata=d.grid, Intervall="Vorhersage", Ebene=0,9))
d.grid[,"lower95p"] = exp(d.grid[,"lwr"])
Betreff: Äquivalent der Funktion „expand.grid“ in RDanke Malcolm,
Die Diskussion geht in diese Richtung. Tatsächlich ist der Kunde auch für andere Optionen offen.
