由买买提看人间百态

boards

本页内容为未名空间相应帖子的节选和存档,一周内的贴子最多显示50字,超过一周显示500字 访问原贴
Statistics版 - 问一个R问题
相关主题
在线急等:做maximum likelihood estimation,用optimization怎么都的不出正确的数值请问proc genmod, sas问题
any one did EM to find MLE of mixed effects model in Rres log pseudo-likelihood function
MCMC算法的Posterior Std. 一定是减少的么?[合集] what is the simplest way to generate 2 unif RVs of give
请问如何在R中写recursive function?parameter estimation
一般的统计仿真实验和monte carlo simulation是一回事吗?[合集] EM 算法与Kalman filter有什么区别?
救命啊,一个很简单的matlab问题estimate inverse gaussian in R
在R里怎么跳到指定的地方?论马尔可夫链和谈恋爱(原创)
请教matlab non-linear optimization的问题[合集] multivariate garch parameter estimation
相关话题的讨论汇总
话题: maxt话题: 参数话题: time话题: estpara话题: startt
进入Statistics版参与讨论
1 (共1页)
w*******9
发帖数: 1433
1
我要循环地generate random data然后估计参数,估计参数的时候会极其偶尔地遇到收
敛很慢的情,我想如果这次估计参数超过一分钟的话就跳过这次运算,进入下一次
iteration。有没有类似tryCatch的东西能处理这种情况?
for(i in 1:100000){
# generate data
# Estimate Parameter <-- 我想做的是 if(time)>1min, then quit and goto i=
i+1.
}
t******g
发帖数: 372
2
EstPara<-function(x, time=as.integer(Sys.time()), maxt=60){
# say the time consuming part is in a for loop
for (j in 1:2){
# do something with x
Sys.sleep(j)
if(as.integer(Sys.time())-time>maxt){
stop('runing too long!')
}
}
return(j)
}
for(i in 1:10){
startt<-as.integer(Sys.time())
res<-tryCatch(EstPara(i, time=startt, maxt=i), error=function(e){
return('not converged')
})
print(paste(i, res))
}
大概如此,只要你可以在你估计参数的函数里经常check下用时然后处理下error/
warning之类的

i=

【在 w*******9 的大作中提到】
: 我要循环地generate random data然后估计参数,估计参数的时候会极其偶尔地遇到收
: 敛很慢的情,我想如果这次估计参数超过一分钟的话就跳过这次运算,进入下一次
: iteration。有没有类似tryCatch的东西能处理这种情况?
: for(i in 1:100000){
: # generate data
: # Estimate Parameter <-- 我想做的是 if(time)>1min, then quit and goto i=
: i+1.
: }
:

w*******9
发帖数: 1433
3
TooYoung, I really appreciate it. It probably meets my needs. I will try it
later. Thanks again.
1 (共1页)
进入Statistics版参与讨论
相关主题
[合集] multivariate garch parameter estimation一般的统计仿真实验和monte carlo simulation是一回事吗?
Question on SAS estimation result救命啊,一个很简单的matlab问题
I have a question about logistic regression model在R里怎么跳到指定的地方?
R LOOP自动退出的问题请教matlab non-linear optimization的问题
在线急等:做maximum likelihood estimation,用optimization怎么都的不出正确的数值请问proc genmod, sas问题
any one did EM to find MLE of mixed effects model in Rres log pseudo-likelihood function
MCMC算法的Posterior Std. 一定是减少的么?[合集] what is the simplest way to generate 2 unif RVs of give
请问如何在R中写recursive function?parameter estimation
相关话题的讨论汇总
话题: maxt话题: 参数话题: time话题: estpara话题: startt