Commit ad15b69b authored by Gervaise H. Henry's avatar Gervaise H. Henry 🤠

Prep for muPr_DWS.CVctrl (lower res for lin id to 0.1)

parent f1bf8579
......@@ -28,7 +28,7 @@ module load R/3.6.1-gccmkl
if [[ "$1" == "huPr_PdPgb" ]] || [[ "$1" == "huPr_PdPb" ]] || [[ "$1" == "huBl" ]] || [[ "$1" == "huPr_PdPgsb" ]]
then
Rscript ../r.scripts/SingleR.R --p "$1" --s "$2" --o "$3"
elif [[ "$1" == "muPrUr" ]]
elif [[ "$1" == "muPrUr" ]] || [[ "$1" == "muPr_DWS.CVctrl" ]]
then
Rscript ../r.scripts/huPr_muPr.R --p "$1" --r "$3"
fi
......
......@@ -76,8 +76,8 @@ lin.se <- lin.se[common,]
sc10x.se <- sc10x.se[common,]
rm(common)
singler.lin <- SingleR(sc10x.se,ref=lin.se,method="cluster",clusters=sc10x.se$integrated_snn_res.1,labels=lin.se$label.fine,BPPARAM=MulticoreParam(workers=10))
sc10x$lin <- singler.lin$labels[match(sc10x.se$integrated_snn_res.1,singler.lin@rownames)]
singler.lin <- SingleR(sc10x.se,ref=lin.se,method="cluster",clusters=sc10x.se$integrated_snn_res.0.1,labels=lin.se$label.fine,BPPARAM=MulticoreParam(workers=10))
sc10x$lin <- singler.lin$labels[match(sc10x.se$integrated_snn_res.0.1,singler.lin@rownames)]
#singler.lin <- SingleR(sc10x.se,ref=lin.se,method="single",labels=lin.se$label.main,BPPARAM=MulticoreParam(workers=10))
#sc10x$lin <- singler.lin$labels
sc10x$lin[sc10x$lin %in% unique(lin.se@colData[lin.se@colData[,"label.main"] %in% c(
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment