r - Rコードで関数を最適化するにはどうすればよいですか

Question

最適化問題：

    x <- c(1,2,3,4,5,8,4,5)
    y <- c(-2,-5,-3,1,5,3,-17,1)
    alpha0 <- 0.01
    Q <- function(alpha, x, y){

        n <- length(x)
        mu <- sum(x)
        xi <- sum(y)
        L1 <- (alpha-1)*xi-alpha*mu
              -n*log(gamma(alpha))+n*alpha*log(alpha)
        return(-L1)
    }

      Qaf <- nlm(Q, alpha0, x, y, hessian=T)
      Qaf <- nlminb(alpha0, Q, x=x, y=y, hessian=T)
      Qaf <- optim(alpha0, Q, x=x, y=y, hessian=T, method="CG")
      Qaf <- optim(alpha0, Q, x=x, y=y, hessian=T, method="BFGS")
      Qaf <- optim(alpha0, Q, x=x, y=y, hessian=T, method="SANN")
      Qaf <- optim(alpha0, Q, x=x, y=y, hessian=T)

すべての最適化により、NaNのパラメーターの結果とwarnings（）の数が異なります。この問題を解決して、パラメーター空間を再パラメーター化できますか？パラメータ空間を再パラメータ化し、Rコードを使用してこれを解決するにはどうすればよいですか？