Linear regression model implementation to the expression levels of mRNA-miRNA

Question

[alptaciroglu]

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.

[alptaciroglu]

"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

[Sean Davis]

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

[alptaciroglu]

4 samples for mRNA 2 for control and 2 for treatment

2 samples for miRNA 1 for control and 1 for treatment