Найдите лучший показатель модели наблюдаемых данных R

Я хотел бы знать, как подогнать эту модель/уравнение

y~x^(-p)

к моим данным, чтобы найти лучшее значение показателя степени p

Мои данные:

y=c(1.1178329,1.0871448,1.0897010,1.0759255,1.0535190,0.8725332)
x=c(6,5,4,3,2,1)

Я пробовал следующим образом

mod <- nlsLM( sigmasG ~ x^(-p),
 start = c(p = 0.01) , 
trace = TRUE, lower=c(0.01) , upper=c(1))

а результаты не меняются...

It.    0, RSS =  0.0671647, Par. =       0.01
It.    1, RSS =  0.0671647, Par. =       0.01

Спасибо всем, кто может помочь с этой настройкой, чтобы найти p


r model-fitting
person david clarck    schedule 21.08.2017    source источник
comment
Вы пробовали использовать журнал? результат ‹- lm(log(y) ~ x); p ‹- - результирующие $коэффициенты[2]   -  person csgroen    schedule 22.08.2017
comment
@csgroen, спасибо за подсказку, в вашем ответе значение p определяется как p ‹- -result$coefficients[2] ? С минусовым синусом перед результирующими $коэффициентами[2] ?   -  person david clarck    schedule 22.08.2017
comment
Пожалуйста, убедитесь, что размещенные вопросы являются полными и автономными. Оператор библиотеки отсутствует на этом. См. минимальный воспроизводимый пример   -  person G. Grothendieck    schedule 22.08.2017
comment
@davidclarck Да. Но если подумать, на самом деле это должно быть lm(log(y) ~ log(x)). Но p по-прежнему будет -result$coefficients[2]. Это похоже на применение log к обеим сторонам, чтобы p мог стать линейным коэффициентом. Однако обратите внимание, что R-квадрат для этой модели составляет около ~ 0,78.   -  person csgroen    schedule 22.08.2017
comment
@csgroen, при создании отношения lm (log (y) ~ log (x)), мое p было: p (log (x)) = -0,1247224. Это имеет смысл? спасибо   -  person david clarck    schedule 22.08.2017
comment
@davidclarck, это не идеальная модель. Попробуйте выполнить x^-p и посмотрите, насколько далеки значения от y и приемлемо ли это.   -  person csgroen    schedule 22.08.2017

Ответы (1)


arrow_upward
3
arrow_downward

Если вы попробуете его без верхней и нижней границ, то он найдет отрицательное значение для p, поэтому неудивительно, что он приводит значение к наименьшему возможному значению в допустимой области.

Кроме того, в качестве отступления не рекомендуется иметь начальное значение на границе допустимой области, хотя в данном случае кажется, что это не имеет значения.

person G. Grothendieck    schedule 21.08.2017