原创

MADlib——基于SQL的数据挖掘解决方案(12)——回归之广义线性模型

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://wxy0327.blog.csdn.net/article/details/79205296

一、广义线性模型简介

        在一些实际问题中,变量间的关系并不都是线性的,这种情况就应该用曲线去进行拟合。用曲线拟合数据首先要解决的问题是回归方程中的参数如何估计。下面以一元非线性回归为例,讨论解决这一问题的基本思路。

        对于曲线回归建模的非线性目标函数,通过某种数学变换,使之“线性化”为一元线性函数的形式,继而利用线性最小二乘估计的方法估计出参数ab,用一元线性回归方程来描述vu间的统计规律性,然后再用逆变换还原为目标函数形式的非线性回归方程。

        比如,对于指数函数,令,则。通过这样的形式,就可以将一些非线性函数转化为线性函数,这样就可以利用线性回归方法进行回归。因此这种转换模型又称为广义线性模型(Generalized Linear Models)。

        广义线性模型是一般线性模型的直接扩展,它使因变量的总体均值通过一个非线性连接函数(link function,如上例中的ln),而依赖于线性预测值,同时还允许响应概率分布为指数分布族中的任何一员。

        广义线性模型在两个方面对普通线性模型进行了扩展:

  • 一般线性模型中要求因变量是连续的且服从正态分布。在广义线性模型中,因变量的分布可扩展到非连续的,如二项分布、泊松分布、负二项分布等。
  • 一般线性模型中,自变量的线性预测值就是因变量的估计值,而广义线性模型中,自变量的线性预测值是因变量的连接函数估计值。

 

        一个广义线性模型包括以下三个组成部分:

  • 线性成分(linear component):
  • 随机成分(random component):
  • 连接函数(link function):

 

        MADlib1.10.0实现的分布族及其相应的连接函数如表1所示。

分布族

连接函数

