User: hmgeiger

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

Posts by hmgeiger

<prev • 13 results • page 1 of 2 • next >
1
vote
1
answer
99
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 6 months ago by hmgeiger10 • updated 6 months ago by Michael Love23k
0
votes
1
answer
488
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 19 months ago by hmgeiger10
0
votes
1
answer
488
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 20 months ago by hmgeiger10
0
votes
1
answer
488
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 20 months ago by hmgeiger10
0
votes
1
answer
488
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 20 months ago by hmgeiger10
0
votes
1
answer
488
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 20 months ago by hmgeiger10
0
votes
1
answer
488
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 20 months ago by hmgeiger10
0
votes
1
answer
488
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 20 months ago by hmgeiger10
0
votes
4
answers
944
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.2 years ago by hmgeiger10 • updated 3.2 years ago by Michael Love23k
0
votes
4
answers
944
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.2 years ago by hmgeiger10

Latest awards to hmgeiger

No awards yet. Soon to come :-)

Help
Access

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