Question

Saving output from glmTreat to a csv file?

0

Entering edit mode

mictadlo ▴ 10

@mictadlo-10885

Last seen 4.9 years ago

Hi,

How is it possible save the output from glmTreat to a csv file?

> tr <- glmTreat(fit, contrast=B.LvsP, lfc=log2(1.5))
> topTags(tr)
Coefficient:  -1Leaves.2 1Leaves.3

          logFC unshrunk.logFC logCPM  PValue     FDR
sp0090975   2.6            2.6    2.8 5.7e-14 2.7e-09
sp0037632  -3.0           -3.0    3.4 1.7e-13 2.7e-09
sp0074153  -3.9           -3.9    3.8 1.8e-13 2.7e-09
sp0008306   3.2            3.2    2.9 1.8e-13 2.7e-09
sp0073530  -4.5           -4.5    3.4 2.3e-12 2.3e-08
sp0025713  -3.9           -3.9    4.4 2.6e-12 2.3e-08
sp0037721   7.8            8.0    2.4 2.7e-12 2.3e-08
sp0083660   2.0            2.0    4.4 3.2e-12 2.3e-08
sp0052245  -2.9           -2.9    4.9 3.9e-12 2.6e-08
sp0071520  -3.3           -3.3    2.7 5.8e-12 3.4e-08

Thank you in advance.

edger rna-seq • 1.1k views

ADD COMMENT • link updated 7.3 years ago by Aaron Lun ★ 28k • written 7.3 years ago by mictadlo ▴ 10

Gordon Smyth · Accepted Answer · 2017-09-18

3

Entering edit mode

Aaron Lun ★ 28k

@alun

Last seen 9 hours ago

The city by the bay

Try using the write.csv function.

library(edgeR)
example(glmTreat)
res <- topTags(tr, n=Inf)
write.csv(file="output.csv", res$table)

ADD COMMENT • link updated 7.3 years ago by Gordon Smyth 52k • written 7.3 years ago by Aaron Lun ★ 28k

0

Entering edit mode

By any chance, do you know how it would be possible to create an R named vector where the names are genes and the values are p-values from `tr` without save `tr` to a file?

ADD REPLY • link 7.3 years ago mictadlo ▴ 10

0

Entering edit mode

Use names<- or setNames on res$table$PValue.

ADD REPLY • link 7.3 years ago Aaron Lun ★ 28k