I have 4 .csv files containing gene level (transcripts collapsed) counts coming from infected (2 biological replicates) and non-infected (2 biological replicates) generated via RSEM. I dont have access to raw data files. Their contents are as follows;
gene_id transcript_id(s) length expected_count FPKM AT1G01010 AT1G01010.1 1688 88 2.68
However, number of rows in each file is different (My guess is 0 counts were removed). Can someone please help me by advising (or share a R script)me on how to get these data combined to match the input of tximport? so I can do the differential expression analysis using DESeq2. Thank you!