nn=length(timeserial)
    n=nn*(nn-1)/2
    s=numeric(n)

   k=1
      for(i in 2:n)
       for(j in 1:(i-1))
        {
        s[k]=(timeserial[i]-timeserial[j])/(i-j)
         k=k+1
         }
    q=quantile(s,0.5,na.rm=T)   
不知道哪里出了问题，计算出来的s值，后面多了很多分量NA