For COX proportional model ONLY. This function is to derive the deviance distribution based on the permutation method;This function is not to be used independently but will be called by the function selectnew()
update_dev_cox(data, vslist, C, ncores, weights)the variable 'data' within seqcutoff()
a vector of selected variables
the number of permutation times
=TRUE: if weighted version is desired, =FALSE, otherwise (binomial model); weights: =vector of weights of the same size as the sample size N: if weighted version is desired, =FALSE, otherwise (other generalized model)
the number of cores to use for parallel computation
dev: a vector of deviance after C permutations (length OF this vector is C)