Date: Wed, 28 Dec 2011 07:15:19 -0800 David Marso "SPSSX(r) Discussion" David Marso Re: Fitting logistic curve to growth data <1325070409338-5105041.post@n5.nabble.com> text/plain; charset=UTF-8

See the NLR (Non Linear Regression) Procedure. Following is an example of simulated data and the set up and results of NLR. Note that the parameter estimates used to generate the data are recovered quite nicely! You will wish to consult the docs on how to set up and run NLR from the dialog boxes. HTH, David --- INPUT PROGRAM. COMPUTE #K=1.5 . COMPUTE #A=4.0. COMPUTE #i= 1.7 . LOOP ID=1 TO 1000. COMPUTE d= Uniform(1)+4. COMPUTE Y=#A/(1+EXP(-#K * (d-#i)))+ UNIFORM(.01). END CASE. END LOOP. END FILE. END INPUT PROGRAM. EXE.

GRAPH /SCATTERPLOT(BIVAR)=d WITH y /MISSING=LISTWISE .

* NonLinear Regression. MODEL PROGRAM A=3 K=1 I=1 . COMPUTE PRED_ = A/(1+EXP(-K * (d-I))). NLR y /PRED PRED_ /CRITERIA SSCONVERGENCE 1E-8 PCON 1E-8 .

All the derivatives will be calculated numerically.

Iteration Residual SS A K I

1 1064.692791 3.00000000 1.00000000 1.00000000 1.1 1.348304567 3.99063649 1.49442721 1.88384875 2 1.348304567 3.99063649 1.49442721 1.88384875 2.1 .0146986193 4.00583423 1.47908955 1.69332655 3 .0146986193 4.00583423 1.47908955 1.69332655 3.1 .0080481338 4.00572338 1.47960600 1.66912002 4 .0080481338 4.00572338 1.47960600 1.66912002 4.1 .0080459976 4.00572123 1.47971278 1.66884117 5 .0080459976 4.00572123 1.47971278 1.66884117 5.1 .0080459976 4.00572119 1.47971372 1.66884238

Run stopped after 10 model evaluations and 5 derivative evaluations. Iterations have been stopped because the relative reduction between successive residual sums of squares is at most SSCON = 1.000E-08

Nonlinear Regression Summary Statistics Dependent Variable Y

Source DF Sum of Squares Mean Square

Regression 3 15524.41480 5174.80493 Residual 997 8.045998E-03 8.070208E-06 Uncorrected Total 1000 15524.42285

(Corrected Total) 999 .74860

R squared = 1 - Residual SS / Corrected SS = .98925

Asymptotic 95 % Asymptotic Confidence Interval Parameter Estimate Std. Error Lower Upper

A 4.005721193 .001171063 4.003423162 4.008019224 K 1.479713721 .026491205 1.427728805 1.531698638 I 1.668842378 .036748380 1.596729333 1.740955424

Asymptotic Correlation Matrix of the Parameter Estimates

A K I

A 1.0000 -.9791 -.9656 K -.9791 1.0000 .9981 I -.9656 .9981 1.0000

DP_Sydney wrote > > Hi all, > > I need to fit a curve to data on mass gain, which follows a sigmoidal > pattern and is best described using a logistic equation (based on similar > published data, with which I wish to compare my data) of the form: > > Y = A/(1 + e^-K(d-i)) > where Y = mass > A = asymptotic mass > K = growth rate > d = age (time) > i = inflection point > (note ^ indicates to the power) > > The data I have are body mass recorded daily for several individuals > (nestling birds), which will be analysed separately. > > I need to output the parameters (A, K, i) rather than just produce a > curve. I've searched the archives and Google but haven't found > instructions for how to do this (I have instructions for running this in > GENSTAT but of course do not have GENSTAT). > > Any help would be great...but beware I am a drop-down menus and buttons > person when it comes to software :-) > > Regards, > Dean >

-- View this message in context: http://spssx-discussion.1045642.n5.nabble.com/Fitting-logistic-curve-to-growth-data-tp5105041p5105450.html Sent from the SPSSX Discussion mailing list archive at Nabble.com.

===================== To manage your subscription to SPSSX-L, send a message to LISTSERV@LISTSERV.UGA.EDU (not to SPSSX-L), with no body text except the command. To leave the list, send the command SIGNOFF SPSSX-L For a list of commands to manage subscriptions, send the command INFO REFCARD

Back to: Top of message | Previous page | Main SPSSX-L page