User: Ahdee

gravatar for Ahdee
Ahdee40
Reputation:
40
Status:
New User
Location:
United States
Last seen:
2 weeks, 2 days ago
Joined:
3 years, 9 months ago
Email:
s********@gmail.com

Posts by Ahdee

<prev • 38 results • page 1 of 4 • next >
0
votes
1
answer
85
views
1
answers
Comment: C: Limma model design: is it better to just build a separate design for each specif
... Thanks James and Gordon! This is very reassuring - I been using single models ( that is analyzing everything together and using contrast to parse out the different comparison ) however started to freak out mostly because the logfc were different. I though that logFC was simply calculated from the j ...
written 6 weeks ago by Ahdee40
0
votes
1
answer
85
views
1
answers
Comment: C: Limma model design: is it better to just build a separate design for each specif
... Thanks James, I should be more specific. This is for an RNA-seq that has been TMM normalized and voom. I'm using voom as the input. So in general what you are saying is that leaving everything as a single model and defining the contrast is "better" but only if the variance between the samples are n ...
written 6 weeks ago by Ahdee40
3
votes
1
answer
85
views
1
answer
Limma model design: is it better to just build a separate design for each specific comparison?
... Hi suppose I have these 3 groups, a b and c. data= c ( rep("a", 3), rep("b",3), rep("c",3) ) mm = model.matrix(~0 + data ) colnames( mm ) = c("a","b","c") makeContrasts( avsb = a - b, levels=mm) a b c 1 1 0 0 2 1 0 0 3 1 0 0 4 0 1 0 5 0 1 0 ...
limma written 6 weeks ago by Ahdee40 • updated 6 weeks ago by James W. MacDonald50k
0
votes
1
answer
102
views
1
answers
Comment: C: Why does limma voom E value give different result than the voom object
... @thokall oh yes that right. I read it more carefully now whereby I thought previously that the E value were not only normalized but weighted as well, my bad, thanks! ...
written 8 weeks ago by Ahdee40
0
votes
1
answer
194
views
1
answers
Comment: C: Seurat 2.3.4 installation
... Hi you could also install directly from github. However first you need to `remove.packages("Seurat")` then make sure that you have install.packages("devtools") library(devtools); devtools::install_github(repo = 'satijalab/seurat', ref = 'release/2.4', dependencies=T) make sure to rest ...
written 8 weeks ago by Ahdee40
1
vote
1
answer
102
views
1
answer
Why does limma voom E value give different result than the voom object
... Hi so I notice something interesting. data <- voom(x, plot=T, design = design2) vfit <- lmFit(data) vfit <- eBayes(vfit) topTable(vfit,coef=2,sort.by="P") vfit <- lmFit(data$E, design = design2) vfit <- eBayes(vfit) topTable(vfit,coef=2,sort.by="P" ...
limma voom rna-seq written 9 weeks ago by Ahdee40 • updated 8 weeks ago by thokall160
0
votes
0
answers
75
views
0
answers
Is there a package to convert DNA coordinates to Cytogenetic location?
... Hi I was wondering if there is package to draw chromosome and/or list out cytogenetic locations based on DNA location. For example here is hg19, chr15:69,373,210-69,373,239. thanks! A ...
annotation coordinates cytogenetic written 3 months ago by Ahdee40
0
votes
0
answers
84
views
0
answers
Help with mapping Kegg hsa01100 using pathview
... Dear all, recently I did an enrichment analysis showing that my RNA.seq data is enriched for the global pathway hsa01100, however on trying to plot this with pathview it does not seem to highlight the genes. I found a similar post previously but its not exactly what I wanted since I'm interested i ...
pathview kegg pathway written 4 months ago by Ahdee40
1
vote
1
answer
169
views
1
answer
benefits of pearson distance for hierarchical clustering
... Dear all, I notice recently that  the heatmap3 package uses pearson distance instead of the default euclidean/manhattan. as.dist(1 - cor(df, use = "pa"))   Is there a benefit of using correlation instead of euclidean when it comes to calculating distance? The reason why I ask is because the hea ...
heatmap heatmap.3 written 7 months ago by Ahdee40 • updated 7 months ago by James W. MacDonald50k
0
votes
0
answers
167
views
0
answers
(Closed) Help with Censoring using the survival package.
... Hi all I would like to do some survival analysis using the survival package. The most difficult part for me is setting up the censoring data correctly and I would appreciated if someone can help confirm if I'm doing this correctly. According to the ?Surv function. "For interval censored data, the s ...
survival kaplan-meier survival analysis hazard ratio written 11 months ago by Ahdee40

Latest awards to Ahdee

Popular Question 16 months ago, created a question with more than 1,000 views. For Combat generating negative values, should I log2 transform first?
Popular Question 2.1 years ago, created a question with more than 1,000 views. For Combat generating negative values, should I log2 transform first?
Popular Question 3.3 years ago, created a question with more than 1,000 views. For Illumina MethylationEPIC vs. 450K

Help
Access

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