*<PRE>;
%macro decile(y=,w=1,dsin=_last_,dout=,gout=,suffix=,print=1);                  
%* -------------------------------------------------------------------          
   Calculates decile shares etc for variable y with weight w.                   
   (gini coefficients are also calculated).                                     
   Special Note:  cases with weights which place them on the decile             
                  boundaries are split in order to provide exact decile         
                  shares.  However this does not deal with tied cases.          
                  If you want to allocate these you should add a small          
                  random variable to the y variable.                            
   Input:   y     variable whose distribution is to be described.               
            w     weight variable (default = 1).                                
            dsin  data set containing these variables (program is most          
                  efficient if no other variables are present)                  
                  (default = _last_).                                           
            dout  name of output dataset for deciles                            
            gout  name of output dataset for gini                               
            suffix Specify two letters (or numbers) to be appended to           
                  the output variable names.  This is useful if you             
                  want to combine output from separate runs of this             
                  program.                                                      
            print Specify 0 if you do not want printed output (def = 1)         
            Note: The program uses the library format decile.                   
    Output:       Decile boundaries, mean incomes, shares, cumulative           
                  shares and cumulative shares*mean.                            
                  Also gini coefficient.                                        
                  The dataset dsout will contain 11 cases representing          
                  the 10 deciles plus a total (placed last) with the            
                  following variables (where zz represents the suffix           
                  if specified),                                                
                    decile - decile 1...10 and 99 for the total.                
                    wcaseszz-number of weighted cases in decile.                
                    lbzz  -  lower bound of decile,                             
                    ubzz  -  upper bound                                        
                    meanzz - mean income in decile,                             
                    sharezz- income share of decile (in %)                      
                    csharezz-cumulative share (in %)                            
                    cmzz -  cumulative mean.                                    
                  In addition a data set with name gout will be                 
                  created with one case with the following variables,           
                    ginizz - gini coefficient                                   
                    meanzz, wcaseszz - as defined above,                        
                    nzz    - unweighted cases.                                  
    Author:       BWB Sept 90.                                                  
    Modifications BWB May 91, CM rather than CSM output.                        
  --------------------------------------------------------------------;         
                                                                                
* sort file                                                                     
  ---------;                                                                    
proc sort data=&dsin out=sorttmp;by &y;                                         
                                                                                
* Get the grand totals of the weight variable                                   
  -------------------------------------------;                                  
proc means noprint;                                                             
  var &w;                                                                       
  output out=gtotals sum=gtwgts;                                                
                                                                                
* Calculate decile rankings                                                     
  -------------------------;                                                    
