R 求方程的根
我想解下面给出的方程式wrt c 其中n=k=30,x如下所示 我使用了R 求方程的根,r,equation,R,Equation,我想解下面给出的方程式wrt c 其中n=k=30,x如下所示 我使用了rootSolve包中的uniroot命令和下面的代码,但它给了我一个错误 model <- function(c){ sx <- apply(x,1, function(x) exp(c*x)) numerator <- sum(sx[2:30])-29*sx[1] denumerator <- (30*29*(sx[1]-1)) numerator/denumerator-0.2
rootSolve
包中的uniroot
命令和下面的代码,但它给了我一个错误
model <- function(c){
sx <- apply(x,1, function(x) exp(c*x))
numerator <- sum(sx[2:30])-29*sx[1]
denumerator <- (30*29*(sx[1]-1))
numerator/denumerator-0.2542
}
ss <- uniroot(model, c(0,1))
model未定义c=0的分母。尝试uniroot(型号c(1e-10,1))