I've calculated as follows:
matusita<-function(f1,f2,lw,up){
fmatu<-function(x) ( sqrt(f1(x)) - sqrt(f2(x)) )^2
sqrt (integrate( fmatu, lower=lw, upper=up, subdivisions=1000000)$value)
}
Josué Almansa
Unitat de Recerca en Serveis Sanitaris
Institut Municipal d'Investigació Mèdica, IMIM
c/ Doctor Aiguader, 80
08003 Barcelona
e-mail: jalmansa@imim.es
[[alternative HTML version deleted]]