data frank;                                                                     
  set sorttmp;                                                                  
  if _n_=1 then set gtotals;                                                    
                                                                                
  * check for missing cases;                                                    
  if &y=. | &w=. then error "missing values found";                             
                                                                                
  * initialise some variables;                                                  
  retain cs 0;                                                                  
  length  decile 2;                                                             
  retain cy 0;                                                                  
                                                                                
  cs = cs+&w;   /* cumulative sum of case weights */                            
  frank = cs/gtwgts;  /* fractional rank */                                     
  * determine the decile of the present case;                                   
  select;                                                                       
      when (frank<=0.10) decile = 1;  /* bottom decile */                       
      when (frank<=0.20) decile = 2;  /* next decile */                         
      when (frank<=0.30) decile = 3;                                            
      when (frank<=0.40) decile = 4;                                            
      when (frank<=0.50) decile = 5;                                            
      when (frank<=0.60) decile = 6;                                            
      when (frank<=0.70) decile = 7;                                            
      when (frank<=0.80) decile = 8;                                            
      when (frank<=0.90) decile = 9;                                            
      otherwise          decile = 10;  /* top decile    */                      
  end;                                                                          
                                                                                
  * Gini calculations;                                                          
  wy = &y * &w;    /* total income represented by present case */               
  cy = cy + wy;   /* cumulative income */                                       
  area = (cy - wy/2) * &w;    /* un-normalised lorenz area */                   
                                                                                
  original = 1;   /* flag for original cases */                                 
  length original 2;                                                            
                                                                                
  keep &y &w decile area wy original;                                           
                                                                                
  * adjustments for cases which straddle the quantile boundary;                 
  ldecile = lag(decile);                                                        
  if ldecile^=. & ldecile^=decile then do;                                      
    * i.e. just for the first case over the boundary;                           
    * find the fractional rank of the boundary just passed;                     
    select;                                                                     
        when (decile = 2) bound = .1         ;                                  
        when (decile = 3) bound = .2         ;                                  
        when (decile = 4) bound = .3;                                           
        when (decile = 5) bound = .4;                                           
        when (decile = 6) bound = .5;                                           
        when (decile = 7) bound = .6;                                           
        when (decile = 8) bound = .7;                                           
        when (decile = 9) bound = .8;                                           
        when (decile = 10) bound = .9;                                          
    end;                                                                        
    * create temporary variables;                                               
    tempw = &w;                                                                 
    tdecile = decile;                                                           
                                                                                
    * output a case below boundary;                                             
    decile = ldecile;  /* decile is value below boundary */                     
    &w = &w-(frank-bound)*gtwgts;                                               
    * e.g. if the boundary is 0.1, and the present cases fractional             
      rank is 0.1004, and there are a total of 1,000,000 weighted cases         
      and the weight of the present case is 1000, then this will split          
      into 1000-(0.1004-0.1)*1,000,000 = 400 cases below the boundary           
      and 600 above.;                                                           
    output;                                                                     
                                                                                
    * output case above quantile boundary;                                      
    decile = tdecile;  /* i.e. the original decile variable */                  
    &w = (frank-bound)*gtwgts;                                                  
    area = 0;  /* not needed a second time for ginis */                         
    wy = 0;                                                                     
    original = 0;  /* used in counting no of unweighted cases */                
    output;                                                                     
  end;  /* of cases straddling the boundary */                                  
                                                                                
  else output;   /* output other cases */                                       
                                                                                
run;                                                                            
                                                                                
* stats for deciles                                                             
  -----------------;                                                            
proc summary data=frank ;    /* this outputs one case per decile + */           
  weight &w;                 /* a case for the whole file          */           
  class decile;                                                                 
  var &y ;                                                                      
  output sumwgt(&y)=wcases&suffix mean(&y)=mean&suffix sum(&y)=totaly           
         min(&y)=lb&suffix max(&y)=ub&suffix                                    
        ;                                                                       
                                                                                
data ;                                                                          
  set _last_;                                                                   
  if _type_=0 then do;   /* ie grand total */                                   
    gtw=wcases&suffix;  /* grand total weighted cases */                        
    gty=totaly;  /* grand total weighted y */                                   
    gtm=mean&suffix;   /* overall mean y */                                     
    decile=99;                                                                  
  end;                                                                          
  retain gtw gty gtm;      /* grand totals to be used */                        
  share&suffix= 100*totaly/gty; /* income share of decile in %*/                
  retain cshare&suffix 0;                                                       
  cshare&suffix = share&suffix + cshare&suffix;                                 
                                       /* cumulative income share */            
  if _type_=0 then cshare&suffix = 0;/* reinitialise for totals case */         
  cm&suffix = (cshare&suffix*gtm/100)/(decile/10);                              
                                                                                
  keep _type_ decile share&suffix mean&suffix wcases&suffix                     
     cshare&suffix cm&suffix lb&suffix ub&suffix;                               
  * these are the variables to be found in the dataset dout;                    
                                                                                
* Put the totals last;                                                          
proc sort;                                                                      
  by decile;                                                                    
                                                                                
data &dout;                                                                     
  set _last_;                                                                   
  if decile=99 then do;                                                         
     cshare&suffix=.;                                                           
     cm&suffix=.;                                                               
     end;                                                                       
run;                                                                            
                                                                                
