Update R/analysis.R

update sigtest() routine