Desculpem a falha, consegui rodar, não tinha visto a mensagem anterior!
É 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 -
É 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 -
Cleber,
Sua pergunta extrapola o objetivo deste grupo (não CMR, etc.).
Ademais ela enseja muitas questões antes que possa ser iniciada uma
discussão com assertividade.
Minha recomendação é que você a refaça nos foruns mais voltados para
questões de Estatística (p.ex.:
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
Bom dia listeiros (@)
Estou realizando um experimento e de uma única amostra inicial dividi em
várias subamostras e em cada uma apliquei um tratamento e ao longo dos dias
fiz a análise dos parâmetros. Gostaria de verificar o efeito dos
tratamentos sob as alterações dos parâmetros em relação ao
