2014-02-18 41 views
3

我有timeserie:- [R找到兩條線之間的角度,具有斜率和截距係數時

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) 

我該如何管理?

這是我的圖表:Plot with two redlines

+1

提示:您可以忽略攔截。 –

+1

'angle = atan(abs((m2-m1)/(1 + m1 * m2)))' – Nishanth

回答

1

由於您的第二斜率是0,這條線和所述一個與21.27031斜率之間的角度的計算非常簡單:

atan(21.27031) * 180/pi 

# [1] 87.30828 

角度線與x軸之間的角度約爲87°。

+0

謝謝你的回答。我已經計算過了。 不過看看圖片: http://gyazo.com/85aa9963769fb9f60a4545a328162ed2 有俊男是角87度 – Marta

+3

@ user3227114沒有x軸標籤,你不能看到從圖形的角度。請注意,您可以更改繪圖的高度和寬度,並且儘管角度始終相同,但角度可能會有所不同。 –

+0

非常感謝! 但我仍然有問題。也許你可以幫我解決它? 我還有一個時間序列 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 – Marta