Segmented regression or broken stick regression

10 次查看（过去 30 天）
Wendy Cameron
Wendy Cameron 2018-6-8
评论： Wendy Cameron 2018-6-8
An example of my data is as follows: xdata=[1150 1200 1263.3 1333.3 1429 1497.3 1567.8 1608.3 1643.5 1668.8]; ydata=[6.2 7.2 8.6 10.4 11.2 12.9 13.3 13.3 13.5 13.55];
A scatter plot shows there are arguably 2 linear sections to the graph, the first steeper than the second. I would like to find out the slope of the first linear section (always positive) and the break point (can vary). the slope the second section may be zero, positive or negative (depending on the weather).
I tried the on line suggested solution but get an error message as shown. xdata=[1150 1200 1263.3 1333.3 1429 1497.3 1567.8 1608.3 1643.5 1668.8]; ydata=[6.2 7.2 8.6 10.4 11.2 12.9 13.3 13.3 13.5 13.55]; >> coefs = @(xbreak) [ones(numel(xdata),1),min(xdata,xbreak)-xbreak]\ydata; >> errfun = @(xbreak) norm([ones(numel(xdata),1),min(xdata,xbreak)-xbreak]*coefs(xbreak) - ydata); >> xbreak = fminbnd(errfun,0,1) Error using horzcat Dimensions of matrices being concatenated are not consistent.
Error in @(xbreak)norm([ones(numel(xdata),1),min(xdata,xbreak)-xbreak]*coefs(xbreak)-ydata)
Error in fminbnd (line 228) x= xf; fx = funfcn(x,varargin{:});
Hoping someone can help me?
Kind regards, Wendy

请先登录，再回答此问题。

采纳的回答

Image Analyst
Image Analyst 2018-6-8
I've already given you an answer here in this answer.
I'm attaching it again where I'm using your data in the m-file above.
  1 个评论
Wendy Cameron
Wendy Cameron 2018-6-8
Thank you. This seems to work well and thank you particularly for the many explanations of what each step was doing which really meant I learnt a lot from this as well as solving my problem.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Linear Regression 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by