-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathnothingofnote.Rmd
140 lines (110 loc) · 3.51 KB
/
nothingofnote.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
---
title: "Yuh"
author: "Nick Gembs"
date: "5/11/2021"
output: html_document
---
```{r}
HistObs <- read.csv("C:/Users/Nick/Downloads/48453-annual-hist_obs-tmin.csv")
HistMod <- read.csv("C:/Users/Nick/Downloads/48453-annual-hist_mod-tmin.csv")
ProjMod <- read.csv("C:/Users/Nick/Downloads/48453-annual-proj_mod-tmin.csv", header=T)
```
```{r}
mean1 <- mean(HistObs$tmin)
mean1
mean2 <- mean(HistMod$weighted_mean)
mean2
mean3 <- mean(HistObs$tmin[33:64]) - mean(HistObs$tmin[1:33])
mean3
mean4 <- mean(HistMod$weighted_mean[33:57]) - mean(HistMod$weighted_mean[1:33])
mean4
mean5 <- mean(ProjMod$rcp45_weighted_mean[15:54]) - mean2
mean5
mean6 <- mean(ProjMod$rcp85_weighted.mean[15:54]) - mean2
mean6
mean7 <- mean(ProjMod$rcp45_weighted_mean[55:94]) - mean2
mean7
mean8 <- mean(ProjMod$rcp85_weighted.mean[55:94]) - mean2
mean8
mean9 <- (mean(ProjMod$rcp85_max[55:94]) - mean2)- (mean(ProjMod$rcp85_min[55:94]) - mean2)
mean9
mean10 <- (mean(ProjMod$rcp45_max[55:94]))- (mean(ProjMod$rcp45_min[55:94]))
mean10
mean9-mean10
mean(ProjMod$rcp85_max)
mean(ProjMod$rcp85_min)
```
```{r}
HistObs <- read.csv("C:/Users/Nick/Downloads/48453-annual-hist_obs-days_tmax_gt_90f (1).csv")
HistMod <- read.csv("C:/Users/Nick/Downloads/48453-annual-hist_mod-days_tmax_gt_90f.csv")
ProjMod <- read.csv("C:/Users/Nick/Downloads/48453-annual-proj_mod-days_tmax_gt_90f.csv")
```
```{r}
mean1 <- mean(HistObs$days_tmax_gt_90f)
mean1
mean2 <- mean(HistMod$weighted_mean)
mean2
mean3 <- mean(HistObs$days_tmax_gt_90f[33:64]) - mean(HistObs$days_tmax_gt_90f[1:33])
mean3
mean4 <- mean(HistMod$weighted_mean[33:57]) - mean(HistMod$weighted_mean[1:33])
mean4
mean5 <- mean(ProjMod$rcp45_weighted_mean[15:54]) - mean2
mean5
mean6 <- mean(ProjMod$rcp85_weighted.mean[15:54]) - mean2
mean6
mean7 <- mean(ProjMod$rcp45_weighted_mean[55:94]) - mean2
mean7
mean8 <- mean(ProjMod$rcp85_weighted.mean[55:94]) - mean2
mean8
```
```{r}
HistObs <- read.csv("C:/Users/Nick/Downloads/48453-annual-hist_obs-pcpn.csv")
HistMod <- read.csv("C:/Users/Nick/Downloads/48453-annual-hist_mod-pcpn.csv")
ProjMod <- read.csv("C:/Users/Nick/Downloads/48453-annual-proj_mod-pcpn.csv")
```
```{r}
mean1 <- mean(HistObs$pcpn)
mean1
mean2 <- mean(HistMod$weighted_mean)
mean2
mean3 <- mean(HistObs$pcpn[33:64]) - mean(HistObs$pcpn[1:33])
mean3
mean4 <- mean(HistMod$weighted_mean[33:57]) - mean(HistMod$weighted_mean[1:33])
mean4
mean5 <- mean(ProjMod$rcp45_weighted_mean[15:54]) - mean2
mean5
mean6 <- mean(ProjMod$rcp85_weighted.mean[15:54]) - mean2
mean6
mean7 <- mean(ProjMod$rcp45_weighted_mean[55:94]) - mean2
mean7
mean8 <- mean(ProjMod$rcp85_weighted.mean[55:94]) - mean2
mean8
mean9 <- mean(ProjMod$rcp85_max) - mean(ProjMod$rcp85_min)
mean9
```
```{r}
HistObs <- read.csv("C:/Users/Nick/Downloads/48453-annual-hist_obs-days_pcpn_gt_1in.csv")
HistMod <- read.csv("C:/Users/Nick/Downloads/48453-annual-hist_mod-days_pcpn_gt_1in.csv")
ProjMod <- read.csv("C:/Users/Nick/Downloads/48453-annual-proj_mod-days_pcpn_gt_1in.csv")
```
```{r}
mean1 <- mean(HistObs$days_pcpn_gt_1in)
mean1
mean2 <- mean(HistMod$weighted_mean)
mean2
mean3 <- mean(HistObs$days_pcpn_gt_1in[33:64]) - mean(HistObs$days_pcpn_gt_1in[1:33])
mean3
mean4 <- mean(HistMod$weighted_mean[33:57]) - mean(HistMod$weighted_mean[1:33])
mean4
mean5 <- mean(ProjMod$rcp45_weighted_mean[15:54]) - mean2
mean5
mean6 <- mean(ProjMod$rcp85_weighted.mean[15:54]) - mean2
mean6
mean7 <- mean(ProjMod$rcp45_weighted_mean[55:94]) - mean2
mean7
mean8 <- mean(ProjMod$rcp85_weighted.mean[55:94]) - mean2
mean8
```
```{r}
pnorm(24,25,1, lower.tail=T)
```