我具有以下分布函数(F):
F<-function(t) {
if(any(t>=0)) (1+ (-1)/(1+t^2)) else (0)
}
绘图时,其外观如下:
plot(Vectorize(F), xlim=c(-100,100), main="Distribution function", xlab="t", ylab="F(t)")
我正在尝试为此F分布计算0.25、0.75和0.50分位数,但出现以下错误消息:
quantile(F, probs = seq(0.25, 0.5,0.75))
分位数中的错误。default(F,概率= seq(0.25,0.5,0.75)): anyNA()应用于非“闭包”类型的(列表或向量)
此外,以下代码也不起作用:
quantile(F, probs = seq(0.25, 0.50, 0.75), na.rm = TRUE)
x [!is.na(x)]中的错误:'closure'类型的对象不可子集化 另外:警告消息: 在is.na(x)中:is.na()应用于类型为“ closure”的非(列表或向量)
我无法解决问题。 谁能告诉我该如何解决?
谢谢