gbm包的示例代码在calibrate部分报错如何解决？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

1.     if（！is.null（offset））{
   
2.       适合<-适合+ offset [model.mask]
   
3.     }
   
4.     if（family ==“ bernoulli”）{
   
5.       适合< -exp （适合）/（1 + exp（适合））
   
6.     }}
   
7.     否则if（family ==“ poisson”）{
   
8.       fits <     -exp （fits）
   
9.     }
   
10. fit.matrix [model.mask，i] <-fits
    
11.     fits < -gbm :: predict.gbm（model.list [[i]] ，x.data [pred.mask，，
    
12.       drop = FALSE]，n.trees = target.trees）
    
13.     如果（！is.null（offset））
    
14.       适合<-fits + offset [pred.mask]
    
15.     fold.fit [pred。掩码] <-适合
    
16.     if（family ==“ bernoulli”）{
    
17.       适合< -exp （fits // 1 + exp（fits））
    
18.     }
    
19.     else if（family ==“泊松”）{
    
20.       适合<-exp（适合）
    
21.     }
    
22.     fit.matrix [pred.mask，i] <-适合
    
23.     y_i <-y.data [pred.mask]
    
24.     u_i <-fit.matrix [pred.mask，i]
    
25.     weight.preds <-site.weights [pred.mask] ]
    
26.     cv.deviance.stats [i] <-calc.deviance（y_i，u_i，weight.preds，
    
27.       family = family）
    
28.     cv.cor.stats [i] <-cor（y_i，u_i）
    
29.     if（family ==“ bernoulli “）{
    
30.       cv.roc.stats [i] <-gbm :: gbm.roc.area（y_i，u_i）
    
31.       cv.calibration.stats [i，] <-gbm :: calibrate.plot（y_i，u_i，
    
32.         ”二项式“）
    
33.       threshold.stats [i] <-大约（ppoints（u_i），sort（u_i，
    
34.         减少= T），患病率）$ y
    
35.     }
    
36.     if（family ==” poisson“）{
    
37.       cv.calibration.stats [i，] <-gbm :: calibrate.plot（y_i，u_i，
    
38.         “泊松”）
    
39.     }
    
40.   }

复制代码

cv.calibration.stats [i，] <-gbm :: calibrate.plot（y_i，u_i，“ binomial”）中的错误：
  被替换的项目不是替换值长度的倍数

报错命令在倒数行左右的位置上，一开始源代码中命令是.calibrate和.roc，修改后更改gbm中相关的命令，还是报错。
有没有用过gbm包的大佬指点一下。

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：calibrate Rate bra ATE ALI 回归树 R语言 回归树 R语言

已解决，library(dismo)即可。