This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Linear regression model implementation to the expression levels of mRNA-miRNA

Dear All,

I need to implement a regression model for the expression relation of miRNAs - mRNAs. In the article http://www.sciencedirect.com/science/article/pii/S0888754314001815 the regression model of y=a0+a1x1+a2x2+...+anxn is provided; y being expression level of mRNAs, x being expression level of miRNAs, the parameters a0, a1 ... being the regression constant and coefficient. I need to apply this formula to my data in R however, I don't even understand how to start doing this.

In my data, I have log fold change of mRNAs and log fold change of miRNAs which should be used as an input to linear regression model. Any help is greatly appreciated.

mrna mirna regression expression
"MiRBase.ID"    "Gene.Symbol"    "miR.Family"    "Gene.ID"    "Transcript.ID"    "PCT"    "affy_hg_u133_plus_2"    "logFCmRNAs"    "adj.P.ValmRNAs"    "logFCmiRNAs"
"1"    "hsa-miR-106a-5p"    "DYRK2"    "miR-17-5p/20-5p/93-5p/106-5p/519-3p/526-3p"    "ENSG00000127334.10"    "ENST00000344096.3"    "0.25"    "202970_at"    0.74088249    0.002534602    -0.2
"2"    "hsa-miR-106a-5p"    "DYRK2"    "miR-17-5p/20-5p/93-5p/106-5p/519-3p/526-3p"    "ENSG00000127334.10"    "ENST00000344096.3"    "0.25"    "202971_s_at"    0.539734916    0.002805065    -0.2
"3"    "hsa-miR-106a-5p"    "PIK3R1"    "miR-17-5p/20-5p/93-5p/106-5p/519-3p/526-3p"    "ENSG00000145675.10"    "ENST00000521381.1"    "0.65"    "212240_s_at"    0.479017416    0.006262084    -0.2
"4"    "hsa-miR-106a-5p"    "GID4"    "miR-17-5p/20-5p/93-5p/106-5p/519-3p/526-3p"    "ENSG00000141034.5"    "ENST00000268719.4"    "0.74"    "228452_at"    1.040642812    1.94e-05    -0.2

On how many samples do you have the mRNA and miRNA measurements?

4 samples for mRNA 2 for control and 2 for treatment

2 samples for miRNA 1 for control and 1 for treatment

0 answers

No answers yet.

Log in to answer this question.