Я предполагаю, что вы имеете в виду, что у вас есть векторы y
и x
, и вы пытаетесь выполнить функцию y(x)=Alog(x)
.
Прежде всего, сборка журнала - плохая идея, потому что она не ведет себя хорошо. К счастью, у нас есть x(y)=exp(y/A)
, поэтому мы можем подобрать экспоненциальную функцию, которая намного удобнее. Мы можем сделать это с помощью нелинейного метода наименьших квадратов:
nls(x~exp(y/A),start=list(A=1.),algorithm="port")
где start
является начальное предположение для A
. Этот подход представляет собой численную оптимизацию, поэтому он может потерпеть неудачу.
Чем более стабильный способ, чтобы превратить его в линейной функции, log(x(y))=y/A
и установить прямую линию с помощью lm
:
lm(log(x)~y)
Больше информации. Что значит «точки данных, представляющие логарифмическую функцию»? –
Вы ищете логистическую регрессию? http://en.wikipedia.org/wiki/Logistic_regression – Amro