diff --git a/bash.scripts/sc_FullAnalysis-DPrF-Lin.sh b/bash.scripts/sc_FullAnalysis-DPrF-Lin.sh new file mode 100644 index 0000000000000000000000000000000000000000..d9b089baf2fd93fadce5bd8bb778a2e7ef9e2036 --- /dev/null +++ b/bash.scripts/sc_FullAnalysis-DPrF-Lin.sh @@ -0,0 +1,37 @@ +#!/bin/bash +# +# CREATED USING THE BIOHPC PORTAL on Thu Nov 30 2017 16:16:26 GMT-0600 (Central Standard Time) +# +# This file is batch script used to run commands on the BioHPC cluster. +# The script is submitted to the cluster using the SLURM `sbatch` command. +# Lines starting with # are comments, and will not be run. +# Lines starting with #SBATCH specify options for the scheduler. +# Lines that do not start with # or #SBATCH are commands that will run. + +# Name for the job that will be visible in the job queue and accounting tools. +#SBATCH --job-name R_FullAnalysis + +# Name of the SLURM partition that this job should run on. +#SBATCH -p 256GB,256GBv1 # partition (queue) +# Number of nodes required to run this job +#SBATCH -N 1 + +# Time limit for the job in the format Days-H:M:S +# A job that reaches its time limit will be cancelled. +# Specify an accurate time limit for efficient scheduling so your job runs promptly. +#SBATCH -t 7-0:0:0 + +# The standard output and errors from commands will be written to these files. +# %j in the filename will be replace with the job number when it is submitted. +#SBATCH -o job_%j.out +#SBATCH -e job_%j.err + +# Send an email when the job status changes, to the specfied address. +#SBATCH --mail-type ALL +#SBATCH --mail-user gervaise.henry@utsouthwestern.edu + +module load R/3.4.1-gccmkl + +# COMMAND GROUP 1 +Rscript ../r.scripts/sc_QuSAGE.Lineage.R --p="DPrF" +# END OF SCRIPT diff --git a/bash.scripts/sc_FullAnalysis-DPrF.sh b/bash.scripts/sc_FullAnalysis-DPrF.sh new file mode 100755 index 0000000000000000000000000000000000000000..136c5fce9745ed8718ecc248c413fbc29fc42806 --- /dev/null +++ b/bash.scripts/sc_FullAnalysis-DPrF.sh @@ -0,0 +1,50 @@ +#!/bin/bash +# +# CREATED USING THE BIOHPC PORTAL on Thu Nov 30 2017 16:16:26 GMT-0600 (Central Standard Time) +# +# This file is batch script used to run commands on the BioHPC cluster. +# The script is submitted to the cluster using the SLURM `sbatch` command. +# Lines starting with # are comments, and will not be run. +# Lines starting with #SBATCH specify options for the scheduler. +# Lines that do not start with # or #SBATCH are commands that will run. + +# Name for the job that will be visible in the job queue and accounting tools. +#SBATCH --job-name R_FullAnalysis + +# Name of the SLURM partition that this job should run on. +#SBATCH -p 256GB,256GBv1 # partition (queue) +# Number of nodes required to run this job +#SBATCH -N 1 + +# Time limit for the job in the format Days-H:M:S +# A job that reaches its time limit will be cancelled. +# Specify an accurate time limit for efficient scheduling so your job runs promptly. +#SBATCH -t 7-0:0:0 + +# The standard output and errors from commands will be written to these files. +# %j in the filename will be replace with the job number when it is submitted. +#SBATCH -o job_%j.out +#SBATCH -e job_%j.err + +# Send an email when the job status changes, to the specfied address. +#SBATCH --mail-type ALL +#SBATCH --mail-user gervaise.henry@utsouthwestern.edu + +module load R/3.4.1-gccmkl + +# COMMAND GROUP 1 +Rscript ../r.scripts/sc_Demultiplex.R --p="DPrF" --d=7 +#Rscript ../r.scritps/sc_D-SampleReorder.R +#Rscript ../r.scripts/sc_Seurat.Score.CellCycle.R --p="DPrF" +Rscript ../r.scripts/sc_QC.R --p="DPrF" --cc=FALSE +Rscript ../r.scripts/sc_Cluster.R --p="DPrF" +Rscript ../r.scripts/sc_PC.Score.Stress.R --p="DPrF" +#Rscript ../r.scripts/sc_QuSAGE.Lineage.R --p="DPrF" +#Rscript ../r.scripts/sc_LineageSubClust.R --p="DPrF" +#Rscript ../r.scripts/sc_QuSAGE_EpiSubClust.R --p="DPrF" +#Rscript ../r.scripts/sc_QuSAGE_StSubClust.R --p="DPrF" +#Rscript ../r.scripts/sc_MergeSubClust.R --p="DPrF" +#Rscript ../r.scripts/sc_PC.Score.NE.R --p="DPrF" +#Rscript ../r.scripts/sc_DEG.R --p="DPrF" +#Rscript ../r.scripts/sc_Tables.R --p="DPrF" +# END OF SCRIPT diff --git a/r.scripts/sc_QuSAGE.Lineage.R b/r.scripts/sc_QuSAGE.Lineage.R index 9c287948776f7ae6fd00dcde854f570c75a94534..076e212d0fae5cc1bedc576d6f4d817e3b7da2cc 100755 --- a/r.scripts/sc_QuSAGE.Lineage.R +++ b/r.scripts/sc_QuSAGE.Lineage.R @@ -55,8 +55,8 @@ if (!is.na(opt$s)){ } else { rnd <- 1:ncol(sc10x.Group@data) } -eset <- as.data.frame(as.matrix(sc10x.Group@data)) -eset <- eset[,rnd] +eset <- sc10x.Group@data[,rnd] +eset <- as.data.frame(as.matrix(eset)) labels <- paste0("Cluster_",as.vector(factor(sc10x.Group@ident))) labels <- labels[rnd] rm(rnd) diff --git a/r.scripts/sc_QuSAGE_EpiSubClust.R b/r.scripts/sc_QuSAGE_EpiSubClust.R index 91fd7efd7c6b68be08045eada4af58a0643cf923..be4edef44c00f2c9110b656f6486373db3658239 100755 --- a/r.scripts/sc_QuSAGE_EpiSubClust.R +++ b/r.scripts/sc_QuSAGE_EpiSubClust.R @@ -62,8 +62,8 @@ if (!is.na(opt$s)){ } else { rnd <- 1:ncol(sc10x.Group.Epi@data) } -eset <- as.data.frame(as.matrix(sc10x.Group.Epi@data)) -eset <- eset[,rnd] +eset <- sc10x.Group@data[,rnd] +eset <- as.data.frame(as.matrix(eset)) labels <- paste0("Cluster_",as.vector(factor(sc10x.Group.Epi@ident))) labels <- labels[rnd] rm(rnd) diff --git a/r.scripts/sc_QuSAGE_StSubClust.R b/r.scripts/sc_QuSAGE_StSubClust.R index 586f2bf79732269365e0dff15be742f1ecbfac4f..6409726b5104640e96fab6153f2f6ea3399b65d0 100755 --- a/r.scripts/sc_QuSAGE_StSubClust.R +++ b/r.scripts/sc_QuSAGE_StSubClust.R @@ -60,8 +60,8 @@ if (!is.na(opt$s)){ } else { rnd <- 1:ncol(sc10x.Group.St@data) } -eset <- as.data.frame(as.matrix(sc10x.Group.St@data)) -eset <- eset[,rnd] +eset <- sc10x.Group@data[,rnd] +eset <- as.data.frame(as.matrix(eset)) labels <- paste0("Cluster_",as.vector(factor(sc10x.Group.St@ident))) labels <- labels[rnd] rm(rnd)