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@usp.br -----Mensagem original----- De: R-br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de ASANTOS Enviada em: quarta-feira, 2 de julho de 2014 16:37 Para: r-br@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@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== _______________________________________________ R-br mailing list R-br@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.