gpt4 book ai didi

r - 如何使用非线性函数拟合数据并绘制数据并使用 ggplot() 拟合

转载 作者:行者123 更新时间:2023-12-04 19:09:14 30 4
gpt4 key购买 nike

测量显示一个信号,其形式类似于具有偏移量和因子的平方根函数。如何找到系数并在一个图中绘制原始数据和拟合曲线?

require(ggplot2)
require(nlmrt) # may be this will help later..?

# generate simulated measurement data

time <- seq(-10,20,0.2)

signal <- sqrt(time + 2) # generate sqrt signal; includes some NA
signal[is.na(signal)] <- 0 # set all NA to zero
signal <- signal + rnorm(length(time)) * 0.1 # add noise

df <- data.frame(x=time, y=signal)

# find coefficiants for y ~ b * sqrt(x - a)
# no idea how...

# plot raw data and fitted curve in one ggplot diagram

ggplot()+
geom_point(data=df, aes(x=x, y=y))

enter image description here

最佳答案

如果您知道切割点在哪里并且切割点之前的值为零:

sfun <- function(x,brk,a=1) {
ifelse(x<brk,0,suppressWarnings(a*sqrt(x-brk)))
}

( suppressWarnings() 存在是因为 ifelsex 的所有值都计算 if 和 else 情况,我们不希望出现关于取负数平方根的警告)

测试(未显示):
curve(sfun(x,1,1),from=0,to=10) ## test (not shown)

模拟一些数据:
x <- seq(0,10,length=101)
set.seed(1)
y <- rnorm(length(x),sfun(x,1,1),sd=0.25)
DF <- data.frame(x,y)

因为我们只需要弄清楚平方根函数是如何缩放的,我们可以通过原点回归来做到这一点(如果你想让切割点下方的值不为零,则去掉 -1):
library("ggplot2")
theme_set(theme_bw())
ggplot(DF,aes(x,y))+geom_point()+
geom_smooth(method="lm",
formula=y~sfun(x,brk=1)-1)
ggsave("truncsqrt.png")

enter image description here

关于r - 如何使用非线性函数拟合数据并绘制数据并使用 ggplot() 拟合,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/17114009/

30 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com