I want to do
RNA seq analysis for multiple files together without any replication and control. I do not want to perform a comparison with control, I want to basically run a linear regression. I have input files (RNA paired end) for 50 varieties. There are 15 chemical compounds concentration each for all 50 varieties. I am using
salmon to create
TPM values and then I want to test for regression relationship between independent variable (each chemical compound which is a vector of all 50 varieties) and dependent variable as
TPM values from all the varieties. I have generated the
TPM values in
quant.sf files using
salmon and I need some help about the steps to move forward to perform linear regression analysis in
R as this is my first time doing it.
Thank you in advance for any suggestion/comment/guidance!