Hello everyone,  I am new to statistical programming for a multilevel/mixed modeling. I am conducting a study where the outcome is (YES/NO) binary indicating the provision of HIV treatment at the facility level. My primary independent variable is the facility payment type i.e. accepts private, public insurance or cash payment; covariates include the census region, state policy on HIV programs, and state HIV prevalence rate.  I am conducting logistic regression analysis as follows:  Provision of HIV treatment= B0 +B1 (payment type) + B2 (census region) + B3 (state HIV policy) + B4 ( tertile of state HIV prevalence rate) My goal is to assess the association between payment type and provision of HIV treatment. My question is that the unit of analysis for this study is the FACILITY whereas the information for state HIV policy and prevalence rate is pulled from the state level data and then extended to facility level depending on which state the facility belongs to.   For eg: If California has highest HIV prevalance rate, according to the tertile category, it falls in the highest tertile group with value 2; because the tertile has values 0,1,2. Now because California has the value of prevalence rate=2, all the facilities belonging to the state of California has value of 2 assigned to them for the purpose of the analysis. Similarly facilities belonging to other states have been assigned values depending on which category the state is assigned.   I am not sure which is the correct model to use here? does simple proc logistic will do and if the methodology of extending state level information to facility level is correct because there might be states with a large number of facilities or vice versa?   Your help is appreciated!   ... View more

Hi everyone,  I am trying to create a map using the gmap procedure. I was able to successfully create the map as shown in the attached images. However, the block option allows us to create a 3-D map and the journal I want to publish this in does not accept 3-D graphs. How can I show the counts of the facilities that is shown using the block option in a 2-D format? Basically, how can I incorporate the information coming from two response variables in a choropleth map? Help to change this block into a 2-D map   ... View more

Coordinates Not Being Displayed in proc gmap  PROC GMAP color change by level  CHORO and ANNOTATE combination I have a USA map ready with the default PROC GMAP settings categorized based on the tertile of a service need in the US.   I used the following code: proc gmap data = data2019 map=maps.us; id state; choro service_need_tertile; run; quit;   Now, I have the count of facilities in each of the US that is present to fulfil that service need. I have to add that information on the map; where one dot should represent one facility in that state (similar to the second attached graph with red and green color) Any help in accomplishing this task will be helpful. The name to the variable that gives the count is SERVICE_COUNT       MY current map What I want to achieve ... View more