我是R的新人
我正在完成我的硕士论文的一部分,并希望在一段时间内制作一些例如甲烷生产的昼夜节。
现在我想在同一时间看到它的时间变化及其与另一个因素的相关性。然后我有两个问题。
第一:
如何定义xlim
和ylim
增加2小时。它有自己的默认值,当我举例说:
xlim = c(0, 23)
然后它从0开始,并在5小时内上升。我想让它在2小时内上升。
第二:
如何在同一时间段内放置另一个可能与我的第一个变量相关的变量。假设23小时内的甲烷产量可能与氧气消耗有关,就像一个例子。如何将氧气和甲烷放在同一轴(y)与时间(x)之间?
如果你能帮助我,我将非常感激。
种类,法哈德
at
函数调用中使用labels
和axis
参数来自定义标签和刻度位置。axis
的side = 4
函数在图形右侧创建自定义y轴。请参阅以下代码,说明上述要点:
set.seed(123)
x <- 0:23
df<- data.frame(
x,
ch4 = 1000 - x ^ 2,
o2 = 2000 - 2 * (x - 10) ^ 2
)
par(mar = c(5, 5, 2, 5))
with(df, plot(x, ch4,
type = "l", col = "red3",
ylab = "CH4 emission",
lwd = 3,
xlim = c(0, 23),
xlab = "",
xaxt = "n"))
axis(1, at = seq(0, 23, 2), labels = seq(0, 23, 2))
par(new = TRUE)
with(df, plot(x, o2,
pch = 16, axes = FALSE,
xlab = NA, ylab = NA, cex = 1.2))
axis(side = 4)
mtext(side = 4, line = 3, "O2 consumption")
legend("topright",
legend = c("O2", "CH4"),
lty = c(1, 0),
lwd = c(3, NA),
pch = c(NA, 16),
col = c("red3", "black"))