PLV=function(x,y,f)
{
    I=ifreq(x,f,phase=TRUE,plot=FALSE)
    J=ifreq(y,f,phase=TRUE,plot=FALSE)
    z=an((I$p[,2]-J$p[,2])*180/pi)
    return(abs(sum(z)/length(z)))
}