R中的正弦曲线拟合答案

【问题标题】：Sinusoidal curve fitting in RR中的正弦曲线拟合
【发布时间】：2021-09-01 16:32:38
【问题描述】：

我正在尝试使用以下代码将正弦曲线拟合到数据中。我不知道为什么它变得时髦，线条似乎从 17 到 9 的观察然后从 9 到 20 再回到 12。我还希望正弦曲线延伸到整个 24 小时。如果你能在 ggplot2 中做到这一点，那就太棒了。谢谢！

Test_case <- structure(list(TEST_RESULT = c(1.31, 65.44, 8.59, 9.91), Time = c(17, 
9, 20, 12)), row.names = c(NA, -4L), class = c("tbl_df", "tbl", 
"data.frame"))

Time <- test_case$Time
test_result <- test_case$TEST_RESULT
  
xc<-cos(2*pi*Time/24)
xs<-sin(2*pi*Time/24)
fit.lm <- lm(test_result ~ xc+xs)

fit <- fitted(fit.lm)  

# find predictions for original time series
pred <- predict(fit.lm, newdata=data.frame(Time=Time))    

plot(test_result ~ Time, data= test_case, xlim=c(1, 24), ylim=c(0,500))
lines(Time, pred, col="blue")

【问题讨论】：

标签： r curve

【解决方案1】：

确保时间变量出现在公式中并且没有被 xc 和 xs 变量隐藏。还要确保预测中使用的时间是按顺序排列的，即 0:23 是按升序排列的。另一方面，输入数据不需要按顺序排列，因为我们只是将其绘制为点，因此绘制输入数据的顺序无关紧要。

因此从我们的问题中使用test_case：

fit.lm <- lm(test_result ~ cos(2*pi*Time/24) + sin(2*pi*Time/24), test_case)
 
Time2 <- 0:23
pred2 <- predict(fit.lm, newdata = list(Time = Time2))    
 
plot(test_result ~ Time, data = test_case, xlim = c(1, 24), ylim = c(0, 500))
lines(pred2 ~ Time2, col = "blue")

【讨论】：

这看起来不错，你能把正弦曲线也延长到整个 24 小时吗？
我得到的结果与您的图像不同。曲线倒过来的样子。
我确实多次更改了答案。请您在新的 R 实例上再试一次。首先复制并粘贴问题中的Test_case，然后将代码复制并粘贴到答案中。
我有一个问题，当我运行 summary(fit.lm) 时，截距为 43.2，但在图中明显高于 50。知道为什么会有差异吗？
只有 4 个点和 3 个参数，所以我认为您不能期望太多。建议您模拟一些点数大得多的数据，并将结果与用于生成数据的模型进行比较。