Dear Madam/Sir,   I would like to have standard deviations (aqresstd) on five values (variable name: aqres) for the past four years (t-4,t-3,t-2,t-1) and the current year (t). I ran the data using the suggested program from this community and there is no error message as follows:   data s9; do n=1 by 1 until(last.gvkey); set s8; by gvkey fyear; array t(50); if n>=5 then aqresstd=std(of t(*)); t(n)=aqres; if n(of t(*))>4 then do;k=n-4; call missing(t(k));end; output; end; drop n t:; run; NOTE: There were 206263 observations read from the data set WORK.S8. NOTE: The data set WORK.S9 has 206263 observations and 70 variables. NOTE: DATA statement used (Total process time): real time 0.54 seconds cpu time 0.54 seconds   However, I do not have correct standard deviation values as follows;                                                                                    gvkey        fyear 34567891011121314151617181920212223242526272829303132 001004 1989 29.51 . 001004 1990 -14.73 . 001004 1991 10.10 . 001004 1992 1.12 . 001004 1993 2.64 18.456 001004 1994 -0.51 10.439 001004 1995 -6.02 4.688 001004 1996 26.90 3.774 001004 1997 22.23 14.544 001004 1998 9.38 16.334 001004 1999 25.65 14.752 001004 2000 1.07 8.017 001004 2001 -67.39 11.407 001004 2002 -24.66 41.002 001004 2003 -3.85 39.755 001004 2004 0.23 31.183 001004 2005 0.65 30.961 001004 2006 3.30 12.008 001004 2007 56.96 2.953 001004 2008 81.35 27.817 001004 2009 40.46 40.058 001004 2010 -38.11 32.775 001004 2011 48.90 51.657 001004 2012 -8.55 50.672 001004 2013 53.31 41.223 001004 2014 -235.42 44.676 001004 2015 11.76 136.268 001004 2016 -25.15 129.714 001004 2017 -7.02 128.429 001004 2018 -52.38 115.299   It will be greatly appreciative if you can advise me how to fix the program. Thanks.   Sincerely, Joon ... View more

Hi Madam/Sir,   In the SAS 9.2, the output of "proc surveylogistic" generates wald chi-square statistics. But in the SAS 9.4, it generates t-statistics.  In the SAS 9.4, proc logistic still generates wald ch-square statistics. But I need to use proc surveylogstic to do firm-level clustering in the standard variation of variables. Is there a way to report wald chi-square statistics? Your help will be highly appreciated. Joon.   proc surveylogistic data=s19; class gvkey; model gc(descending) = sdd1 starmeasure lev chlev lnage return beta sd_residual zscore lnasset ploss investment ffinancing ocf reportlag big4 y2001-y2016 id1-id65/stb; run; Analysis of Maximum Likelihood EstimatesParameter Estimate StandardError t Value Pr > |t| StandardizedEstimateInterceptsdd1starmeasurelevchlev -6.6960 0.9303 -7.20 <.0001   1.3863 0.2235 6.20 <.0001 0.2345 1.0422 0.2931 3.56 0.0004 0.2327 0.6787 0.2085 3.26 0.0011 0.1338 0.5872 0.3805 1.54 0.1229 0.0683 ... View more

Hi Madam/Sir,   I would like to choose the last observation out of two observations in the same firm-year. For instance, I would like to choose observations # 2, 4, and 6 and drop observations #1,2, and 3. It will be highly appreciative if you can provide the programming code. Thank you. Joon1 Obs CIKCODE cyear ANNUALREPORTDATE dirind123456 0000002024 2016 20160101 10 0000002024 2016 20161201 11 0000010254 2014 20140301 3 0000010254 2014 20141201 5 0000018169 2011 20110101 4 0000018169 2011 20111201 5 ... View more

Hi Madam/Sir,   I would like to count the number of observations when ind_dir=1 in the same company (cikcode) and in the same year (annualreportdate). The following code was tried and printed. However, as shown in the excel sheet containing data, the whole observation per company-year was counted regardless whether ind_dir=1 or ind_dir=0. How can I fix the program to count the number of observations only when ind_dir=1 per company-year pair? 163 proc sql; 163! create table count1 as select*, count(ind_dir=1) as dirindsum from os5 group by 163! cikcode, annualreportdate; run; NOTE: The query requires remerging summary statistics back with the original data. NOTE: Table WORK.COUNT1 created, with 2474560 rows and 40 columns. NOTE: PROC SQL statements are executed immediately; The RUN statement has no effect. 164 NOTE: PROCEDURE SQL used (Total process time): real time 8:43.94 cpu time 1:09.04 165 proc print data=count1(obs=1000); var cikcode annualreportdate ind_dir dirindsum; run; NOTE: There were 1000 observations read from the data set WORK.COUNT1. NOTE: PROCEDURE PRINT used (Total process time): real time 0.32 seconds cpu time 0.31 seconds   Thank you very much for your help in advance   SIncerely, Joon1   ... View more

          Dear Madam/Sir,   I would like to construct standard deviation of ROA with a rolling window of the past five years. After referring several comments in this forum, I made the following command and had the error messages (attached).   proc sql; create table m24 as select *, (std(a.roa) from have where a.fyear between a.fyear-6 and a.fyear-1 and a.gvkey) from m23 as a; quit; Obs gvkey FYEAR roa 1 1004 1999 0.04745 2 1004 2000 0.02640 3 1004 2001 -0.08299 4 1004 2002 -0.01807 5 1004 2003 0.00494 6 1004 2004 0.02110 7 1004 2005 0.03592 8 1004 2006 0.05494 9 1004 2007 0.05517 10 1004 2008 0.05710 11 1004 2009 0.02973 12 1004 2010 0.04098 13 1004 2011 0.03084 14 1004 2012 0.02574 15 1004 2013 0.03314 16 1004 2014 0.00673 17 1004 2015 0.03308 18 1004 2016 0.03756 19 1013 2000 0.21864 20 1013 2001 -0.51514 21 1013 2002 -1.00070 22 1013 2003 -0.05914 23 1013 2004 0.01148 24 1013 2005 0.07212 25 1013 2006 0.04077 26 1013 2007 0.06023 27 1013 2008 -0.02181 28 1013 2009 -0.35301 29 1013 2010 0.04205        The error logs are attached. Any of your help will be highly appreciated.   Sincerely, Joon ... View more