Skip to main content
V
ViktorKorobov1-Visitor
1-Visitor
October 7, 2017
Question

Parameter Estimation Problem

  • October 7, 2017
  • 3 replies
  • 3591 views

Dear Experts,

I want to find confidence intervals for estimated parameters using LeastSquaresFit. Without units, there is no problem. When using units, an error occurs if the fitting function defined as a program block. How to solve this problem?

With units:

ParEst_withUnits.png

Without units:

ParEst_withoutUnits.png

 

3 replies

L
LucMeekes23-Emerald IV
23-Emerald IV
October 7, 2017

With this kind of problems it's usually wise to check the labels of variables.

I found that most k2's are either labelled as 'variable' or as an 'automatic' 'variable', except for the one indicated:

LM_20171007_ParEst.png

Success!
Luc

V
ViktorKorobov1-VisitorAuthor
1-Visitor
October 7, 2017
@LucMeekes wrote:

With this kind of problems it's usually wise to check the labels of variables.

I found that most k2's are either labelled as 'variable' or as an 'automatic' 'variable', except for the one indicated:

 

 

Thanks for the reply, Luc. Unfortunately, it does not solve the problem.

T
21-Topaz II
terryhendicott21-Topaz II
21-Topaz II
October 9, 2017

Hi,

 

The values used in the second phase without units are k1=1 and k2=2 as start guesses.

The values used in the first phase with units are k1=0.278 and k2=5.556*10^-4 as start guesses.  see the sheet.

That is why 0.046 and 1.278x10^-4 are found by Minerr. 

 

The values used in the second phase without units are [1 2 3 4 6] T for t and [0.135 0.217 etc]T for Ce.

The values used in the first phase with units are [3600 7200 etc]T for t and [135 217 etc]T for Ce. see the sheet

 

These two the first and second phases thus are not the same numerical problem.  I only have express version of Mathcad Prime 4.0 so it will not do the couple of lines of programming to show you how to correct your sheet.  Essentially with units you must force the value by dividing by the units to use to get the value you want. 

T
21-Topaz II
terryhendicott21-Topaz II
21-Topaz II
October 9, 2017

Hi,

 

The values used in the second phase without units are k1=1 and k2=2 as start guesses.

The values used in the first phase with units are k1=0.278 and k2=5.556*10^-4 as start guesses.  see the sheet.

That is why 0.046 and 1.278x10^-4 are found by Minerr. 

 

The values used in the second phase without units are [1 2 3 4 6] T for t and [0.135 0.217 etc]T for Ce.

The values used in the first phase with units are [3600 7200 etc]T for t and [135 217 etc]T for Ce. see the sheet

 

These two the first and second phases thus are not the same numerical problem.  I only have express version of Mathcad Prime 4.0 so it will not do the couple of lines of programming to show you how to correct your sheet.  Essentially with units you must force the value by dividing by the units to use to get the value you want. 

V
ViktorKorobov1-VisitorAuthor
1-Visitor
October 9, 2017

Thank you, Terry.

The problem is not in which units the guess values are given. Functions Minerr and LeastSquaresFit return the result in the base units of the system SI. The problem is that the construction of Sol1 does not work!
Unfortunately I can not open files Prime 4.0.

 

T
21-Topaz II
terryhendicott21-Topaz II
21-Topaz II
October 9, 2017

Hi,

Captured graphic

Terms & ConditionsCookie settingsAccessibility statement