We examine a sample of detected and undetected Isolated Neutron Stars (INSs), selected on the basis of their energy loss and distance, in order to maximize their detection probability. Since the sample we consider encompasses more upper limits than detections, we make use, for the first time in this field of astronomy, of survival analysis procedures through the ASURV software package (Feigelson \& Nelson 1985, Isobe et al. 1986). We show that these techniques lead to a better understanding of the physical processes at work in high-energy emission from INSs. In particular, the recent detection of PSR 1951+32 and upper limits from pulsars pointed but not detected by EGRET show that the $\gamma $-ray efficiency of ISN is not correlated to any simple pulsar parameter.

After reviewing the multifrequency behaviour of the Isolated Neutron Stars detected so far, we analyze for each object the efficiency of conversion of the star's rotational energy loss into optical, X and gamma radiation. Although the number of pulsars detected at different wavelengths is rather limited, a pattern is seen to emerge from our analysis pointing towards the period derivative as the leading parameter to describe the multifrequency emission of INS. One object in particular, PSR 1509-58, stands out as archetypal for the Pdot-dependence of its luminosity at different wavelengths.