自定义博客皮肤VIP专享

*博客头图:

格式为PNG、JPG,宽度*高度大于1920*100像素,不超过2MB,主视觉建议放在右侧,请参照线上博客头图

请上传大于1920*100像素的图片!

博客底图:

图片格式为PNG、JPG,不超过1MB,可上下左右平铺至整个背景

栏目图:

图片格式为PNG、JPG,图片宽度*高度为300*38像素,不超过0.5MB

主标题颜色:

RGB颜色,例如:#AFAFAF

Hover:

RGB颜色,例如:#AFAFAF

副标题颜色:

RGB颜色,例如:#AFAFAF

自定义博客皮肤

-+
  • 博客(4)
  • 收藏
  • 关注

原创 R:多元正态密度函数计算

R中计算多元正态密度函数可以使用mvtnorm包中的dmvnorm函数。dvnorm(x,mu,sigma,log=T/F)x为服从多元正态分布的随机向量mu和sigma分别为多元正态均值和协方差阵(可以是对角阵或者非对角阵)log=T时,对密度取自然对数当x为维时,输出结果与dnorm相同可以同时输出多个多元正态密度:当二维随机向量服从以mu为均值,sigma为协方差阵的多元正态...

2020-03-18 17:00:21 5444

原创 正态逆威沙特分布NIW

在研究变点模型时,假设每个区组内的队列观测都独立同分布于该区组的正态分布(每个区组内的正态参数可能不同,用omega存储所有区组的正态参数)。如:共有10个队列,分为3个区组,变点位置Q1=3,Q2=8,则队列1,2属于区组1,队列3-7属于区组2,队列8-10属于区组3——x1,x2 iid服从于N(mu1,var1),x3-x7 iid服从于N(mu2,var2), x8-x10 iid服从...

2020-03-08 16:51:40 2623

原创 R:变点模型中区组划分及参数匹配

getblock=function(c,Q){#Q是分割点向量 blk=sum(c>Q)+1 return(blk)}

2020-03-05 12:00:49 431

原创 R选择国内镜像安装程序包

R使用install.packages 安装程序包时,由于默认选择国外CRAN镜像,会导致下载程序包特别慢或者由于无法连接国外镜像而安装失败Warning in install.packages : unable to access index for repository https://cran.rstudio.com/src/contrib: cannot open URL 'ht...

2020-03-04 12:00:01 7305 1

空空如也

空空如也

TA创建的收藏夹 TA关注的收藏夹

TA关注的人

提示
确定要删除当前文章?
取消 删除