使用rms包进行OLS模型校准

时间:2016-07-28 16:47:20

标签: r calibration statistics-bootstrap

我正在尝试使用1000次复制(rms包)的引导来执行OLS模型的校准。有趣的是,我使用默认的40次复制没有问题,但是当我进入例如100,我得到以下文字:

  

lowess(x,y,iter = 0)中的错误:'delta'必须是有限的且> 0

如果我反复运行该命令(大约10次),我可以解决这个问题,但是如果我想去,例如500次重复,我必须重新运行代码大约100次。不幸的是,我必须做1000次重复,无论我重新运行多少次,我都不会得到结果。有什么想法吗?

这是脚本:

fit <- ols( stiffness ~ cap + cd4 + cd8 + age + BMI + cd48 + diabetes + 
                         hypertonia + triglicerid + cholesterin + 
                         lipodystrophia + betegsegfennalla + darunavir + 
                         atanazavir + raltegravir + etravirine + 
                         nevirapine + efavirenz + tenofovir + abacavir + 
                         zidovudine + lamivudine + lopinavir,
             data = S[!is.na( S$cholesterin ) & !is.na( S$betegsegfennalla ) & 
                      !is.na(S$triglicerid ) & !is.na( S$hypertonia ) & 
                      !is.na( S$BMI ), ], 
             x = TRUE, y = TRUE ) 

plot( calibrate( fit, B = 1000 ) )
  

lowess(x,y,iter = 0)中的错误:'delta'必须是有限的且&gt; 0

0 个答案:

没有答案