About barchan

barchan · ‎05-28-2015

Just a hint: Are the system options XSYNC and XWAIT turned on? options xwait xsync; proc options option=(xwait xsync); run;

barchan · ‎12-28-2014

Thanks Rick, I like the graphical presentations and the simple series subroutine. The PDF can also be calculated by replacing the CDF of the gamma distribution in the module by the corresponding PDF. The numerical differentiation that you supplied was a valuable check. Regards

barchan · ‎12-15-2014

An IML module is provided for calculating the distribution function and the density of a positive linear combination of independent central chi-square variates. The cumulative distribution function (CDF) is obtained by inverting the moment generating function as described in the paper by Moschopoulos and Canada: "The distribution function of a linear combination of chi-squares" (Comp. and Maths. with Appls., Vol.10, Nº 4/5, pp 383-386, 1984) . The distribution function is expressed as an infinite series of gamma distribution functions. Similarly, the density (PDF) can be expressed as an infinite series of gamma densities. The number of terms of the series used in calculating the CDF and the PDF should be specified by the user. An upper limit of the truncation error is provided for both, PDF and CDF. Optionally, the CDF may be estimated by simulation.

barchan · ‎07-18-2014

Here is another proposal: data have; input char $char20.; cards; 1234567890 123456 987654321 8888888 999999999 4242424 run; data want; set have; length char_new $5; char_new=left(reverse(char)); char_new=reverse(char_new); run; proc print data=want; run;

barchan · ‎04-16-2014

Hi, what exactly is the syntax with PROC TEMPLATE to change the number of decimals displayed in the output. For example, in the FitStatistics table of PROC REG (template Stat.REG.FitStatistics) display the value of R-square (which is cValue2 in the table) with 2 decimals instead of the usual 4 decimals. Thanks!

barchan · ‎02-20-2014

My previous code was wrong because the output in the datastep CHANGE stopped prematurely. The following is better. data hold; length portfol $4 year qtr 8 stock $4 amt 8; input portfol--amt; cards; p1 2000 1 s1 100 p1 2000 1 s2 200 p1 2000 1 s3 50 p1 2000 1 s4 100 p1 2000 2 s1 150 p1 2000 2 s3 100 p1 2000 2 s4 100 p1 2000 3 s1 200 p1 2000 3 s3 100 p1 2000 3 s4 150 p1 2000 3 s5 100 p1 2000 4 s3 100 p1 2000 4 s4 150 p1 2000 4 s5 100 p2 2000 1 s1 100 p2 2000 1 s2 200 p2 2000 1 s3 50 p2 2000 1 s4 100 p2 2000 2 s3 100 p2 2000 2 s4 100 p2 2000 3 s1 200 p2 2000 3 s3 100 p2 2000 3 s4 150 p2 2000 3 s5 100 run; data prev(rename=(amt=prev_amt)); set hold; qtr=qtr+1; year=year; if qtr=5 then do; qtr=1; year=year+1; end; run; proc sort data=hold(keep=portfol year qtr) out=qtrs nodupkey; by portfol year qtr; run; data change; merge hold prev; by portfol year qtr stock; change=max(0,amt)-max(0,prev_amt); run; data out; merge qtrs(in=_qtrs) change; by portfol year qtr; if _qtrs; run; options missing="-"; proc print data=out; by portfol year qtr; id portfol year qtr; run;

barchan · ‎02-20-2014

Try this one: data hold; length portfol $4 year qtr 8 stock $4 amt 8; input portfol--amt; cards; p1 2000 1 s1 100 p1 2000 1 s2 200 p1 2000 1 s3 50 p1 2000 1 s4 100 p1 2000 2 s1 150 p1 2000 2 s3 100 p1 2000 2 s4 100 p1 2000 3 s1 200 p1 2000 3 s3 100 p1 2000 3 s4 150 p1 2000 3 s5 100 run; data prev(rename=(amt=prev_amt)); set hold; qtr=qtr+1; year=year; if qtr=5 then do; qtr=1; year=year+1; end; run; data change(drop=flag); merge hold(in=_hold) prev; by portfol year qtr stock; if first.qtr & ^_hold then stop; if _N_>1 & qtr^=lag(qtr) then flag+1; if flag then do; prev_amt=max(0,prev_amt); amt=max(0,amt); change=amt-prev_amt; end; run; proc print data=change; by portfol year qtr; id portfol year qtr; run;

barchan · ‎12-12-2013

The following code extracts the very first record and the first record after change of CODES: data want(drop=xcode flag); length xcode $4; do until(last.id); set have; by id; if first.id then do; output; xcode=codes; end; if codes^=xcode then do; if ^flag then output; flag=1; end; end; run;

barchan · ‎12-12-2013

Here is another: proc sort data=have; by group descending xnum; run; data want(keep=group--xmax); length group $1 xxx $100 xmax 8; do until(last.group); set have; by group; if first.group then xmax=xnum; if xnum=xmax then xxx=catx("/",xxx,x); end; run;

Online Status	Offline
Date Last Visited	‎09-01-2015 07:12 AM

Re: X command and call system not working correctly

Re: Distribution function of a linear combination of chi-squares

Distribution function of a linear combination of chi-squares

Re: Last 5 characters from string(irrespective of the length)

Changing the number of decimals displayed

Re: Creating a change-variable with unique data set

Re: Creating a change-variable with unique data set

Re: SUBSETTING

Re: Hard question: Selecting observation(s) conditionally by group

Re: Hard question: Selecting observation(s) conditionally by group

Re: Creating a change-variable with unique data set

Re: SUBSETTING

Distribution function of a linear combination of chi-squares

Re: Distribution function of a linear combination of chi-squares

Re: X command and call system not working correctly

Re: Distribution function of a linear combination of chi-squares

Distribution function of a linear combination of chi-squares

Re: Last 5 characters from string(irrespective of the length)

Changing the number of decimals displayed

Re: Creating a change-variable with unique data set

Re: Creating a change-variable with unique data set

Re: SUBSETTING

Re: Hard question: Selecting observation(s) conditionally by group