Tente mod<-nls(log(var.H.termes)~log(a)+b*log(p.H.termes*(1-p.H.termes)),start=list(a=1,b=1)) Você errou na forma de utilizar a função log() log[p.H.termes(1-p.H.termes)] : ERRADO log(p.H.termes*(1-p.H.termes)): CORRETO (parenteses e não colchetes att Em Qua 02 Jul 2014 16:37:26 BRT, ASANTOS escreveu:
Boa tarde Pessoal,
Estou tentando ajustar os coeficientes de um modelo não linear e uma das variáveis esta sendo confundida com uma função, mudei o nome e não resolveu, segue CRM:
var.H.termes<-c(0.07692308,0.20923521,0.24819625,0.03059211,0.03995434,0.11122031,0.25221018,
0.1343785,0.09014085,0.07893114,0.23298351,0.10613884,0.01098901,0.10188088,0.25029985)## Vetor da variância
p.H.termes<-c(0.08333333,0.29292929,0.56565657,0.03125,0.04109589,0.12612613,
0.50909091,0.15789474,0.09859155,0.08571429,0.36206897,0.11904762,0.01098901,0.11363636,0.45689655)## Vetor de probabilidade
## Lei de potência de Taylor binomial - log(v)=log(a)+b*log[p(1-p)]
mod<-nls(log(var.H.termes)~log(a)+b*log[p.H.termes(1-p.H.termes)],start=list(a=1,b=1))
Alguém poderia me ajudar?