The log fold change from DE_timepoints is found from topTable in limma, but the p-values and adj pvalues are pulled from the results of the eBayes function (the fit object). However, topTable resorts the results that they provide, so the log-fold change are not getting matched up with the right genes. We need to reorder the lfc from topTable so that they match the right genes from fit.
The log fold change from DE_timepoints is found from
topTable
in limma, but the p-values and adj pvalues are pulled from the results of theeBayes
function (thefit
object). However,topTable
resorts the results that they provide, so the log-fold change are not getting matched up with the right genes. We need to reorder the lfc fromtopTable
so that they match the right genes fromfit
.