您可以使用以下代碼 library(minpack.lm)k_0 = 0.21k_1 = 0.21m = 52alpha = 0.05beta = 0.2pi_0 = 0.669power <- 1-betacz <- 20z_alpha <- qnorm(p= alpha/2, lower.tail=FALSE)Z_beta <- qnorm(p= beta, lower.tail=FALSE)fun <- as.formula(cz ~ 1 + ((z_alpha + Z_beta)^2)*((pi_0*(1-pi_0)/m + pi_1*(1 - pi_1)/m + (((k_0)^2)*((pi_0)^2) + ((k_1)^2)*((pi_1)^2)))/((pi_0 - pi_1)^2)))df <- data.frame(cz, z_alpha, Z_beta, alpha, m, beta, pi_0, k_0, k_1)#Fitting model using minpack.lm packagenls.out <- nlsLM(fun, data = df, start=list(pi_1=1), algorithm = "LM", control = nls.lm.control(maxiter = 500))summ