About Ksharp

Ksharp · ‎03-26-2026

data have; infile cards expandtabs; input ID $ visit VISDATE :$20. PAY; cards; A 10 2/18/2022 1 A 11 4/18/2022 1 A 12 5/6/2023 1 A 13 10/6/2023 2 A 14 7/8/2024 2 A 15 12/8/2024 1 A 16 2/28/2025 2 A 17 9/25/2025 1 A 18 12/7/2025 1 ; data want; set have; by id pay notsorted; retain new_visitdate; if first.pay and pay=1 then new_visitdate=VISDATE; if pay ne 1 then new_visitdate=VISDATE; run;

Ksharp · ‎03-25-2026

Use left-side SUBSTR() function: data news; val='06MAY98'; substr(val, 6, 2) = 99; proc print data=news; run;

Ksharp · ‎03-25-2026

I would prefer to use this : %macro ds_rowct(ds); %let dsid=%sysfunc(open(&ds.)); %let nobs=%sysfunc(attrn(&dsid.,nlobs)); %let dsid=%sysfunc(close(&dsid.)); %put &ds. have &nobs. obs.; %mend; %ds_rowct(sashelp.class) 1 %macro ds_rowct(ds); 2 %let dsid=%sysfunc(open(&ds.)); 3 %let nobs=%sysfunc(attrn(&dsid.,nlobs)); 4 %let dsid=%sysfunc(close(&dsid.)); 5 6 %put &ds. have &nobs. obs.; 7 %mend; 8 9 10 11 %ds_rowct(sashelp.class) sashelp.class have 19 obs.

Ksharp · ‎03-25-2026

Yes. I think so. Since proc bglimm is the Bayes version of proc glimmix, they should get the similar parameter estimators. And I suggest you to post this question to Stat Forum to let statistician to notice it. https://communities.sas.com/t5/Statistical-Procedures/bd-p/statistical_procedures

Ksharp · ‎03-24-2026

1) You have too many strata/cluster in your model due to two RANDOM statement. I would suggest to get rid of the second RANDOM statement and say whether it is working or not. 2) Or using Bayes version of GLIMMIX: proc bglimm

Ksharp · ‎03-24-2026

It really looks weird. But if you change this, the issue is gone. /* ---> / * Take Tom 's code for example: data test; string1 = 'this has / * a comment */ in it'; string2 = resolve(string1); string3 = resolve(tranwrd(string1,'*','%str(*)')); put (string1-string3) (=/); run; 55 data test; 56 string1 = 'this has / * a comment */ in it'; 57 string2 = resolve(string1); 58 string3 = resolve(tranwrd(string1,'*','%str(*)')); 59 put (string1-string3) (=/); 60 run; string1=this has / * a comment */ in it string2=this has / * a comment */ in it string3=this has / * a comment */ in it

Ksharp · ‎03-21-2026

This graph would be better for Olympic medal count. https://blogs.sas.com/content/graphicallyspeaking/2014/02/12/sochi-medal-graphs/

Ksharp · ‎03-21-2026

Yeah. The most difficult thing is the Y axis which is non-linear scale . Otherwise , I would like to use PROC SGPLOT to replicate this graph .

Ksharp · ‎03-20-2026

You want this ? data have; input Obs Bimodal_Var; cards; 1 11.18227 2 8.349733 3 14.40938 4 24.66227 5 9.279151 6 7.238792 7 15.46339 8 22.47347 9 15.132 10 21.10275 11 8.501284 12 15.77826 13 16.20075 14 17.67517 15 7.934397 16 18.39657 17 11.40699 18 8.399611 19 8.490487 20 22.70707 21 11.09486 22 11.00888 23 22.85953 24 6.783415 25 11.31967 26 25.19409 27 12.12261 28 10.3787 29 22.13491 30 20.40574 31 9.802019 32 15.45317 33 17.87584 34 10.8164 35 10.8377 36 11.4462 37 6.82155 38 10.80242 39 19.35702 40 17.50303 41 16.63369 42 21.50525 43 23.79174 44 10.77436 45 9.775864 46 10.62838 47 9.338187 48 8.149213 49 17.10799 50 26.5875 ; /*Get N Mean Std */ proc summary data=have ; var Bimodal_Var; output out=mean_std_n n=n mean=mean std=std; run; data _null_; set mean_std_n; call symputx('mean',put(mean,8.2 -l)); call symputx('std' ,put(std,8.3 -l)); call symputx('n' ,n); run; %put &=mean &=std; /*Get AD statistic and P-Value*/ ods select none; ods output TestsForNormality= TestsForNormality; proc univariate data=have normal ; var Bimodal_Var; run; ods select all; data _null_; set TestsForNormality(where=(Test='Anderson-Darling')); call symputx('AD',put(Stat,8.3 -l)); call symputx('pvalue',cats(pSign,vvalue(pValue))); run; %put &=ad &=pvalue; /*Get the normal probability table*/ ods graphics /reset=index noborder; ods listing gpath="%sysfunc(pathname(work))" style=htmlblue;; *Save this plot into a path; ods select ProbabilityPlot; proc reliability data=have ; probplot Bimodal_Var/NOINSET; run; ods select all; %sganno data sganno; %SGIMAGE(IMAGE="%sysfunc(pathname(work))\ProbabilityPlot1.png",ANCHOR="topleft",BORDER="FALSE",DRAWSPACE="LAYOUTPERCENT" ,x1=-1,y1=101) %SGTEXT(LABEL="Mean", BORDER= "FALSE",DRAWSPACE= "LAYOUTPERCENT" , LAYER= "FRONT",X1=10,Y1=90) %SGTEXT(LABEL="&Mean", BORDER= "FALSE",DRAWSPACE= "LAYOUTPERCENT" , LAYER= "FRONT",X1=20,Y1=90) %SGTEXT(LABEL="StDev", BORDER= "FALSE",DRAWSPACE= "LAYOUTPERCENT" , LAYER= "FRONT",X1=10,Y1=85) %SGTEXT(LABEL="&std", BORDER= "FALSE",DRAWSPACE= "LAYOUTPERCENT" , LAYER= "FRONT",X1=20,Y1=85) %SGTEXT(LABEL="N", BORDER= "FALSE",DRAWSPACE= "LAYOUTPERCENT" , LAYER= "FRONT",X1=10,Y1=80) %SGTEXT(LABEL="&n", BORDER= "FALSE",DRAWSPACE= "LAYOUTPERCENT" , LAYER= "FRONT",X1=20,Y1=80) %SGTEXT(LABEL="AD", BORDER= "FALSE",DRAWSPACE= "LAYOUTPERCENT" , LAYER= "FRONT",X1=10,Y1=75) %SGTEXT(LABEL="&ad", BORDER= "FALSE",DRAWSPACE= "LAYOUTPERCENT" , LAYER= "FRONT",X1=20,Y1=75) %SGTEXT(LABEL="P-Value", BORDER= "FALSE",DRAWSPACE= "LAYOUTPERCENT" , LAYER= "FRONT",X1=10,Y1=70) %SGTEXT(LABEL="&pvalue", BORDER= "FALSE",DRAWSPACE= "LAYOUTPERCENT" , LAYER= "FRONT",X1=20,Y1=70) run; ods graphics/reset noborder; proc sgplot data=sashelp.class sganno=sganno; scatter x=weight y=height; run;

