É reproduzível este código? Caso afirmativo, necessita de alguma bibliotéca?
Delcio,
Basta fazer assim:durner <- as.formula(theta ~ w0 * (thetar + ((thetas - thetar)/(1 + (alpha0 * psi)^n0)^(1 - 1/n0))) + (1 - w0) * (thetar + ((thetas -
Delcio,
Basta fazer assim:
durner <- as.formula(theta ~
w0 * (thetar + ((thetas - thetar)/(1 + (alpha0 *
psi)^n0)^(1 - 1/n0))) +
(1 - w0) * (thetar + ((thetas - thetar)/(1 +
(alpha1 * psi)^n1)^(1 - 1/n1
dur <- nlsLM(durner,
Vejam o seguinte script:
rm(list=ls(all=TRUE))
library(minpack.lm)
dados<-read.table(header = TRUE, dec=",",text=
"psi theta
0,01 0,5
0,5 0,48
1 0,42
3 0,4
10 0,39
30 0,35
100 0,3
300 0,28
1000 0,27
")
dados
theta=dados$theta
psi=dados$psi