我在 R 中使用 BayesTree 包。使用作者的示例:
##simulate data (example from Friedman MARS paper)
f = function(x){
10*sin(pi*x[,1]*x[,2]) + 20*(x[,3]-.5)^2+10*x[,4]+5*x[,5]
}
sigma = 1.0 #y = f(x) + sigma*z , z~N(0,1)
n = 100 #number of observations
set.seed(99)
x=matrix(runif(n*10),n,10) #10 variables, only first 5 matter
Ey = f(x)
y=Ey+sigma*rnorm(n)
##run BART
set.seed(99)
bartFit = bart(x,y)
plot(bartFit) # plot bart fit
执行 plot(bartFit) 时,它会生成两个图“绘图方法将 mfrow 设置为 c(1,2) 并制作两个图。” 我希望只使用此命令生成的第二个图。有没有一种简单的方法可以只提取第二个情节?
这也提供:
## S3 method for class 'bart':
plot((
x,
plquants=c(.05,.95), cols =c('blue','black'),
...))
我对 R 比较陌生,一般来说,如果一个命令给出的情节超过 1 个,你怎么能选择你想要的特定情节?例如,在线性模型 (lm) 绘图命令中,您还可以获得多个绘图。