Problem with differential expression analysis about sns HOT 14 CLOSED

If there is missing output, then something went wrong. There may not be a specific error. Have you checked the full sbatch log file?

from sns.

The sbatch log file has 13 slurm.out files in comparison to my original 12 fastq files. The other log folders look intact as well.

from sns.

That makes sense. There should be an rna-star log for every sample and one for rna-star-groups-dge.

from sns.

Yes, I think the main issue is that I am not seeing the deseq2.dds.rds file as I want to make graphs with my data.

from sns.

If it's missing, something went wrong. Check the full sbatch log file, especially where it ended.

from sns.

Unless you are talking about the logs-sbatch file, I do not have a full sbatch log file. However, the logs files that were outputted (logs-trimmomatic, logs-sbatch, logs-featureCounts, logs-align-star), all seem to have checked out. Was there supposed to be another logs file folder that somehow did not go through?

from sns.

In the logs-sbatch directory, there is a log file for every sbatch job. The last one should be for the rna-star-groups-dge job.

from sns.

Hey Igor,
upon looking in the logs-sbatch directory, there seems to be a rna-star-groups-dge log file, as the date modified is about the same time I performed the rna-star-groups-dge command. In addition, there were 13 log files in the logs-sbatch directory (1 rna-star log for each of the 12 samples, and 1 rna-star-groups-dge log file).

from sns.

Hey Igor,
here attached I have a screenshot of my folder. Maybe that will help?

from sns.

I can see you have 13 sbatch log files (as previously mentioned). The most recent one was likely generated by the rna-star-groups-dge job. You can view it in the terminal with cat logs-sbatch/slurm-23702903.out. You can also download it and open it with any text editor. If you are not sure what to look for, you can copy the contents here (or attach the whole file).

from sns.

I am looking at the file right now and there does seem to be an error. Here is the message:

Error in checkForExperimentalReplicates(object, modelMatrix) :

The design matrix has the same number of samples and coefficients to fit,
so estimation of dispersion is not possible. Treating samples
as replicates was deprecated in v1.20 and no longer supported since v1.22.

Calls: DESeq ... estimateDispersions -> .local -> checkForExperimentalReplicates
Execution halted

from sns.

According to the message, you do not have replicates. You cannot perform differential expression without replicates.

from sns.

Hey Igor,
I realized that my group names were different for each sample, resulting in no replicates. After altering the names, the rdata folder was populated and my problem was fixed. Thank you so much for your help!

from sns.

Happy to hear everything is okay now.

from sns.