Deutsches R-Forum

Code: Alles auswählen

# Ordnerpfad auswählen
setwd("...")

# Tabelle einlesen
dateiname <- ".....dat"
table <- read.table(dateiname, header = TRUE)

# Parameterauswahl: T = Temperatur ...
Parameter <- "T"
von <- 285
bis <- 290
# Achsen erstellen
a <- seq(from = 25, to = 5975, by = 50)
b <- c( 0.0150, table$Z[2:110])

Breite <- table$X
Hoehe <- table$Z
Temp <- table$T

# Schriftgröße 
font_size <- 0.7

# Matrix erstellen
mtrx1 <- acast(table, Breite~Hoehe, value.var=Parameter) 
# falsch rum mtrx2 <- acast(table28, Hoehe~Breite, value.var = "Parameter")

# Maßeinheiten
inMeter <- function(x) {
  paste0(x, " m")
}
inCelsius <- function(x) {
  paste0(x, " K")
}

# contour anzeigen
filled.contour(x= a, y=b, mtrx1,
               xlim = range(a, finite = TRUE),
               ylim = range(0, 10, finite = TRUE),
               zlim = range(von, bis, finite = TRUE), 
               plot.axes = { 
                 #axis(1, labels = inMeter(a), at = as.numeric(a), las =1 , cex.axis = 0.7) # x-Achse
                 #axis(2, labels = inMeter(b), at = as.numeric(b), las =1, cex.axis = 0.5) # y-Achse
                 axis(1, at = NULL, cex.axis = font_size) # x-Achse
                 axis(2, at = NULL, cex.axis = font_size) # y-Achse
               },
               key.axes = {
                 axis(4, at = pretty(Temp), labels = inCelsius(pretty(Temp))) 
               })

Code: Alles auswählen

pdf(...) # Datei anlegen/öffnen
... # Grafik erzeugen
dev.off() # Datei schließen

Code: Alles auswählen

help("png")