[R-br] Problema para ajustar modelo nls

Fernando Souza nandodesouza em gmail.com
Quarta Julho 2 16:53:22 BRT 2014


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?
>


Mais detalhes sobre a lista de discussão R-br