如何在R中生成带线性分量的余弦模型？

Question

您好StackOverflow社区， 我不是程序员，而是自学成才的R用户，用于统计和数据可视化。这是我的第一个问题，因为我总是发现其他成员的帖子回答了我的问题，所以谢谢你，请原谅我的问题中的任何礼仪失误。

我正在尝试在R中创建一个带有线性分量的余弦模型。

这样，= {2（（ - ∅）/τ）} ++。其中A是正弦波的幅度，τ是设定为24小时的周期，t是时间，Ø是正弦曲线的acrophase，C和D分别是线性分量的斜率和y轴截距。我查看过cosinor.lm {cosinor}和cosinor {psych}的文档，但看不到这个功能。

是否可以将线性分量添加到余弦模型中，还是可以输出另一个可以输出幅度，acrophase，斜率和y截距的显着性和值的模型？

以下是我的示例数据：

library(cosinor)
library(ggplot2)

#data frame
time<- c(2, 6, 10, 18, 22, 26, 30, 34, 38)
resp<- c(2.54, 0.13, -0.38, -0.57,  0.11, -0.20, -0.26, -0.62, -0.73)
df<-data.frame(time,resp)

# data fits both linear and cosinor models (but not very well)
cos.model <- cosinor.lm(resp ~ time(time), period=24, data=df)
summary(cos.model)
l.model <- lm(resp ~ time, data=df)
summary(l.model)

#plot with a loess smoother looks like a combination of the models 
#plot to see fit

ggplot(df, aes(x = time, y = resp))+
  geom_point()+
  geom_smooth(method = "lm", se = FALSE)

ggplot.cosinor.lm(cos.model)+
  geom_point(data=df, aes(y=resp, x=time))

ggplot(df, aes(x = time, y = resp))+
  geom_point()+
  geom_smooth(method ="loess", se = FALSE)

如果有人能帮我创建一个包含余弦和线性分量的模型，我将非常感激！

提前致谢。