Dear SAS community,   Simple question: How do I create a bar chart?   I would like to graph a bar chart instead of the line graph below for the following data:   data graph4; input buckets change; datalines;   1 4.5 2 1.2 3 2.3 4 3.2 5 0.4 ; run;   title "Graph 4: Change in 2012"; PROC SGPLOT = WORK.graph4; series x=buckets y=change; xaxis grid values=(0 to 5 by 1); RUN;   I would like to change the number of buckets, hence the automation. I would like to draw up buckets in the x-axis and change in the y-axis.   Thanks for your help. ... View more

Hi All,   I would like to calculate the average of costs at each percentile based on contracts. I got as far as calculating the percentile limits, but this gives me a dataset with varilables; P_1, P-2, ... P_100. I need to somehow merge these into my original dataset and then use proc means I guess to calculate the averege cost per percentile.   What's the most efficient way of doing this?   data claims; input client year cost contract; datalines;   18 2011 589.92 10.2 19 2010 1629.8 11.5 20 2011 1813.29 2.3 21 2012 412.06 3.2 22 2012 219.82 0.4 23 2010 3669.98 5.9 24 2012 4879.63 6.8 ; run;   proc univariate data=claims noprint; var contract; output pctlpre=P_ pctlpts= 1 to 100; run;   ... View more

  Dear Members,   I have a large dataset with healthcare products (over 4000 products). Each product is made up of activities. The required activities change a little bit each year, and some of the products are repackaged (or renamed a different product name) without the underlying activities changing (see product 2 in 2014). I would like to link these repackaged products together by giving them a new id.    What is the most efficient way of doing this for a large dataset?    data Have;    input product_id      year activity; cards; 1     2012    25      1     2012    32 1     2013  25 1     2013    32 1     2013    11 2     2014    25      2     2014    32 ; run;   Desired output:   Product_id year activity new_id 1 2012 25 1 1 2012 32 1 1 2013 25 2 1 2013 32 2 1 2013 11 2 2 2014 25 1 2 2014 32 1     ... View more

    data clients; input client year cost; datalines; 18 2011 589.92 18 2010 1629.8 18 2011 1813.29 18 2012 412.06 19 2012 219.82 20 2010 3669.98 20 2012 4879.63   ; run;   data work.client2; set work.client; dif_cost=dif(cost); if first.client then dif_cost = . ; run;   I get an error message that first.client is uninitialized. Why? My data is sorted by client and year.    Thanks ... View more

Hello,    The code below is from the following link (https://stats.idre.ucla.edu/sas/faq/how-can-i-bootstrap-estimates-in-sas/) for bootstrapping. I am trying to run this code, but I get an error msg in the beginning. "Label2" seems to be uninitialized.    Basic statistics are saved into dataset t0, and we are trying to save the value where the  label says RSquare.   My code is the following:   ods output FitStatistics = t0; proc glm data=data_predict; class year (ref="2010"); model cost=contract|year|age|female /solution; run; quit;   *store the estimated r-square; data _null_; set t0; if label2="RSquare" then call symput('r2bar', cvalue2); run; ... View more

Hello,    I'm trying to run a Proc GLM with categorical variables (year and age-groups). I would like to see the interaction of each year and each age-group (plus female, etc.) explicitly WITHOUT having to create dummy variables for each year, as this would make my model very cumbersome.    Indicating year and age_gr in a CLASS statement works well with PROC logistic on the same model. I get individual coefficients for year=2010  to year=2013, but this does not work with Proc GLM and I get one coefficient for the year variable only.   What am I doing wrong?     ********************************************************   proc logistic data=data_predict descending; class year age_gr; model treated=contract|year|age|female; run;   proc glm data=data_predict; class year age_gr; model prob_score=contract|year|age_gr|female; run; ... View more

Hello,    I'm still struggling with the EffectPlot procedure. I would like to graph the probability of getting treated (y-axis) per year (x-axis) with 95% confidence interval.    ODS GRAPHICS ON; proc logistic data=data descending; class year; model treated=year; effectplot slicefit(x=year)/noobs; run; ODS GRAPHICS Off;   I get the following error message:   ERROR: The X= variable must be a continuous variable used in the model.   Thanks ... View more

Hello,    I would like to create a plot with the probability of being treated on the y-axis and the contract-type (continuous variable between 0-1) on the x-axis. I would like to compare the effect of an increasing contract-type between the years, so I want to see this line graph per year and stacked in the same graph.    This is what I have, but I am not sure how to show this with a different graph per year.   proc logistic data=data descending;     model treated = contract y2011 y2012 y2013 y2011*contract y2012*contract y2013*contract; store logit; run;   title "Change in the probability of being treated by year and contract-type";   ODS GRAPHICS ON; proc plm source=logit;    effectplot slicefit(x=contract plotby=...); run; ODS GRAPHICS Off; ... View more