Skip to content

Commit fe958b4

Browse files
Add files via upload
1 parent a9b523f commit fe958b4

File tree

1 file changed

+18
-5
lines changed

1 file changed

+18
-5
lines changed

Sea-Level-Predictor/sea_level_predictor.py

Lines changed: 18 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -4,19 +4,32 @@
44

55
def draw_plot():
66
# Read data from file
7-
7+
df = pd.read_csv('epa-sea-level.csv')
8+
x = df['Year']
9+
y = df['CSIRO Adjusted Sea Level']
810

911
# Create scatter plot
10-
12+
plt.scatter(x, y)
1113

1214
# Create first line of best fit
13-
15+
line1 = linregress(x, y)
16+
x_pred1 = pd.Series([i for i in range(1880,2051)])
17+
y_pred1 = line1.slope * x_pred1 + line1.intercept
18+
plt.plot(x_pred1, y_pred1, 'r')
1419

1520
# Create second line of best fit
16-
21+
newdf = df.loc[df["Year"] >= 2000]
22+
x_new = newdf['Year']
23+
y_new = newdf['CSIRO Adjusted Sea Level']
24+
line2 = linregress(x_new, y_new)
25+
x_pred2 = pd.Series([i for i in range(2000,2051)])
26+
y_pred2 = line2.slope * x_pred2 + line2.intercept
27+
plt.plot(x_pred2, y_pred2, 'b')
1728

1829
# Add labels and title
19-
30+
plt.title('Rise in Sea Level')
31+
plt.xlabel('Year')
32+
plt.ylabel('Sea Level (inches)')
2033

2134
# Save plot and return data for testing (DO NOT MODIFY)
2235
plt.savefig('sea_level_plot.png')

0 commit comments

Comments
 (0)