User: wewolski

gravatar for wewolski
wewolski10
Reputation:
10
Status:
New User
Location:
Zurich
Last seen:
2 months, 1 week ago
Joined:
2 years, 7 months ago
Email:
w*******@gmail.com

Posts by wewolski

<prev • 7 results • page 1 of 1 • next >
0
votes
1
answers
118
views
1
answers
Comment: C: limma - contrast.fit returns coefficients without names
... Thank you! [ ] default behaviour strikes again! "Omitting drop = FALSE when subsetting matrices and data frames is one of the most common sources of programming errors. (It will work for your test cases, but then someone will pass in a single column data frame and it will fail in an unexpected and ...
written 10 weeks ago by wewolski10
1
vote
1
answer
118
views
1
answer
limma - contrast.fit returns coefficients without names
... I have the following code: fit <- lmFit(intmat , designMatrix) tmpfit <- fit[,-1] lmfit.cont <- contrasts.fit(fit[,-1], cont[-1,2]) lmfitebayes <- eBayes(lmfit.cont)​ topTable(lmfitebayes, coef=name, number=Inf) Which fails when executing topTable sometimes. Thats because contrast. ...
limma written 12 weeks ago by wewolski10 • updated 12 weeks ago by Aaron Lun18k
0
votes
1
answers
456
views
1
answers
Comment: C: Why does limma contrast.fit produces NA's for all contrast although missing data
... "Is it possible that the expression values have simply failed to meet some detection threshold for this condition?" No, it is not. Let me assure you, I did try to understand and model where the missingness is coming from.  If I used the imputation you suggest I would introduce a heavy bias. Bias is ...
written 5 months ago by wewolski10 • updated 5 months ago by Gordon Smyth32k
0
votes
1
answers
456
views
1
answers
Comment: A: Why does limma contrast.fit produces NA's for all contrast although missing data
... Thanks, Gordon, I have now a solution based on the post:  https://support.bioconductor.org/p/19879/ but it is a hack, workaround, and features usually do not need workarounds.  Still, if you think that contrast.fit not handling NA is a feature why not just add contrast.fit.NA to limma which is a ...
written 6 months ago by wewolski10 • updated 6 months ago by Gordon Smyth32k
2
votes
1
answer
456
views
1
answer
Why does limma contrast.fit produces NA's for all contrast although missing data only in one condition?
... I will try to describe briefly what I am doing: my design matrix designMatrix <- model.matrix( ~ 0 + Condition + Plant, data=grp2$annotation_) > unique(grp2$annotation_$Condition) [1] "PC_16h_sys"  "PC_16h_test" "PC_38h_sys"  "PC_38h_test" "PC_96h_sys"  "PC_96h_test" Then I define the con ...
limma contrast lmfit missing data makecontrasts written 6 months ago by wewolski10 • updated 5 months ago by Gordon Smyth32k
3
votes
2
answers
241
views
2
answers
limma for Paired Samples vs. lme
... Hi, Just comparing the output of limma and lme. The model that I fit with lme is : le <- lme(lmIntensity ~ Condition, random = ~1|Donor, data = x) Since I have repeated measures (Condition) on the same Donor and I am not interested in the Donor effect I model it as random effect. In lim ...
limma linear model paired samples written 9 months ago by wewolski10 • updated 9 months ago by Aaron Lun18k
0
votes
6
answers
484
views
6
answers
Comment: C: matching of AAStringSet vs. another AAStringSet
... Hi Tobi, Here you have a performant and IMHO easier to use solution based on pure CRAN packages seqinr and AhoCorasickTrie. use : seqinr::read.fasta(file = file, as.string = TRUE, seqtype="AA") to read your fasta. then perform the matching with: system.time(res <- AhoCorasickTrie::AhoCorasi ...
written 14 months ago by wewolski10

Latest awards to wewolski

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 2.2.0
Traffic: 124 users visited in the last hour