Warning: file_get_contents(/data/phpspider/zhask/data//catemap/4/r/64.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181

Warning: file_get_contents(/data/phpspider/zhask/data//catemap/6/jenkins/5.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181
为什么R的物流增长产出与另一个存在差异?_R_Statistics_Nls - Fatal编程技术网

为什么R的物流增长产出与另一个存在差异?

为什么R的物流增长产出与另一个存在差异?,r,statistics,nls,R,Statistics,Nls,但是给了我们: 我知道有些语言有不同的输出。这很正常,但我想知道你的想法是什么 提前谢谢。问题是不同的型号适合 nls具有自启动功能的SSlogis适合该型号(请参见help('SSlogis')) 或者,用你的符号 Asym/(1+exp((xmid-input)/scal)) 笔和纸显示以下转换在网页中给出结果 phi1/(1 + exp((phi2 - input)/phi3)) fit谢谢。尽管公式不同,为什么我们要使用这个呢?所以,正如我所知,公式是y=c/(1+a*(e^-b*

但是给了我们:

我知道有些语言有不同的输出。这很正常,但我想知道你的想法是什么


提前谢谢。

问题是不同的型号适合

nls
具有自启动功能的
SSlogis
适合该型号(请参见
help('SSlogis')

或者,用你的符号

Asym/(1+exp((xmid-input)/scal))
笔和纸显示以下转换在网页中给出结果

phi1/(1 + exp((phi2 - input)/phi3))

fit谢谢。尽管公式不同,为什么我们要使用这个呢?所以,正如我所知,公式是y=c/(1+a*(e^-b*x))。@nerdakgul就像我在文章中说的,
SSlogis
适合不同的模型,它适合公式
Asym/(1+exp((xmid-x)/scal))
。这等于Asym/(1+exp(xmid/scal)*exp(-1/scal*x))
。然后,
c=Asym
a=exp(xmid/scal)
b=1/scal
Asym=coef[1]
xmid=coef[2]
scal=coef[3]
。在函数
transf
中有转换。
Asym/(1+exp((xmid-input)/scal))
phi1/(1 + exp((phi2 - input)/phi3))
fit <- nls(var ~ SSlogis(time, phi1, phi2, phi3), data = df)

kappa <- coef(fit)[1]
alpha <- exp(coef(fit)[2]/coef(fit)[3])
beta <- 1/coef(fit)[3]
c(kappa = unname(kappa), alpha = unname(alpha), beta = unname(beta))
#      kappa       alpha        beta 
#105.7373679   6.8832991   0.2595471 
transf <- function(x){
  kappa <- coef(x)[1]
  alpha <- exp(coef(x)[2]/coef(x)[3])
  beta <- 1/coef(x)[3]
  c(kappa = unname(kappa), alpha = unname(alpha), beta = unname(beta))
}

transf(fit)