oracle 求解多元线性回归方程

Question

求回归系数

Answer 1

REGR_ (Linear Regression) Functions 
功能描述：这些线性回归函数适合最小二乘法回归线，有9个不同的回归函数可使用。 
REGR_SLOPE：返回斜率，等于COVAR_POP(expr1, expr2) / VAR_POP(expr2) 
REGR_INTERCEPT：返回回归线的y截距，等于 
AVG(expr1) - REGR_SLOPE(expr1, expr2) * AVG(expr2) 
REGR_COUNT：返回用于填充回归线的非空数字对的数目 
REGR_R2：返回回归线的决定系数，计算式为： 
If VAR_POP(expr2) = 0 then return NULL 
If VAR_POP(expr1) = 0 and VAR_POP(expr2) != 0 then return 1 
If VAR_POP(expr1) > 0 and VAR_POP(expr2 != 0 then 
return POWER(CORR(expr1,expr),2) 
REGR_AVGX：计算回归线的自变量(expr2)的平均值，去掉了空对(expr1, expr2)后，等于AVG(expr2) 
REGR_AVGY：计算回归线的应变量(expr1)的平均值，去掉了空对(expr1, expr2)后，等于AVG(expr1) 
REGR_SXX： 返回值等于REGR_COUNT(expr1, expr2) * VAR_POP(expr2) 
REGR_SYY： 返回值等于REGR_COUNT(expr1, expr2) * VAR_POP(expr1) 
REGR_SXY: 返回值等于REGR_COUNT(expr1, expr2) * COVAR_POP(expr1, expr2) 
（下面的例子都是在SH用户下完成的） 
SAMPLE 1：下例计算1998年最后三个星期中两种产品（260和270）在周末的销售量中已开发票数量和总数量的累积斜率和回归线的截距 
SELECT t.fiscal_month_number "Month", t.day_number_in_month "Day", 
REGR_SLOPE(s.amount_sold, s.quantity_sold) 
OVER (ORDER BY t.fiscal_month_desc, t.day_number_in_month) AS CUM_SLOPE, 
REGR_INTERCEPT(s.amount_sold, s.quantity_sold) 
OVER (ORDER BY t.fiscal_month_desc, t.day_number_in_month) AS CUM_ICPT 
FROM sales s, times t 
WHERE s.time_id = t.time_id 
AND s.prod_id IN (270, 260) 
AND t.fiscal_year=1998 
AND t.fiscal_week_number IN (50, 51, 52) 
AND t.day_number_in_week IN (6,7) 
ORDER BY t.fiscal_month_desc, t.day_number_in_month; 
Month Day CUM_SLOPE CUM_ICPT 
---------- ---------- ---------- ---------- 
12 12 -68 1872 
12 12 -68 1872 
12 13 -20.244898 1254.36735 
12 13 -20.244898 1254.36735 
12 19 -18.826087 1287 
12 20 62.4561404 125.28655 
12 20 62.4561404 125.28655 
12 20 62.4561404 125.28655 
12 20 62.4561404 125.28655 
12 26 67.2658228 58.9712313 
12 26 67.2658228 58.9712313 
12 27 37.5245541 284.958221 
12 27 37.5245541 284.958221 
12 27 37.5245541 284.958221 
SAMPLE 2：下例计算1998年4月每天的累积交易数量 
SELECT UNIQUE t.day_number_in_month, 
REGR_COUNT(s.amount_sold, s.quantity_sold) 
OVER (PARTITION BY t.fiscal_month_number ORDER BY t.day_number_in_month) 
"Regr_Count" 
FROM sales s, times t 
WHERE s.time_id = t.time_id 
AND t.fiscal_year = 1998 AND t.fiscal_month_number = 4; 
DAY_NUMBER_IN_MONTH Regr_Count 
------------------- ---------- 
1 825 
2 1650 
3 2475 
4 3300