二项分布(Binomial

logit, probit

伽马分布(Gamma)

inverse, identity, log

高斯分布(Gaussian)

identity, inverse, log

逆高斯分布(Inverse Gaussian)

inverse of square, inverse, identity, log

泊松分布(Poisson)

log, identity, square-root

表1 MADlib 1.10.0支持的连接函数

 

二、MADlib广义线性模型相关函数

1.  培训函数

(1)  语法

glm(source_table,  
    model_table,  
    dependent_varname,  
    independent_varname,  
    family_params,  
    grouping_col,  
    optim_params,  
    verbose  
    )

(2)  参数

参数名称

数据类型

描述

source_table

VARCHAR

包含训练数据的源表名。

model_table

VARCHAR

包含模型的输出表名。主输出表列和概要输出表列如表3、4所示。

dependent_varname

VARCHAR

训练数据中因变量列的名称。

independent_varname

VARCHAR

评估使用的自变量的表达式列表,一般显式地由包括一个常数1项的自变量列表提供。

family_params(可选)

VARCHAR

当前支持如下分布族参数:

family=poisson并且link=[log|identity|sqrt]。

family=gaussian并且link=[identity|log|inverse]。当family=gaussian并且link=identity时,GLM模型等同于线性回归。

family=gamma并且link=[inverse|identity|log]。

family=inverse_gaussian并且link=[sqr_inverse|log|identity|inverse]。

family=binomial并且link=[probit|logit]。

grouping_col(可选)

VARCHAR

缺省值为NULL。和SQL中的“GROUP BY”类似,是一个将输入数据集分成离散组的表达式,每个组运行一个回归。此值为NULL时,将不使用分组,并产生一个单一的结果模型。

optim_params(可选)

VARCHAR

优化相关参数,缺省值为‘max_iter=100,optimizer=irls,tolerance=1e-6’。

verbose(可选)

BOOLEAN

缺省值为FALSE,指定是否输出训练结果的详细信息。

表2 glm函数参数说明

 

列名

数据类型

描述

<...>

TEXT

分组列,取决于grouping_col输入,可能是多个列。

coef

FLOAT8[]

线性预测的回归系数向量。

log_likelihood

FLOAT8

对数似然值l(β)。训练函数使用离散度参数的极大似然估算值计算对数似然值。

std_err

FLOAT8[]

系数的标准方差向量。

z_stat或t_stats

FLOAT8[]

系数的z-统计量向量(泊松分布或二项分布)或t-统计量向量(其它分布)。

p_values

FLOAT8[]

系数的P值向量。

dispersion

FLOAT8

Pearson估计的离散度。当family=poisson或family=binomial时,离散度总是1。

num_rows_processed

BIGINT

处理的数据行数。

num_missing_rows_skipped

BIGINT

训练时由于缺少值或失败跳过的行数。

num_iterations

INTEGER

实际完成的迭代次数。

表3 glm函数主输出表列说明

 

        训练函数在产生输出表的同时,还会创建一个名为<out_table>_summary的概要表,具有以下列:

列名

数据类型

描述

Method

VARCHAR

'glm'

source_table

VARCHAR

源数据表名称。

model_table

VARCHAR

模型表名。

dependent_varname

VARCHAR

因变量名。

independent_varname

VARCHAR

自变量名。

family_params

TEXT

分布族相关参数,形式为'family=..., link=...'

grouping_col

TEXT

分组列名称。

optimizer_params

TEXT

包含优化参数的字符串,形式为'optimizer=..., max_iter=..., tolerance=...'。

num_all_groups

INTEGER

分组数。

num_failed_groups

INTEGER

失败的分组数。

total_rows_processed

BIGINT

所有分组处理的总行数。

total_rows_skipped

BIGINT

由于缺失值或失败跳过的总行数。

表4 glm函数概要输出表列说明

 

2.  预测函数

(1)  语法

glm_predict(coef,
            col_ind_var
            link)

(2)  参数

  • coef:FLOAT8[]类型,训练模型获得的回归系数向量。
  • col_ind_var:FLOAT8[]类型,包含自变量列名索引的数组,应该与训练函数中‘independent_varname’参数得到的数组具有相同的数组长度。
  • link:TEXT类型,连接函数字符串,应该与训练函数使用相同的连接函数。

 

三、示例

1.      问题提出

        为了了解百货商店销售额x与流通率(这是反映商业活动的一个质量指标,指每元商品流转额所分摊的流通费用)y之间的关系,收集了九个商店的有关数据(见表5)。假设数据服从伽马分布,请建立销售额与流通率之间的回归模型。

样本点

x-销售额(万元)

y-流通费率(%)

1

1.5

7.0

2

4.5

4.8

3

7.5

3.6

4

10.5

3.1

5

13.5

2.7

6

16.5

2.5

7

19.5

2.4

8

22.5

2.3

9

25.5

2.2

表5 销售额与流通费率数据

 

2. 建立表并添加数据

drop table if exists t1;    
create table t1 (x float, y float);    
insert into t1 values    
(1.5,7.0), (4.5,4.8), (7.5,3.6), (10.5,3.1), (13.5,2.7),
(16.5,2.5), (19.5,2.4), (22.5,2.3), (25.5,2.2);

 

3.      使用identity连接函数

(1)  训练生成模型

drop table if exists t1_glm, t1_glm_summary;    
select madlib.glm( 't1',    
                   't1_glm',    
                   'y',    
                   'array[1, x]',  
                   'family=gamma, link=identity'  
                  ); 

(2)  查询结果

\x on  
select * from t1_glm; 

        结果:

-[ RECORD 1 ]------+-------------------------------------------
coef               | {5.032482978893,-0.124291210170803}
log_likelihood     | -8.09384975775831
std_err            | {0.594336150129185,0.0304615709894519}
t_stats            | {8.46740178567152,-4.08026264350718}
p_values           | {6.33195030885822e-05,0.00468813243359916}
dispersion         | 0.0500321476619188
num_rows_processed | 9
num_rows_skipped   | 0
num_iterations     | 11

(3)  利用预测函数估计残差

\x off  
select x,y,predict,y-predict residual   
  from (select t1.*,  
               madlib.glm_predict(coef,array[1, x]::float8[],'identity') as predict  
          from t1, t1_glm m) t;  
        结果:

  x   |  y  |     predict      |      residual      
------+-----+------------------+--------------------
  1.5 |   7 | 4.84604616363679 |   2.15395383636321
  4.5 | 4.8 | 4.47317253312439 |  0.326827466875614
  7.5 | 3.6 | 4.10029890261198 | -0.500298902611978
 10.5 | 3.1 | 3.72742527209957 |  -0.62742527209957
 13.5 | 2.7 | 3.35455164158716 | -0.654551641587162
 16.5 | 2.5 | 2.98167801107475 | -0.481678011074754
 19.5 | 2.4 | 2.60880438056235 | -0.208804380562346
 22.5 | 2.3 | 2.23593075004994 | 0.0640692499500619
 25.5 | 2.2 | 1.86305711953753 |   0.33694288046247
(9 rows)

 

4.      使用inverse连接函数

(1)  训练生成模型

drop table if exists t1_glm, t1_glm_summary;    
select madlib.glm( 't1',    
                   't1_glm',    
                   'y',    
                   'array[1, x]',  
                   'family=gamma, link=inverse'  
                  ); 

(2)  查询结果

\x on  
select * from t1_glm;

        结果:

-[ RECORD 1 ]------+--------------------------------------------
coef               | {0.137708424903124,0.0147719820924827}
log_likelihood     | -1.12823356728084
std_err            | {0.0137829792192191,0.0013221473412605}
t_stats            | {9.9911944081801,11.1727200376923}
p_values           | {2.15188055884168e-05,1.02542310032808e-05}
dispersion         | 0.00982968278042609
num_rows_processed | 9
num_rows_skipped   | 0
num_iterations     | 6

(3)  利用预测函数估计残差

\x off  
select x,y,predict,y-predict residual   
  from (select t1.*,  
               madlib.glm_predict(coef,array[1, x]::float8[],'inverse') as predict  
          from t1, t1_glm m) t; 

        结果:

  x   |  y  |     predict      |      residual       
------+-----+------------------+---------------------
  1.5 |   7 | 6.25522318791614 |   0.744776812083859
  4.5 | 4.8 | 4.89758310560006 | -0.0975831056000596
  7.5 | 3.6 | 4.02417255104088 |  -0.424172551040882
 10.5 | 3.1 | 3.41513449166629 |  -0.315134491666294
 13.5 | 2.7 |  2.9662132018307 |  -0.266213201830703
 16.5 | 2.5 | 2.62160216829753 |  -0.121602168297529
 19.5 | 2.4 | 2.34872962402895 |  0.0512703759710513
 22.5 | 2.3 | 2.12730643262039 |   0.172693567379613
 25.5 | 2.2 | 1.94403523699905 |   0.255964763000946
(9 rows)

 

5.      使用log连接函数

(1)  训练生成模型

drop table if exists t1_glm, t1_glm_summary;    
select madlib.glm( 't1',    
                   't1_glm',    
                   'y',    
                   'array[1, x]',  
                   'family=gamma, link=log'  
                  );

(2)  查询结果

\x on  
select * from t1_glm; 

        结果:

-[ RECORD 1 ]------+--------------------------------------------
coef               | {1.75614873826444,-0.0442718284575659}
log_likelihood     | -5.67230519309458
std_err            | {0.111953106344754,0.00719290385898647}
t_stats            | {15.6864672683264,-6.15493121074527}
p_values           | {1.03559068682089e-06,0.000465286458020498}
dispersion         | 0.0279384475992961
num_rows_processed | 9
num_rows_skipped   | 0
num_iterations     | 7

(3)  利用预测函数估计残差

\x off  
select x,y,predict,y-predict residual   
  from (select t1.*,  
               madlib.glm_predict(coef,array[1, x]::float8[],'log') as predict  
          from t1, t1_glm m) t; 

        结果:

  x   |  y  |     predict      |      residual       
------+-----+------------------+---------------------
  1.5 |   7 | 5.41807721742676 |    1.58192278257324
  4.5 | 4.8 |  4.7442127661024 |  0.0557872338976022
  7.5 | 3.6 | 4.15415909866612 |   -0.55415909866612
 10.5 | 3.1 | 3.63749238658367 |  -0.537492386583666
 13.5 | 2.7 | 3.18508524786705 |  -0.485085247867048
 16.5 | 2.5 | 2.78894550366558 |  -0.288945503665581
 19.5 | 2.4 | 2.44207498924096 | -0.0420749892409602
 22.5 | 2.3 | 2.13834592509533 |   0.161654074904671
 25.5 | 2.2 | 1.87239266423715 |   0.327607335762848
(9 rows)

 

        从自变量x的标准差、t统计值和P值来看,在符合伽马分布的假设下,inverse连接函数的拟合程度最高,其次是log连接函数,而误差最大的是identity连接函数。与madlib.linregr_train线性回归训练函数不同,madlib.glm不返回R2决定系数,而是用对数似然值评估模型的拟合程度。统计学中,似然函数是一种关于统计模型参数的函数。给定输出x时,关于参数θ的似然函数L(θ|x)(在数值上)等于给定参数θ后变量X的概率:L(θ|x)=P(X=x|θ)。涉及到似然函数的许多应用中,更方便的是使用似然函数的自然对数形式,即“对数似然函数”。求解一个函数的极大化往往需要求解该函数的关于未知参数的偏导数。由于对数函数是单调递增的,而且对数似然函数在极大化求解时较为方便,所以对数似然函数常用在最大似然估计及相关领域中。

        在该例子中,伽马分布的三种连接函数identity、inverse和log对应的似然函数值分别为-8.0938、-1.1282和-5.6723。对数似然函数绝对值越小,表示拟合程度越好。因此从结果来看,可以认为inverse连接函数形式更符合因变量与自变量之间的关系。

 


文章最后发布于: 2018-01-30 14:38:44
展开阅读全文
0 个人打赏
私信求帮助

没有更多推荐了,返回首页

©️2019 CSDN 皮肤主题: 大白 设计师: CSDN官方博客

分享到微信朋友圈

×

扫一扫,手机浏览