epi_econ_code Tracker

Mixed model PROC LOGISTIC

epi_econ_code — Wed, 08 May 2024 18:25:57 GMT

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!

Re: PROC GMAP

epi_econ_code — Mon, 29 Apr 2024 19:02:06 GMT

@GraphGuy I came across your solution for other graph-related questions. Is this something you can help with? I really would appreciate it.

Re: PROC GMAP

epi_econ_code — Mon, 29 Apr 2024 16:43:25 GMT

Not sure how to code for adding the count as bubble (I am new to using PROC GMAP)
I used the following code to get the current map: legend1 label=(f=swissb j=c '% of households who need meals') across=4 down=2 frame position=(top center
inside);
legend2 label=(f=swissb j=c 'Number of school meal programs') across=3 down=1 frame position=(bottom center
inside);
title h=4 f=swissb 'Availability of meal programs and percentage of households in need';
proc gmap data = ccs2019 map=maps.us ;

id state;
area les_tert/legend=legend1;
block ccs_count/levels=3 legend=legend2 ;
format les_tert les_tert.;
run;
quit;

PROC GMAP

epi_econ_code — Mon, 29 Apr 2024 14:48:42 GMT

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

Annotating using dot PROC GMAP

epi_econ_code — Sun, 28 Apr 2024 06:53:16 GMT

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 mapWhat I want to achieve