* print decile shares;                                                          
%if &print %then %do;                                                           
  proc print split='|';                                                         
    Title3 "Decile Shares for &y";                                              
    title4 "from data set &dsin with weight variable &w";                       
    format decile decile. /* this format is in the sasutil library */           
           lb&suffix ub&suffix comma9.0                                         
           mean&suffix comma8.0 share&suffix cshare&suffix 6.2                  
           cm&suffix comma8.0 wcases&suffix comma10.0;                          
    id decile;                                                                  
    var lb&suffix ub&suffix mean&suffix share&suffix cshare&suffix              
        cm&suffix wcases&suffix ;                                               
    label lb&suffix = 'Lower|Bound';                                            
    label ub&suffix = 'Upper|Bound';                                            
    label mean&suffix = 'Mean';                                                 
    label share&suffix = 'Share';                                               
    label cshare&suffix = 'Cumulative|Share';                                   
    label cm&suffix = 'Cumulative|Mean';                                        
    label wcases&suffix = 'Weighted|Cases';                                     
%end;                                                                           
                                                                                
                                                                                
* do gini calculation                                                           
  -------------------;                                                          
proc summary data=frank;   /* produces just one case */                         
  var area wy &w original;                                                      
  output sum(area)=tarea sum(wy)=ty                                             
         sum(&w)=sumwgt sum(original)=n&suffix;                                 
                                                                                
data &gout;                                                                     
  set _last_;                                                                   
  wcases&suffix = sumwgt;  /* grand total weighted cases */                     
  mean&suffix = ty/sumwgt;  /* mean  y */                                       
  gini&suffix = 1 - 2*(tarea/(ty*sumwgt));                                      
  keep gini&suffix mean&suffix wcases&suffix n&suffix ;                         
run;                                                                            
* print gini statistics;                                                        
%if &print %then %do;                                                           
  proc print;                                                                   
    title3 "Gini Coefficient for variable &y weighted by &w";                   
    title4 "from data set &dsin";                                               
    format gini&suffix 5.3 ;                                                    
    format mean&suffix comma10.0  ;                                             
    format wcases&suffix comma20.0;                                             
    format n&suffix comma8.0;                                                   
  run;                                                                          
%end;                                                                           
* reset the titles;                                                             
title3;                                                                         
                                                                                
%mend decile;
                                                                   
*</PRE>;

<!--'"</title></head>-->

<script type="text/javascript">
//OwnerIQ
var __oiq_pct = 50;
if( __oiq_pct>=100 || Math.floor(Math.random()*100/(100-__oiq_pct)) > 0 ) {
var _oiqq = _oiqq || [];
_oiqq.push(['oiq_addPageBrand','Lycos']);
_oiqq.push(['oiq_addPageCat','Internet > Websites']);
_oiqq.push(['oiq_addPageLifecycle','Intend']);
_oiqq.push(['oiq_doTag']);
(function() {
var oiq = document.createElement('script'); oiq.type = 'text/javascript'; oiq.async = true;
oiq.src = document.location.protocol + '//px.owneriq.net/stas/s/lycosn.js';
var s = document.getElementsByTagName('script')[0]; s.parentNode.insertBefore(oiq, s);
})();
}
//Google Analytics
var _gaq = _gaq || [];
_gaq.push(['_setAccount','UA-21402695-19']);
_gaq.push(['_setDomainName','tripod.com']);
_gaq.push(['_setCustomVar',1,'member_name','schick',3]);
_gaq.push(['_trackPageview']);
(function() {
var ga = document.createElement('script'); ga.type = 'text/javascript'; ga.async = true;
ga.src = ('https:' == document.location.protocol ? 'https://ssl' : 'http://www') + '.google-analytics.com/ga.js';
var s = document.getElementsByTagName('script')[0]; s.parentNode.insertBefore(ga, s);
})();
//Lycos Init
function getReferrer() {
var all= this.document.cookie;
if (all== '') return false;
var cookie_name = 'REFERRER=';
var start = all.lastIndexOf(cookie_name);
if (start == -1) return false;
start += cookie_name.length;
var end = all.indexOf(';', start);
if (end == -1) end = all.length;
return all.substring(start, end);
}
function getQuery() {
var rfr = getReferrer();
if (rfr == '') return false;
var q = extractQuery(rfr, 'yahoo.com', 'p=');
if (q) return q;
q = extractQuery(rfr, '', 'q=');
return q ? q : "";
}
function extractQuery(full, site, q_param) {
var start = full.lastIndexOf(site);
if (start == -1) return false;
start = full.lastIndexOf(q_param);
if (start == -1) return false;
start += q_param.length;
var end = full.indexOf('&', start);
if (end == -1) end = full.length;
return unescape(full.substring(start, end)).split(" ").join("+");
}
function generateHref(atag, template){
atag.href=template.replace('_MYURL_', window.location.href.replace('http://', '')).replace('_MYTITLE_','Check%20out%20this%20Tripod%20Member%20site!'); 
}
var lycos_ad = Array();
var lycos_onload_timer;
var cm_role = "live";
var cm_host = "tripod.lycos.com";
var cm_taxid = "/memberembedded";
var tripod_member_name = "schick";
var tripod_member_page = "schick/decile.sas";
var tripod_ratings_hash = "1714073027:3fa79e90a0bc1a5ada0e75fc2d3e22ff";

