2016-01-26 25 views
1

我已經開始使用R中的樣條線了。我想創建一個樣條線,它被迫通過數據中的數據點。下面是一些示例代碼:強制通過數據點的樣條線

x <- 1:12 
y <- c(2172205, 5332408, 10453117, 16038564, 21109922, 
     27966857, 33563819, 39948003, 42629618, 45226902, 
     46172149, 46805245) 
z <- spline(x, y, 36) 

> plot(z) 
> print(z) 
$x 
[1] 1.000000 1.314286 1.628571 1.942857 2.257143 2.571429 2.885714 3.200000 3.514286 3.828571 4.142857 
[12] 4.457143 4.771429 5.085714 5.400000 5.714286 6.028571 6.342857 6.657143 6.971429 7.285714 7.600000 
[23] 7.914286 8.228571 8.542857 8.857143 9.171429 9.485714 9.800000 10.114286 10.428571 10.742857 11.057143 
[34] 11.371429 11.685714 12.000000 

$y 
[1] 2172205 2869024 3855824 5086171 6514487 8101175 9809181 11599445 13406684 15147279 16741530 18233441 
[13] 19806294 21649308 23806067 26054346 28143697 29930437 31596698 33386543 35470344 37630794 39523543 40845316 
[25] 41679456 42321140 43055048 43928388 44777465 45433385 45831051 46053800 46196328 46338922 46528029 46805245 

我想什麼是如下所示,但是通過強迫每一個人x和y分量最初給出的行擁有在花鍵的各個點的值。因此,5332408應顯示在y輸出中,2應顯示在x中。

我意識到這可能會給我一個更糟,更不平滑的曲線,但對我的目的來說沒關係。任何幫助非常感謝,並提前感謝您!

回答

2

使用splinefun代替spline看到你實際上得到相交的曲線與點:

zf <- splinefun(x, y) 

print(list(seq(36), y2=zf(seq(36))) ) 
#---------- 
[[1]] 
[1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 
[25] 25 26 27 28 29 30 31 32 33 34 35 36 

$y2 
[1] 2172205 5332408 10453117 16038564 21109922 27966857 
[7] 33563819 39948003 42629618 45226902 46172149 46805245 
[13] 48751454 53350661 61942754 75867617 96465137 125075200 
[19] 163037692 211692499 272379506 346438600 435209667 540032593 
[25] 662247264 803193565 964211384 1146640605 1351821115 1581092800 
[31] 1835795545 2117269238 2426853763 2765889007 3135714856 3537671196 

y %in% zf(seq(36)) 
[1] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE 
+0

真棒,我就知道有事情'splinefun',但並沒有將其連接。謝謝。 – giraffehere

相關問題