User: wewolski

gravatar for wewolski
wewolski10
Reputation:
10
Status:
New User
Location:
Zurich
Last seen:
2 days, 12 hours ago
Joined:
2 years, 11 months ago
Email:
w*******@gmail.com

Posts by wewolski

<prev • 8 results • page 1 of 1 • next >
0
votes
1
answer
55
views
1
answer
Getting Significance for second factor (bioreplicate/subject) when using limma.
... Daer Forum, I want to get an estimate of the pairing variable BioReplicate What I would do when using linear model is: x1<-rnorm(10) x2<-1+rnorm(10) # Now create a dataframe for lme myDat <- data.frame(c(x1,x2), c(rep("x1", 10), rep("x2", 10)), rep(paste("S", seq(1,10), sep=""), 2)) na ...
limma anova written 24 days ago by wewolski10 • updated 23 days ago by Aaron Lun19k
0
votes
1
answers
161
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 6 months ago by wewolski10
1
vote
1
answer
161
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 6 months ago by wewolski10 • updated 6 months ago by Aaron Lun19k
0
votes
1
answers
514
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 9 months ago by wewolski10 • updated 9 months ago by Gordon Smyth33k
0
votes
1
answers
514
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 10 months ago by wewolski10 • updated 10 months ago by Gordon Smyth33k
2
votes
1
answer
514
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 10 months ago by wewolski10 • updated 9 months ago by Gordon Smyth33k
3
votes
2
answers
326
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 13 months ago by wewolski10 • updated 13 months ago by Aaron Lun19k
0
votes
6
answers
594
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 18 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: 141 users visited in the last hour