【问题标题】:R find angle between two lines, when have slope and intercept coefficientsR找到两条线之间的角度,当有斜率和截距系数时
【发布时间】:2014-02-18 07:25:55
【问题描述】:

我有时间:

x

4557  9940  9855  9894 10142  9501  9532  9229  9169  9214  9347  9176  8951  9344  9873  9970
9139  9420  9476  9205  9271  8632  8730  9336  9150  9601 10012  9841  9951  9222  8799  9316
10087  9677  9154  9019 10549  9629  9131  9560 10246 10982 11748  9054  8690  9923

tt<-1:length(x)

plot(x, xaxt = "n", type = "l", xlab = NULL, ylab = NULL, col = "royalblue2", lwd = 2.3)

我要添加线性趋势:

fit <- lm(x ~ tt)
co <- coef(fit)
co
    (Intercept)             tt 

    8940.23478           21.27031 

然后我需要找到两条红线之间的角度:

abline(8940.23478, 21.27031, col = "red", lwd = 2)
abline(8940.23478, 0, col = "red", lwd = 2)

我该如何管理?

这是我的图表:

【问题讨论】:

  • 提示:你可以忽略拦截。
  • angle = atan(abs((m2-m1)/(1+m1*m2)))

标签: r angle linear-regression intercept


【解决方案1】:

因为你的第二个斜率是0,所以计算这条线与斜率为21.27031的线之间的角度很容易:

atan(21.27031) * 180 / pi

# [1] 87.30828

直线与x轴的夹角约为87°。

【讨论】:

  • 感谢您的回答。我已经计算过了。但是看图:gyazo.com/85aa9963769fb9f60a4545a328162ed2不可能有87度角
  • @user3227114 如果没有 x 轴标签,您将无法从图形中看到角度。请注意,您可以更改绘图的高度和宽度,尽管角度始终相同,但角度可能会有所不同。
  • 非常感谢!但我仍然有一个问题。也许你可以帮我解决它?我还有一个时间序列X 1.133421 1.309256 1.322476 1.351425 1.319168 1.334070 1.323227 1.328096 1.321736 1.340243 1.348989 1.325523 1.321640 1.360659 1.366150 1.344534 1.327169 1.326433 1.372520 1.395546 1.365872 1.335496 1.380183 1.425236 1.330601 1.334444 1.314822 1.395691 1.393629 1.356322 1.285601 1.333727 1.358382 1.356309 1.339961 1.307129 1.331690 1.301070 1.290220 1.302197 1.497267 1.531415 1.549200 1.393307 1.324051 1.323692
  • @SvenHohenstein ,图表看起来相同,但斜率为 0.001841031。 gyazo.com/ae2de73dcec3d95941b7c63746b2c68d这是为什么呢?
  • @user3227114 使用这些命令绘制具有相同 x 轴和 y 轴限制的图:fit &lt;- lm(x ~ seq_along(x)); co &lt;- coef(fit); plot(x, type = "l", xlim = c(0, length(x)), ylim = c(0, length(x))); abline(co); abline(c(co[1], 0))。现在,您可以看到斜率非常小。
猜你喜欢
  • 1970-01-01
  • 2022-11-07
  • 2011-03-22
  • 2014-08-26
  • 1970-01-01
  • 1970-01-01
  • 2013-02-13
  • 1970-01-01
  • 2018-07-04
相关资源
最近更新 更多