User: hmgeiger

gravatar for hmgeiger
hmgeiger10
Reputation:
10
Status:
New User
Location:
Last seen:
11 months ago
Joined:
3 years, 7 months ago
Email:
h*******@nygenome.org

Posts by hmgeiger

<prev • 13 results • page 1 of 2 • next >
1
vote
1
answer
138
views
1
answer
Differential expression where only certain samples are paired
... I have the following experimental design (n=9). Trying to compare the two cell types. Celltype 1 from individuals A, B, and C Celltype 2 from individuals A, B, D, E, F, and G Is there any way to include in the design the fact that 2/3 cell type 1 samples have a corresponding paired sample in cell ...
deseq2 written 11 months ago by hmgeiger10 • updated 11 months ago by Michael Love25k
0
votes
1
answer
544
views
1
answers
Comment: C: Looking at interaction term in DESeq2 while also controlling for individual vari
... I also realized just now that I did not supply the model matrix to estimateDispersions. The code should be estimateDispersions(dds,modelMatrix = design_table_model_matrix). I imagine this could have affected the results as well. ...
written 2.0 years ago by hmgeiger10
0
votes
1
answer
544
views
1
answers
Comment: C: Looking at interaction term in DESeq2 while also controlling for individual vari
... Update: Also tried a super high number of iterations (maxit=10000) with only a few additional rows converging. ...
written 2.1 years ago by hmgeiger10
0
votes
1
answer
544
views
1
answers
Comment: C: Looking at interaction term in DESeq2 while also controlling for individual vari
... There are a few genes that are DE if you do not remove these rows, and for these the counts are often high in quite a few samples. For example, here are normalized counts for one gene that did not converge, for group1 and group3. Note, in the real data there are 4 individuals for group 1 and 9 indi ...
written 2.1 years ago by hmgeiger10
0
votes
1
answer
544
views
1
answers
Comment: A: Looking at interaction term in DESeq2 while also controlling for individual vari
... Finally got results that make sense. This is my code. Only issue left is I am still getting 206 rows that do not converge in beta even after strict expression filtering and increasing maxit (following the tips here: https://support.bioconductor.org/p/65091/). Any other tips on how to get them all t ...
written 2.1 years ago by hmgeiger10
0
votes
1
answer
544
views
1
answers
Comment: A: Looking at interaction term in DESeq2 while also controlling for individual vari
... Update: Tried running the following code afterwards, but get the same p-values no matter which two groups are being compared. Also there are way too many DEGs given the expected strength of the signal. Don't think the above design is correct, but not sure how to fix it. ddsClean <- dds[which(mc ...
written 2.1 years ago by hmgeiger10
0
votes
1
answer
544
views
1
answers
Comment: A: Looking at interaction term in DESeq2 while also controlling for individual vari
... That makes sense, although for the last part ("now this matrix m1 can be provided to the full argument of DESeq") I am still unclear on the syntax. Can you take a look at my code and let me know? First, I add the patient indices to the design table. patient_indices <- rep("Index",times=nrow(de ...
written 2.1 years ago by hmgeiger10
0
votes
1
answer
544
views
1
answer
Looking at interaction term in DESeq2 while also controlling for individual variable
... I am running an experiment where we have samples from a number of individuals, pre and post treatment. These individuals were then also grouped according to features found in the DNA. So, overall the design looks like this. I named pre treatment "before" so that DESeq2 will automatically take it as ...
deseq2 written 2.1 years ago by hmgeiger10
0
votes
4
answers
1.1k
views
4
answers
Comment: C: DESeq2 only runs interaction factor x continuous variable with old not new versi
... Got it! Thanks. This is how I ended up setting up my script. A is normal, B and C are my two test conditions. dds <- DESeq(DESeqDataSetFromMatrix(countData=featureCounts[,paste(1:24,"D",sep="")], colData=design,design=~Age + Diagnosis + Diagnosis:Age)) B_vs_control_plus_interaction_age <- ...
written 3.6 years ago by hmgeiger10 • updated 3.6 years ago by Michael Love25k
0
votes
4
answers
1.1k
views
4
answers
Answer: A: DESeq2 only runs interaction factor x continuous variable with old not new versi
... The investigator wants to look at the interaction between age and diagnosis. They wouldn't be opposed to grouping age into factors as appropriate I don't think, though I also wanted to just see if it could be done keeping age as a continuous variable. ...
written 3.6 years ago by hmgeiger10

Latest awards to hmgeiger

Popular Question 11 months ago, created a question with more than 1,000 views. For DESeq2 only runs interaction factor x continuous variable with old not new version

Help
Access

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