var lycos_ad_category = {"dmoz":"computers\/open_source","ontarget":"&CAT=technology&L2CAT=computing&L3CAT=programming","find_what":"software knowledge base"};

var lycos_ad_remote_addr = "209.202.244.9";
var lycos_ad_www_server = "www.tripod.lycos.com";
var lycos_ad_track_small = "";
var lycos_ad_track_served = "";
var lycos_search_query = getQuery();
</script>

<script type="text/javascript" src="https://scripts.lycos.com/catman/init.js"></script>

<script type='text/javascript'>
 var googletag = googletag || {};
 googletag.cmd = googletag.cmd || [];
 (function() {
   var gads = document.createElement('script');
   gads.async = true;
   gads.type = 'text/javascript';
   var useSSL = 'https:' == document.location.protocol;
   gads.src = (useSSL ? 'https:' : 'http:') +
     '//www.googletagservices.com/tag/js/gpt.js';
   var node = document.getElementsByTagName('script')[0];
   node.parentNode.insertBefore(gads, node);
 })();
</script>

<script type='text/javascript'>
 googletag.cmd.push(function() {
   googletag.defineSlot('/95963596/TRI_300X250_dfp', [300, 250], 'div-gpt-ad-1450204159126-0').addService(googletag.pubads());
   googletag.defineSlot('/95963596/TRI_above_728x90_dfp', [728, 90], 'div-gpt-ad-1450204159126-1').addService(googletag.pubads());
   googletag.defineSlot('/95963596/TRI_below_728x90_dfp', [728, 90], 'div-gpt-ad-1450204159126-2').addService(googletag.pubads());
   googletag.pubads().enableSingleRequest();
   googletag.enableServices();
 });
</script>


<script type="text/javascript"> 
(function(isV)
{
    if( !isV )
    {
        return;
    }
    var adMgr = new AdManager();
    var lycos_prod_set = adMgr.chooseProductSet();
    var slots = ["leaderboard", "leaderboard2", "toolbar_image", "toolbar_text", "smallbox", "top_promo", "footer2", "slider"];
    var adCat = this.lycos_ad_category;
    adMgr.setForcedParam('page', (adCat && adCat.dmoz) ? adCat.dmoz : 'member');
    if (this.lycos_search_query)
    {
        adMgr.setForcedParam("keyword", this.lycos_search_query);
    } 
    else if(adCat && adCat.find_what)
    {
        adMgr.setForcedParam('keyword', adCat.find_what);
    }
    
    for (var s in slots)
    {
        var slot = slots[s];
        if (adMgr.isSlotAvailable(slot))
        {
            this.lycos_ad[slot] = adMgr.getSlot(slot);
        }
    }

    adMgr.renderHeader();
    adMgr.renderFooter();
}((function() {

var w = 0, h = 0, minimumThreshold = 300;

if (top == self)
{
    return true;
}
if (typeof(window.innerWidth) == 'number' )
{
    w = window.innerWidth;
    h = window.innerHeight;
}
else if (document.documentElement && (document.documentElement.clientWidth || document.documentElement.clientHeight))
{
    w = document.documentElement.clientWidth;
    h = document.documentElement.clientHeight;
}
else if (document.body && (document.body.clientWidth || document.body.clientHeight))
{
    w = document.body.clientWidth;
    h = document.body.clientHeight;
}
return ((w > minimumThreshold) && (h > minimumThreshold));
}())));




