如何利用极大似然函数提取r中的极大似然值

我的概率分布函数是

f(x;r)=(1+r)/Γ(1/(1+r))exp(-x^(1+r)) , 0<x<Inf, r>-1

但我无法获得该值，因为它会给出以下错误

fit1$系数中出错：$运算符未为此S4类定义

现在我得到了答案

nLL <- function(r)-sum(log(1+r)-log(gamma(1/(1+r)))-x1^(1+r))
fit0 <- mle(nLL, start = list(r = 0.1), nobs = NROW(x1))
fit1 <- mle(nLL, start = list(r = 0.1), nobs = NROW(x1),
        method = "Brent", lower = -1, upper = 200)
stopifnot(nobs(fit0) == length(x1))
class(fit1)
str(fit1)
fiteg[k] = fit1@coef

---

nLL您到底尝试了什么，哪些不起作用？你确定那是PDF文件吗？例如，对于
r=.5
，它似乎没有集成到1。你是想用数值方法还是解析方法来解决这个问题？这是一个pdf，因为它是伽马分布的一个特例，c=1，比例参数，b=1/（1+r），形状参数。我很抱歉打字。我正在学习安装MathJax。它真的被称为“非正规”发行版吗？我很抱歉。我只是证明它是伽马分布的一个特例。我将改变我的热带。我为我的错误道歉。@Sania，谢谢你的帮助-。
fitv = fit1$Coefficients 

nLL <- function(r)-sum(log(1+r)-log(gamma(1/(1+r)))-x1^(1+r))
fit0 <- mle(nLL, start = list(r = 0.1), nobs = NROW(x1))
fit1 <- mle(nLL, start = list(r = 0.1), nobs = NROW(x1),
        method = "Brent", lower = -1, upper = 200)
stopifnot(nobs(fit0) == length(x1))
class(fit1)
str(fit1)
fiteg[k] = fit1@coef