r中的对数线性拟合绘图

Question

我有一个数据帧对应每个变量的频率 （例如：变量1出现1984次，变量2出现974 ......）

dff<-data.frame(Var1=c(1 ,2 ,3 ,4 ,5 ,6 ,7 ,8 ,9 ,10 ,11, 12 ,13 ,14 ,15 ,16 ,27, 30 ,35 ,36 ,38 ,39 ,40 ,41), Freq=c(1984,  974,  464 , 251 , 127 ,  83 ,  45 ,  26  , 16  , 12 ,   9   , 5 ,   5  ,  2    ,3  ,  1   , 1  ,  1 ,   1   , 2  ,  5,4,2,1))

plot(dff$Var1,log(dff$Freq))

日志情节

正如我们在图片中看到的那样，有一个线性回归，

我只是想找到一种方法，它只能适合线性拟合与x轴相交的回归部分。

因为，我需要提取线性拟合中不适合的其他点（dff$Var1=20之后的点），即平均值不是噪点

Answer 1

这个怎么样：

plot(dff$Var1, log(dff$Freq))
lr <- lm(log(Freq) ~ Var1, data = dff[dff$Var1 < 20, ])
abline(lr)

截止点是20.但你可以根据你的行为改变它。

如果你想计算

  

其中线性拟合与x轴相交。

获取系数：

coef(lr)
(Intercept)        Var1 
  7.4636699  -0.4741615 

并解决等式7.4636699 + Var1*(-0.4741615) = 0.