推 wwfc: 謝謝c大,我看太快了,沒有多try幾次 08/03 09:47
※ 引述《wwfc (月老工讀生)》之銘言:
: [軟體程式類別]:
: R
: [程式問題]:
: 建立折線模型
: [軟體熟悉度]:
: 高(1年以上)
: [問題敘述]:
: 最近接了一個分析數據的case,
: 其中一個項目是要針對資料建立折線模型,
: 我虛擬了一筆資料並且用小畫家把想要的折線畫上去。
: http://imgur.com/Ri844kS
: 目前是利用segments畫出客戶要的,
: 作法:目測不同線段區間,再用lm, predict.lm,去把節點算出來,
: 最後用segments畫線,有時候需要一點人工微調,把線接在一起。
: 可是這種作法太曠日費時,
: 請問有沒有package可以直接畫折線模型,
: google "R 折線圖",可是出來的結果不是我要的。
: 請有經驗的版友分享一下,謝謝。
R code:
x <- rnorm(50)
y <- sin(x*pi)
plot(x, y)
lines(spline(x, y, 7, "natural"), col = 2)
lines(spline(x, y, 8, "natural"), col = 3)
lines(lowess(x, y, 2/3), col = 4)
lines(lowess(x, y, 0.5), col = 5)
legend("topright", c("spline (n = 7)","spline (n = 8)",
"lowess (f = 2/3)", "lowess (f = 0.5)"), col = 2:5, lty = 1)
http://i.imgur.com/WCipLSJ.png
lowess是取附近點以距離作為weight做回歸分析,得到特定x的估計值
spline是根據你給的切割點個數去最小化每一個區間的MSE,幾次可微通常是可以設定,
這裡給的是0次可微,只有連續
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 123.205.27.107
※ 文章網址: https://www.ptt.cc/bbs/Statistics/M.1438504718.A.923.html
※ 編輯: celestialgod (111.82.229.42), 08/02/2015 16:55:44