在非参数回归中存储和检索带宽对象

Question

我使用R s np包（Racine＆amp; Hayfield）非参数估计条件密度。

估计带宽是非常耗时的（我的实际例子约20分钟，而不是下面的可再现的例子）。因此，在说明类中的最终结果之前，我想预先计算带宽对象，存储它并在课堂上检索它。

我想到这样的事情，没有成功：

library(np)

# some artificial data
n <- 100
y <- rnorm(n)
X1 <- rnorm(n)
X2 <- rnorm(n)

# nonparametric regression
bw <- npcdensbw(y ~ X1+X2, bwmethod="cv.ls")
fhat <- npcdens(bws=bw)

save(bw,file="bandwidth.Rdata") # saves something, but...

# ...this does not work
rm(list=ls())
bw <- load("bandwidth.Rdata")
fhat <- npcdens(bws=bw)