求助R语言回归分析 y=ksin(x/a) 如何利用lm()求解a的值

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

在做回归分析时 需要对参数进行二次处理 如sin（x/a）形式 lm()回归只能求解sinx前的系数 无法求解a的值 应该如何解决

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：求助R语言 回归分析 R语言

如果想要表达的公式为y=ksinx，则仅需要输入lm（y~sinx）即可，但是如果需要表达的公式为y=ksin(x/a)，存在两个系数k与a，其中k可以求出，a无法再公式中进行表达，遇到这种情况应该如何处理呢？

这从本质上来讲，就是一个拟合的问题。如果不知道用R来做，可以用其他软件来做吧。

nkunku 发表于 2022-4-11 08:14
这从本质上来讲，就是一个拟合的问题。如果不知道用R来做，可以用其他软件来做吧。

因为后续的分析都是基于R做的 现在只有拟合这一部分出现了问题 所以需要解决这个问题

可以百度一下sin(x)的泰勒展开，变为多项式再回归试试

这不是一个线性问题，lm没法搞定啊。
e = y - k * sin (a * x)
min (sum {e_i ^ 2})非线性函数，自己写个最优求解的目标函数吧，R中也有现成的求解器

上面有人提到用泰勒展开，如果a*x的绝对值较小（泰勒展开需要一定效率，绝对值越小越好呗），可以这么干，如果不是的，那还是按照上面的方式弄吧。

abelus 发表于 2022-4-12 11:25
这不是一个线性问题，lm没法搞定啊。
e = y - k * sin (a * x)
min (sum {e_i ^ 2})非线性函数，自己写个 ...

可以具体说一下么？

如果只是周期数据的拟合，可以考虑用TSA包的harmonic()函数，谐波函数

llb_321 发表于 2022-4-13 10:07
如果只是周期数据的拟合，可以考虑用TSA包的harmonic()函数，谐波函数

我的思路是先对数据进行拟合，基于可视化对不同模型进行对比分析，再对模型进行改进。
周期数据的拟合是为了模型的可视化做准备