识别曲线中的点
Identifying points in a curve
我觉得这个问题很简单...
如何识别图形中的坐标?我绘制了一些数据,使用 unireg(uniReg 包)制作样条曲线,并希望从一个点提取数据。
library(uniReg)
P0mM <- read.table(text="
Time FeuM
0.04 138.8181818
7 1258.636364
14 1320.545455
21 2110.37037
28 13730.37037
35 1550.909091",header=TRUE)
z=seq(min(P0mM$Time),max(P0mM$Time),length=201)
uf=with(P0mM,unireg(Time,FeuM,g=5,sigma=1))
plot(FeuM~Time,P0mM,ylim=c(0,16000),ylab="Fe2+ uM", xlab="Time", main="0mM P")
lines(z,uf$unimod.func(z))
我能够找到曲线的最大 y 值(即 14444)
max((uf$unimod.func(z)))
我想确定在 x 轴上发生这种情况的位置。 (应该在 30 左右,但我想准确一点)。
你是怎么做到的?
谢谢!
似乎是 optimise
或 optimize
的情况(取决于您对英式或美式英语的熟悉程度):
optimise(uf$unimod.func, maximum=TRUE, interval=range(P0mM$Time))
#$maximum
#[1] 29.27168
#
#$objective
# [,1]
#[1,] 14444.85
我觉得这个问题很简单...
如何识别图形中的坐标?我绘制了一些数据,使用 unireg(uniReg 包)制作样条曲线,并希望从一个点提取数据。
library(uniReg)
P0mM <- read.table(text="
Time FeuM
0.04 138.8181818
7 1258.636364
14 1320.545455
21 2110.37037
28 13730.37037
35 1550.909091",header=TRUE)
z=seq(min(P0mM$Time),max(P0mM$Time),length=201)
uf=with(P0mM,unireg(Time,FeuM,g=5,sigma=1))
plot(FeuM~Time,P0mM,ylim=c(0,16000),ylab="Fe2+ uM", xlab="Time", main="0mM P")
lines(z,uf$unimod.func(z))
我能够找到曲线的最大 y 值(即 14444)
max((uf$unimod.func(z)))
我想确定在 x 轴上发生这种情况的位置。 (应该在 30 左右,但我想准确一点)。
你是怎么做到的?
谢谢!
似乎是 optimise
或 optimize
的情况(取决于您对英式或美式英语的熟悉程度):
optimise(uf$unimod.func, maximum=TRUE, interval=range(P0mM$Time))
#$maximum
#[1] 29.27168
#
#$objective
# [,1]
#[1,] 14444.85