Displaying 2 results from an estimated 2 matches for "nzw_smk".
Did you mean:
new_sk
2004 Dec 20
2
problems with limma
...as.factor(ts)
> levels(ts) <- c("nzwC","nzwT","akrC","akrT")
> design <- model.matrix(~0+ts)
> colnames(design) <- levels(ts)
> fit4 <- lmFit(all,design)
> cont.matrix <- makeContrasts(
+ Baseline = akrC - nzwC,
+ NZW_Smk = nzwT - nzwC,
+ AKR_Smk = akrT - akrC,
+ Diff = (akrT - akrC) - (nzwT - nzwC),
+ levels=design)
> fit42 <- contrasts.fit(fit4,cont.matrix)
> fit42 <- eBayes(fit42)
> #
> topTable(fit42,coef="Baseline",adjust="fdr",number=5)...
2004 Dec 21
0
Fwd: problems with limma
...s(ts) <- c("nzwC","nzwT","akrC","akrT")
>> design <- model.matrix(~0+ts)
>> colnames(design) <- levels(ts)
>> fit4 <- lmFit(all,design)
>> cont.matrix <- makeContrasts(
> + Baseline = akrC - nzwC,
> + NZW_Smk = nzwT - nzwC,
> + AKR_Smk = akrT - akrC,
> + Diff = (akrT - akrC) - (nzwT - nzwC),
> + levels=design)
>> fit42 <- contrasts.fit(fit4,cont.matrix)
>> fit42 <- eBayes(fit42)
>> #
>> topTable(fit42,coef="Baseline",adjust="fdr&...