From 9458dac842f77252744c65f6e82deb81bf94203e Mon Sep 17 00:00:00 2001
From: Mahsa DEG analysis
function in logarithm scale
{r} fit <- lmFit(Ex, design = design) contrast <- makeContrasts(non_V600E-v600e, levels = design) fit2 <- contrasts.fit(fit = fit, contrasts = contrast) fit2 <- eBayes(fit2) deg <- topTable(fit2, number = Inf)
logFC AveExpr t P.Value adj.P.Val B
+ENSG00000114302 2.059983 9.016187 9.425437 5.230565e-18 7.353333e-14 30.26257
+ENSG00000135709 -1.866375 2.860313 -9.361123 8.118053e-18 7.353333e-14 29.83572
+ENSG00000106733 -1.558516 2.832475 -9.152146 3.356666e-17 1.843978e-13 28.45755
+ENSG00000157954 -2.057048 2.989810 -9.123564 4.071489e-17 1.843978e-13 28.27013
+ENSG00000164944 1.905614 7.507759 9.077754 5.544899e-17 2.009028e-13 27.97029
+ENSG00000187098 -1.353065 2.716279 -9.004784 9.056078e-17 2.533801e-13 27.49412
diff --git a/README.md b/README.md
index 6471062..586a274 100644
--- a/README.md
+++ b/README.md
@@ -133,10 +133,18 @@ fit2 <- eBayes(fit2)
deg <- topTable(fit2, number = Inf)
```
+```
+ logFC AveExpr t P.Value adj.P.Val B
+ENSG00000114302 2.059983 9.016187 9.425437 5.230565e-18 7.353333e-14 30.26257
+ENSG00000135709 -1.866375 2.860313 -9.361123 8.118053e-18 7.353333e-14 29.83572
+ENSG00000106733 -1.558516 2.832475 -9.152146 3.356666e-17 1.843978e-13 28.45755
+ENSG00000157954 -2.057048 2.989810 -9.123564 4.071489e-17 1.843978e-13 28.27013
+ENSG00000164944 1.905614 7.507759 9.077754 5.544899e-17 2.009028e-13 27.97029
+ENSG00000187098 -1.353065 2.716279 -9.004784 9.056078e-17 2.533801e-13 27.49412
+```
-![](https://github.com/Mahsa-Ehsanifard/GEO_microarray/blob/main/voom.png?raw=true)