List lsf job groups bjgroup -s /username
Make a new one: bgadd -L <job_limit> <name>
Alter the limit later: bgmod
Cancer_type | Lifetime_cancer_incidence | Total_cells_tissue | Total_Stem_Cells | Stem_cell_divisions_per_year | Stem_cell_divisions_per_lifetime | LSCD | |
---|---|---|---|---|---|---|---|
ALL | 0.0041 | 3000000000000 | 135000000 | 12 | 960 | 129900000000 | |
BCC | 0.3 | 180000000000 | 5820000000 | 7.6 | 608 | 3550000000000 | |
CLL | 0.0052 | 3000000000000 | 135000000 | 12 | 960 | 129900000000 | |
Colorectal | 0.048 | 30000000000 | 200000000 | 73 | 5840 | 1168000000000 | |
Colorectal_FAP | 1 | 30000000000 | 200000000 | 73 | 5840 | 1168000000000 | |
Colorectal_Lynch | 0.5 | 30000000000 | 200000000 | 73 | 5840 | 1168000000000 | |
Duodenum_adenocarcinoma | 0.0003 | 680000000 | 4000000 | 24 | 1947 | 7796000000 | |
Duodenum_adenocarcinoma_with_FAP | 0.035 | 680000000 | 4000000 | 24 | 1947 | 7796000000 | |
Esophageal_squamous_cell_carcinoma | 0.001938 | 3240000000 | 846000 | 17.4 | 1390 | 1203000000 |
List lsf job groups bjgroup -s /username
Make a new one: bgadd -L <job_limit> <name>
Alter the limit later: bgmod
df $path and df -i $path to check space and inodes |
rm -f simulated_hets.bam && python simulate_bam_reads.py ~/dat/ref/hs37d5/all_sequences.fa 0.5 1:100000 1:100200 1:10040 0 1:100600 1:100800 | samtools view -bh > simulated_hets.bam && samtools index simulated_hets.bam | |
rm -f simulated_hets.bam && python simulate_bam_reads.py ~/dat/ref/hs37d5/all_sequences.fa 1.0 1:100000 1:100200 1:10040 0 1:100600 1:100800 | samtools view -bh > simulated_homs.bam && samtools index simulated_homs.bam | |
rm -f cis_ase_tumor_rna.bam && python simulate_bam_reads.py ~/src/regtools/tests/integration-test/data/fa/test_chr22.fa 0 22:1050 22:1550 22:1950 22:2050 22:2550 22:2950 22:3050 22:3550 22:3950 | samtools view -bh > cis_ase_tumor_rna.bam && samtools index cis_ase_tumor_rna.bam | |
rm -f cis_ase_tumor_dna.bam && python simulate_bam_reads.py ~/src/regtools/tests/integration-test/data/fa/test_chr22.fa 0.5 22:1050 22:1550 22:1950 22:2050 22:2550 22:2950 22:3050 22:3550 22:3950 | samtools view -bh > cis_ase_tumor_dna.bam && samtools index cis_ase_tumor_dn |
sacct -u aru
JobID JobName Partition Account AllocCPUS State ExitCode
------------ ---------- ---------- ---------- ---------- ---------- --------
5718685 interacti+ interacti+ dcl 1 COMPLETED 0:0
5718686 interacti+ interacti+ dcl 1 COMPLETED 0:0
sacct -j 5719109 --format=JobID,JobName,MaxRSS,Elapsed
JobID JobName MaxRSS Elapsed
The default soft (warning) and hard (error) size limits are set to 500,000 and 10,000,000 bytes but can be overriden through the hooks.filesizesoftlimit and hooks.filesizehardlimit settings respectively:
$ git config hooks.filesizesoftlimit 100000
$ git config hooks.filesizehardlimit 4000000
> counts <- rnbinom(1000, size = 10, p = 0.1) | |
> mu <- mean(counts) | |
> p <- mu/var(counts) | |
> r <- mu ^ 2/(var(counts) - mu) | |
> hist(counts, freq = F) | |
> curve(dnbinom(x, size = r, prob = p), from = 0, to = max(counts), n = max(counts) + 1) |
first_time <- function() { | |
t1 <- read.table("../ALL.chrY.phase3_integrated_v1b.20130502.genotypes.vcf.gz") | |
print(nrow(t1)) | |
nucs <- c('A', 'C', 'G', 'T') | |
t1 <- t1[t1$V4 %in% nucs & t1$V5 %in% nucs, ] | |
print(nrow(t1)) | |
t2 <- t1[, 10:ncol(t1)] | |
t2 <- sapply(t2, as.character) | |
save(t2, file = "t2") | |
} |