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.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? -- ====================================================================== 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 ======================================================================
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?
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.
O colchetes é um operador reservado do R para acessar índices, ou seja, fazer seleção nos objetos do R. x[1], y[2:4], X[2,3:8], L[[5]], m0[["residuals"]]. Nas linguagens de programação em geral não se usa os operadores (, [ e { como aprendemos na matemática para separar/priorizar as operações. No R apenas os parenteses devem ser usados em operações matemáticas, colchetes são para índices e chaves para blocos de código ou dentro da expression(), dentro outras excessões. Por exemplo, cada um deles ainda tem um significado reservado quando o assunto é expressão regular. À disposição. Walmes.
participantes (4)
-
ASANTOS -
Davi Butturi-Gomes -
Fernando Souza -
walmes .