I've been using the sm.denstiy.compare function and would like to know some specifics: a) What kind of bootstrap sampling is used, with or without replacement? b) A p-value is returned after running the following script (for example): library(sm) y <- rnorm(100) g <- rep(1:2, rep(50,2)) sm.density.compare(y, g, model="equal") to what statistic does this p-value correspond? c) Also, in the help file this description appears for "Value" of sm.density.compareomparison of univariate density estimates When model is set to "none", nothing is returned. When "model" is set to "equal", a list containing the smoothing parameter and the p-value of the test is returned. When band takes the value TRUE, and there are only two groups to compare, the list contains in addition the upper and lower end-points of the reference band for equality. however, the mentioned "list" is not returned when I run the analysis... only the p-value is returned. Am I doing something wrong? Any help is appreciated. Thank you, Claudia [[alternative HTML version deleted]]