У меня есть следующая функция распределения (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)")
Я пытаюсь вычислить квантили 0,25,0,75 и 0,50 для этого F-распределения, но получаю следующее сообщение об ошибке:
quantile(F, probs = seq(0.25, 0.5,0.75))
Ошибка в quantile.default (F, probs= seq (0.25, 0.5, 0.75)): anyNA (), примененная к non (списку или вектору) типа 'closure'
Кроме того, следующий код также не работает:
quantile(F, probs = seq(0.25, 0.50, 0.75), na.rm = TRUE)
Ошибка в x [! Is.na (x)]: объект типа «замыкание» не может быть поднабором. Дополнительно: Предупреждение: In is.na (x): is.na () применяется к не- (список или вектор) типа 'замыкание'
Я не могу решить проблему.Может кто-нибудь сказать мне, как это исправить?
Спасибо