a1 <- proc.time() source("/seppdata/sepp/linkage/release/split/run_all.R") a2 <- proc.time() Tdiff <- a2-a1 user system elapsed mm <- (c(42310.496, 24.732, 42406.355) + c(41457.897, 27.561, 41556.689) + c(42588.23, 34.35, 42696.52) + c(41029.422, 26.162, 41125.169) + c(41812.917, 25.188, 41912.710))/5 > mm user system elapsed [1] 41839.7924 27.5986 41939.4886 > mm/3600 [1] 11.622164556 0.007666278 11.649857944 11.65 hours > ((mm/3600)/100)*639 [1] 74.26563151 0.04898752 74.44259226 74 hours for whole analysis (3 days two hours) frameshift deletion frameshift insertion 3724 3430 frameshift substitution nonframeshift deletion 49 2126 nonframeshift insertion nonframeshift substitution 505 17 nonsynonymous SNV stopgain SNV 275837 5458 stoploss SNV synonymous SNV 353 194307 unknown 4565 c("frameshift deletion","frameshift insertion","frameshift substitution","nonframeshift deletion","nonframeshift insertion","nonframeshift substitution","nonsynonymous SNV","stopgain SNV","stoploss SNV","synonymous SNV","unknown") > table(variant$V1) downstream exonic exonic;splicing 248906 482226 8241 intergenic intronic ncRNA_exonic 23271824 15019781 75702 ncRNA_intronic ncRNA_splicing ncRNA_UTR3 928810 374 5385 ncRNA_UTR5 ncRNA_UTR5;ncRNA_UTR3 splicing 833 2 4823 upstream upstream;downstream UTR3 215458 8386 333952 UTR5 UTR5;UTR3 53666 166