Mathematica 9 is now available
Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
2004
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*December
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 2004

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

Search the Archive

Re: Fitting question

  • To: mathgroup at smc.vnet.net
  • Subject: [mg51605] Re: Fitting question
  • From: D Herring <dherring at at.uiuc.dot.edu>
  • Date: Wed, 27 Oct 2004 01:54:08 -0400 (EDT)
  • References: <cldov3$4a8$1@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com

(* How's this? *)

(* Given data *)
runtotpergen = {1, 23, 136, 568, 1735, 4382, 9099, 16384, 25993, 
36694, 	  47752, 58044, 66761, 73534, 78429, 81753, 83750, 84873, 85482, 
85779,
   85937, 86025, 86069, 86081, 86087};

totpergen = {1, 22, 113, 432, 1167, 2647, 4717, 7285, 9609, 10701,
   11058, 10292, 8717, 6773, 4895, 3324, 1997, 1123, 609, 297, 158, 88,
   44, 12, 6};

ListPlot[totpergen];

(* Do the fit *)
<< Statistics`ContinuousDistributions`
soln = FindFit[totpergen,
     a*PDF[NormalDistribution[b, c], x], {a, {b, 12}, c}, x]
Plot[a*PDF[NormalDistribution[b, c], x] /. soln, {x, 0, 25}];

(* Check the fit error *)
err = Table[
       totpergen[[i]] - (a*PDF[NormalDistribution[b, c], x] /. soln) /.
         x -> i, {i, 25}];
ListPlot[err];

relerr = err/totpergen;
ListPlot[relerr];


(* Hope that helps.  -- Daniel *)


  • Prev by Date: Bug in Graphics Text under Linux
  • Next by Date: Re: Re: Step by step answer ?
  • Previous by thread: Fitting question
  • Next by thread: Re: Fitting question