') + ')', 'gi'); if (regex.test(text)) { found = true; var frag = document.createDocumentFragment(); var parts = text.split(regex); parts.forEach(function(part, i) { if (i % 2 === 0) { frag.appendChild(document.createTextNode(part)); } else { var span = document.createElement('span'); span.className = 'userscript-highlight'; span.textContent = part; frag.appendChild(span); } }); node.parentNode.replaceChild(frag, node); } }); } else if (node.nodeType === 1 && node.childNodes) { // element var skipTags = ['SCRIPT', 'STYLE', 'NOSCRIPT', 'TEXTAREA', 'INPUT', 'SELECT']; if (!skipTags.includes(node.tagName)) { Array.from(node.childNodes).forEach(highlight); } } } highlight(document.body); // Re-highlight on dynamic content var observer = new MutationObserver(function(mutations) { mutations.forEach(function(m) { m.addedNodes.forEach(function(node) { if (node.nodeType === 1 || node.nodeType === 3) highlight(node); }); }); }); observer.observe(document.body, { childList: true, subtree: true }); })(); } } catch(__e) { console.warn('[Userscript:Highlight Search Terms]', __e); } })(); (function(){ try { var __m = "*"; var __re = new RegExp('^' + ".*" + ', 'i'); if (__m === '*' || __re.test(location.href)) { // Strip utm_, fbclid, gclid, etc. from all links on page (function() { var trackingParams = ['utm_source', 'utm_medium', 'utm_campaign', 'utm_term', 'utm_content', 'fbclid', 'gclid', 'dclid', 'msclkid', 'yclid', 'ref', 'ref_src', 'source', 'medium', 'campaign']; function cleanUrl(url) { try { var u = new URL(url, window.location.origin); var changed = false; trackingParams.forEach(function(p) { if (u.searchParams.has(p)) { u.searchParams.delete(p); changed = true; } }); return changed ? u.toString() : url; } catch (e) { return url; } } function cleanLinks() { document.querySelectorAll('a[href]').forEach(function(a) { var clean = cleanUrl(a.href); if (clean !== a.href) a.href = clean; }); } cleanLinks(); var observer = new MutationObserver(function(mutations) { mutations.forEach(function(m) { m.addedNodes.forEach(function(node) { if (node.nodeType === 1) { if (node.tagName === 'A') cleanLinks(); node.querySelectorAll('a[href]').forEach(function(a) { var clean = cleanUrl(a.href); if (clean !== a.href) a.href = clean; }); } }); }); }); observer.observe(document.body, { childList: true, subtree: true }); })(); } } catch(__e) { console.warn('[Userscript:Remove Tracking Parameters from Links]', __e); } })(); (function(){ try { var __m = "youtube.com"; var __re = new RegExp('^' + "youtube\\.com" + ', 'i'); if (__m === '*' || __re.test(location.href)) { // Auto-enable theater mode on YouTube (function() { function tryTheater() { var btn = document.querySelector('button[aria-label="Theater mode"], ytd-player #player button[title="Theater mode"]'); if (btn && !btn.classList.contains('activated')) { btn.click(); } } // Try immediately tryTheater(); // Try after navigation (SPA) var lastUrl = location.href; setInterval(function() { if (location.href !== lastUrl) { lastUrl = location.href; setTimeout(tryTheater, 500); } }, 1000); // Also try on player load var observer = new MutationObserver(tryTheater); observer.observe(document.body, { childList: true, subtree: true }); })(); } } catch(__e) { console.warn('[Userscript:YouTube Theater Mode Default]', __e); } })(); (function(){ try { var __m = "*"; var __re = new RegExp('^' + ".*" + ', 'i'); if (__m === '*' || __re.test(location.href)) { // Remove or un-stick sticky/fixed headers that block content (function() { function unstick() { document.querySelectorAll('header, nav, [role="banner"], .header, .navbar, .sticky, .fixed-top, [style*="position: fixed"], [style*="position:sticky"]').forEach(function(el) { if (el.style.position === 'fixed' || el.style.position === 'sticky' || getComputedStyle(el).position === 'fixed' || getComputedStyle(el).position === 'sticky') { el.style.position = 'static'; el.style.top = 'auto'; el.style.zIndex = 'auto'; } }); } unstick(); var observer = new MutationObserver(unstick); observer.observe(document.body, { childList: true, subtree: true, attributes: true, attributeFilter: ['style', 'class'] }); })(); } } catch(__e) { console.warn('[Userscript:Kill Sticky Headers]', __e); } })(); })(); GitHub - zcslab/SSMD: A semi-supervised approach for a robust identification of cell types and deconvolution of mouse transcriptomics data https://zcslab.github.io/SSMD/ · GitHub
Skip to content
BranchesTags

Repository files navigation

SSMD: A semi-supervised approach for a robust identification of cell types and deconvolution of mouse transcriptomics data

Description

Multiple deconvolution methods have been developed for investigating the heterogeneous immune and stromal (I/S) cell types in human cancer tissue to estimate their relative abundances using transcriptomic data. However, there is a lack of a robust method and user-friendly software for mouse transcriptomic data deconvolution. Here, we developed a novel semi-supervised approach, namely SSMD, by (i) deriving potential I/S cell signature genes from a large collection of mouse data sets to form a marker labeling matrix; (ii) implementing a rank-1 sub matrix identification method to test the presence of I/S cell types and identify data set specific I/S cell markers; and (iii) utilizing a constrained non-negative matrix factorization (NMF) based framework to account for diversity of mouse models. The new method was validated on single cell RNA-seq simulated bulk tissue data and independent immuno-assay data. The method is applied to mouse prostate cancer data sets to infer the level of anti-cancer immune cell populations.

Installation

install.packages("devtools")
devtools::install_github("zy26/SSMD")

Usage

estimate.proportion <- function(data, lambda = lambda)

Arguments

  • data input gene expression matrix. MGI gene symbol should be as their row names
  • parameter threshold of mean correlation to define rank-1 co-expression module

Value

An object of class is also invisibly returned. This is a list containing the following components:

  • Stat_all statistics for all rank-1 co-expression module. CT: cell type; mean: mean correlation inside the module; Core_overlap_number: Overlap number with core marker list; Core_overlap_rate: overlap rate with core marker list; BCV_rank: bcv rank of the first base
  • module_keep modules with the high overlap number with core marker list for each cell type
  • proportion estimated proportion for each cell type

Examples

#load your own gene expression data
load('example_bulk.RData')
estimate.proportion(data, lambda = 0.8)

About

A semi-supervised approach for a robust identification of cell types and deconvolution of mouse transcriptomics data https://zcslab.github.io/SSMD/

Resources

Readme
Activity
Custom properties

Stars

4 stars

Watchers

0 watching

Forks

0 forks
Report repository

Releases

Packages

Contributors

Languages