为什么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)