[R-br] Problema para ajustar modelo nls

[Davi Butturi-Gomes]

Não testei, mas o problema não está nesse colchete
log[p.H.termes(1-p.H.termes)] e na ausência da multiplicação?
mod<-nls(log(var.H.termes)~log(a)+b*log[p.H.termes(1-p.H.termes)],start=list
(a=1,b=1))

Talvez:
mod<-nls(log(var.H.termes)~log(a)+b*log(p.H.termes*(1-p.H.termes)),start=lis
t(a=1,b=1))

____________________
MSc. Davi Butturi-Gomes
ESALQ-USP
davibg em usp.br

-----Mensagem original-----
De: R-br [mailto:r-br-bounces em listas.c3sl.ufpr.br] Em nome de ASANTOS
Enviada em: quarta-feira, 2 de julho de 2014 16:37
Para: r-br em listas.c3sl.ufpr.br
Assunto: [R-br] Problema para ajustar modelo nls

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.111
22031,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?

-- 
======================================================================
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial
Cáceres - MT                      CEP: 78.200-000
Fone: (+55) 65 8132-8112 (TIM)   (+55) 65 9686-6970 (VIVO)
e-mails:alexandresantosbr em yahoo.com.br
         alexandre.santos em cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================

_______________________________________________
R-br mailing list
R-br em listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.