R: Используя уравнение с натуральным логарифмом в НСГ

Question

Добрый день,

Я борюсь с R и натуральный логарифм (LN). Во-первых, я не могу найти функцию ln (x) в R. Я заметил, что log (x) совпадает с ln (x) (при использовании ln (x) с калькулятором).

В R:

log(5) = 1.609438 

И с калькулятором:

ln(5) = 1.609438 
log(5) = 0.69897 

Я пытаюсь подогнать уравнение R (это точно, как я нашел в литературе 3 ссылки) :

у = а + Ь (^х /) + с (^х /) + D пер (/ _х) + F пер (/ _х)

Правильно ли использовать следующее синтаксис в R для использования уравнения?

y ~ a + b*(x/305) + c*((x/305)^2) + d*log(305/x) + f*(log(305/x))^2 

Идея заключается в том, чтобы использовать эту функцию с NLS() в R. Заранее спасибо!

Answer 1

In R, log - натуральный логарифм. В калькуляторах журнал обычно означает логарифм базы 10. Для этого в R вы можете использовать функцию log10.

log(5) 
## [1] 1.609438 
log10 
## [1] 0.69897(5) 

Что касается вашей формулы, представляется правильным, так как log натуральный логарифм.

Answer 2

Кроме того, я укажу, что ваша модель

y ~ a + b*(x/305) + c*((x/305)^2) + d*log(305/x) + f*(log(305/x))^2 

линейна в статистическом смысле является линейным в коэффициентах; он не должен быть линейным по x.

Вам не нужны nls для этой модели, вы можете использовать lm().

Но не забудьте взглянуть на функцию I(), чтобы выразить такие выражения, как (x/305)^2.

ЕТА пример:

aDF <- data.frame(x=abs(rnorm(100)), y=rnorm(100)) 
lm(y ~ 1 + I(x/305) + I((x/305)^2) + log(305/x) + I(log(305/x)^2), data=aDF)