关闭

基于R语言的模型组合

标签: R语言的模型组合
990人阅读 评论(0) 收藏 举报
分类:
组合预测模型---基于R语言的模型组合 
算术平均法、 最优权数法、 方差倒数法

模型中各参数的 t 值均显著 ,且通过 F 检验
和拟合忧度检验

c=c(1:20)
b=c-0.1
a=c-0.3

方差倒数法
e1=sum((c-b)^2)
e2=sum((c-a)^2)

w1=(1/e1)/(1/e1+1/e2)
w2=(1/e2)/(1/e1+1/e2)

x=w1*a+w2*b

构造损失函数
e1t=sum(c-b)
e2t=sum(c-a)
------------------推导
ft=w1*a+w2*b
et=c-w1*b-w2*a=w1*e1t+w2*e2t
j=et^2=(w1*e1t+w2*e2t)*(w1*e1t+w2*e2t)=w1^2*e1t^2+2*w1*w2*e1t*e2t+w2^2*e2t
=WT*En*W
WT=(w1,w2)--1xn矩阵
En=[ei*ej]---nxn矩阵
  e1t^2     e1t*e2t
  e1t*e2t   e2t^2
W=[w1  ---nx1矩阵
   w2]  
Rn=t(t(rep(1,2)))  nx1的全是1的矩阵   
拉格朗日乘数法
W*=(En^-1)*Rn/Rn^T*(En^-1)*Rn  En的逆矩阵*Rn/Rn的转置*En的逆矩阵*Rn
J*=1/Rn^T*(En^-1)*Rn   1//Rn的转置*En的逆矩阵*Rn
---------------------
En=matrix(rep(0,4),nrow=2,ncol=2)
En[1,1]=e1t^2
En[1,2]=e1t*e2t
En[2,1]=e1t*e2t
En[2,2]=e2t^2
Rn=t(t(rep(1,2)))
#solve(En) #矩阵的逆
library("MASS")
#ginv(En)  #矩阵的广义逆
#En%*%Rn等效crossprod(En,Rn) #矩阵的乘法
w=(ginv(En)%*%Rn)/(t(Rn)%*%ginv(En)%*%Rn)[1,1]
0
0
查看评论

LR模型常见问题小议

经常说的2/8原则,LR肯定就是能解决80%问题中那20%的工具。所以LR还是值得好好研究的。发现以前对LR重视不够,总想着赶紧把其他算法也学了,才能拉小跟同事之间机器学习的gap。其实LR用得还是挺多的,而且效果还是不错的。一些高大上的算法,在公司这种大数据面前不一定跑得动,即使跑得动,效果也不一...
  • Dinosoft
  • Dinosoft
  • 2016-01-10 23:33
  • 3209

组合预测模型

组合预测模型---基于R语言的模型组合  组合预测模型的普遍形式为各个单项预测模型的加权平均, 因此组合预测模型的重点在于加权系数的确定。如果对各个单项预测模型的加权系数赋值合理, 那么整个组合预测模型的预测精度也会相应提高。 目前常用的方法有算术平均法、 最优权数法、 方差倒数法等 ...
  • q383700092
  • q383700092
  • 2016-05-20 15:10
  • 1532

R语言-快速提取所有的组合

以 1, 2,3为例,要提取他们的所有组合的可能性,那么一共有C16+C26+C36C^1_6+C^2_6 + C^3_6 种情况,计算结果为:considerNum <- sapply(1:3, function(x) choose(3,x)) # 3 3 1 sum(considerNu...
  • qq_27755195
  • qq_27755195
  • 2016-12-08 16:35
  • 1051

基于opencpu构建R语言的RESTful API

http://xccds1977.blogspot.jp/2015/01/opencpurrestful-api.html 一个R语言模型的落地应用有几种方式:可能是一个简单的分析报告,可能是将预测结果输出到数据库,即离线模型,也可能需要实时读入数据进行实时预测输出。第三种情况是最有用,也最难处理...
  • u014032673
  • u014032673
  • 2015-03-13 11:06
  • 753

基于R语言构建的电影评分预测模型

电影评分系统是一种常见的推荐系统。现在使用R语言基于协同过滤算法来构建一个电影评分预测模型。 一,前提准备         1.R语言包:ggplot2包(绘图),recommenderlab包,reshape包(数据处理)     ...
  • wzgl__wh
  • wzgl__wh
  • 2016-08-05 02:25
  • 6002

基于R语言的用户征信行为分类预测模型搭建总结

数据集:15万行数据量,11个数据维度。目标变量为是否有重大违约倾向(定义为超过90天没有还贷) 建模流程: 1. 缺失值处理 发现数据集的数据缺失主要集中在亲属数量和月收入两个属性。其中,亲属数量却缺失属性为3000多条,且这3000多条记录中,月收入属性也同时缺失,故直接删除。 对于有1万多...
  • droprainbow
  • droprainbow
  • 2016-07-25 18:41
  • 3257

R语言实践:排列与组合

组合数:choose(n, k) 注意:从n个中选出k个 阶乘(k!):factorial(k) 排列数:choose(n, k) * factorial(k) 幂 : ^ 余数: %% 整数商:%/% == combn(x,n) -- 列出所有组合数矩阵 t(combn(x,n)) -- 将结果...
  • flyfrommath
  • flyfrommath
  • 2017-08-05 17:34
  • 979

基于R语言的多元线性回归--我国经济增长的定量研究

基于R语言的多元线性回归--我国经济增长的定量研究 一、实验介绍 1.1 实验内容 经济增长一直以来都是我国宏观经济政策的目标之一,研究影响经济增长的因素对促进我国经济快速发展有着重要意义。本实验运用 R 语言编写代码拟合多元线性回归模型,对模型拟合结果进行诊断,即对假设前提进行检验,并选择最...
  • oxuzhenyi
  • oxuzhenyi
  • 2017-08-10 12:53
  • 420

多因子模型之组合构建与优化器(上)

根据多因子模型,或者说alpha策略的开发顺序,我们应当是按照:因子--》alpha 模型--》风险模型--》组合构建 这样几个模块来的。今天来说说组合构建这个事。        组合构建是在你有了alpha模型和风险模型之后,也...
  • qtlyx
  • qtlyx
  • 2017-10-21 21:26
  • 792

R语言决策树分类模型

rm(list=ls()) gc() memory.limit(4000) library(corrplot) library(rpart) data_health<-read.csv("D:/smart_data0608/smart_data_section_good_15.tx...
  • DF_XIAO
  • DF_XIAO
  • 2015-06-10 15:23
  • 1683
    个人资料
    • 访问:198258次
    • 积分:3168
    • 等级:
    • 排名:第12865名
    • 原创:157篇
    • 转载:17篇
    • 译文:0篇
    • 评论:43条
    github地址
    https://github.com/lytforgood
    最新评论