imagine o tanto de aproximacoes nao aconteceram quando vc chamou integrate()... experimente all.equal(integrate(f=fdp_expweibull,par=c(1,1,1),0,Inf)$value, 1) b Em 5 de julho de 2013 18:23, Pedro Rafael <pedro.rafael.marinho@gmail.com> escreveu:
Acho que percebi. Convertendo o resultado para carácter o valor não é 1.
as.character(integrate(f=fdp_expweibull,par=c(1,1,1),0,Inf)$value)
[ ], Pedro Rafael Diniz Marinho.
Em 5 de julho de 2013 18:00, Pedro Rafael <pedro.rafael.marinho@gmail.com> escreveu:
Pessoal eu tenho uma dúvida. Considere o código abaixo:
function(par,x){ beta = par[1] c = par[2] a = par[3] a * beta * c * exp(-(beta*x)^c) * (beta*x)^(c-1) * (1 - exp(-(beta*x)^c))^(a-1) }
integrate(f=fdp_expweibull,par=c(1,1,1),0,Inf)$value == 1
Porque o resultado é false?
[ ], Pedro Rafael Diniz Marinho.
_______________________________________________ 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.