2016-11-08 14 views
3

У меня возникли проблемы с изменением меток x и y на частичном графике для модели gbm. Мне нужно переименовать их в статью журнала.Изменение меток x и y на частичном графике gbm

Я прочитал это и создать сюжет следующим образом:

library(gbm) 
final<- readRDS(final_gbm_model) 
summary(final, n.trees=final$n.trees) 

Вот резюме выход:

             var rel.inf 
ProbMn50ppb          ProbMn50ppb 11.042750 
ProbDOpt5ppm          ProbDOpt5ppm 7.585275 
Ngw_1975            Ngw_1975 6.314080 
PrecipMinusETin_1971_2000_GWRP PrecipMinusETin_1971_2000_GWRP 4.988598 
N_total            N_total 4.776950 
DTW60YrJurgens         DTW60YrJurgens 4.415016 
CVHM_TextZone         CVHM_TextZone 4.225048 
RiverDist_NEAR         RiverDist_NEAR 4.165035 
LateralPosition        LateralPosition 4.036406 
CAML1990_natural_water     CAML1990_natural_water 3.720303 
PctCoarseMFUpActLayer     PctCoarseMFUpActLayer 3.668184 
BioClim_BIO12         BioClim_BIO12 3.561071 
MFDTWSpr2000Faunt       MFDTWSpr2000Faunt 3.383900 
PBot_krig           PBot_krig 3.362289 
WaterUse2           WaterUse2 3.291040 
AVG_CLAY            AVG_CLAY 3.280454 
Age_yrs            Age_yrs 3.144734 
MFVelSept2000         MFVelSept2000 3.064030 
AVG_SILT            AVG_SILT 2.882709 
ScreenLength          ScreenLength 2.683542 
HydGrp_C            HydGrp_C 2.666106 
AVG_POR            AVG_POR 2.563147 
MFVelFeb2000          MFVelFeb2000 2.505106 
HiWatTabDepMin         HiWatTabDepMin 2.421521 
RechargeAnnualmmWolock     RechargeAnnualmmWolock 2.252706 

Я могу создать частичную зависимость участок следующим образом:

plot(final,"ProbMn50ppb",n.trees=final$n.trees) 

enter image description here Но если я попытаюсь установить ярлык argume nts Я получаю следующую ошибку:

plot(final,"ProbMn50ppb",n.trees=final$n.trees,ylab="LNNO3") 

Error in plot.default(X$X1, X$y, type = "l", xlab = x$var.names[i.var], : 
    formal argument "ylab" matched by multiple actual arguments 

Как изменить метки оси y и x?

+0

Спасибо за попытку предоставить ваши данные, но мы не позволяем личным облачным ссылкам из-за вредоносного ПО и потому, что они ломаются со временем. Также, как наилучшая практика и избежать повреждения файлов, я бы рекомендовал вам предоставить вашим RDS-файлам расширение «.RDS». –

+0

Хорошо, спасибо! Есть ли лучший способ предоставить такие данные? – user29609

+0

Как правило, лучше всего воссоздать то, что вы можете в своем коде, и 'dput()' что-нибудь еще. Таким образом, он легко воспроизводится и содержит все данные класса/мета. –

ответ

4

Функция plot.gbm передает свое собственное имя в общую функцию графика, так что они сталкиваются друг с другом. Таким образом, вы не сможете настроить сюжет так, как вы хотите в этом режиме. Но авторы предоставили альтернативу, если вы установили return.grid=TRUE. Вместо того, чтобы строить сюжет, он выдаст сами данные. Затем вы можете использовать это для любого сюжета, включая ggplot2.

plotdata <- plot(gbm1, return.grid=TRUE) 
plot(plotdata, type="l", ylab="ylab", xlab="xlab") 

enter image description here

Пример данные помощи (GBM)

+0

Это намного лучше, чем решение, о котором я думал, что должно было бы изменить исходный код (потому что 'plot.gbm' включает в себя как лаборатории по умолчанию, так и' ... ', из-за чего возникает проблема). –

0

Вы также можете изменить сам объект GBM перед построением (или функции):

your_gbm_obj$var.names[index] = "axis label"