Ваш псевдоним должен быть в вашем файле .profile не в вашем скрипте, если вы вызываете его в приглашении.
Если вы поместите псевдоним в свой скрипт, тогда вы должны называть его в своем скрипте .
Источник файла является правильным ответом при попытке запустить скрипт, внутри которого есть псевдоним.
source yourscript.sh
x и y меняются местами, и e ^ (...) должно быть exp (...). Также я обнаружил, что установка d в 0 помогла.
d <- 0 # fix d at 0
st <- list(a = mean(y), b = 1/sd(x), c = mean(x))
fm <- nls(y ~ a/(1+exp(-b*(x-c))) + d, start = st)
fm
дает:
Nonlinear regression model
model: y ~ a/(1 + exp(-b * (x - c)))
data: parent.frame()
a b c
19.96517 0.02623 99.73842
residual sum-of-squares: 1.82
Number of iterations to convergence: 9
Achieved convergence tolerance: 9.023e-06
При составлении графика это, кажется, подходит визуально:
plot(y ~ x)
lines(fitted(fm) ~ x, col = "red")