fractional curve fitting of the function $y=a+bx^{alpha}$
$begingroup$
Assume I have a set of data $(x_{i},y_{i})$, $i=1,...,m$. How can we find the best values of the parameters $a$ and $b$ and $alpha$ such that the curve $y=a+bx^{alpha}$ best fits the data.
This question arises in many problems that use fractional calculus in modelling. I see it is almost always done numerically using a certain routine in Matlab or Maple. Can we do it analytically using the least squares minimization method ?
If $a=0$, we can linearize the model by taking the logarithm of both sides. The difficulty here is that, in practice, the parameter $a$ is unknown and needs to be determined.
Another question: Can Wolfram Alpha solve such a problem ??
Thanks a lot.
statistics regression least-squares weighted-least-squares
asked Nov 26 '18 at 3:43
MedoMedo
617214
$endgroup$
add a comment |
$begingroup$
Assume I have a set of data $(x_{i},y_{i})$, $i=1,...,m$. How can we find the best values of the parameters $a$ and $b$ and $alpha$ such that the curve $y=a+bx^{alpha}$ best fits the data.
This question arises in many problems that use fractional calculus in modelling. I see it is almost always done numerically using a certain routine in Matlab or Maple. Can we do it analytically using the least squares minimization method ?
If $a=0$, we can linearize the model by taking the logarithm of both sides. The difficulty here is that, in practice, the parameter $a$ is unknown and needs to be determined.
Another question: Can Wolfram Alpha solve such a problem ??
Thanks a lot.
statistics regression least-squares weighted-least-squares
asked Nov 26 '18 at 3:43
MedoMedo
617214
$endgroup$
add a comment |
$begingroup$
Assume I have a set of data $(x_{i},y_{i})$, $i=1,...,m$. How can we find the best values of the parameters $a$ and $b$ and $alpha$ such that the curve $y=a+bx^{alpha}$ best fits the data.
This question arises in many problems that use fractional calculus in modelling. I see it is almost always done numerically using a certain routine in Matlab or Maple. Can we do it analytically using the least squares minimization method ?
If $a=0$, we can linearize the model by taking the logarithm of both sides. The difficulty here is that, in practice, the parameter $a$ is unknown and needs to be determined.
Another question: Can Wolfram Alpha solve such a problem ??
Thanks a lot.
statistics regression least-squares weighted-least-squares
asked Nov 26 '18 at 3:43
MedoMedo
617214
$endgroup$
Assume I have a set of data $(x_{i},y_{i})$, $i=1,...,m$. How can we find the best values of the parameters $a$ and $b$ and $alpha$ such that the curve $y=a+bx^{alpha}$ best fits the data.
This question arises in many problems that use fractional calculus in modelling. I see it is almost always done numerically using a certain routine in Matlab or Maple. Can we do it analytically using the least squares minimization method ?
If $a=0$, we can linearize the model by taking the logarithm of both sides. The difficulty here is that, in practice, the parameter $a$ is unknown and needs to be determined.
Another question: Can Wolfram Alpha solve such a problem ??
Thanks a lot.
statistics regression least-squares weighted-least-squares
statistics regression least-squares weighted-least-squares
asked Nov 26 '18 at 3:43
MedoMedo
617214
asked Nov 26 '18 at 3:43
MedoMedo
617214
asked Nov 26 '18 at 3:43
MedoMedo
617214
asked Nov 26 '18 at 3:43
MedoMedo
617214
asked Nov 26 '18 at 3:43
MedoMedo
617214
617214
add a comment |
add a comment |
1 Answer
1
active
oldest
votes
$begingroup$
Based on $n$ data points $(x_i,y_i)$, you want to fit the model $$y=a+bx^{alpha}$$ which is nonlinear which means that you need estimates to start the nonlinear regression and the question is then : how to get these estimates ?
If you think about it, the model is nonlinear just because of $alpha$. So, suppose you give a value to $alpha$, define $t_i=x_i^{alpha}$ and you just face a linear model $y=a+b t$. So, for this arbitrary value of $alpha$, you van get $a(alpha$), $b(alpha$) and $SSQ(alpha)$.
So, an idea is to try different values of $alpha$ and, plotting $SSQ(alpha)$, to look where it goes more or less through a minimum. Let us call $alpha_*$ this point. To it correspond $a_*$ and $b_*$. These are your estimates for the nonlinear regression.
If you do not have such a tool, then contine the process zooming more and more around the minimum.
It this clear for you ? If not, I could elaborate using an example you could add to the post.
What you can also do, hoping that the errors are not too large is to take three points (say for example points $1$, $n$ and one somewhere in the middle (say point $m$) and write the three equations
$$y_1=a+bx_1^{alpha}qquad
y_m=a+bx_m^{alpha}qquad
y_n=a+bx_n^{alpha}$$ So
$$y_m-y_1=b (x_m^{alpha}-x_1^{alpha})qquad
y_n-y_m=b (x_n^{alpha}-x_m^{alpha})$$
$$frac {y_m-y_1 } {y_n-y_m }=frac {x_m^{alpha}-x_1^{alpha} } { x_n^{alpha}-x_m^{alpha}}$$ which is a single equation in $alpha$ that you can solve using a plot. This gives you $alpha_*$ and, going backward in he three equations, then $a_*$ and $b_*$.
Now, you are ready !
answered Nov 26 '18 at 5:05
Claude LeiboviciClaude Leibovici
120k1157132
$endgroup$
$begingroup$
Thank you for the insight.
$endgroup$
– Medo
Nov 26 '18 at 23:52
add a comment |
Your Answer
StackExchange.ifUsing("editor", function () {
return StackExchange.using("mathjaxEditing", function () {
StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
});
});
}, "mathjax-editing");
StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "69"
};
initTagRenderer("".split(" "), "".split(" "), channelOptions);
StackExchange.using("externalEditor", function() {
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled) {
StackExchange.using("snippets", function() {
createEditor();
});
}
else {
createEditor();
}
});
function createEditor() {
StackExchange.prepareEditor({
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader: {
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
},
noCode: true, onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});
}
});
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3013801%2ffractional-curve-fitting-of-the-function-y-abx-alpha%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
1 Answer
1
active
oldest
votes
1 Answer
1
active
oldest
votes
active
oldest
votes
active
oldest
votes
$begingroup$
Based on $n$ data points $(x_i,y_i)$, you want to fit the model $$y=a+bx^{alpha}$$ which is nonlinear which means that you need estimates to start the nonlinear regression and the question is then : how to get these estimates ?
If you think about it, the model is nonlinear just because of $alpha$. So, suppose you give a value to $alpha$, define $t_i=x_i^{alpha}$ and you just face a linear model $y=a+b t$. So, for this arbitrary value of $alpha$, you van get $a(alpha$), $b(alpha$) and $SSQ(alpha)$.
So, an idea is to try different values of $alpha$ and, plotting $SSQ(alpha)$, to look where it goes more or less through a minimum. Let us call $alpha_*$ this point. To it correspond $a_*$ and $b_*$. These are your estimates for the nonlinear regression.
If you do not have such a tool, then contine the process zooming more and more around the minimum.
It this clear for you ? If not, I could elaborate using an example you could add to the post.
What you can also do, hoping that the errors are not too large is to take three points (say for example points $1$, $n$ and one somewhere in the middle (say point $m$) and write the three equations
$$y_1=a+bx_1^{alpha}qquad
y_m=a+bx_m^{alpha}qquad
y_n=a+bx_n^{alpha}$$ So
$$y_m-y_1=b (x_m^{alpha}-x_1^{alpha})qquad
y_n-y_m=b (x_n^{alpha}-x_m^{alpha})$$
$$frac {y_m-y_1 } {y_n-y_m }=frac {x_m^{alpha}-x_1^{alpha} } { x_n^{alpha}-x_m^{alpha}}$$ which is a single equation in $alpha$ that you can solve using a plot. This gives you $alpha_*$ and, going backward in he three equations, then $a_*$ and $b_*$.
Now, you are ready !
answered Nov 26 '18 at 5:05
Claude LeiboviciClaude Leibovici
120k1157132
$endgroup$
$begingroup$
Thank you for the insight.
$endgroup$
– Medo
Nov 26 '18 at 23:52
add a comment |
$begingroup$
Based on $n$ data points $(x_i,y_i)$, you want to fit the model $$y=a+bx^{alpha}$$ which is nonlinear which means that you need estimates to start the nonlinear regression and the question is then : how to get these estimates ?
If you think about it, the model is nonlinear just because of $alpha$. So, suppose you give a value to $alpha$, define $t_i=x_i^{alpha}$ and you just face a linear model $y=a+b t$. So, for this arbitrary value of $alpha$, you van get $a(alpha$), $b(alpha$) and $SSQ(alpha)$.
So, an idea is to try different values of $alpha$ and, plotting $SSQ(alpha)$, to look where it goes more or less through a minimum. Let us call $alpha_*$ this point. To it correspond $a_*$ and $b_*$. These are your estimates for the nonlinear regression.
If you do not have such a tool, then contine the process zooming more and more around the minimum.
It this clear for you ? If not, I could elaborate using an example you could add to the post.
What you can also do, hoping that the errors are not too large is to take three points (say for example points $1$, $n$ and one somewhere in the middle (say point $m$) and write the three equations
$$y_1=a+bx_1^{alpha}qquad
y_m=a+bx_m^{alpha}qquad
y_n=a+bx_n^{alpha}$$ So
$$y_m-y_1=b (x_m^{alpha}-x_1^{alpha})qquad
y_n-y_m=b (x_n^{alpha}-x_m^{alpha})$$
$$frac {y_m-y_1 } {y_n-y_m }=frac {x_m^{alpha}-x_1^{alpha} } { x_n^{alpha}-x_m^{alpha}}$$ which is a single equation in $alpha$ that you can solve using a plot. This gives you $alpha_*$ and, going backward in he three equations, then $a_*$ and $b_*$.
Now, you are ready !
answered Nov 26 '18 at 5:05
Claude LeiboviciClaude Leibovici
120k1157132
$endgroup$
$begingroup$
Thank you for the insight.
$endgroup$
– Medo
Nov 26 '18 at 23:52
add a comment |
$begingroup$
Based on $n$ data points $(x_i,y_i)$, you want to fit the model $$y=a+bx^{alpha}$$ which is nonlinear which means that you need estimates to start the nonlinear regression and the question is then : how to get these estimates ?
If you think about it, the model is nonlinear just because of $alpha$. So, suppose you give a value to $alpha$, define $t_i=x_i^{alpha}$ and you just face a linear model $y=a+b t$. So, for this arbitrary value of $alpha$, you van get $a(alpha$), $b(alpha$) and $SSQ(alpha)$.
So, an idea is to try different values of $alpha$ and, plotting $SSQ(alpha)$, to look where it goes more or less through a minimum. Let us call $alpha_*$ this point. To it correspond $a_*$ and $b_*$. These are your estimates for the nonlinear regression.
If you do not have such a tool, then contine the process zooming more and more around the minimum.
It this clear for you ? If not, I could elaborate using an example you could add to the post.
What you can also do, hoping that the errors are not too large is to take three points (say for example points $1$, $n$ and one somewhere in the middle (say point $m$) and write the three equations
$$y_1=a+bx_1^{alpha}qquad
y_m=a+bx_m^{alpha}qquad
y_n=a+bx_n^{alpha}$$ So
$$y_m-y_1=b (x_m^{alpha}-x_1^{alpha})qquad
y_n-y_m=b (x_n^{alpha}-x_m^{alpha})$$
$$frac {y_m-y_1 } {y_n-y_m }=frac {x_m^{alpha}-x_1^{alpha} } { x_n^{alpha}-x_m^{alpha}}$$ which is a single equation in $alpha$ that you can solve using a plot. This gives you $alpha_*$ and, going backward in he three equations, then $a_*$ and $b_*$.
Now, you are ready !
answered Nov 26 '18 at 5:05
Claude LeiboviciClaude Leibovici
120k1157132
$endgroup$
Based on $n$ data points $(x_i,y_i)$, you want to fit the model $$y=a+bx^{alpha}$$ which is nonlinear which means that you need estimates to start the nonlinear regression and the question is then : how to get these estimates ?
If you think about it, the model is nonlinear just because of $alpha$. So, suppose you give a value to $alpha$, define $t_i=x_i^{alpha}$ and you just face a linear model $y=a+b t$. So, for this arbitrary value of $alpha$, you van get $a(alpha$), $b(alpha$) and $SSQ(alpha)$.
So, an idea is to try different values of $alpha$ and, plotting $SSQ(alpha)$, to look where it goes more or less through a minimum. Let us call $alpha_*$ this point. To it correspond $a_*$ and $b_*$. These are your estimates for the nonlinear regression.
If you do not have such a tool, then contine the process zooming more and more around the minimum.
It this clear for you ? If not, I could elaborate using an example you could add to the post.
What you can also do, hoping that the errors are not too large is to take three points (say for example points $1$, $n$ and one somewhere in the middle (say point $m$) and write the three equations
$$y_1=a+bx_1^{alpha}qquad
y_m=a+bx_m^{alpha}qquad
y_n=a+bx_n^{alpha}$$ So
$$y_m-y_1=b (x_m^{alpha}-x_1^{alpha})qquad
y_n-y_m=b (x_n^{alpha}-x_m^{alpha})$$
$$frac {y_m-y_1 } {y_n-y_m }=frac {x_m^{alpha}-x_1^{alpha} } { x_n^{alpha}-x_m^{alpha}}$$ which is a single equation in $alpha$ that you can solve using a plot. This gives you $alpha_*$ and, going backward in he three equations, then $a_*$ and $b_*$.
Now, you are ready !
answered Nov 26 '18 at 5:05
Claude LeiboviciClaude Leibovici
120k1157132
answered Nov 26 '18 at 5:05
Claude LeiboviciClaude Leibovici
120k1157132
answered Nov 26 '18 at 5:05
Claude LeiboviciClaude Leibovici
120k1157132
answered Nov 26 '18 at 5:05
Claude LeiboviciClaude Leibovici
120k1157132
120k1157132
$begingroup$
Thank you for the insight.
$endgroup$
– Medo
Nov 26 '18 at 23:52
add a comment |
$begingroup$
Thank you for the insight.
$endgroup$
– Medo
Nov 26 '18 at 23:52
$begingroup$
Thank you for the insight.
$endgroup$
– Medo
Nov 26 '18 at 23:52
$begingroup$
Thank you for the insight.
$endgroup$
– Medo
Nov 26 '18 at 23:52
add a comment |
Thanks for contributing an answer to Mathematics Stack Exchange!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
Use MathJax to format equations. MathJax reference.
To learn more, see our tips on writing great answers.
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3013801%2ffractional-curve-fitting-of-the-function-y-abx-alpha%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
