Re: Dose-response association

Posted by Hosam on
URL: http://spssx-discussion.165.s1.nabble.com/Dose-response-association-tp5740322p5740343.html

*Generate Data.
set seed 98765432.
new file.
input program.
loop ID = 1 to 100.
     compute ModVig_minsweek = rv.normal(0,1).
     compute b0 = 0.5.
     compute b1 = 1.0.
     compute b2 = -1.5.
     compute b3 = 2.0.
     compute error = rv.normal(0,1).
     compute HRQoL = b0 + b1*ModVig_minsweek + b2*ModVig_minsweek**2 +
b3*ModVig_minsweek**3 + error.
     end case.
   end loop.
end file.
end input program.
execute.

delete variables b0 b1 b2 b3 error.

GRAPH
  /SCATTERPLOT(BIVAR)=ModVig_minsweek WITH HRQoL
  /MISSING=LISTWISE.

compute ModVig_minsweek_2nd = ModVig_minsweek**2.
compute ModVig_minsweek_3rd = ModVig_minsweek**3.
execute.

REGRESSION
  /MISSING LISTWISE
  /STATISTICS COEFF OUTS R ANOVA
  /CRITERIA=PIN(.05) POUT(.10)
  /NOORIGIN
  /DEPENDENT HRQoL
  /METHOD=ENTER ModVig_minsweek ModVig_minsweek_2nd ModVig_minsweek_3rd.

* Curve Estimation.
CURVEFIT
  /VARIABLES=HRQoL WITH ModVig_minsweek
  /CONSTANT
  /MODEL=CUBIC
  /PLOT FIT.



--
Sent from: http://spssx-discussion.1045642.n5.nabble.com/

=====================
To manage your subscription to SPSSX-L, send a message to
[hidden email] (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