modified: dodseq.r