Ksharp · ‎03-20-2026

Are they exactly one-to-one matched ? One parent might match multiple child ? And could you post some data and desired output to explain your problem ?

Ksharp · ‎03-19-2026

Interesting . Although there is an ERROR in your LOG. 5 proc http 6 method="get" 7 url='https://api.census.gov/data/timeseries/intltrade/imports/hs?get=CON_QY1_MO,CTY_CODE,CTY_NAME,I_COMMODITY,I_COMMODITY_SDESC,UNIT_QY1,GEN_VA 7 ! L_MO&YEAR=2025&MONTH=01&COMM_LVL=HS10&I_COMMODITY=76*' 8 verbose 9 out=out 10 headerout=hdrout; URL = https://api.census.gov/data/timeseries/intltrade/imports/hs?get=CON_QY1_MO,CTY_CODE,CTY_NAME,I_COMMODITY,I_COMMODITY_SDESC,UNIT_QY1,GEN_VAL_MO&YEAR=20 25&MONTH=01&COMM_LVL=HS10&I_COMMODITY=76* METHOD = get Out = c:\temp\Output.txt Header Out = c:\temp\Response.txt 11 debug level=3; 12 run; ERROR: 连接至“148.129.75.191:443”时出错。(连接超时。) ERROR: 无法连接 Web 服务器。 But the json file you want to download is already in there , check this: filename out "c:\temp\Output.txt"; filename hdrout "c:\temp\Response.txt"; proc http method="get" url='https://api.census.gov/data/timeseries/intltrade/imports/hs?get=CON_QY1_MO,CTY_CODE,CTY_NAME,I_COMMODITY,I_COMMODITY_SDESC,UNIT_QY1,GEN_VAL_MO&YEAR=2025&MONTH=01&COMM_LVL=HS10&I_COMMODITY=76*' verbose out=out headerout=hdrout; debug level=3; run; Therefore you could ignore this ERROR info I guess.

Ksharp · ‎03-19-2026

You need to resort to the function VVALUE(). Here is an example: proc format; value a 1='a' 2='b' ; value b 1='xxx' 2='yyy' ; run; data have; input a b; format a a. b b.; cards; 1 2 2 2 1 1 ; data want; set have; array x{*} a b; array y{*} $ 40 v1 v2; do i=1 to dim(x); y{i}=vvalue(x{i}); end; run;

Ksharp · ‎03-18-2026

That issue is due to the encoding of "DAT.UPDATED_DIRECTORY" is different with SAS encoding. Make SAS encoding be same with your dataset. You could use the following code to check your dataset encoding: proc contents data=DAT.UPDATED_DIRECTORY; run;

Ksharp · ‎03-17-2026

The way you are using "proc mianalyze" is wrong. "proc mianalyze" only can pool estimated parameter like t-value , p-value. Since your "proc freq" does not include these statistics, you can not pool this results. And better post it at Statistical Forum: https://communities.sas.com/t5/Statistical-Procedures/bd-p/statistical_procedures

Ksharp · ‎03-16-2026

You gotta need the help of SAS/IML to use matrix multiply operator. Calling @Rick_SAS https://blogs.sas.com/content/iml/2014/02/19/scoring-a-regression-model-in-sas.html

Online Status	Online
Date Last Visited	2 hours ago