- 博客(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 5530
原创 正态逆威沙特分布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 2695
原创 R:变点模型中区组划分及参数匹配
getblock=function(c,Q){#Q是分割点向量 blk=sum(c>Q)+1 return(blk)}
2020-03-05 12:00:49 443
原创 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 7500 1
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人