proc surveylogistic data =newnsch3a; class sex (ref = '0') race(Ref = '0') FPL (ref= '0') smokes(Ref= '0') adulteducate (ref='0') family (ref='0') sleep (ref='0') comp (ref='0') age3_1718 (ref='2') MOLD (ref='2') exercise (ref='0') pest (ref='0')/param=glm ; strata Fipsst; cluster hhid; weight fwc; Model ast (event ='1') = age3_1718 sex race FPL smokes comp adulteducate MOLD exercise pest family sleep smokes*comp; lsmeans smokes*comp/oddsratio cl diff; slice smokes*comp/ sliceby=comp diff=control('0' '0') oddsratio cl; run; Hello statisticians, I am interested in finding the P for trend for the Odd Ratios (OR) that i got from the interaction between smoking and a composite variable that i created. Composite is a variable that has 4 level (0,1,2,3) while smokes have 2 levels (0,1). The interaction between composite and smokes was significant. So i have 4 ORs. My problem is that i would like to get a P for trend on these 4 ORs. Is there any other syntax that i should add to my code (above) to enable me get a P-value for the trend?   ... View more

Hello programmers,   I am very familiar with using the SELECTION= option in sas in order to perform forward, backward or stepwise model selection. However, is anyone aware if i can do this in proc survey logistic? If so, any referral to SAS documentation and syntax will be helpful. This document did not address model selection in proc survey logistic....and i could not find any yet. https://support.sas.com/documentation/onlinedoc/stat/131/surveylogistic.pdf Thanks   ... View more

Hello statisticians, Please i'll be glad to get any input on this as mixed models are not my strong suit. It's a clinical trial data comparing 2 treatments.  I want to know 1. if the two treatments differ in their effects on length (outcome) 2. If there's a  difference in the pattern of  change between subjects receiving the two treatments and if  one treatment show results more quickly?  Attached is my code proc mixed data= new1 COVTEST method=ml; Class ID treat monthcat; MODEL lenght= month treat month*treat /solution; RANDOM intercept month /SUB=ID TYPE=UN G V; repeated monthcat/subject=id type=toep r ; run;     My thought is that Number 1 is asking for treatment effects on outcome, so i will use -0.1278. But what does the negative effect mean? Then what does "pattern of change" imply in number 2. Is it the treatment*time effect?     ... View more

Dear programmers,   I want to use the estimate statement in proc genmod but my code is not giving me what i want. My outcome is status (Yes vs No) . My predictors are dose ( 1 vs 0) and Time measurements (1, 2, 3 and 4). Simply put, the status of the subjects were taken at these 4 time points. I also want to include the interaction term between dose and time.  The problem i have is specifying the estimate statement so that i can get the Odds Ratios. Attached are my codes. proc genmod data = newdata descending; class ID dose time/PARAM=EFFECT DESCENDING ; model status=dose time dose*time/dist=bin link=logit; repeated subject=ID/type=ar(1) corrw covb modelse within=time; estimate "dose" dose -1 1 / exp; estimate "time4" time 1 0 0/exp; estimate "time3" time 0 1 0/ exp; estimate "time2" time 0 0 1/exp; ESTIMATE "time1" Time -1 -1 -1/ exp; run; ... View more