User: Andrew_McDavid

Reputation:
80
Status:
Trusted
Location:
Website:
https://amcdavid.githu...
Twitter:
EquivMeasures
Last seen:
2 weeks, 4 days ago
Joined:
1 year, 2 months ago
Email:
A*************@urmc.rochester.edu

Posts by Andrew_McDavid

<prev • 15 results • page 1 of 2 • next >
0
votes
1
answers
106
views
1
answers
Comment: C: Differential Expression analysis between distinct scRNA-seq datasets
... This all depends on how you define "batch."  As I understand your design, statistics won't be able to tell you if any differences you see should be ascribed to "treatment" or "plate".  Unless you run some replicates to estimate, or bound, the plate-to-plate variability, these two factors are confoun ...
written 18 days ago by Andrew_McDavid80
0
votes
2
answers
85
views
2
answers
Answer: A: What is the most useful distance measure to use for cell clustering based on MAS
... The L2-norm makes some sense to me, since the residuals may be approximately N(0,1) distributed under a null distribution of exchangeable cells, given the covariates included in the model. The clustering in this case could seen as searching for latent Gaussian structure among the residuals?  The opt ...
written 6 weeks ago by Andrew_McDavid80
0
votes
1
answers
159
views
1
answers
Comment: C: How do you test a contrast in MAST?
... Can you be more specific when you say "interface for contrasts" The contrasts(colData(x)$strain) = "contr.sum" business is base R and determines the coding for your factors, hence interpretation of your coefficients.  Venables & Ripley is the canonical reference for that.   Specifying linear co ...
written 9 weeks ago by Andrew_McDavid80
0
votes
1
answers
159
views
1
answers
Answer: A: How do you test a contrast in MAST?
... You can fit such a model a "cellmean model" (zlm_cellmean = zlm(~ 0 + group, x)) though it is possibly a little less statistically efficient [1] than using contrasts that explicitly includes an intercept: ## assuming you manually crossed strain and condition to get group zlm_dummy = zlm(~stra ...
written 9 weeks ago by Andrew_McDavid80
0
votes
1
answers
375
views
1
answers
Answer: A: NAs in continuous coefficient and fold change in `zlm`
... TL;DR One of the classes you are comparing has no detectable expression for a gene. The log-fold change is trying to get an estimate of effect size, but that is not well-defined when we never observe expression in a class. P-values from lrTest or waldTest are still well-defined. Fine print From t ...
written 3 months ago by Andrew_McDavid80
0
votes
1
answer
375
views
1
answer
NAs in continuous coefficient and fold change in `zlm`
... I noticed that MAST sometimes has NA for the continuous coefficient. When this happens, it also has NA for the estimated fold change. I couldn't find much in the documentation about this, but I was wondering if someone can explain it? This is a problem for me because I'm trying to select the most d ...
differential gene expression scrnaseq mast written 3 months ago by Andrew_McDavid80
0
votes
3
answers
330
views
3
answers
Answer: A: Estimating cell size factors for single cell RNA-seq data
... You might also try Bacher R, Chu L, Leng N, Gasch AP, Thomson JA, Stewart RM, Newton MA and Kendziorski C (2017). “SCnorm: robust normalization of single-cell RNA-seq data.” Nature Methods. https://www.nature.com/nmeth/journal/vaop/ncurrent/full/nmeth.4263.html. Software is available in Bioconduct ...
written 3 months ago by Andrew_McDavid80
0
votes
1
answers
576
views
1
answers
Answer: A: Intepretation of Z scores for gene set enrichment analysis in MAST
... The model coefficients that are being tested for enrichment are the coefficients from the zlmfit objects. With A as a ref, then positive coefficients in the ZlmFit object means that the mean of B > mean of A. The GSEA is calculating the average ZlmFit coefficient inside the set. So a positive Z-s ...
written 3 months ago by Andrew_McDavid80
0
votes
1
answers
116
views
1
answers
Comment: C: Modeling Dose-response data in MAST
... Our goal is to support any model that you can run with `lm`.  Using the numeric values for your ordered factor corresponds to a very particular model: it will estimate the linear effect of dose.  That is a restriction from the factor model, which allows the mean function to follow *any* effect of do ...
written 3 months ago by Andrew_McDavid80
0
votes
1
answers
116
views
1
answers
Answer: A: Modeling Dose-response data in MAST
... First off, in base R, I have only dabbled with the ordered factor contrasts so don't know them well enough to give reliable advice about how to interpret them. That being said, there are two ways to apply them with MAST. First, you can set the contrasts attribute to the column of interest: librar ...
written 3 months ago by Andrew_McDavid80

Latest awards to Andrew_McDavid

Scholar 7 months ago, created an answer that has been accepted. For A: MAST reported warning "Coefficients ... are never estimible and will be dropped.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 198 users visited in the last hour