Hi I am currently building a dashboard on viya for learners as part of my masters program and for one part I need to display values on a bar chart but I also would like to add a line on each bar that is the cut-off value but I can't find a object that allows me to do it in such a way. Is there any way I can do this ?   This is how my Bar chart looks I have values for each year for the size but I still need to display the cut-off value on this bar chart in a line or something for each year can this be done in VA? ... View more

Hello, everyone   I managed to estimate the parameters for a Weibull accelerated failure time model from the PROC LIFEREG procedure. I have the point and interval estimations for predictor coefficients. In addition, I have the estimations for the Weibull shape parameter.   After I get these estimations, I would like to simulate the cumulative incidence curve/survival curve. I know I should use the point estimation of these coefficients and the shape parameter in the prediction model. But I am not sure whether to use the interval estimation of coefficients, or the shape parameter.   You answer is much appreciated! ... View more

I've been looking for a good, efficient method to create "Table 1", which has a more or less standard format in research journals: Variable Name Category 1 Category 2 P-value Variable 1       1 N(%) N(%) Chi-sq or Fisher's 2 N(%) N(%) 3 N(%) N(%) Variable 2 Mean (SD) Mean (SD) t-test   Median (IQR) Median (IQR)   I found this PDF paper. It seems to work pretty well, but I'm having 2 issues. The first is that the chi-square results are saved into a file and then merged into the main file. But the value of the "variable" variable is the var label in 1 file and the var name in the other, so they don't merge.   The second issue is that Proc Report does not produce a table and instead gives me this: NOTE: Groups are not created because the usage of levels is DISPLAY. To avoid this note, change all GROUP variables to ORDER variables. WARNING: A GROUP, ORDER, or ACROSS variable is missing on every observation.   data have; infile datalines dsd dlm=',' truncover; input DEM_AGE DEM_SEX cohort_flag TM_group; datalines; 3,1,1,0 2,1,1,1 3,2,1,1 3,2,1,1 3,2,1,0 2,2,1,1 2,1,1,1 3,1,1,1 2,1,1,1 3,2,1,0 2,1,1,0 2,2,1,1 3,2,1,0 2,2,0, 3,2,1,1 3,2,1,1 3,1,1,0 3,2,1,0 2,1,1,0 3,1,1,1 3,2,1,1 3,2,1,0 3,2,1,1 3,2,1,1 3,2,1,1 ; RUN; /*Load formats from existing file in temp folder*/ options fmtsearch=(temp.formats); /*DEM_SEX sex 1 Male 2 Female*/ /*DEM_AGE AGE2GRP 1:Age Group <65 2:Age Group [65,75) 3:Age Group >=75*/ /*TM_group yesno 1 Yes 2 No*/ /*Generate descriptive statistics*/ proc means data = temp.have noprint n sum mean; class DEM_AGE DEM_SEX; var TM_group /*MA_group*/; ways 1; output out = temp.expl_PreTable n = sum = mean = / autoname; WHERE cohort_flag = 1; run; /*Format descriptive stats*/ data temp.expl_Table (keep = variable levels TM_group_N TM_group_sum TM_group_mean pct ExpPct indexvar); set temp.expl_PreTable; length variable $ 20; /* These four variables */; length levels $ 20; /* will describe the first */; length pct $ 8; /* four columns of the table */; length ExpPct $ 15; if DEM_AGE ne . then do; /*Building "variable" and "Levels" columns for "DEM_AGE"*/; variable = 'Age category'; levels = put(DEM_AGE, age2grp.); IndexVar = 1; /*This index is included just in case the order of data presentation needs to be changed*/; end; if DEM_SEX ne . then do; /*Building "variable" and "Levels" columns for "DEM_SEX"*/; variable = 'Sex'; levels = put(DEM_SEX,sex.); IndexVar = 2; end; pct = put(TM_group_mean*100,4.1); /*Calculate % exposed */; ExpPct = compress(put(TM_group_sum,comma4.),' ') ||' '||'('||compress(pct,' ')||')'; /*creating data in the form of "count (%)" */; run; /*Run chi-square significance tests and use ODS to create a dataset of these results*/ ods trace on; ods output chisq = temp.expl_ChiData; proc freq data = temp.have; table TM_group*DEM_AGE / chisq; table TM_group*DEM_SEX / chisq; WHERE cohort_flag = 1; run; ods trace off; /*Rearrange chi-square dataset so it can be merged with descriptive stats table*/ data temp.expl_ChiData2 (keep = variable prob); set temp.expl_ChiData (where = (statistic = 'Chi-Square')); length variable $ 20; variable = scan(table,-1,' '); /* Returns the last word in a character value from the "table" variable*/ run; /*Sort both tables so they will merge*/ PROC SORT data=temp.expl_Table OUT=temp.expl_Table_sort; BY variable; RUN; PROC SORT data=temp.expl_ChiData2 OUT=temp.expl_ChiData2_sort; BY variable; RUN; /**************MERGE DOES NOT WORK BECAUSE The value of the "variable" variable is the var label in 1 file and the var name in the other. DUE TO THE 'IF a' STATEMENT, NOTHING FROM THE CHIDATA2 FILE IS MERGED IN****/ /*Merge descriptive stats table with chi-square table*/ DATA temp.expl_TableData; MERGE temp.expl_Table_sort (in = a) temp.expl_ChiData2_sort (in = b); BY variable; IF a; RUN; /*Use PROC REPORT to create final output table*/ proc report data = temp.expl_TableData nowd; column variable levels TM_group_N ExpPct prob; define variable / "Variable" group format = $variable.; define levels / " " ; define TM_group_N / "TM" /*format = comma5.*/; define ExpPct / "TM Group/n (%)"; define prob / "p-value" group format = pvalue6.4; Title "Table 1. Descriptive characteristics of individuals in the sample"; RUN; ... View more