window.onload = function()
{
    var f = document.getElementById("FooterAd");
    var b = document.getElementsByTagName("body")[0];
    b.appendChild(f);
    f.style.display = "block";
    document.getElementById('lycosFooterAdiFrame').src = '/adm/ad/footerAd.iframe.html';
    


    
    // DOM Inj Ad
    (function(isTrellix)
    {
        var e = document.createElement('iframe');
        e.style.border = '0';
        e.style.margin = 0;
        e.style.display = 'block';
        e.style.cssFloat = 'right';
        e.style.height = '254px';
        e.style.overflow = 'hidden';
        e.style.padding = 0;
        e.style.width = '300px';


        var isBlokedByDomain = function( href )
        {
            var blockedDomains = [
                "ananyaporn13000.tripod.com",
                "xxxpornxxx.tripod.com"
            ];
            var flag = false;
            
            for( var i=0; i<blockedDomains.length; i++ )
            {
                if( href.search( blockedDomains[ i ] ) >= 0 )
                {
                    flag = true;
                }
            }
            return flag;
        }

        var getMetaContent = function( metaName )
        {
            var metas = document.getElementsByTagName('meta');
            for (i=0; i<metas.length; i++)
            { 
                if( metas[i].getAttribute("name") == metaName )
                { 
                    return metas[i].getAttribute("content"); 
                } 
            }
            return false;
        }
        
        var getCommentNodes = function(regexPattern)
        {
            var nodes = {};
            var nodesA = [];
            var preferredNodesList = ['a', 'c', 'b'];
        
            (function getNodesThatHaveComments(n, pattern)
            {
                if (n.hasChildNodes())
                {
                    if (n.tagName === 'IFRAME')
                    {
                        return false;
                    }
                    for (var i = 0; i < n.childNodes.length; i++)
                    {
                        if ((n.childNodes[i].nodeType === 8) && (pattern.test(n.childNodes[i].nodeValue)))
                        {
                            var areaName = pattern.exec(n.childNodes[i].nodeValue)[1];
                            nodes[areaName] = n;
                        }
                        else if (n.childNodes[i].nodeType === 1)
                        {
                            getNodesThatHaveComments(n.childNodes[i], pattern);
                        }
                    }
                }
            }(document.body, regexPattern));

            for (var i in preferredNodesList)
            {
                if (nodes[preferredNodesList[i]])
                {
                    if( isTrellix && nodes[preferredNodesList[i]].parentNode.parentNode.parentNode.parentNode )
                    {
                        nodesA.push(nodes[preferredNodesList[i]].parentNode.parentNode.parentNode.parentNode);
                    }
                    else
                    {
                        nodesA.push( nodes[preferredNodesList[i]] );
                    }
                }
            }
            return nodesA;
        }
        
        
        var properNode = null;
        var areaNodes = getCommentNodes( new RegExp( '^area Type="area_(\\w+)"' ) );

        for (var i = 0; i < areaNodes.length; i++)
        {
            var a = parseInt(getComputedStyle(areaNodes[i]).width);
            if ((a >= 300) && (a <= 400))
            {
                properNode = areaNodes[i];
                break;
            }
        }


        var propertyName = getMetaContent("property") || false;
        if( isTrellix && (properNode) )
        {
            e.src = '/adm/ad/injectAd.iframe.html';
            properNode.insertBefore(e, properNode.firstChild);
        }
        else if( isTrellix && !( properNode ) ) // Slap the ad eventhought there is no alocated slot
        {
            e.src = '/adm/ad/injectAd.iframe.html';
            e.style.cssFloat = 'none';
            var cdiv = document.createElement('div');
            cdiv.style = "width:300px;margin:10px auto;";
            cdiv.appendChild( e );
            b.insertBefore(cdiv, b.lastChild);
        }
        else if( !isBlokedByDomain( location.href ) )
        {
            var injF = document.createElement('iframe');
            injF.style.border = '0';
            injF.style.margin = 0;
            injF.style.display = 'block';
            injF.style.cssFloat = 'none';
            injF.style.height = '254px';
            injF.style.overflow = 'hidden';
            injF.style.padding = 0;
            injF.style.width = '300px';
            injF.src = '/adm/ad/injectAd.iframe.html';

            if( b && ( !isTrellix || ( typeof isTrellix == "undefined" ) ) ) // All other tripod props
            {
                var cdiv = document.createElement('div');
                cdiv.style = "width:300px;margin:10px auto;";
                cdiv.appendChild( injF );
                b.insertBefore(cdiv, b.lastChild);
            } 
        }
  }( document.isTrellix ));
}

