-
Notifications
You must be signed in to change notification settings - Fork 0
/
POT with GPD (JPFA).R
83 lines (76 loc) · 1.42 KB
/
POT with GPD (JPFA).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
library(readxl)
POT=read_excel("POT.xlsx")
# untuk l_1
l_1 <- function(df) {
n <- length(df)
df <- sort(df)
hasil <- 0
for (i in 1:n) {
hasil <- hasil + df[i]/n
}
return(hasil)
}
# Untuk l_2
l_2 <- function(df) {
n <- length(df)
df <- sort(df) #Mengurutukan df
hasil <- 0
for (i in 1:(n-1)) {
for (j in (i+1):n) {
hasil <- hasil + (df[j] - df[i])
}
}
hasil <- hasil / (2 * choose(n, 2))
return(hasil)
}
# Untuk l_3
l_3 <- function(df) {
n <- length(df)
df <- sort(df) # Mengurutkan df
hasil <- 0
for (i in 1:(n-2)) {
for (j in (i+1):(n-1)) {
for (k in (j+1):n) {
hasil <- hasil + (df[k] - 2 * df[j] + df[i])
}
}
}
hasil <- hasil / (3 * choose(n, 3))
return(hasil)
}
#untuk l_4
# Fungsi untuk menghitung persamaan l_4
l_4 <- function(df) {
n <- length(df)
df <- sort(df) # Mengurutkan df
hasil <- 0
for (i in 1:(n-3)) {
for (j in (i+1):(n-2)) {
for (k in (j+1):(n-1)) {
for (l in (k+1):n) {
hasil <- hasil + (df[l] - 3 * df[k] + 3 * df[j] - df[i])
}
}
}
}
hasil <- hasil / (4 * choose(n, 4))
return(hasil)
}
##2. Untuk df JPFA
df=POT$JPFA
#Untuk l_1
l_1=l_1(df)
#Untuk l_2
l_2=l_2(df)
#Untuk l_3
l_3=l_3(df)
#Untuk l_4
l_4=l_4(df)
#tau_3
tau_3=l_3/l_2
#tau_4
tau_4=l_4/l_2
##GEV untuk df JPFA
LM3=data.frame(l_1=l_1,l_2=l_2,l_3=l_3,
l_4=l_4,tau_3=tau_3,tau_4=tau_4)
LM3