作者celestialgod (攸蓝)
看板Statistics
标题Re: [程式] R 折线模型
时间Sun Aug 2 16:38:35 2015
※ 引述《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://webptt.com/cn.aspx?n=bbs/Statistics/M.1438504718.A.923.html
※ 编辑: celestialgod (111.82.229.42), 08/02/2015 16:55:44
1F:推 wwfc: 谢谢c大,我看太快了,没有多try几次 08/03 09:47