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.