Hi all, I am a new SAS user and am working on a sensitivity analysis where I am predicting covid infection by water, sanitation and hygiene data. I am using proc genmod with link=log to produce risk instead of odds due to my study design. I am also testing an interaction between district and my primary predictor variable, a wash composite score. I have three districts and would like to produce risk for each of my districts (chiradzulu, chikawa, and blantyre) for a composite score of 1, 5 and 10. I also have 2 covariates, disability_num(binary 0/1) and num_members (count) Can someone help me? Here is what I have so far. I know I have to use an estimate statement but am struggling with what to enter!    PROC GENMOD DATA = CAPSTONE.FINAL DESCENDING; CLASS DISTRICT (REF = 'Chikwawa')/ PARAM = REF ; MODEL SENSITIVITY (EVENT='1')= SUM_COMP DISABILITY_NUM NUM_MEMBERS DISTRICT DISTRICT|SUM_COMP/ DIST = BIN LINK = LOG; estimate 'chiradzulu 1' sum_comp 1 district 1/exp; estimate 'chiradzulu 10' sum_comp 10 district 1/exp; estimate 'blantyre 1' sum_comp 1 district -1/exp; estimate 'blantyre 10' sum_comp 10 district -1/exp; estimate 'chikwawa 1' sum_comp 1 district 0/exp; estimate 'chikwawa10' sum_comp 10 district 0/exp; RUN;   my design matrix:    Class Level Information Class Value Design Variables District Blantyre 1 0   Chikwawa 0 0   Chiradzulu 0 1   I would appreciate some help! It is my first time doing an interaction with a three level categorical variable so please be kind. thank you   ... View more

Hi all,   I am a new sas user trying to clean up a dataset. I am interested in coding some categorical variables into a composite variable of water and sanitation quality.    I have a dataset from nicaragua where responses have accents. It was a previous SPSS file, and when I import to SAS all accent characters convert to unknown characters and prevent me from running a proc freq. I simply need 9 categories to code as 0/1, so the variable could simply be converted to numeric if I know what the values mean. Can someone please advise on how to get rid of the unknown character? The variable in question is S1P15   dataset is attached. Current proc import below:  *import spss dataset and convert; proc import datafile = "/folders/myfolders/sasuser.v94/WFP/datasets/EMNV14-02 DATOS DE LA VIVIENDA Y EL HOGAR (1).SAV" out= work.nicaragua dbms=sav replace; run;   Thank you!  ... View more

Hi everyone,    I am a newer SAS user and am working on a project that requires me to create macro that creates a histogram with a variety of macros inside. I have created the following code, but although I get no errors or warnings in my log, I get no results. Can anyone help?   Thanks so much.    my code:    %MACRO mymacro ( DSNm = , Vbl = , Stats = N MEAN STDDEV , Ndec = 1); TITLE1 "title1"; proc sgplot data = work.&DSNm; histogram &Vbl/N MEAN STDDEV; run; %MEND title2; %( DSNm = , Vbl = , Stats = N MEAN STDDEV , Ndec = 1); ... View more

Hi all,   I am working on creating a graph that includes two lines of xaxis labels, one bolded and the second not. I am unsure how to selectively bold the first line in proc sgplot. I would appreciate any help. Thanks!    graph I am attempting to produce:      current code:    TITLE 'title1'; PROC SGPLOT DATA=work.test; VBAR State / RESPONSE=Col FILLATTRS=(color=blue); XAXIS LABEL = 'State ' LABELATTRS=(WEIGHT=bold); YAXIS LABEL='Percent' LABELATTRS=(WEIGHT=bold) GRID VALUES=(0 to 4 by 1); REFLINE 1.8 / AXIS=y LINEATTRS= (COLOR=indigo pattern=dash) LABEL='NA = 1.8' LABELLOC=inside; RUN; TITLE;   Thank you! ... View more

Hi all,   I have been working on this problem for quite a while and cannot seem to figure out how to break up my proc report statements so as to display mean and std simultaneously for a numeric variable.    The table I wish to produce is:   Average Pressures   Systolic Diastolic Race Code N MEAN STD MEAN STD Asian 11 117.1 11.1 69.1 9.0 Black 83 123.3 15.3 74.0 7.0 Other Race 30 115.2 15.5 69.8 11.0 White 449 120.8 15.5 70.0 7.6   573 120.8 15.4 70.6 7.8   My current code will not allow me to simultaneously display MEAN and STD: PROC REPORT DATA = hypanal.hypanalysis1 Split="/" MISSING; COLUMN RaceCd N ("Average Pressures" ("Systolic" SBP SBP) ("Diastolic" DBP DBP)); DEFINE RaceCd /group; Define N / "N"; DEFINE SBP / ANALYSIS MEAN STD format=4.1 "MEAN" "STD"; DEFINE DBP / ANALYSIS MEAN STD format=4.1 "MEAN" "STD"; RBREAK AFTER / SUMMARIZE; RUN;   the result is close, but only displays STD instead of STDs and MEANS:   Average Pressures  Systolic DiastolicRace Code N MEANSTD MEANSTD MEANSTD MEANSTD A 11 11.1 11.1 9.0 9.0 B 83 15.3 15.3 7.0 7.0 O 30 15.5 15.5 11.0 11.0 W 449 15.5 15.5 7.6 7.6   573 15.4 15.4 7.8 7.8   I would appreciate any help!!    Thanks so much ... View more