全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 经管代码库
494 0
2023-09-13
求助,该怎么处理呢# Initialization> rm(list=ls())> library(stochvol)> set.seed(1000) # for replication> options(digits=17)> vt   = read.table('vft.txt',sep = '\t')> T    = dim(vt)[1]> N    = dim(vt)[2]> for (i in 1:N){+         if(min(log(vt[,i]^2))== -Inf){+           vt[vt < 1e-10] <- 1e-5 + vt[vt < 1e-10] +         }+ }> # Run MCMC algorithm and store draws> S    = 50000> burn = 50000> m    = matrix(0,T+3,N)> g    = matrix(0,3,N)> for (i in 1:N){+         draws  = svsample(vt[,i],draws=S,burnin=burn,quiet=TRUE,thinpara=10,thinlatent=10)+         all    = cbind(draws$para,draws$latent)+         all    = cbind(draws$para[[1]][,-(4:5)],draws$latent[[1]])+         m[,i]  = colMeans(all)+         g[,i]  = geweke.diag(draws$para)$z+         name   = sprintf('svfdraws%d.txt',i)+ #        write(t(all),file=name,ncolumn=dim(all)[2])+ }Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) :   NA/NaN/Inf in 'y'

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群