LINEAR-REGRESSION.R

#load data
data=read.csv("ADMIX.csv", header = TRUE)

#analyze
admix.time=lm(ADMIXTURE~INITIAL, data=data)
summary(admix.time)

## 
## Call:
## lm(formula = ADMIXTURE ~ INITIAL, data = data)
## 
## Residuals:
##          1          2          3          4          5 
##  1.689e-02 -2.133e-02 -8.674e-18 -8.444e-03  1.289e-02 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -0.6916533  0.0322271  -21.46 0.000221 ***
## INITIAL      0.0059556  0.0001702   34.99 5.13e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01805 on 3 degrees of freedom
## Multiple R-squared:  0.9976, Adjusted R-squared:  0.9967 
## F-statistic:  1224 on 1 and 3 DF,  p-value: 5.133e-05