</script>

<div id="tb_container" style="background:#DFDCCF; border-bottom:1px solid #393939; position:relative; z-index:999999999!important">
<!--form name="search" onSubmit="return searchit()" id='header_search' >
<input type="text" placeholder="Search" size=30 name="search2" value="">
<input type="button" value="Go!" onClick="searchit()">
</form>
<style>
form#header_search {
    width: 916px;
    margin: 0 auto 8px;
    position: relative;
}


form#header_search input {
    height: 40px;
    font-size: 14px;
    line-height: 40px;
    padding: 0 8px;
    box-sizing: border-box;
    background: #F4F2E9;
    border: 1px solid #BBB8B8;
    transition: background-color 300ms ease-out,
                color 300ms ease;
}

form#header_search input[type="text"] {
    width: 100%;
}
form#header_search input[type="text"]:focus {
    border-color: #A2D054;
    background-color: #fff;
    box-shadow: 0 0px 12px -4px #A2D054;
}



form#header_search input[type="button"] {
    position: absolute;
    top: 1px;
    right: 1px;
    opacity: 1;
    background: #DFDCCF;
    color: #463734;
    width: 125px;
    cursor: pointer;
    height: 38px;
    border: none;
}
form#header_search input[type="text"]:focus ~ input[type='button']:hover,
form#header_search input[type='button']:hover {
    background-color: #A5CE56;
    color: #fff;
}
form#header_search input[type="text"]:focus ~ input[type='button'] {
    background-color: #52AEDF;
    color: #fff;
}

</style>

<script>
function searchit(){
    
    // determine environment 
    var search_env 
    if (lycos_ad_www_server.indexOf(".pd.") > -1) {
       search_env = 'http://search51.pd.lycos.com/a/';
    } else if (lycos_ad_www_server.indexOf(".qa.") > -1) {
       search_env = 'http://search51.qa.lycos.com/a/';
    } else {
       search_env = 'http://search51.lycos.com/a/';
    }

var search_term = encodeURIComponent(document.search.search2.value)
var search_url = search_env+search_term;
window.open(search_url);

return false
}
</script-->
<style>
    .adCenterClass{margin:0 auto}
</style>
<div id="tb_ad" class="adCenterClass" style="display:block!important; overflow:hidden; width:916px;">

<div id="ad_container" style="display:block!important; float:left; width:728px ">
<script type="text/javascript">
if (typeof lycos_ad !== "undefined" && "leaderboard" in lycos_ad) {
  document.write(lycos_ad['leaderboard']);
}
</script>
</div>
</div>
</div>
<script type="text/javascript">
if (typeof lycos_ad !== "undefined" && "slider" in lycos_ad) {
  document.write(lycos_ad['slider']);
}
</script> <!-- added 7/22 -->
<div id="FooterAd" style="background:#DFDCCF; border-top:1px solid #393939; clear:both; display:none; width:100%!important; position:relative; z-index:999999!important; height:90px!important"> 
<div class="adCenterClass" style="display:block!important; overflow:hidden; width:916px;">
<div id="footerAd_container" style="display:block!important; float:left; width:728px">
<iframe id="lycosFooterAdiFrame" style="border:0; display:block; float:left; height:96px; overflow:hidden; padding:0; width:750px"></iframe>
</div>
</div>
</div>