2016-01-07 7 views
8

Wie kann ich die X-Achsen-Grenzwerte für jedes Panel in einem stratigraphischen Diagramm manuell anpassen?Ändern der X-Achsen-Grenzwerte für stratigraphische Diagramme (z. B. Multi-Panel-Diagramme)

Zum Beispiel, hier ist Stratiplot von analogue:

library(analogue) 
data(V12.122) 
Depths <- as.numeric(rownames(V12.122)) 

(plt <- Stratiplot(Depths ~ O.univ + G.ruber + G.tenel + G.pacR, 
        data = V12.122, type = c("h","l","g","smooth"))) 

enter image description here

Wie könnte ich zum Beispiel der Xlim von G.ruber bis c verändern (0,3, 0,9) und G.pacR zu c (0, 0.75)?

Oder eine andere Möglichkeit, hier ist strat.plot von rioja:

library(rioja) 
library(vegan) ## decorana 
data(RLGH) 
# remove less abundant taxa 
mx <- apply(RLGH$spec, 2, max) 
spec <- RLGH$spec[, mx > 3] 
depth <- RLGH$depths$Depth 
#basic stratigraphic plot 
strat.plot(spec, y.rev=TRUE) 

enter image description here

Wie könnte ich zum Beispiel die Xlim von TA004A bis c ändern (0, 20)?

Ich denke, dass ich etwas liefern muss, um den zugrunde liegenden Gitter/Basis-Plot-Code zu adressieren, aber ich bin mir nicht sicher, wie ich damit anfangen soll.

+0

Vielen Dank für einen Blick. Ich hoffte auf etwas "Magie", um einige Argumente zu überbringen, aber das scheint nicht der Fall zu sein. – Ben

Antwort

5

Dies ist eine kurze Teil Antwort durch die Grenzen des Grundstücks Objekt zwicken und das Gitter Extra-Paket mit, um die Platten die Größe

# your data 
library(analogue) 
data(V12.122) 
Depths <- as.numeric(rownames(V12.122)) 

(plt <- Stratiplot(Depths ~ O.univ + G.ruber + G.tenel + G.pacR, 
        data = V12.122, type = c("h","l","g","smooth"))) 

# manually change the limits of second panel 
# this auto updates ticks and labels 
plt$x.limits[[2]] <- c(0.25,0.95) ; 

# resize the panels 
latticeExtra::resizePanels(plt, w=c(5,5,5,5)) 

die

enter image description here

nicht Dies aber tut gibt Geben Sie die kleine Leerstelle am Anfang jedes Segments ein.

+1

Toll, danke, Ihre Antwort hat mir wirklich die Tür geöffnet, um diese Plots weiter anzupassen, sehr hilfreich. – Ben

+0

sind Sie sehr willkommen. (Wenn Sie sich den Bearbeitungsverlauf anschauen, sehen Sie eine Möglichkeit, die Grobs direkt zu bearbeiten und ein bisschen mehr Kontrolle zu bekommen - Sie erhalten die Grobnamen, indem Sie entweder 'getNames()' oder 'grid.ls (print = FALSE) aufrufen $ name' nach der Handlung) – user20650

+0

Ich sehe es, danke nochmal! – Ben

Verwandte Themen