-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathFig4_Sginificance_Testing.R
93 lines (73 loc) · 4.26 KB
/
Fig4_Sginificance_Testing.R
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
#install.packages("exactRankTests")
library(exactRankTests)
#Load in Screened and Unscreened Esoil data
#Mosaic
Mosaic_og=read.table('/Volumes/REESEN/SMAP/Validation_Data/NLDAS2/Mosaic/Esoil_ET_ratio/Mosaic_esoil_noScreening.csv',header=FALSE,sep=",");
Mosaic_og=as.matrix(Mosaic_og)
Mosaic_og=as.vector(Mosaic_og)
Screened_Mosaic=read.table('/Volumes/REESEN/SMAP/Validation_Data/NLDAS2/Mosaic/Esoil_ET_ratio/Mosaic_esoil_Screened.csv',header=FALSE,sep=",");
Screened_Mosaic=as.matrix(Screened_Mosaic)
Screened_Mosaic=as.vector(Screened_Mosaic)
Screened_Mosaic<-Screened_Mosaic[!is.na(Mosaic_og)]
Mosaic_og<-Mosaic_og[!is.na(Mosaic_og)]
#determine if the mean is significantly different with wilcox rank sum test
wilcox.exact(x=Screened_Mosaic,y=Mosaic_og,alternative="greater",paired=TRUE)
#Noah
Noah_og=read.table('/Volumes/REESEN/SMAP/Validation_Data/NLDAS2/Noah/Esoil_ET_ratio/Noah_esoil_noScreening.csv',header=FALSE,sep=",");
Noah_og=as.matrix(Noah_og)
Noah_og=as.vector(Noah_og)
Screened_Noah=read.table('/Volumes/REESEN/SMAP/Validation_Data/NLDAS2/Noah/Esoil_ET_ratio/Noah_esoil_Screened.csv',header=FALSE,sep=",");
Screened_Noah=as.matrix(Screened_Noah)
Screened_Noah=as.vector(Screened_Noah)
Screened_Noah<-Screened_Noah[!is.na(Noah_og)]
Noah_og<-Noah_og[!is.na(Noah_og)]
#determine if the mean is significantly different with wilcox rank sum test
wilcox.exact(x=Screened_Noah,y=Noah_og,alternative="greater",paired=TRUE)
#GLEAM
GLEAM_og=read.table('/Volumes/REESEN/SMAP/Validation_Data/GLEAM_9km/Esoil_ET_ratio/GLEAM_esoil_noScreening.csv',header=FALSE,sep=",");
GLEAM_og=as.matrix(GLEAM_og)
GLEAM_og=as.vector(GLEAM_og)
Screened_GLEAM=read.table('/Volumes/REESEN/SMAP/Validation_Data/GLEAM_9km/Esoil_ET_ratio/GLEAM_esoil_Screened.csv',header=FALSE,sep=",");
Screened_GLEAM=as.matrix(Screened_GLEAM)
Screened_GLEAM=as.vector(Screened_GLEAM)
Screened_GLEAM<-Screened_GLEAM[!is.na(GLEAM_og)]
GLEAM_og<-GLEAM_og[!is.na(GLEAM_og)]
#determine if the mean is significantly different with wilcox rank sum test
wilcox.exact(x=Screened_GLEAM,y=GLEAM_og,alternative="greater",paired=TRUE)
##===============================================================================================
##Now test for Esoil/ET ratio
##===============================================================================================
#Load in Screened and Unscreened Esoil data
#Mosaic
Mosaic_og=read.table('/Volumes/REESEN/SMAP/Validation_Data/NLDAS2/Mosaic/Esoil_ET_ratio/Mosaic_ratio_noScreening.csv',header=FALSE,sep=",");
Mosaic_og=as.matrix(Mosaic_og)
Mosaic_og=as.vector(Mosaic_og)
Screened_Mosaic=read.table('/Volumes/REESEN/SMAP/Validation_Data/NLDAS2/Mosaic/Esoil_ET_ratio/Mosaic_ratio_Screened.csv',header=FALSE,sep=",");
Screened_Mosaic=as.matrix(Screened_Mosaic)
Screened_Mosaic=as.vector(Screened_Mosaic)
Screened_Mosaic<-Screened_Mosaic[!is.na(Mosaic_og)]
Mosaic_og<-Mosaic_og[!is.na(Mosaic_og)]
#determine if the mean is significantly different with wilcox rank sum test
wilcox.exact(x=Screened_Mosaic,y=Mosaic_og,alternative="greater",paired=TRUE)
#Noah
Noah_og=read.table('/Volumes/REESEN/SMAP/Validation_Data/NLDAS2/Noah/Esoil_ET_ratio/Noah_ratio_noScreening.csv',header=FALSE,sep=",");
Noah_og=as.matrix(Noah_og)
Noah_og=as.vector(Noah_og)
Screened_Noah=read.table('/Volumes/REESEN/SMAP/Validation_Data/NLDAS2/Noah/Esoil_ET_ratio/Noah_ratio_Screened.csv',header=FALSE,sep=",");
Screened_Noah=as.matrix(Screened_Noah)
Screened_Noah=as.vector(Screened_Noah)
Screened_Noah<-Screened_Noah[!is.na(Noah_og)]
Noah_og<-Noah_og[!is.na(Noah_og)]
#determine if the mean is significantly different with wilcox rank sum test
wilcox.exact(x=Screened_Noah,y=Noah_og,alternative="greater",paired=TRUE)
#GLEAM
GLEAM_og=read.table('/Volumes/REESEN/SMAP/Validation_Data/GLEAM_9km/Esoil_ET_ratio/GLEAM_ratio_noScreening.csv',header=FALSE,sep=",");
GLEAM_og=as.matrix(GLEAM_og)
GLEAM_og=as.vector(GLEAM_og)
Screened_GLEAM=read.table('/Volumes/REESEN/SMAP/Validation_Data/GLEAM_9km/Esoil_ET_ratio/GLEAM_ratio_Screened.csv',header=FALSE,sep=",");
Screened_GLEAM=as.matrix(Screened_GLEAM)
Screened_GLEAM=as.vector(Screened_GLEAM)
Screened_GLEAM<-Screened_GLEAM[!is.na(GLEAM_og)]
GLEAM_og<-GLEAM_og[!is.na(GLEAM_og)]
#determine if the mean is significantly different with wilcox rank sum test
wilcox.exact(x=Screened_GLEAM,y=GLEAM_og,alternative="greater",paired=TRUE)