library(tidyverse)
library(modelr)
library(readxl)
Agronomic_Yields <- read_excel("~/Agronomic Yields.xlsx")
View(Agronomic_Yields)
Agronomic_Yields_Model <- lm(formula = fertilizer_rate_kg_ha ~ yield_kg_ha, data = Agronomic_Yields)
summary(Agronomic_Yields_Model)

Call:
lm(formula = fertilizer_rate_kg_ha ~ yield_kg_ha, data = Agronomic_Yields)

Residuals:
    Min      1Q  Median      3Q     Max 
-170.78  -64.38  -17.98   53.87  238.59 

Coefficients:
             Estimate Std. Error t value Pr(>|t|)    
(Intercept) 1.462e+01  5.203e+00    2.81  0.00505 ** 
yield_kg_ha 1.410e-02  6.975e-04   20.21  < 2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 78.2 on 1035 degrees of freedom
  (159 observations deleted due to missingness)
Multiple R-squared:  0.2829,    Adjusted R-squared:  0.2823 
F-statistic: 408.4 on 1 and 1035 DF,  p-value: < 2.2e-16
LS0tDQp0aXRsZTogIkFncm9ub21pYyBHcm93dGgiDQpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sNCi0tLQ0KDQpgYGB7cn0NCmxpYnJhcnkodGlkeXZlcnNlKQ0KbGlicmFyeShtb2RlbHIpDQpgYGANCg0KYGBge3J9DQpsaWJyYXJ5KHJlYWR4bCkNCkFncm9ub21pY19ZaWVsZHMgPC0gcmVhZF9leGNlbCgifi9BZ3Jvbm9taWMgWWllbGRzLnhsc3giKQ0KVmlldyhBZ3Jvbm9taWNfWWllbGRzKQ0KYGBgDQpgYGB7cn0NCkFncm9ub21pY19ZaWVsZHNfTW9kZWwgPC0gbG0oZm9ybXVsYSA9IGZlcnRpbGl6ZXJfcmF0ZV9rZ19oYSB+IHlpZWxkX2tnX2hhLCBkYXRhID0gQWdyb25vbWljX1lpZWxkcykNCmBgYA0KDQpgYGB7cn0NCnN1bW1hcnkoQWdyb25vbWljX1lpZWxkc19Nb2RlbCkNCmBgYA0KDQo=