User: PJ

gravatar for PJ
PJ0
Reputation:
0
Status:
New User
Last seen:
2 years, 7 months ago
Joined:
4 years, 3 months ago
Email:
p******@colostate.edu

Profile information, website and location are not shown for new users.

This helps us discourage the inappropriate use of our site.

Posts by PJ

<prev • 9 results • page 1 of 1 • next >
0
votes
1
answer
612
views
1
answers
Comment: C: metagenomeSeq log transformation
... Yes, you are right. I added a second normalization step based on reference gene length and 16S bacterial reads so I am not using the read counts. Could this affect the statistical output? Where can I see the actual formula to calculate average expression? Thanks again! ...
written 2.7 years ago by PJ0
0
votes
1
answer
612
views
1
answers
Comment: C: metagenomeSeq log transformation
... Thanks Joseph for the quick answer. What about the negative tansformed log values? Is the program taking the absolute values for statistics?  ...
written 2.7 years ago by PJ0
0
votes
1
answer
612
views
1
answers
Comment: C: metagenomeSeq log transformation
... In addition, how is it possible that features with negative log2 values in the transformed count matrix come up with positive average expression values in the statistical output? Thanks! ...
written 2.7 years ago by PJ0
0
votes
1
answer
612
views
1
answer
metagenomeSeq log transformation
... Hi, Normalized counts from the matrix table are transformed to log2 values by metagenomeseq, aren´t they? My question is about the zeros in the count matrix. How are they transformed to allow statistical analysis? (as log2 of zero gives an error). Thanks! ...
metagenomeseq logtransformation written 2.7 years ago by PJ0 • updated 2.7 years ago by Joseph Nathaniel Paulson270
1
vote
1
answer
692
views
1
answer
limma interaction term
... Hi everybody, I have a model with 2 main factors and its interaction: designmodel = model.matrix(~0 + Factor1 + Factor2 + Factor1:Factor2) Factor 1 has two levels and factor 2 has 4 levels. My question is if when I test the interaction term I can add all 8 groups as follow:  contrastInteraction ...
limma interaction model written 3.2 years ago by PJ0 • updated 3.2 years ago by Gordon Smyth38k
0
votes
2
answers
780
views
2
answers
Random effects metagenomeSeq
... Does metagenomeSeq handle mixed models including random effects? If yes, how can I incorporate those random effects in the model?Thanks! ...
model.matrix written 3.7 years ago by PJ0 • updated 3.6 years ago by noelle.noyes30
0
votes
0
answers
495
views
0
answers
Model in metagenomeSeq
... Hi, My experimental model is repeated measures with 2 fixed effects and 1 random effect. How can I specify that in the model statement? I am having a hard time trying to figure it out. Thanks!  ...
repeated measures written 3.9 years ago by PJ0
0
votes
0
answers
592
views
0
answers
metagenomeSeq heatmap error
... Hi,  I am creating a heat map based on the following script: trials=pData(object)$treatment heatmapColColors=brewer.pal(12, "Set3")[as.integer(factor(trials))] heatmapCols=colorRampPalette(brewer.pal(9, "RdBu"))(50) plotMRheatmap(obj=object, n=6, cexRow=0.4, cexCol=0.4, trace="none", col=heatmap ...
heatmap written 4.3 years ago by PJ0
0
votes
1
answer
857
views
1
answer
model.matrix in metagenomeSeq
... Hi, I am running fitZig for differential abundance testing. When I place: treatment=pData(obj)$Treatment Before=pData(obj)$Time normFactor=normFactors(obj) normFactor=log2(normFactor/median(normFactor)+1) mod = model.matrix (Treatment + Before + normFactor), the answer is: Error: $ operator i ...
fitzig written 4.3 years ago by PJ0 • updated 4.3 years ago by Hector Corrada Bravo60

Latest awards to PJ

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: 297 users visited in the last hour