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]]
Maybe Matching Threads
- Cross-checking a custom function for separability indices
- Curso de Introducción a R (6ª Edición) - Barcelona, 19-22 Septiembre 2011
- ROC curve with survival data
- How to disable output messages (prints or cats) from functions in R?
- get data from Perl arrays
Wisdom of the Ancients
