lm.fit求助

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

求助，该怎么处理呢# 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'

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：fit Replication Algorithm singular options