stata 计量数据标准化(min-max, z-scores)_24Mar2020

重要性:数据标准化是消除量纲影响的重要手段,构造新指标或应用某些算法前需要使用到。

写作理由:在 R 或者 Python 中,即使不使用第三方包,实现标准化也是非常容易的,但在 stata 中就显得不太顺手(仅我个人手感而言,并不代表我否认stata是优秀的软件),
本人也没多去寻找现存的包,顺手写了一个小工具,分享给需要的人。

可实现:z-score标准化,minmax标准化

* by Rujiang
* 19Dec2019
capture program drop stand
program stand
    version 14.1
    syntax varlist [if] [in] [, prefix(name) methods(name)]
    if "`methods'" == "" {
    
        if "`prefix'" == "" {
            local prefix "norm_"
        }
        
        foreach var in `varlist' {
            qui: summarize `var' `if' `in'
            qui: generate `prefix'`var' = (`var'- r(mean)) / r(sd) `if' `in'
            qui: label variable `prefix'`var' "normalized `var'"
            di "====> `prefix'`var' generated"
        }
            
        di "-----------xxxxxx-----------"
        di "all the variables you specified are **normalized**"
    }
    
    
    if "`methods'" == "minmax" {
    
            if "`prefix'" == "" {
                local prefix "minmax_"
            }
        
        foreach var in `varlist' {
            qui: summarize `var' `if' `in'
            qui: generate `prefix'`var' = (`var'- r(min)) / (r(max) - r(min)) `if' `in'
            qui: label variable `prefix'`var' "minmax standardized `var'"
            di "====> `prefix'`var' generated"
        }
            
        di "-----------xxxxxx-----------"
        di "all the variables you specified are **minmax standardized**"
    }

    
    if "`methods'" != "" {
        if "`methods'" != "minmax" {
            di as err `"nothing happened, only methods("minmax") allowed"'
        }
    }
    
end


使用举例


. clear

. sysuse auto
(1978 Automobile Data)

. des

Contains data from D:\来自华硕_RJL_25Apr2018\stata\StataCorp Stata 14.1 (Revision
>  20.06.2016)\ado\base/a/auto.dta
  obs:            74                          1978 Automobile Data
 vars:            12                          13 Apr 2014 17:45
 size:         3,182                          (_dta has notes)
---------------------------------------------------------------------------------
              storage   display    value
variable name   type    format     label      variable label
---------------------------------------------------------------------------------
make            str18   %-18s                 Make and Model
price           int     %8.0gc                Price
mpg             int     %8.0g                 Mileage (mpg)
rep78           int     %8.0g                 Repair Record 1978
headroom        float   %6.1f                 Headroom (in.)
trunk           int     %8.0g                 Trunk space (cu. ft.)
weight          int     %8.0gc                Weight (lbs.)
length          int     %8.0g                 Length (in.)
turn            int     %8.0g                 Turn Circle (ft.)
displacement    int     %8.0g                 Displacement (cu. in.)
gear_ratio      float   %6.2f                 Gear Ratio
foreign         byte    %8.0g      origin     Car type
---------------------------------------------------------------------------------
Sorted by: foreign


. stand price
====> norm_price generated
-----------xxxxxx-----------
all the variables you specified are **normalized**

. stand price, methods("minmax")
====> minmax_price generated
-----------xxxxxx-----------
all the variables you specified are **minmax standardized**


. stand price-head, methods("minmax") prefix("stand_")
====> stand_price generated
====> stand_mpg generated
====> stand_rep78 generated
====> stand_headroom generated
-----------xxxxxx-----------
all the variables you specified are **minmax standardized**

. su stand*

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
 stand_price |         74    .2278444    .2338086          0          1
   stand_mpg |         74    .3205965    .1995001          0          1
 stand_rep78 |         69    .6014493    .2474831          0          1
stand_head~m |         74    .4266409    .2417128          0          1


最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 216,240评论 6 498
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 92,328评论 3 392
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 162,182评论 0 353
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,121评论 1 292
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,135评论 6 388
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,093评论 1 295
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,013评论 3 417
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,854评论 0 273
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,295评论 1 310
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,513评论 2 332
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,678评论 1 348
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,398评论 5 343
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,989评论 3 325
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,636评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,801评论 1 268
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,657评论 2 368
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,558评论 2 352

推荐阅读更多精彩内容