在我的程序中,我需要有一个函数,它可以从任意样本中进行连续密度估计(在实数上定义)。所以我选择library(ks)
并发现它有时会产生错误的对象,而其他函数(例如plot
)在访问时会崩溃R会话。
请,我希望其他人检查是否是包中的错误,我的R版本,(甚至我做错了)。
所以,重现崩溃的代码:
library(ks)
set.seed(8192)
density_generator<-function(s)
{
# this function returns kernel density estimate built on sample 's'
hpi1 <- hpi(x=s) # calculating h parameter for kernel estimation
fhat.pi1 <- kde(x=s, H=hpi1) # generating density object
fhat.pi1
}
## testing the density generator
conditional_density_object<-density_generator(c(1,2,3,4,5))
foo<-function(z){predict(conditional_density_object,x=z)}
y<-seq(from=-7,to=11, by=0.01) # R session fails for some of the parameters
plot(y,sapply(y,foo),pch=".")
其他参数,例如y<-seq(from=-7,to=5, by=0.01)
不会导致R会话崩溃:
如果一切正常(-7,11) - 请检查其他(大)数字,这可能是系统特定的。