R:库(ks) - bug还是不行?

时间:2016-05-21 16:09:23

标签: r crash

在我的程序中,我需要有一个函数,它可以从任意样本中进行连续密度估计(在实数上定义)。所以我选择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会话崩溃: Rplot

如果一切正常(-7,11) - 请检查其他(大)数字,这可能是系统特定的。

0 个答案:

没有答案