Mathematica 9 is now available
Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2012

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Fit data

  • To: mathgroup at smc.vnet.net
  • Subject: [mg124710] Re: Fit data
  • From: Darren Glosemeyer <darreng at wolfram.com>
  • Date: Wed, 1 Feb 2012 03:52:55 -0500 (EST)
  • Delivered-to: l-mathgroup@mail-archive0.wolfram.com
  • References: <201201301009.FAA17575@smc.vnet.net>

On 1/30/2012 4:09 AM, Miguel Camino wrote:
> Hi, I'm new in this forum and with Mathematica..I've got 2 questions;Is too much?
> 1. I've to regress some data with a formula that I've implemented.This is an Example of my data and this is my formula:sigma=a*pc^(0.618)*tc*(1+omega)^c*(1-tr)^dtr,                omega,         tc,                   pc             sigma0.703784817, 0.310237474, 361.0989249, 34.97280104, 12.410.710106871, 0.310283818, 361.1126025, 34.94104003, 12.140.715649046, 0.310330225, 361.1262985, 34.90923611, 11.770.720912444, 0.310376694, 361.1400131, 34.87738919, 11.520.726476477, 0.31043254, 361.1564951, 34.83911601, 11.30.731293672, 0.310479148, 361.1702505, 34.80717419, 11I used this formulafit = NonlinearModelFit[data,   a*pc^(0.618)*tc^b*(1 + omega)^c*(1 - Tr)^d, {a, b, c, d}, {Tr,    omega, tc, pc}]and I try to calculate best {a,b,c,d}.
> 2. I don't know if 0.618 is the best value. Is thes a way to find the best number that Fit my data?Thanks and sorry for two questions!M.
>
>

For question 1, your data will need to be a numeric matrix with 5 
columns. The first 4 columns will be the values of tr, omega, tc, and 
pc, and the last will be the sigma values. Note that Tr is a built-in 
Mathematica function, so I would recommend using tr or something else 
that doesn't start with a capital letter. You may also need to give 
starting values for the parameters a, b, c, and d to get a good fit; 
this is just the nature of nonlinear optimization.

For question2, you could replace 0.618 in the formula by some symbol 
(say p) and add that symbol to the list of parameters to estimate. If 
.618 seems like a reasonable value for your data, you could use that as 
a starting value for p.

Darren Glosemeyer
Wolfram Research



  • Prev by Date: Re: Trying to close some loose plotting
  • Next by Date: Re: large lists with variable and if condition (vrs. 8.0.4)
  • Previous by thread: Re: Trying to close some loose plotting
  • Next by thread: Re: large lists with variable and if condition (vrs. 8.0.4)