R 正在尝试匹配我的函数
我尝试将这个函数循环1000次,并记录模型的系数 但有时当迭代超过50次时,函数停止。如何使函数继续运行并记录时间失败?谢谢 我正在使用tryCatch,但我不知道如何使它工作。 我如何计算1000次中有多少次错误? 请帮忙R 正在尝试匹配我的函数,r,function,for-loop,data-fitting,R,Function,For Loop,Data Fitting,我尝试将这个函数循环1000次,并记录模型的系数 但有时当迭代超过50次时,函数停止。如何使函数继续运行并记录时间失败?谢谢 我正在使用tryCatch,但我不知道如何使它工作。 我如何计算1000次中有多少次错误? 请帮忙 e0<-1:1000 emax<-1:1000 e50<-1:1000 for (i in 1:1000) { rho=0.4 #5g pre.5mg<-rnorm(n=50,mean=-0.5,sd=0.2) x=0.4*rho-0.5
e0<-1:1000
emax<-1:1000
e50<-1:1000
for (i in 1:1000)
{
rho=0.4
#5g
pre.5mg<-rnorm(n=50,mean=-0.5,sd=0.2)
x=0.4*rho-0.5
post.5mg<-rnorm(n=50,mean=pre.5mg*rho+x,sd=sqrt((0.3*0.3*(1-rho*rho))))
#25g
pre.25mg<-rnorm(n=50,mean=-0.5,sd=0.2)
x25=0.4*rho-1
post.25mg<-rnorm(n=50,mean=pre.25mg*rho+x25,sd=sqrt((0.3*0.3*(1-rho*rho))))
#100g
pre.100mg<-rnorm(n=50,mean=-0.5,sd=0.2)
x100<-0.4*rho-1.2
post.100mg<-rnorm(n=50,mean=pre.100mg*rho+x100,sd=sqrt((0.3*0.3*(1-rho*rho))))
#250g
pre.250mg<-rnorm(n=50,mean=-0.5,sd=0.2)
x250<-0.4*rho-1.4
post.250mg<-rnorm(n=50,mean=pre.250mg*rho+x250,sd=sqrt((0.3*0.3*(1-rho*rho))))
#response
Y<-c(post.5mg,post.25mg,post.100mg,post.250mg)
#dose
dose <- c(rep(10,50),rep(30,50),rep(120,50),rep(240,50))
w<-data.frame(Y,dose)
#MLE hyperbolic EMAX
abc=nls(formula=Y~e0+((dose*emax)/(dose+ed50)),start=c(e0=1000,emax=10000,ed50=20)
,data=w,algorithm ="default")
bcd<-summary(abc)
e0[i]<-bcd$coef[1]
emax[i]<-bcd$coef[2]
e50[i]<-bcd$coef[3]
tryCatch(1, finally = print("failed"))
}
附上您的代码如下:
try( {
# put your code here
# ...
},silent=FALSE)
上面写着“错误:意外”,“输入:好的,我知道了,我可以运行1000次,但是我如何才能从1000次失败中获得数字”在这个链接上读到,其他的则从谷歌上读到: