南有嘉鱼 | 南有樛木

Posted by Leung ZhengHua on 2017-11-26

多元统计

[文献阅读] Bayesian adaptive Lasso quantile regression

是的，我又来读文献了。找了两篇贝叶斯Lasso相关的文章，是一个大佬Rahim Alhamzawi写的，他还贡献了一个R包Brq，对应他论文所用的方法。三位作者分别是Rahim Alhamzawi， Keming Yu and Dries F Benoit，我发誓我只打一次他们的名字。理论分位数回归分位数回归模型可以表示如下： y_i=\beta_0+x_i'\beta+\epsil......

Posted by Leung ZhengHua on 2017-11-26

文献阅读

[R语言] 利用R合成卡农和最炫民族风

这是一次严肃的电子音乐实验

程序代码参考网上的MATLAB代码，由于流传已久，找不到原作者，仅引出两篇相对整洁的资源帖另外，可以上github搜matlab的钢琴gui程序。我想在R里面自动播放音乐，那么在等运行结果的时候就可以听到提示音然而调出来的是MediaPlayer，我不是很想看到一个新的窗口出现我想隐式在后台自动响起还有，R调用MediaPlayer的时候会一直等待MediaPlayer进程的结束才......

Posted by Leung ZhengHua on 2017-11-18

R

[Manual] Rcpp使用手册2

这是一篇关于RcppArmadillo的基础性技术文档

如果你是第一次阅读这篇文章，以下资源贴可能是你感兴趣的： RcppArmadillo简明手册Rcpp相关知识整理RcppArmadillo官方文档Armadillo C++ linear algebra library 学习笔记（4）——矩阵的运算Armadillo矩阵库的使用（二）之API接口 .cpp文件的储存位置当我们在Rstudio中新建C++文件时，它会直接新建一个文档，方便我们......

Posted by Leung ZhengHua on 2017-11-16

Rcpp R

[R语言] 一种特殊的lasso求解

这是一篇关于提高计算速度的技术文档

最近开始忙的要死，一大堆东西堆在一起但又不想做的感觉。今晚研究了Rcpp发现又有新的收获，Rcpp用来做一些简单的循环和判断还可以，涉及矩阵运算时难道还要逐行逐列相乘相加吗？想想都要放弃了，直到我遇上了RcppArmadillo包。它是C++里面专门做科学运算的，用法类似Matlab。有了这个包，我终于找到了一个更有效更容易提高计算速度的办法，在尊贵的阿苏斯上的测试对比发现，Rcpp:Mat......

Posted by Leung ZhengHua on 2017-11-16

R

[R语言] apply函数

这是一篇关于R语言科学计算的技术文档

最近在使用apply函数，想要往函数FUN里面输入多个参数再调用apply，下面是实验的结果： apply1234567ryan=function(x,y,msg) { return(paste0(x,y,msg))}ryan(5,6,"-hello")apply(matrix(1:6,ncol = 1), 1, ryan,y=1:8,msg=c("-GOOD MORI......

Posted by Leung ZhengHua on 2017-11-15

R

[文献阅读] Multi-style rotation strategies using augmented Black Litterman factor model

对Black Litterman模型一无所知的可以参考这篇文档文章作者是Ginanjar Dewandaru、Rumi Masih、Obiyathulla Ismath Bacha和A. Mansur. M. Masih，多风格轮动策略采用了3种风格： momentum value quality 看了一会发现出现很多Islamic相关的宗教词汇，我还以为这篇文章是要宣扬伊斯兰教义。查......

Posted by Leung ZhengHua on 2017-11-12

文献阅读

[R语言] RStudio美化

这是一篇关于RStudio美化的小技巧

help page 输出文档美化这个美化文档的最初想法来自Raining_wcc，后来在他的CSS上大刀阔斧地又调出了符合本人审美的CSS，一个特点就是字体比较大。将下列文件保存为R.css文件，替换到Rstudio安装文件夹下的RStudio\resources\R.css。123456789101112131415161718192021222324252627282930313233......

Posted by Leung ZhengHua on 2017-11-11

R RStudio

[计量] 异方差性

这是一篇关于异方差性的介绍

定义同方差意为方差具有相同的分布，或者说相同的方差。 $E(u_i^2)=\sigma^2,i=1,2,…,n$ 不管变量$X_i$取什么值，$Y_i$的条件方差都保持不变。异方差表明$Y_i$的方差不再保持不变，会随着$X_i$的变化而变化,符号上记为： $E(u_i^2)=\sigma_i^2,i=1,2,…,n$ OLS估计记双变量模型为：$Y_i=\beta_1+\beta......

Posted by Leung ZhengHua on 2017-11-09

计量

[Statistics] 多元正态分布

这是一篇关于多元正态分布的重要提醒

一元正态分布f(x)=\frac{1}{\sqrt{2\pi\sigma^2}}\exp\{-\frac{(x-\mu)^2}{2\sigma^2}\},-\infty...