R: en Utilisant l'équation avec le logarithme naturel de la nls

Bonne journée,

J'ai du mal avec le R et le logarithme naturel (ln). Tout d'abord, je ne peux pas trouver un ln(x) en fonction de R. j'ai remarqué que log(x) est le même que ln(x) (lors de l'utilisation de ln(x) avec une calculatrice).

Dans R:

log(5) = 1.609438

Et avec une calculatrice:

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

Je suis en essayant d'intégrer une équation dans R (c'est exactement comment je l'ai trouvé dans la littérature de 3 références):

y = a + b(x/305) + c(x/305)2 + d ln(305/x) + f ln2(305/x)

Est-il correct d'utiliser la syntaxe suivante dans la R utilisation de l'équation?

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

L'idée est d'utiliser cette fonction avec nls() dans R.
Merci à l'avance!

  • log dans R désigne le logarithme népérien. C'est la convention de mathématiciens, depuis le "commun" des logarithmes ont aucun intérêt mathématique. Le "ln" abréviation est quelque chose qui a été introduit pour rendre les choses moins à confusion pour les étudiants.
  • Vous n'avez pas besoin d'une calculatrice pour savoir ce log n'. Vous avez seulement besoin de lire la documentation: help("log"). Votre formule semble correct.
InformationsquelleAutor wernerfeuer | 2014-06-19
  1. 17

    Dans R, log est le logarithme naturel. Dans les calculatrices, les journaux signifie généralement logarithme de base 10. Pour l'accomplir dans R, vous pouvez utiliser le log10 fonction.

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

    Comme pour votre formule, il semble correct, puisque log est le logarithme naturel.

    • Merci pour la réponse! J'en attendais beaucoup et a juste besoin de conformation. Merci encore!
    InformationsquelleAutor alko989
  2. 4

    En outre, je ferai remarquer que votre modèle

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

    est linéaire dans le sens statistique du terme de linéaire à coefficients; il n'a pas besoin d'être linéaire en x.

    Vous n'avez pas besoin nls pour s'adapter à ce modèle, vous pouvez utiliser lm().

    Mais n'oubliez pas de regarder le I() la fonction d'exprimer des termes tels que (x/305)^2.

    ETA exemple:

    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)
    InformationsquelleAutor user20637