Tahapan Pengestimasian Direct dan Traditional Indirect Estimation:
library(readxl)
Jatim_pengeluaran <- read_excel("D:/Kuliah/tkt IV/The Journey Start/SAE/Data SAE/Jatim_pengeluaran.xlsx",
col_types = c("text", "text", "numeric",
"text", "text", "numeric", "numeric",
"numeric", "numeric"))
Xnew_ya <- read_excel("D:/Kuliah/tkt IV/The Journey Start/SAE/Data SAE/Xnew ya.xlsx",
sheet = "Sheet2", col_types = c("text", "text", "text", "text", "text", "text",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric"))
#variabel_X_signifikan
x7 <- aggregate(Xnew_ya$X7, list(Xnew_ya$R102),sum,na.rm=TRUE)
x13 <- aggregate(Xnew_ya$X13, list(Xnew_ya$R102),sum,na.rm=TRUE)
x14 <- aggregate(Xnew_ya$X14, list(Xnew_ya$R102),sum,na.rm=TRUE)
x18 <- aggregate(Xnew_ya$X18, list(Xnew_ya$R102),sum,na.rm=TRUE)
x19 <- aggregate(Xnew_ya$X19, list(Xnew_ya$R102),sum,na.rm=TRUE)
x20 <- aggregate(Xnew_ya$X20, list(Xnew_ya$R102),sum,na.rm=TRUE)
x21 <- aggregate(Xnew_ya$X21, list(Xnew_ya$R102),sum,na.rm=TRUE)
x24 <- aggregate(Xnew_ya$X24, list(Xnew_ya$R102),sum,na.rm=TRUE)
x25 <- aggregate(Xnew_ya$X25, list(Xnew_ya$R102),sum,na.rm=TRUE)
x34 <- aggregate(Xnew_ya$X34, list(Xnew_ya$R102),sum,na.rm=TRUE)
dt2 <- data.frame(cbind("bs_o"=Jatim_pengeluaran$`No Kode Sampel`,"pop"=Jatim_pengeluaran$Exp_krb,"kab"=Jatim_pengeluaran$R102))
mat_x <- matrix(cbind(1,x7$x,x13$x,x14$x,x18$x,x19$x,x20$x,x21$x,x24$x,x25$x,x34$x),ncol = 11)
# Fungsi Menghitung Jumlah desa
jml_desa <- function(data){
r102 = unique(Jatim_pengeluaran$R102)
de = data.frame(r102)
Peng_jatim <- data.frame(Jatim_pengeluaran)
ds <- list()
for (i in r102){
ds[[i]] <- subset(Peng_jatim, R102 %in% i)
}
bsz <- list()
for(j in r102){
bsz[[j]] <- unique(ds[[j]]$No.Kode.tiga)
}
Ni <- c()
for(k in r102){
Ni[k] <- length(bsz[[k]])
}
Ni <- Ni[-c(30:70)]
N_desa <- data.frame("desa"=Ni)
return(N_desa)
}
Ni_desa <- jml_desa(Jatim_pengeluaran)Jumlah desa setiap kabupaten/kota:
Ni_desa## desa
## 1 34
## 2 43
## 3 42
## 4 51
## 5 52
## 6 64
## 7 80
## 8 50
## 9 87
## 10 71
## 11 45
## 12 44
## 13 55
## 14 63
## 15 72
## 16 50
## 17 52
## 18 47
## 19 41
## 20 32
## 21 46
## 22 54
## 23 53
## 24 53
## 25 52
## 26 42
## 27 43
## 28 40
## 29 52Fungsi Estimasi untuk metode Two stage cluster SRSWOR-SRSWOR
resampling_2st_2 = function(dataY, dataX, f1, f2, iterasi ){
est_2st_cl <- function(dataY, f1,f2){
sp_2_ct <- function(dataY,f1,f2){
samp_th1 <- function(mydata,f1){ #pengambilan sampel tahap 1
bs_u <-unique(mydata$bs_o)
bs_u <- data.frame(bs_u)
samp_bs <- sample(1:nrow(bs_u), size = f1*nrow(bs_u), replace = FALSE)
sampi <- bs_u[samp_bs,]
sampi2 <- as.data.frame(paste(sampi))
colnames(sampi2) <- c("sampel_bs")
return(sampi2)
}
hs_sp_t1 <- samp_th1(dataY,f1)
samp_th2 <- function(dt_kb_1,hs_sp_1,f2){
index <- c(1:nrow(hs_sp_t1))
hs_sp_1 <- cbind(index,hs_sp_t1)
dt_2 <- list();dt_3 <- list() ; samp_rt <- list()
hs_samp_srs <- list() ; hs_smp_srs <- list()
for(i in hs_sp_1$index){ #Pengambilan sampel ssu scr srs wor
dt_2[[i]] <- subset(hs_sp_1, index %in% i)
dt_3[[i]] <- subset(dt_kb_1, bs_o %in% dt_2[[i]]$sampel_bs)
samp_rt[[i]] <- sample(1:length(dt_3[[i]]$pop), size = f2*length(dt_3[[i]]$pop), replace = FALSE)
hs_samp_srs[[i]] <- dt_3[[i]][samp_rt[[i]],]
}
gb_samp = data.frame()
for (j in hs_sp_1$index) {
gb_samp <- rbind(gb_samp,hs_samp_srs[[j]])
}
return(gb_samp)
}
res_sp <- samp_th2(dataY,hs_sp_t1,f2)
res_sp$pop <- as.numeric(paste(res_sp$pop))
res_sp$bs_o <- as.character(paste(res_sp$bs_o))
dataY$bs_o <- as.character(paste(dataY$bs_o))
ybar_i. <- aggregate(res_sp$pop, list(res_sp$bs_o), mean)
m_i <- data.frame(table(res_sp$bs_o))
mi <- m_i$Freq
Mi <- round(mi/f2)
n <- length(unique(res_sp$bs_o))
N <- length(unique(dataY$bs_o))
Y_cap_est <- (N/n) * (sum(Mi * ybar_i.$x))
sw_i_2 <- 1/(mi-1) * sum((res_sp$pop - rep(ybar_i.$x,length(res_sp$pop)))^2)
sb_2 <- 1/(n-1) * sum((Mi*ybar_i.$x) - (1/n *(sum(Mi*ybar_i.$x)))^2)
var_cap_est <- (N^2 * (1/n - 1/N) * sb_2 ) + (N/n * sum(Mi^2 * (1/mi - 1/Mi) * sw_i_2))
rel_var_y <- var_cap_est/Y_cap_est
hasil <- cbind(Y_cap_est,var_cap_est,rel_var_y)
return(hasil)
}
r102 = unique(dt2$kab)
datakoeh = list()
for(b in r102){
datakoeh[[b]] <- subset(dt2, kab %in% b)
}
try_ye <- list()
for(c in r102){
try_ye[[c]] <- sp_2_ct(data = datakoeh[[c]],f1,f2)
}
this_it = data.frame()
for (q in r102) {
this_it <- rbind(this_it, try_ye[[q]])
}
this_it = cbind(r102, this_it)
}
#fungsi estimasi synthetic dan composite
Y_synth_comp <- function(x,y,vary){
#synthetic function
xty <- t(x)%*%y
xtx <- t(x)%*%x
bt_xy <- solve(xtx)%*%xty
bt_xy = t(bt_xy)
Y_syn <- c(rep(0,nrow(x)))
for(i in 1: ncol(x)){
Y_syn = Y_syn + x[,i]*bt_xy[,i]
}
Y_syn <- as.data.frame(Y_syn)
mse_syn <- abs(((Y_syn-y)^2)-vary)
rel_mse_syn <- mse_syn/Y_syn
#composite function
tetha_i <- mse_syn/(vary+mse_syn)
tetha_i[tetha_i<0]<-0
tetha_i[tetha_i>1]<-1
Y_comp <- (tetha_i*y) + ((1-tetha_i)*Y_syn)
mse_comp <- (1-tetha_i)*mse_syn
rel_mse_comp <- mse_comp/Y_comp
Syn_Comp <- data.frame(cbind(Y_syn, Y_comp,mse_syn, mse_comp, rel_mse_syn, rel_mse_comp))
colnames(Syn_Comp) = c("Y syn","Y comp", "MSE syn", "MSE comp", "Relatif MSE syn","Relatif MSE comp")
return(Syn_Comp)
}
dataset = list()
est_dataset = list()
for( j in 1:iterasi){
set.seed(j)
dataset[[j]] = est_2st_cl(dataY, f1, f2)
est_dataset [[j]] = Y_synth_comp(dataX,dataset[[j]]$Y_cap_est,dataset[[j]]$var_cap_est)
}
databind = data.frame()
databind_ind = data.frame()
for (b in 1:iterasi) {
databind = rbind(databind, dataset[[b]])
databind_ind = rbind(databind_ind, est_dataset[[b]])
}
datafinal = aggregate(list(databind$Y_cap_est,databind_ind$`Y syn`,databind_ind$`Y comp`,
databind$var_cap_est,databind_ind$`MSE syn`,databind_ind$`MSE comp`
,databind$rel_var_y,databind_ind$`Relatif MSE syn`,databind_ind$`Relatif MSE comp`),
by=list(databind$r102),
FUN=mean)
colnames(datafinal) = c("kab", "Y dir", "Y syn","Y com", "MSE dir","MSE syn","MSE com",
"Relatif MSE dir","Relatif MSE syn","Relatif MSE com")
return(datafinal)
}b. MEnggunakan rumus MSE synthetic yang perbaikanresampling2_2st_2 = function(dataY, dataX, Ni, f1, f2, iterasi ){
est_2st_cl <- function(dataY, f1,f2){
sp_2_ct <- function(dataY,f1,f2){
samp_th1 <- function(mydata,f1){ #pengambilan sampel tahap 1
bs_u <-unique(mydata$bs_o)
bs_u <- data.frame(bs_u)
samp_bs <- sample(1:nrow(bs_u), size = f1*nrow(bs_u), replace = FALSE)
sampi <- bs_u[samp_bs,]
sampi2 <- as.data.frame(paste(sampi))
colnames(sampi2) <- c("sampel_bs")
return(sampi2)
}
hs_sp_t1 <- samp_th1(dataY,f1)
samp_th2 <- function(dt_kb_1,hs_sp_1,f2){
index <- c(1:nrow(hs_sp_t1))
hs_sp_1 <- cbind(index,hs_sp_t1)
dt_2 <- list();dt_3 <- list() ; samp_rt <- list()
hs_samp_srs <- list() ; hs_smp_srs <- list()
for(i in hs_sp_1$index){ #Pengambilan sampel ssu scr srs wor
dt_2[[i]] <- subset(hs_sp_1, index %in% i)
dt_3[[i]] <- subset(dt_kb_1, bs_o %in% dt_2[[i]]$sampel_bs)
samp_rt[[i]] <- sample(1:length(dt_3[[i]]$pop), size = f2*length(dt_3[[i]]$pop), replace = FALSE)
hs_samp_srs[[i]] <- dt_3[[i]][samp_rt[[i]],]
}
gb_samp = data.frame()
for (j in hs_sp_1$index) {
gb_samp <- rbind(gb_samp,hs_samp_srs[[j]])
}
return(gb_samp)
}
res_sp <- samp_th2(dataY,hs_sp_t1,f2)
res_sp$pop <- as.numeric(paste(res_sp$pop))
res_sp$bs_o <- as.character(paste(res_sp$bs_o))
dataY$bs_o <- as.character(paste(dataY$bs_o))
ybar_i. <- aggregate(res_sp$pop, list(res_sp$bs_o), mean)
m_i <- data.frame(table(res_sp$bs_o))
mi <- m_i$Freq
Mi <- round(mi/f2)
n <- length(unique(res_sp$bs_o))
N <- length(unique(dataY$bs_o))
Y_cap_est <- (N/n) * (sum(Mi * ybar_i.$x))
sw_i_2 <- 1/(mi-1) * sum((res_sp$pop - rep(ybar_i.$x,length(res_sp$pop)))^2)
sb_2 <- 1/(n-1) * sum((Mi*ybar_i.$x) - (1/n *(sum(Mi*ybar_i.$x)))^2)
var_cap_est <- (N^2 * (1/n - 1/N) * sb_2 ) + (N/n * sum(Mi^2 * (1/mi - 1/Mi) * sw_i_2))
rel_var_y <- var_cap_est/Y_cap_est
hasil <- cbind(Y_cap_est,var_cap_est,rel_var_y)
return(hasil)
}
r102 = unique(dt2$kab)
datakoeh = list()
for(b in r102){
datakoeh[[b]] <- subset(dt2, kab %in% b)
}
try_ye <- list()
for(c in r102){
try_ye[[c]] <- sp_2_ct(data = datakoeh[[c]],f1,f2)
}
this_it = data.frame()
for (q in r102) {
this_it <- rbind(this_it, try_ye[[q]])
}
this_it = cbind(r102, this_it)
}
#fungsi estimasi synthetic dan composite
Y_synth_comp <- function(x,y,vary,Ni){
#synthetic function
N <- Ni
xty <- t(x)%*%y
xtx <- t(x)%*%x
bt_xy <- solve(xtx)%*%xty
bt_xy = t(bt_xy)
Y_syn <- c(rep(0,nrow(x)))
for(i in 1: ncol(x)){
Y_syn = Y_syn + x[,i]*bt_xy[,i]
}
Y_syn <- as.data.frame(Y_syn)
var_y_cap_is <- matrix(var(Y_syn),nrow = nrow(N))
Ni_2 <- N^2/nrow(N) #Ni kuadrat*1/m
var_Ysyn_y <- matrix(var(Y_syn - y),nrow = nrow(N))
ht_2 <- 1/N^2*var_Ysyn_y
y_sel <- data.frame((Y_syn - y)^2)
ba_2 <- (1/N^2)*(y_sel)-ht_2
mse_syn <- var_y_cap_is + (Ni_2*ba_2)
rel_mse_syn <- mse_syn/Y_syn
#composite function
tetha_i <- mse_syn/(vary+mse_syn)
tetha_i[tetha_i<0]<-0
tetha_i[tetha_i>1]<-1
Y_comp <- (tetha_i*y) + ((1-tetha_i)*Y_syn)
mse_comp <- (1-tetha_i)*mse_syn
rel_mse_comp <- mse_comp/Y_comp
Syn_Comp <- data.frame(cbind(Y_syn, Y_comp,mse_syn, mse_comp, rel_mse_syn, rel_mse_comp))
colnames(Syn_Comp) = c("Y syn","Y comp", "MSE syn", "MSE comp", "Relatif MSE syn","Relatif MSE comp")
return(Syn_Comp)
}
dataset = list()
est_dataset = list()
for( j in 1:iterasi){
set.seed(j)
dataset[[j]] = est_2st_cl(dataY, f1, f2)
est_dataset [[j]] = Y_synth_comp(dataX,dataset[[j]]$Y_cap_est,dataset[[j]]$var_cap_est,Ni)
}
databind = data.frame()
databind_ind = data.frame()
for (b in 1:iterasi) {
databind = rbind(databind, dataset[[b]])
databind_ind = rbind(databind_ind, est_dataset[[b]])
}
datafinal = aggregate(list(databind$Y_cap_est,databind_ind$`Y syn`,databind_ind$`Y comp`,
databind$var_cap_est,databind_ind$`MSE syn`,databind_ind$`MSE comp`
,databind$rel_var_y,databind_ind$`Relatif MSE syn`,databind_ind$`Relatif MSE comp`),
by=list(databind$r102),
FUN=mean)
colnames(datafinal) = c("kab", "Y dir", "Y syn","Y com", "MSE dir","MSE syn","MSE com",
"Relatif MSE dir","Relatif MSE syn","Relatif MSE com")
return(datafinal)
}Fungsi Estimasi untuk metode Two stage cluster PPSWR-SRSWOR
#Pakai rumus awal
resampling_2st_1 = function(dataY, dataX, muatan, f1, f2, iterasi ){
est_2st_cl <- function(dataY, f1,f2){
sp_2_ct <- function(dataY,f1,f2){
kb_1 <- dataY
kb_1$bs_o <- as.character(paste(kb_1$bs_o))
Mi_kb_1 <- data.frame(table(kb_1$bs_o))
colnames(Mi_kb_1) <- c("bs_o","Muatan_x")
samp_th1 <- function(data,f1){ #pengambilan sampel tahap 1
pps_wr <- function(sizes,n){
N <- length(sizes)
cumsizes <- cumsum(sizes)
totsize <- cumsizes[N]
s <- numeric(n)
for (u in 1:n) {
r <- runif(1,0,totsize)
i <- 1
while (cumsizes[i] < r) {
i <- i+1
}
s[u] <- i
}
return(s)
}
n_length <- f1*nrow(data)
data_sa <- data[pps_wr(data$Muatan_x,n_length),]
return(data_sa)
}
hs_sp_t1 <- samp_th1(Mi_kb_1,f1)
samp_th2 <- function(dt_kb_1,hs_sp_1,f2){
index <- c(1:nrow(hs_sp_1))
hs_sp_1 <- cbind(index,hs_sp_1)
dt_2 <- list();dt_3 <- list() ; samp_rt <- list() ; kode <- list() ; mi <- list()
hs_samp_srs <- list() ; hs_smp_srs <- list(); n_samp <- list()
for(i in hs_sp_1$index){ #Pengambilan sampel ssu scr srs wor
dt_2[[i]] <- subset(hs_sp_1, index %in% i)
dt_3[[i]] <- subset(dt_kb_1, bs_o %in% dt_2[[i]]$bs_o)
samp_rt[[i]] <- sample(1:length(dt_3[[i]]$pop), size = f2*length(dt_3[[i]]$pop), replace = FALSE)
hs_samp_srs[[i]] <- dt_3[[i]][samp_rt[[i]],]
n_samp [[i]] <- nrow(hs_samp_srs[[i]])
kode[[i]] <- rep(i,n_samp[[i]])
mi[[i]] <- rep(n_samp[[i]],n_samp[[i]])
hs_smp_srs[[i]] <- cbind("kode"=kode[[i]],hs_samp_srs[[i]], "mi"=mi[[i]] )
}
gb_samp = data.frame()
for (j in hs_sp_1$index) {
gb_samp <- rbind(gb_samp,hs_smp_srs[[j]])
}
return(gb_samp)
}
res_sp <- samp_th2(kb_1,hs_sp_t1,f2) #kab 1
res_sp$pop <- as.numeric(paste(res_sp$pop))
res_sp$mi <- as.numeric(paste(res_sp$mi))
y_i. <- aggregate(res_sp$pop, list(res_sp$kode), mean)
M_nol <- sum(Mi_kb_1$Muatan_x)
Y_cap_est <- 1/ (length(y_i.$x)) * (sum(y_i.$x * M_nol))
pk_1 <- M_nol*y_i.$x
pk_2 <- pk_1 - Y_cap_est
pk_2_x <- pk_2^2
var_cap_est <- (sum(pk_2_x))/(((length(y_i.$x)))*((length(y_i.$x))-1))
rel_var_y <- var_cap_est/Y_cap_est
hasil <- cbind(Y_cap_est,var_cap_est,rel_var_y)
return(hasil)
}
r102 = unique(dt2$kab)
datakoeh = list()
for(b in r102){
datakoeh[[b]] <- subset(dt2, kab %in% b)
}
try_ye <- list()
for(c in r102){
try_ye[[c]] <- sp_2_ct(data = datakoeh[[c]],f1,f2)
}
this_it = data.frame()
for (q in r102) {
this_it <- rbind(this_it, try_ye[[q]])
}
this_it = cbind(r102, this_it)
}
#fungsi estimasi synthetic dan composite
Y_synth_comp <- function(x,y,vary){
#synthetic function
xty <- t(x)%*%y
xtx <- t(x)%*%x
bt_xy <- solve(xtx)%*%xty
bt_xy = t(bt_xy)
Y_syn <- c(rep(0,nrow(x)))
for(i in 1: ncol(x)){
Y_syn = Y_syn + x[,i]*bt_xy[,i]
}
Y_syn <- as.data.frame(Y_syn)
mse_syn <- abs(((Y_syn-y)^2)-vary)
rel_mse_syn <- mse_syn/Y_syn
#composite function
tetha_i <- mse_syn/(vary+mse_syn)
tetha_i[tetha_i<0]<-0
tetha_i[tetha_i>1]<-1
Y_comp <- (tetha_i*y) + ((1-tetha_i)*Y_syn)
mse_comp <- (1-tetha_i)*mse_syn
rel_mse_comp <- mse_comp/Y_comp
Syn_Comp <- data.frame(cbind(Y_syn, Y_comp,mse_syn, mse_comp, rel_mse_syn, rel_mse_comp))
colnames(Syn_Comp) = c("Y syn","Y comp", "MSE syn", "MSE comp", "Relatif MSE syn","Relatif MSE comp")
return(Syn_Comp)
}
dataset = list()
est_dataset = list()
for( j in 1:iterasi){
set.seed(j)
dataset[[j]] = est_2st_cl(dataY, f1, f2)
est_dataset [[j]] = Y_synth_comp(dataX,dataset[[j]]$Y_cap_est,dataset[[j]]$var_cap_est)
}
databind = data.frame()
databind_ind = data.frame()
for (b in 1:iterasi) {
databind = rbind(databind, dataset[[b]])
databind_ind = rbind(databind_ind, est_dataset[[b]])
}
datafinal = aggregate(list(databind$Y_cap_est,databind_ind$`Y syn`,databind_ind$`Y comp`,
databind$var_cap_est,databind_ind$`MSE syn`,databind_ind$`MSE comp`
,databind$rel_var_y,databind_ind$`Relatif MSE syn`,databind_ind$`Relatif MSE comp`),
by=list(databind$r102),
FUN=mean)
colnames(datafinal) = c("kab", "Y dir", "Y syn","Y com", "MSE dir","MSE syn","MSE com",
"Relatif MSE dir","Relatif MSE syn","Relatif MSE com")
return(datafinal)
}b. MEnggunakan rumus MSE synthetic yang perbaikanresampling2_2st_1 = function(dataY, dataX, Ni, f1, f2, iterasi ){
est_2st_cl <- function(dataY, f1,f2){
sp_2_ct <- function(dataY,f1,f2){
kb_1 <- dataY
kb_1$bs_o <- as.character(paste(kb_1$bs_o))
Mi_kb_1 <- data.frame(table(kb_1$bs_o))
colnames(Mi_kb_1) <- c("bs_o","Muatan_x")
samp_th1 <- function(data,f1){ #pengambilan sampel tahap 1
pps_wr <- function(sizes,n){
N <- length(sizes)
cumsizes <- cumsum(sizes)
totsize <- cumsizes[N]
s <- numeric(n)
for (u in 1:n) {
r <- runif(1,0,totsize)
i <- 1
while (cumsizes[i] < r) {
i <- i+1
}
s[u] <- i
}
return(s)
}
n_length <- f1*nrow(data)
data_sa <- data[pps_wr(data$Muatan_x,n_length),]
return(data_sa)
}
hs_sp_t1 <- samp_th1(Mi_kb_1,f1)
samp_th2 <- function(dt_kb_1,hs_sp_1,f2){
index <- c(1:nrow(hs_sp_1))
hs_sp_1 <- cbind(index,hs_sp_1)
dt_2 <- list();dt_3 <- list() ; samp_rt <- list() ; kode <- list() ; mi <- list()
hs_samp_srs <- list() ; hs_smp_srs <- list(); n_samp <- list()
for(i in hs_sp_1$index){ #Pengambilan sampel ssu scr srs wor
dt_2[[i]] <- subset(hs_sp_1, index %in% i)
dt_3[[i]] <- subset(dt_kb_1, bs_o %in% dt_2[[i]]$bs_o)
samp_rt[[i]] <- sample(1:length(dt_3[[i]]$pop), size = f2*length(dt_3[[i]]$pop), replace = FALSE)
hs_samp_srs[[i]] <- dt_3[[i]][samp_rt[[i]],]
n_samp [[i]] <- nrow(hs_samp_srs[[i]])
kode[[i]] <- rep(i,n_samp[[i]])
mi[[i]] <- rep(n_samp[[i]],n_samp[[i]])
hs_smp_srs[[i]] <- cbind("kode"=kode[[i]],hs_samp_srs[[i]], "mi"=mi[[i]] )
}
gb_samp = data.frame()
for (j in hs_sp_1$index) {
gb_samp <- rbind(gb_samp,hs_smp_srs[[j]])
}
return(gb_samp)
}
res_sp <- samp_th2(kb_1,hs_sp_t1,f2) #kab 1
res_sp$pop <- as.numeric(paste(res_sp$pop))
res_sp$mi <- as.numeric(paste(res_sp$mi))
y_i. <- aggregate(res_sp$pop, list(res_sp$kode), mean)
M_nol <- sum(Mi_kb_1$Muatan_x)
Y_cap_est <- 1/ (length(y_i.$x)) * (sum(y_i.$x * M_nol))
pk_1 <- M_nol*y_i.$x
pk_2 <- pk_1 - Y_cap_est
pk_2_x <- pk_2^2
var_cap_est <- (sum(pk_2_x))/(((length(y_i.$x)))*((length(y_i.$x))-1))
rel_var_y <- var_cap_est/Y_cap_est
hasil <- cbind(Y_cap_est,var_cap_est,rel_var_y)
return(hasil)
}
r102 = unique(dt2$kab)
datakoeh = list()
for(b in r102){
datakoeh[[b]] <- subset(dt2, kab %in% b)
}
try_ye <- list()
for(c in r102){
try_ye[[c]] <- sp_2_ct(data = datakoeh[[c]],f1,f2)
}
this_it = data.frame()
for (q in r102) {
this_it <- rbind(this_it, try_ye[[q]])
}
this_it = cbind(r102, this_it)
}
#fungsi estimasi synthetic dan composite
Y_synth_comp <- function(x,y,vary,Ni){
#synthetic function
N <- Ni
xty <- t(x)%*%y
xtx <- t(x)%*%x
bt_xy <- solve(xtx)%*%xty
bt_xy = t(bt_xy)
Y_syn <- c(rep(0,nrow(x)))
for(i in 1: ncol(x)){
Y_syn = Y_syn + x[,i]*bt_xy[,i]
}
Y_syn <- as.data.frame(Y_syn)
var_y_cap_is <- matrix(var(Y_syn),nrow = nrow(N))
Ni_2 <- N^2/nrow(N) #Ni kuadrat*1/m
var_Ysyn_y <- matrix(var(Y_syn - y),nrow = nrow(N))
ht_2 <- 1/N^2*var_Ysyn_y
y_sel <- data.frame((Y_syn - y)^2)
ba_2 <- (1/N^2)*(y_sel)-ht_2
mse_syn <- var_y_cap_is + (Ni_2*ba_2)
rel_mse_syn <- mse_syn/Y_syn
#composite function
tetha_i <- mse_syn/(vary+mse_syn)
tetha_i[tetha_i<0]<-0
tetha_i[tetha_i>1]<-1
Y_comp <- (tetha_i*y) + ((1-tetha_i)*Y_syn)
mse_comp <- (1-tetha_i)*mse_syn
rel_mse_comp <- mse_comp/Y_comp
Syn_Comp <- data.frame(cbind(Y_syn, Y_comp,mse_syn, mse_comp, rel_mse_syn, rel_mse_comp))
colnames(Syn_Comp) = c("Y syn","Y comp", "MSE syn", "MSE comp", "Relatif MSE syn","Relatif MSE comp")
return(Syn_Comp)
}
dataset = list()
est_dataset = list()
for( j in 1:iterasi){
set.seed(j)
dataset[[j]] = est_2st_cl(dataY, f1, f2)
est_dataset [[j]] = Y_synth_comp(dataX,dataset[[j]]$Y_cap_est,dataset[[j]]$var_cap_est,Ni)
}
databind = data.frame()
databind_ind = data.frame()
for (b in 1:iterasi) {
databind = rbind(databind, dataset[[b]])
databind_ind = rbind(databind_ind, est_dataset[[b]])
}
datafinal = aggregate(list(databind$Y_cap_est,databind_ind$`Y syn`,databind_ind$`Y comp`,
databind$var_cap_est,databind_ind$`MSE syn`,databind_ind$`MSE comp`
,databind$rel_var_y,databind_ind$`Relatif MSE syn`,databind_ind$`Relatif MSE comp`),
by=list(databind$r102),
FUN=mean)
colnames(datafinal) = c("kab", "Y dir", "Y syn","Y com", "MSE dir","MSE syn","MSE com",
"Relatif MSE dir","Relatif MSE syn","Relatif MSE com")
return(datafinal)
}res2_2sc_1_1 <- resampling_2st_2(dataY = dt, dataX = mat_x , f1 = 1 , f2 = 1, iterasi = 10)
res2_2sc2_1_1 <- resampling2_2st_2(dataY = dt, dataX = mat_x , Ni = Ni_desa, f1 = 1 , f2 = 1, iterasi = 10)
res2_2sc_0.9_0.9 <- resampling_2st_2(dataY = dt, dataX = mat_x , f1 = 0.9 , f2 = 0.9, iterasi = 10)
res2_2sc2_0.9_0.9 <- resampling2_2st_2(dataY = dt, dataX = mat_x , Ni = Ni_desa, f1 = 0.9 , f2 = 0.9, iterasi = 10)
res2_2sc_0.8_0.8 <- resampling_2st_2(dataY = dt, dataX = mat_x , f1 = 0.8 , f2 = 0.8, iterasi = 10)
res2_2sc2_0.8_0.8 <- resampling2_2st_2(dataY = dt, dataX = mat_x , Ni = Ni_desa, f1 = 0.8 , f2 = 0.8, iterasi = 10)
res2_2sc_0.7_0.7 <- resampling_2st_2(dataY = dt, dataX = mat_x , f1 = 0.7 , f2 = 0.7, iterasi = 10)
res2_2sc2_0.7_0.7 <- resampling2_2st_2(dataY = dt, dataX = mat_x , Ni = Ni_desa, f1 = 0.7 , f2 = 0.7, iterasi = 10)
res2_2sc_0.6_0.6 <- resampling_2st_2(dataY = dt, dataX = mat_x , f1 = 0.6 , f2 = 0.6, iterasi = 10)
res2_2sc2_0.6_0.6 <- resampling2_2st_2(dataY = dt, dataX = mat_x , Ni = Ni_desa, f1 = 0.6 , f2 = 0.6, iterasi = 10)
res2_2sc_0.5_0.5 <- resampling_2st_2(dataY = dt, dataX = mat_x , f1 = 0.5 , f2 = 0.5, iterasi = 10)
res2_2sc2_0.5_0.5 <- resampling2_2st_2(dataY = dt, dataX = mat_x , Ni = Ni_desa, f1 = 0.5 , f2 = 0.5, iterasi = 10)
res2_2sc_0.4_0.4 <- resampling_2st_2(dataY = dt, dataX = mat_x , f1 = 0.4 , f2 = 0.4, iterasi = 10)
res2_2sc2_0.4_0.4 <- resampling2_2st_2(dataY = dt, dataX = mat_x , Ni = Ni_desa, f1 = 0.4 , f2 = 0.4, iterasi = 10)
res2_2sc_0.3_0.3 <- resampling_2st_2(dataY = dt, dataX = mat_x , f1 = 0.3 , f2 = 0.3, iterasi = 10)
res2_2sc2_0.3_0.3 <- resampling2_2st_2(dataY = dt, dataX = mat_x , Ni = Ni_desa, f1 = 0.3 , f2 = 0.3, iterasi = 10)
res2_2sc_0.2_0.2 <- resampling_2st_2(dataY = dt, dataX = mat_x , f1 = 0.2 , f2 = 0.2, iterasi = 10)
res2_2sc2_0.2_0.2 <- resampling2_2st_2(dataY = dt, dataX = mat_x , Ni = Ni_desa, f1 = 0.2 , f2 = 0.2, iterasi = 10)res_2sc_1_1 <- resampling_2st_1(dataY = dt, dataX = mat_x , f1 = 1 , f2 = 1, iterasi = 10)
res_2sc2_1_1 <- resampling2_2st_1(dataY = dt, dataX = mat_x , Ni = Ni_desa, f1 = 1 , f2 = 1, iterasi = 10)
res_2sc_0.9_0.9 <- resampling_2st_1(dataY = dt, dataX = mat_x , f1 = 0.9 , f2 = 0.9, iterasi = 10)
res_2sc2_0.9_0.9 <- resampling2_2st_1(dataY = dt, dataX = mat_x , Ni = Ni_desa, f1 = 0.9 , f2 = 0.9, iterasi = 10)
res_2sc_0.8_0.8 <- resampling_2st_1(dataY = dt, dataX = mat_x , f1 = 0.8 , f2 = 0.8, iterasi = 10)
res_2sc2_0.8_0.8 <- resampling2_2st_1(dataY = dt, dataX = mat_x , Ni = Ni_desa, f1 = 0.8 , f2 = 0.8, iterasi = 10)
res_2sc_0.7_0.7 <- resampling_2st_1(dataY = dt, dataX = mat_x , f1 = 0.7 , f2 = 0.7, iterasi = 10)
res_2sc2_0.7_0.7 <- resampling2_2st_1(dataY = dt, dataX = mat_x , Ni = Ni_desa, f1 = 0.7 , f2 = 0.7, iterasi = 10)
res_2sc_0.6_0.6 <- resampling_2st_1(dataY = dt, dataX = mat_x , f1 = 0.6 , f2 = 0.6, iterasi = 10)
res_2sc2_0.6_0.6 <- resampling2_2st_1(dataY = dt, dataX = mat_x , Ni = Ni_desa, f1 = 0.6 , f2 = 0.6, iterasi = 10)
res_2sc_0.5_0.5 <- resampling_2st_1(dataY = dt, dataX = mat_x , f1 = 0.5 , f2 = 0.5, iterasi = 10)
res_2sc2_0.5_0.5 <- resampling2_2st_1(dataY = dt, dataX = mat_x , Ni = Ni_desa, f1 = 0.5 , f2 = 0.5, iterasi = 10)
res_2sc_0.4_0.4 <- resampling_2st_1(dataY = dt, dataX = mat_x , f1 = 0.4 , f2 = 0.4, iterasi = 10)
res_2sc2_0.4_0.4 <- resampling2_2st_1(dataY = dt, dataX = mat_x , Ni = Ni_desa, f1 = 0.4 , f2 = 0.4, iterasi = 10)
res_2sc_0.3_0.3 <- resampling_2st_1(dataY = dt, dataX = mat_x , f1 = 0.3 , f2 = 0.3, iterasi = 10)
res_2sc2_0.3_0.3 <- resampling2_2st_1(dataY = dt, dataX = mat_x , Ni = Ni_desa, f1 = 0.3 , f2 = 0.3, iterasi = 10)
res_2sc_0.2_0.2 <- resampling_2st_1(dataY = dt, dataX = mat_x , f1 = 0.2 , f2 = 0.2, iterasi = 10)
res_2sc2_0.2_0.2 <- resampling2_2st_1(dataY = dt, dataX = mat_x , Ni = Ni_desa, f1 = 0.2 , f2 = 0.2, iterasi = 10)MEnggunakan rumus MSE synthetic yang awal
res2_2sc_1_1## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1637409 1593637 1637409 0 1915959130 0
## 2 10 2452226 2411973 2452226 0 1620293809 0
## 3 11 1432925 1146047 1432925 0 82299323238 0
## 4 12 1628141 1604098 1628141 0 578061217 0
## 5 13 1660654 1539506 1660654 0 14676947727 0
## 6 14 2677707 2613785 2677707 0 4085950156 0
## 7 15 5059790 4628255 5059790 0 186222104807 0
## 8 16 2702258 2538366 2702258 0 26860745533 0
## 9 17 2063216 2755581 2063216 0 479368784452 0
## 10 18 1802921 2055582 1802921 0 63837560347 0
## 11 19 1893080 2458794 1893080 0 320033300994 0
## 12 2 1631372 1434795 1631372 0 38642737109 0
## 13 20 1933963 1528077 1933963 0 164743829852 0
## 14 21 1543883 2075747 1543883 0 282880026972 0
## 15 22 2024290 1912770 2024290 0 12436499438 0
## 16 23 2145480 1996068 2145480 0 22323787837 0
## 17 24 2721467 2809820 2721467 0 7806205403 0
## 18 25 3923463 3735418 3923463 0 35361064290 0
## 19 26 2049321 1442285 2049321 0 368492681119 0
## 20 27 1850200 1431308 1850200 0 175471128846 0
## 21 28 1632444 1718061 1632444 0 7330251063 0
## 22 29 1945828 1910735 1945828 0 1231521822 0
## 23 3 1794553 2646763 1794553 0 726262701413 0
## 24 4 2297582 2461884 2297582 0 26995005332 0
## 25 5 2083579 2122123 2083579 0 1485661569 0
## 26 6 2526531 2357735 2526531 0 28492221686 0
## 27 7 2716794 3573111 2716794 0 733277373716 0
## 28 71 2297609 2288027 2297609 0 91814618 0
## 29 72 1872873 2049375 1872873 0 31152960873 0
## 30 73 3484448 3408041 3484448 0 5838136073 0
## 31 74 2197198 2010784 2197198 0 34750233064 0
## 32 75 1867766 2087481 1867766 0 48274556977 0
## 33 76 2088732 2024320 2088732 0 4148835200 0
## 34 77 2275392 2420130 2275392 0 20949030793 0
## 35 78 7020112 6437879 7020112 0 338995297458 0
## 36 79 2178314 2013176 2178314 0 27270573971 0
## 37 8 2025630 1700059 2025630 0 105996366605 0
## 38 9 2279804 2477361 2279804 0 39028580079 0
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 0 1202.2555 0
## 2 0 671.7710 0
## 3 0 71811.4957 0
## 4 0 360.3653 0
## 5 0 9533.5455 0
## 6 0 1563.2310 0
## 7 0 40235.9187 0
## 8 0 10581.9061 0
## 9 0 173962.8639 0
## 10 0 31055.7161 0
## 11 0 130158.6221 0
## 12 0 26932.5893 0
## 13 0 107811.2248 0
## 14 0 136278.6447 0
## 15 0 6501.8253 0
## 16 0 11183.8810 0
## 17 0 2778.1870 0
## 18 0 9466.4286 0
## 19 0 255492.3460 0
## 20 0 122594.9822 0
## 21 0 4266.5836 0
## 22 0 644.5279 0
## 23 0 274396.5455 0
## 24 0 10965.1824 0
## 25 0 700.0826 0
## 26 0 12084.5759 0
## 27 0 205221.0163 0
## 28 0 40.1283 0
## 29 0 15201.2016 0
## 30 0 1713.0477 0
## 31 0 17281.9355 0
## 32 0 23125.7447 0
## 33 0 2049.4956 0
## 34 0 8656.1605 0
## 35 0 52656.3659 0
## 36 0 13546.0443 0
## 37 0 62348.6268 0
## 38 0 15754.0973 0 b. f = 90%res2_2sc_0.9_0.9## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1636986 1592285 1614610 2.095245e+12 2.092198e+12 1.046860e+12
## 2 10 2448930 2418047 2433485 6.140818e+12 6.138901e+12 3.069930e+12
## 3 11 1421352 1144085 1280277 2.252747e+12 2.175255e+12 1.106659e+12
## 4 12 1625446 1606477 1615946 2.489491e+12 2.486400e+12 1.243972e+12
## 5 13 1663225 1527483 1595103 2.848169e+12 2.828951e+12 1.419261e+12
## 6 14 2680657 2610325 2645458 5.579538e+12 5.573005e+12 2.788134e+12
## 7 15 5055180 4612941 4833343 3.156969e+13 3.137030e+13 1.573483e+13
## 8 16 2684384 2536543 2610228 3.788269e+12 3.765572e+12 1.888441e+12
## 9 17 2080308 2757703 2445370 3.231330e+12 2.771663e+12 1.491754e+12
## 10 18 1809275 2061249 1937109 2.464443e+12 2.399757e+12 1.215798e+12
## 11 19 1911671 2451197 2197455 2.682271e+12 2.390255e+12 1.263759e+12
## 12 2 1644132 1438254 1540254 2.387403e+12 2.344765e+12 1.182944e+12
## 13 20 1941041 1537992 1735087 3.865736e+12 3.701206e+12 1.890820e+12
## 14 21 1536454 2074466 1822832 2.425896e+12 2.135827e+12 1.135716e+12
## 15 22 2026726 1909133 1967754 2.969411e+12 2.953843e+12 1.480801e+12
## 16 23 2139697 1995460 2067316 3.467823e+12 3.445614e+12 1.728338e+12
## 17 24 2710160 2799417 2754848 5.668209e+12 5.658205e+12 2.831599e+12
## 18 25 3903980 3722677 3813182 1.262227e+13 1.258657e+13 6.302193e+12
## 19 26 2034148 1446449 1715904 2.283182e+12 1.936794e+12 1.047779e+12
## 20 27 1844046 1433882 1630513 2.159767e+12 1.990513e+12 1.035812e+12
## 21 28 1631040 1718602 1674944 2.557203e+12 2.547537e+12 1.276177e+12
## 22 29 1929624 1894598 1912108 3.802542e+12 3.801291e+12 1.900958e+12
## 23 3 1783932 2630155 2246438 4.280966e+12 3.563372e+12 1.944289e+12
## 24 4 2270448 2462961 2367269 5.282075e+12 5.239825e+12 2.630402e+12
## 25 5 2099107 2124879 2111996 4.026698e+12 4.025661e+12 2.013090e+12
## 26 6 2533881 2365748 2449656 8.171273e+12 8.141912e+12 4.078282e+12
## 27 7 2700917 3557805 3158840 5.786313e+12 5.050298e+12 2.696423e+12
## 28 71 2288659 2280277 2284467 4.940105e+12 4.939139e+12 2.469811e+12
## 29 72 1871972 2042432 1957928 2.781307e+12 2.749860e+12 1.382711e+12
## 30 73 3421763 3370910 3396316 1.011121e+13 1.010382e+13 5.053757e+12
## 31 74 2198368 2006751 2101217 1.503782e+12 1.465230e+12 7.421098e+11
## 32 75 1880351 2082021 1982470 2.233723e+12 2.190577e+12 1.105918e+12
## 33 76 2087826 2010660 2049154 2.522888e+12 2.514412e+12 1.259319e+12
## 34 77 2255208 2403001 2329491 4.372781e+12 4.346447e+12 2.179760e+12
## 35 78 6965927 6386184 6675443 8.087929e+13 8.054101e+13 4.035489e+13
## 36 79 2169969 2007866 2088134 1.747177e+12 1.717732e+12 8.661445e+11
## 37 8 2022379 1706131 1861715 3.245834e+12 3.145118e+12 1.597328e+12
## 38 9 2288095 2470224 2379445 6.118439e+12 6.084188e+12 3.050628e+12
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 1276953.8 1313071.1 647326.8
## 2 2505777.4 2538433.6 1260951.9
## 3 1580943.4 1900998.4 863083.9
## 4 1529183.7 1548717.9 769261.0
## 5 1710044.0 1849663.1 888503.2
## 6 2081490.0 2134202.5 1053710.7
## 7 6240168.4 6796802.9 3253218.3
## 8 1410527.4 1484230.3 723206.2
## 9 1552709.6 1005125.5 609973.8
## 10 1360170.6 1164179.1 627155.2
## 11 1402125.9 975177.1 574959.1
## 12 1451074.9 1629468.2 767541.3
## 13 1987388.7 2405683.8 1088218.5
## 14 1578426.3 1029727.4 623009.8
## 15 1463386.2 1547219.2 752008.4
## 16 1620064.9 1725932.0 835609.3
## 17 2090763.7 2021774.0 1027764.7
## 18 3229057.7 3378008.6 1650869.4
## 19 1121090.2 1337578.3 609868.8
## 20 1169786.0 1387162.4 634600.6
## 21 1566176.4 1482953.0 761547.2
## 22 1965360.0 2001008.8 991510.2
## 23 2398289.3 1354705.5 865281.6
## 24 2317903.7 2127102.7 1108807.1
## 25 1917879.5 1894337.6 952996.6
## 26 3220566.8 3438879.5 1663012.9
## 27 2141335.3 1419154.2 853299.9
## 28 2157065.1 2165845.4 1080677.2
## 29 1479264.4 1345092.3 704286.0
## 30 2951257.3 2996663.5 1486754.6
## 31 683655.4 730117.1 353052.3
## 32 1183724.0 1051430.9 556621.2
## 33 1204412.5 1249115.9 613111.9
## 34 1932122.6 1806599.5 933321.3
## 35 11605777.3 12607445.0 6042902.2
## 36 803206.2 855511.8 414236.0
## 37 1604556.7 1842700.0 857684.0
## 38 2670595.0 2461800.9 1280888.4 c. f = 80%res2_2sc_0.8_0.8## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1625221 1582494 1603835 3.324757e+12 3.320831e+12 1.661396e+12
## 2 10 2448295 2419685 2433987 1.118283e+13 1.118051e+13 5.590836e+12
## 3 11 1448698 1146046 1295749 4.488851e+12 4.395904e+12 2.220935e+12
## 4 12 1632550 1606001 1619267 4.405111e+12 4.402680e+12 2.201947e+12
## 5 13 1669021 1555428 1612136 5.245564e+12 5.231226e+12 2.619191e+12
## 6 14 2722098 2644459 2683240 1.117575e+13 1.116472e+13 5.585115e+12
## 7 15 5089721 4670848 4879965 6.186224e+13 6.168152e+13 3.088587e+13
## 8 16 2720981 2529620 2625024 6.896086e+12 6.858517e+12 3.438622e+12
## 9 17 2055816 2754683 2422797 5.697145e+12 5.204675e+12 2.719094e+12
## 10 18 1808414 2048147 1929125 4.566822e+12 4.508019e+12 2.268602e+12
## 11 19 1903053 2454138 2188244 4.810766e+12 4.504220e+12 2.326007e+12
## 12 2 1604995 1420621 1512399 4.155904e+12 4.121022e+12 2.069191e+12
## 13 20 1914344 1512889 1711220 6.955817e+12 6.792817e+12 3.436669e+12
## 14 21 1520171 2072041 1808188 4.083676e+12 3.777049e+12 1.961729e+12
## 15 22 2006261 1904671 1955385 5.282760e+12 5.269499e+12 2.638059e+12
## 16 23 2152017 1998388 2075007 6.518963e+12 6.492174e+12 3.252763e+12
## 17 24 2715240 2815421 2765375 1.029070e+13 1.027944e+13 5.142532e+12
## 18 25 3966398 3761964 3864089 2.444537e+13 2.440280e+13 1.221203e+13
## 19 26 2055269 1425197 1725055 4.373192e+12 3.973835e+12 2.081913e+12
## 20 27 1835848 1431458 1628663 3.569666e+12 3.404916e+12 1.742619e+12
## 21 28 1620923 1710788 1665926 4.453330e+12 4.443578e+12 2.224222e+12
## 22 29 1943464 1911379 1927420 7.221622e+12 7.220510e+12 3.610533e+12
## 23 3 1842051 2670374 2277620 8.109869e+12 7.420943e+12 3.873599e+12
## 24 4 2297315 2459346 2378559 1.078545e+13 1.075383e+13 5.384791e+12
## 25 5 2089675 2126603 2108177 6.712039e+12 6.706089e+12 3.354529e+12
## 26 6 2528929 2368646 2448702 1.409092e+13 1.406323e+13 7.038528e+12
## 27 7 2737518 3617489 3194816 1.070168e+13 9.922251e+12 5.148165e+12
## 28 71 2300918 2272372 2286640 8.865780e+12 8.858944e+12 4.431180e+12
## 29 72 1833508 2030200 1932745 4.686714e+12 4.641032e+12 2.331806e+12
## 30 73 3484788 3420673 3452714 1.936067e+13 1.935021e+13 9.677719e+12
## 31 74 2181289 1991418 2085593 2.602018e+12 2.563900e+12 1.291397e+12
## 32 75 1821274 2068658 1946838 3.444609e+12 3.376938e+12 1.705078e+12
## 33 76 2086997 2012329 2049595 4.360058e+12 4.349213e+12 2.177313e+12
## 34 77 2258396 2403600 2331322 7.738031e+12 7.705019e+12 3.860718e+12
## 35 78 7099338 6509233 6803920 1.463126e+14 1.459581e+14 7.306754e+13
## 36 79 2150902 1993512 2071675 3.117607e+12 3.084941e+12 1.550569e+12
## 37 8 2027749 1692762 1858624 6.022789e+12 5.908984e+12 2.982660e+12
## 38 9 2288287 2474153 2381406 1.173991e+13 1.170213e+13 5.860487e+12
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 2032685 2090019 1030274.2
## 2 4566724 4620312 2296543.3
## 3 3092752 3831573 1711160.9
## 4 2696601 2740654 1359001.7
## 5 3131391 3358287 1620246.2
## 6 4100788 4218799 2079143.4
## 7 12139687 13196771 6322836.6
## 8 2533812 2710343 1309519.4
## 9 2759983 1888555 1120221.0
## 10 2521602 2200023 1174747.7
## 11 2522982 1835214 1061833.2
## 12 2586338 2898199 1366595.8
## 13 3616762 4478723 2000762.0
## 14 2674211 1822407 1082930.5
## 15 2624289 2764529 1346091.4
## 16 3025185 3249114 1566349.5
## 17 3777222 3644732 1854820.2
## 18 6157205 6481083 3157462.8
## 19 2125012 2788829 1205939.7
## 20 1941429 2377392 1068622.0
## 21 2741535 2593764 1332581.1
## 22 3692732 3753689 1861469.1
## 23 4374149 2772263 1695215.9
## 24 4661471 4366002 2253795.5
## 25 3194154 3152241 1586120.4
## 26 5557894 5929686 2868700.1
## 27 3905090 2743189 1610566.3
## 28 3836700 3896715 1932727.3
## 29 2535705 2279286 1199446.4
## 30 5543322 5650903 2797952.5
## 31 1191545 1287263 618745.3
## 32 1877458 1631210 872041.7
## 33 2078207 2160381 1059002.4
## 34 3413150 3206199 1652138.5
## 35 20582448 22399298 10725982.9
## 36 1440611 1548755 746152.1
## 37 2963341 3487098 1601868.6
## 38 5120505 4725111 2457031.7 d. f = 70%res2_2sc_0.7_0.7## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1616307 1578886 1597547 4.456075e+12 4.449383e+12 2.226360e+12
## 2 10 2418382 2405519 2411950 1.495908e+13 1.495504e+13 7.478531e+12
## 3 11 1449058 1169190 1307899 5.753112e+12 5.666065e+12 2.854576e+12
## 4 12 1638281 1600037 1619144 6.218126e+12 6.214070e+12 3.108048e+12
## 5 13 1651059 1537062 1593986 6.845293e+12 6.829498e+12 3.418692e+12
## 6 14 2635718 2594061 2614877 1.436788e+13 1.435801e+13 7.181470e+12
## 7 15 4956678 4570370 4763290 7.702523e+13 7.685770e+13 3.847067e+13
## 8 16 2689935 2548547 2619092 8.829726e+12 8.802099e+12 4.407941e+12
## 9 17 2116274 2749459 2441114 8.831415e+12 8.423534e+12 4.310985e+12
## 10 18 1812652 2048644 1931230 6.352914e+12 6.295497e+12 3.162027e+12
## 11 19 1884646 2449682 2173647 7.330976e+12 7.010324e+12 3.583469e+12
## 12 2 1602153 1433267 1517414 5.594439e+12 5.559969e+12 2.788566e+12
## 13 20 1904186 1540388 1720824 1.039723e+13 1.024520e+13 5.160285e+12
## 14 21 1542188 2062208 1810056 5.825114e+12 5.547794e+12 2.840981e+12
## 15 22 2031159 1896481 1963695 7.387406e+12 7.364446e+12 3.687952e+12
## 16 23 2139800 1971372 2055376 8.400902e+12 8.366207e+12 4.191752e+12
## 17 24 2734362 2815637 2775025 1.378755e+13 1.377748e+13 6.891255e+12
## 18 25 3883256 3707509 3795305 2.826353e+13 2.822327e+13 1.412169e+13
## 19 26 2010666 1451640 1722981 5.658565e+12 5.340924e+12 2.747508e+12
## 20 27 1810387 1412662 1607810 4.600338e+12 4.434958e+12 2.258035e+12
## 21 28 1634796 1714692 1674801 6.505102e+12 6.496529e+12 3.250403e+12
## 22 29 1967574 1933639 1950606 1.059863e+13 1.059735e+13 5.298996e+12
## 23 3 1783744 2579527 2196103 1.064550e+13 1.000564e+13 5.156725e+12
## 24 4 2309992 2454956 2382707 1.585419e+13 1.581458e+13 7.917161e+12
## 25 5 2087482 2109798 2098667 1.085794e+13 1.084984e+13 5.426941e+12
## 26 6 2493951 2339984 2416884 1.882024e+13 1.878858e+13 9.402196e+12
## 27 7 2745183 3539126 3150858 1.521725e+13 1.458276e+13 7.446467e+12
## 28 71 2345345 2281799 2313529 1.245700e+13 1.244093e+13 6.224476e+12
## 29 72 1937661 2044848 1991523 7.676137e+12 7.649751e+12 3.831440e+12
## 30 73 3367689 3338240 3352955 2.451089e+13 2.448213e+13 1.224825e+13
## 31 74 2166577 2010224 2088060 3.593633e+12 3.566266e+12 1.789938e+12
## 32 75 1860847 2087107 1974674 5.748848e+12 5.693340e+12 2.860449e+12
## 33 76 2032231 1994677 2013408 5.312315e+12 5.301073e+12 2.653341e+12
## 34 77 2282816 2405977 2344701 1.120625e+13 1.117676e+13 5.595699e+12
## 35 78 6878159 6308461 6593033 1.787230e+14 1.783903e+14 8.927823e+13
## 36 79 2151334 2010293 2080546 4.157226e+12 4.131819e+12 2.072230e+12
## 37 8 2023873 1705086 1863341 7.942447e+12 7.835658e+12 3.944323e+12
## 38 9 2271166 2466509 2369014 1.488440e+13 1.484203e+13 7.431584e+12
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 2728019 2809475 1383392.4
## 2 6178160 6221204 3099275.1
## 3 3914544 4824928 2158939.1
## 4 3788135 3880354 1916357.8
## 5 4130301 4455796 2142955.3
## 6 5437616 5530388 2741078.7
## 7 15461588 16784154 8046557.0
## 8 3269661 3451080 1678594.4
## 9 4163760 3067338 1764967.3
## 10 3502501 3070874 1635954.4
## 11 3890461 2860972 1648308.8
## 12 3450043 3871823 1823455.0
## 13 5395920 6648665 2975690.7
## 14 3752715 2688636 1564098.3
## 15 3620623 3881276 1872799.9
## 16 3908042 4243791 2033823.9
## 17 5037709 4893566 2482007.2
## 18 7219994 7573378 3695737.7
## 19 2806124 3681426 1592017.0
## 20 2506955 3131410 1391760.7
## 21 3972698 3791038 1939482.7
## 22 5377194 5472057 2712084.8
## 23 5908922 3866200 2334083.5
## 24 6776416 6452746 3302254.6
## 25 5170901 5138855 2576332.3
## 26 7495904 8006326 3870231.9
## 27 5540779 4121088 2362818.1
## 28 5276906 5449484 2679642.4
## 29 3921581 3738709 1911764.8
## 30 7231743 7316166 3634562.7
## 31 1654849 1773763 856015.8
## 32 3086749 2728362 1447672.6
## 33 2593562 2652421 1310552.9
## 34 4871898 4649745 2377508.4
## 35 25888108 28179083 13492107.4
## 36 1921725 2055596 992953.5
## 37 3903917 4591287 2109415.6
## 38 6535503 6007028 3129446.7 e. f = 60%res2_2sc_0.6_0.6## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1587576 1578601 1583102 5.432394e+12 5.423606e+12 2.713995e+12
## 2 10 2427684 2399478 2413577 1.797021e+13 1.796416e+13 8.983594e+12
## 3 11 1401387 1181179 1290766 6.607451e+12 6.554421e+12 3.290399e+12
## 4 12 1605746 1601218 1603481 6.851222e+12 6.840734e+12 3.422985e+12
## 5 13 1644666 1545833 1595167 8.084347e+12 8.070806e+12 4.038780e+12
## 6 14 2701051 2586828 2643883 1.859297e+13 1.857407e+13 9.291754e+12
## 7 15 4825377 4477796 4651376 7.614969e+13 7.599912e+13 3.803715e+13
## 8 16 2696272 2511261 2603575 1.309290e+13 1.305036e+13 6.535789e+12
## 9 17 2028061 2713785 2382801 8.428132e+12 7.949970e+12 4.090130e+12
## 10 18 1826595 2047310 1937587 7.699168e+12 7.643335e+12 3.835530e+12
## 11 19 1840162 2399841 2127498 7.552226e+12 7.233983e+12 3.694330e+12
## 12 2 1588838 1447062 1517760 5.969805e+12 5.943447e+12 2.978292e+12
## 13 20 1986656 1545138 1764048 1.401575e+13 1.380563e+13 6.954856e+12
## 14 21 1506534 2037784 1781546 6.340563e+12 6.043219e+12 3.092976e+12
## 15 22 2027951 1933005 1980430 8.547682e+12 8.534984e+12 4.270663e+12
## 16 23 2142823 2001454 2071975 9.828678e+12 9.798882e+12 4.906869e+12
## 17 24 2728370 2813293 2770946 1.590379e+13 1.588248e+13 7.946551e+12
## 18 25 3850988 3657485 3754127 3.611787e+13 3.605842e+13 1.804405e+13
## 19 26 2077600 1458565 1759050 6.997647e+12 6.609750e+12 3.398971e+12
## 20 27 1831280 1429147 1626982 5.853109e+12 5.681795e+12 2.882983e+12
## 21 28 1612772 1706487 1659770 7.609203e+12 7.588452e+12 3.799399e+12
## 22 29 1905717 1873519 1889616 9.832287e+12 9.831042e+12 4.915832e+12
## 23 3 1773547 2540345 2171883 1.020944e+13 9.611443e+12 4.949100e+12
## 24 4 2333637 2421249 2377563 1.869100e+13 1.866656e+13 9.339373e+12
## 25 5 2080418 2087722 2084103 1.207070e+13 1.206103e+13 6.032928e+12
## 26 6 2414774 2332135 2373435 1.926514e+13 1.925274e+13 9.629468e+12
## 27 7 2695261 3475228 3093646 1.621017e+13 1.559230e+13 7.947156e+12
## 28 71 2249927 2220809 2235364 1.388323e+13 1.386779e+13 6.937751e+12
## 29 72 1849993 1991523 1921148 6.394576e+12 6.364502e+12 3.189721e+12
## 30 73 3349716 3227808 3288709 2.696894e+13 2.693909e+13 1.347700e+13
## 31 74 2182707 1961266 2070970 3.990798e+12 3.931672e+12 1.980463e+12
## 32 75 1805328 2033623 1922817 5.827641e+12 5.749651e+12 2.893462e+12
## 33 76 1990009 1954471 1972260 4.930874e+12 4.919952e+12 2.462698e+12
## 34 77 2122599 2326346 2225660 1.017614e+13 1.011525e+13 5.072580e+12
## 35 78 6605093 6065277 6334971 2.048828e+14 2.045797e+14 1.023656e+14
## 36 79 2100811 1960406 2030333 4.788615e+12 4.762042e+12 2.387630e+12
## 37 8 2030356 1707254 1867714 9.174525e+12 9.062294e+12 4.558996e+12
## 38 9 2257793 2434545 2346304 1.779779e+13 1.776213e+13 8.889965e+12
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 3311768 3398893 1675849.0
## 2 7384362 7476804 3714332.0
## 3 4682943 5534222 2535172.2
## 4 4225859 4255445 2118169.2
## 5 4880990 5188883 2514303.0
## 6 6875284 7160029 3506690.7
## 7 15586780 16838929 8092619.9
## 8 4841793 5198553 2506373.4
## 9 4125485 2928516 1710521.8
## 10 4176776 3727098 1968143.6
## 11 4080624 3011327 1731363.3
## 12 3715849 4096271 1947311.6
## 13 7005101 8916481 3920181.1
## 14 4148414 2970627 1725552.2
## 15 4208514 4416751 2154619.4
## 16 4559335 4900462 2360671.2
## 17 5784689 5654319 2858077.4
## 18 9300757 9828838 4777368.8
## 19 3359730 4530452 1928756.5
## 20 3161816 3966748 1758785.0
## 21 4689484 4447121 2279814.0
## 22 5087226 5171876 2564553.5
## 23 5684633 3779808 2266510.0
## 24 7940934 7706467 3908114.6
## 25 5776207 5787549 2889489.3
## 26 7908411 8237513 4034052.7
## 27 5983669 4481371 2560782.2
## 28 6078456 6220591 3072280.4
## 29 3388443 3178161 1638523.2
## 30 8000599 8323545 4078163.0
## 31 1825532 2009402 956174.7
## 32 3140194 2821165 1482041.6
## 33 2445368 2504074 1236442.0
## 34 4708199 4313367 2248109.1
## 35 30942714 33658694 16120900.4
## 36 2269119 2429638 1172916.5
## 37 4506254 5312690 2437223.8
## 38 7766151 7197747 3734765.8 f. f = 50%res2_2sc_0.5_0.5## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1627076 1628308 1627706 6.813104e+12 6.804606e+12 3.404425e+12
## 2 10 2463634 2403272 2433435 2.060266e+13 2.059169e+13 1.029859e+13
## 3 11 1528357 1173926 1349900 1.050616e+13 1.036958e+13 5.218681e+12
## 4 12 1578191 1597892 1588049 6.840831e+12 6.837279e+12 3.419527e+12
## 5 13 1668720 1535740 1601952 9.661709e+12 9.632895e+12 4.823610e+12
## 6 14 2577312 2589395 2583358 1.758382e+13 1.757337e+13 8.789295e+12
## 7 15 5005324 4578876 4791843 1.091351e+14 1.089180e+14 5.451321e+13
## 8 16 2676179 2499590 2587600 1.367715e+13 1.361721e+13 6.823526e+12
## 9 17 2034597 2703588 2379620 9.974913e+12 9.515098e+12 4.868595e+12
## 10 18 1799060 1995348 1897857 9.088492e+12 9.036188e+12 4.531062e+12
## 11 19 1893053 2440797 2174258 8.655055e+12 8.344195e+12 4.247564e+12
## 12 2 1595747 1433322 1514346 7.713578e+12 7.683416e+12 3.849229e+12
## 13 20 1915551 1501319 1706518 1.236782e+13 1.217568e+13 6.135337e+12
## 14 21 1560879 2051373 1810236 8.992676e+12 8.747094e+12 4.433846e+12
## 15 22 1980581 1868042 1924222 8.104444e+12 8.087443e+12 4.047963e+12
## 16 23 2129981 1964478 2046910 1.106724e+13 1.101353e+13 5.520138e+12
## 17 24 2690581 2832940 2761914 1.276272e+13 1.273461e+13 6.374310e+12
## 18 25 4017313 3757820 3887425 3.993094e+13 3.985470e+13 1.994639e+13
## 19 26 2068217 1422399 1735066 7.404874e+12 6.970641e+12 3.590320e+12
## 20 27 1851414 1410402 1627185 7.038104e+12 6.839315e+12 3.468490e+12
## 21 28 1626755 1700916 1663901 8.253430e+12 8.239676e+12 4.123270e+12
## 22 29 1992458 1962151 1977303 1.360507e+13 1.360402e+13 6.802272e+12
## 23 3 1771144 2615418 2211941 1.396588e+13 1.323656e+13 6.792086e+12
## 24 4 2249290 2429616 2339910 1.962964e+13 1.957567e+13 9.801251e+12
## 25 5 2055809 2082080 2069068 1.164174e+13 1.160754e+13 5.812296e+12
## 26 6 2366850 2256948 2311842 1.723462e+13 1.721327e+13 8.611966e+12
## 27 7 2661085 3536624 3109676 1.801744e+13 1.722860e+13 8.806380e+12
## 28 71 2271207 2279727 2275534 1.597097e+13 1.594209e+13 7.978248e+12
## 29 72 1865281 2041868 1956302 8.992426e+12 8.921028e+12 4.477732e+12
## 30 73 3443270 3348332 3395739 3.472188e+13 3.468734e+13 1.735229e+13
## 31 74 2182305 2006444 2093574 5.092680e+12 5.041630e+12 2.533423e+12
## 32 75 1852538 2084857 1970684 7.901642e+12 7.817190e+12 3.929259e+12
## 33 76 2045557 1991975 2018725 7.323012e+12 7.311767e+12 3.658691e+12
## 34 77 2225991 2408124 2317889 1.164182e+13 1.157186e+13 5.803200e+12
## 35 78 6904907 6313135 6608784 2.437727e+14 2.434028e+14 1.217938e+14
## 36 79 2163777 2005724 2084462 6.107438e+12 6.072865e+12 3.045041e+12
## 37 8 2050571 1718862 1883446 9.938578e+12 9.804626e+12 4.935502e+12
## 38 9 2226086 2445019 2335765 1.715294e+13 1.709872e+13 8.562886e+12
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 4078432 4158173 2056992
## 2 8341444 8549294 4220796
## 3 6712699 8766808 3798360
## 4 4308286 4274309 2144981
## 5 5748686 6169284 2972835
## 6 6824182 6772257 3397514
## 7 21655610 23681538 11308755
## 8 5036517 5484318 2622647
## 9 4845044 3522174 2036029
## 10 4985654 4544847 2374727
## 11 4485846 3413243 1934823
## 12 4791645 5328505 2522067
## 13 6402746 8095852 3571188
## 14 5719407 4253860 2437536
## 15 4087356 4335538 2103373
## 16 5093724 5576440 2658688
## 17 4730269 4495458 2304294
## 18 9913736 10598115 5121643
## 19 3549338 4899891 2056931
## 20 3770060 4818646 2114719
## 21 5028355 4817154 2458359
## 22 6543959 6648019 3297776
## 23 7641053 5007219 3014898
## 24 8371104 8006457 4087539
## 25 5505147 5575101 2764817
## 26 7226514 7618311 3706907
## 27 6736195 4855477 2818810
## 28 6828918 6961321 3442570
## 29 4540154 4347186 2213873
## 30 10000325 10347062 5082468
## 31 2320888 2509594 1204860
## 32 4170085 3759266 1971998
## 33 3553118 3655728 1800779
## 34 5259236 4818354 2508846
## 35 35047495 38296666 18298756
## 36 2798946 3024854 1453150
## 37 4768211 5709588 2594520
## 38 7637790 6932536 3632899 g. f = 40%res2_2sc_0.4_0.4## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1608938 1557611 1583097 5.190718e+12 5.161585e+12 2.588038e+12
## 2 10 2416227 2402915 2409550 2.147376e+13 2.144183e+13 1.072889e+13
## 3 11 1401256 1101907 1249790 8.043403e+12 7.925895e+12 3.991948e+12
## 4 12 1586585 1536892 1561699 8.344274e+12 8.329667e+12 4.168479e+12
## 5 13 1630510 1525620 1577938 9.510893e+12 9.490133e+12 4.750241e+12
## 6 14 2655271 2580683 2617933 2.276281e+13 2.272689e+13 1.137241e+13
## 7 15 4938873 4576759 4757511 9.295250e+13 9.270770e+13 4.641494e+13
## 8 16 2677275 2532628 2604878 1.568204e+13 1.563765e+13 7.829898e+12
## 9 17 2093857 2758605 2444796 1.287228e+13 1.238285e+13 6.305529e+12
## 10 18 1780489 2033975 1909256 7.856432e+12 7.757931e+12 3.903164e+12
## 11 19 1913709 2430644 2177931 9.315265e+12 9.036903e+12 4.586327e+12
## 12 2 1571062 1433270 1501982 7.905083e+12 7.874264e+12 3.944811e+12
## 13 20 1989441 1510772 1747356 1.153874e+13 1.128360e+13 5.704798e+12
## 14 21 1521635 2044751 1790648 7.548474e+12 7.266961e+12 3.701673e+12
## 15 22 2080768 1917682 1998773 1.190246e+13 1.184517e+13 5.936807e+12
## 16 23 2106348 1935569 2020571 9.360562e+12 9.316060e+12 4.669085e+12
## 17 24 2685221 2840781 2763808 1.692074e+13 1.686302e+13 8.445755e+12
## 18 25 4021220 3726820 3873677 5.396468e+13 5.380194e+13 2.694153e+13
## 19 26 2051055 1443206 1737700 8.205829e+12 7.781608e+12 3.992972e+12
## 20 27 1857142 1371470 1609555 7.689773e+12 7.428076e+12 3.778108e+12
## 21 28 1615073 1707875 1662694 7.094971e+12 7.061009e+12 3.538723e+12
## 22 29 1917013 1879589 1898297 1.127190e+13 1.127001e+13 5.635477e+12
## 23 3 1628841 2564287 2152947 7.308381e+12 6.420704e+12 3.403981e+12
## 24 4 2293359 2452534 2374371 1.973957e+13 1.962971e+13 9.841975e+12
## 25 5 2063273 2098413 2081038 1.520771e+13 1.517669e+13 7.596059e+12
## 26 6 2459345 2346421 2402814 2.558885e+13 2.555387e+13 1.278567e+13
## 27 7 2682244 3542710 3122023 2.048295e+13 1.971744e+13 1.004557e+13
## 28 71 2147234 2230905 2189785 1.683686e+13 1.676135e+13 8.399396e+12
## 29 72 1882179 1988176 1936550 1.017501e+13 1.012687e+13 5.075130e+12
## 30 73 3382618 3316332 3349254 3.702574e+13 3.694635e+13 1.849293e+13
## 31 74 2237472 1954930 2094492 5.111396e+12 5.013271e+12 2.530828e+12
## 32 75 1741644 2031275 1889871 4.721770e+12 4.607195e+12 2.331378e+12
## 33 76 1991465 1918001 1954450 7.371362e+12 7.324064e+12 3.673769e+12
## 34 77 2231703 2372365 2303561 1.277247e+13 1.269579e+13 6.366683e+12
## 35 78 7006792 6381056 6693635 2.580246e+14 2.575865e+14 1.289027e+14
## 36 79 2040377 1953384 1996526 3.763868e+12 3.740655e+12 1.876067e+12
## 37 8 1978746 1685276 1831282 1.123883e+13 1.113707e+13 5.593833e+12
## 38 9 2287518 2487688 2387862 2.331137e+13 2.325620e+13 1.164185e+13
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 3153396 3305405 1610077.7
## 2 8669037 8903001 4385860.2
## 3 5602197 7221382 3144957.9
## 4 5226260 5378019 2647447.1
## 5 5708137 6137101 2954877.4
## 6 8469026 8867060 4326480.8
## 7 17897437 19665198 9358970.1
## 8 5767054 6151084 2974510.9
## 9 5905521 4521412 2547718.6
## 10 4250082 3827414 2006724.1
## 11 4816394 3713903 2094288.4
## 12 4926134 5481898 2591957.1
## 13 5672017 7482389 3223154.4
## 14 4854094 3543024 2044111.4
## 15 5687518 6158502 2952139.9
## 16 4437752 4827331 2310070.3
## 17 6078669 5853520 2978828.5
## 18 13014893 14211800 6786495.6
## 19 3953713 5406710 2279209.3
## 20 4014553 5381190 2296374.5
## 21 4286485 4111284 2093962.9
## 22 5753529 5862985 2903776.7
## 23 4262099 2436615 1537332.1
## 24 7968161 7959607 3966028.3
## 25 7195689 7203002 3593540.4
## 26 9976121 10640750 5144614.3
## 27 7560726 5533414 3192014.2
## 28 7432775 7511437 3721771.5
## 29 5051538 4972576 2499079.6
## 30 10737960 10994332 5425413.0
## 31 2267684 2579357 1205975.0
## 32 2694388 2266741 1226902.6
## 33 3544299 3756039 1819954.4
## 34 5599295 5360076 2729553.7
## 35 36287010 39936877 19008900.6
## 36 1834830 1914585 936069.4
## 37 5587364 6610333 3025718.7
## 38 9996805 9209102 4789168.3 h. f = 30%res2_2sc_0.3_0.3## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1557814 1551710 1554777 7.663532e+12 7.653063e+12 3.829144e+12
## 2 10 2418841 2380705 2399763 2.720478e+13 2.717776e+13 1.359563e+13
## 3 11 1434265 1162108 1297300 6.855437e+12 6.772509e+12 3.406848e+12
## 4 12 1684398 1614045 1649052 9.015394e+12 8.969247e+12 4.496109e+12
## 5 13 1518427 1409623 1463594 9.472012e+12 9.420552e+12 4.723030e+12
## 6 14 2663220 2542298 2602623 2.182192e+13 2.177226e+13 1.089852e+13
## 7 15 4968500 4518861 NaN Inf Inf NaN
## 8 16 2733113 2494468 2613102 1.332008e+13 1.323031e+13 6.637448e+12
## 9 17 1921628 2674501 2326564 8.169789e+12 7.579252e+12 3.924254e+12
## 10 18 1751501 2035811 1894954 8.618005e+12 8.520110e+12 4.284271e+12
## 11 19 1876156 2366768 2133329 7.418348e+12 7.141823e+12 3.635780e+12
## 12 2 1685573 1393469 1538185 7.182960e+12 7.087803e+12 3.567443e+12
## 13 20 1902736 1511355 1704861 1.218339e+13 1.196899e+13 6.037441e+12
## 14 21 1470119 2016231 1754656 7.598645e+12 7.276567e+12 3.714944e+12
## 15 22 1945898 1866243 1905930 8.841754e+12 8.815941e+12 4.414402e+12
## 16 23 2183698 2051402 2117247 1.075137e+13 1.070771e+13 5.364718e+12
## 17 24 2572417 2705935 2639454 1.631171e+13 1.627237e+13 8.145979e+12
## 18 25 3771490 3589609 NaN Inf Inf NaN
## 19 26 1943330 1406841 NaN Inf Inf NaN
## 20 27 1708308 1425428 1563803 3.663351e+12 3.561781e+12 1.805543e+12
## 21 28 1586629 1615170 1601043 5.586135e+12 5.556691e+12 2.785668e+12
## 22 29 1973889 1950889 1962387 1.736921e+13 1.736805e+13 8.684315e+12
## 23 3 1758122 2505354 2153269 1.222021e+13 1.163103e+13 5.953524e+12
## 24 4 2347631 2390382 2369614 2.845884e+13 2.835137e+13 1.420241e+13
## 25 5 1940240 2054353 1997634 9.457369e+12 9.424138e+12 4.720319e+12
## 26 6 2453888 2284372 NaN Inf Inf NaN
## 27 7 2630302 3474058 3062679 1.984480e+13 1.910442e+13 9.732379e+12
## 28 71 2042854 2221886 NaN Inf Inf NaN
## 29 72 1795388 1989355 1895935 5.961947e+12 5.876640e+12 2.958812e+12
## 30 73 3128778 3204950 3167099 3.552441e+13 3.545491e+13 1.774477e+13
## 31 74 2209428 1959381 2082362 5.018326e+12 4.912539e+12 2.482198e+12
## 32 75 1767609 2031637 1907242 6.719013e+12 6.617548e+12 3.331896e+12
## 33 76 2249655 2003144 2125137 9.591414e+12 9.465859e+12 4.763922e+12
## 34 77 2338948 2276602 NaN Inf Inf NaN
## 35 78 6572536 6018285 NaN Inf Inf NaN
## 36 79 1992904 1961236 1979320 3.811316e+12 3.744117e+12 1.888099e+12
## 37 8 1968301 1675271 1820513 1.136806e+13 1.123497e+13 5.650450e+12
## 38 9 2296794 2431599 2364370 2.499399e+13 2.495183e+13 1.248643e+13
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 4534153 4673297 2299461.3
## 2 10819093 11206009 5499006.1
## 3 4648710 5820822 2581784.4
## 4 5174193 5544500 2667289.5
## 5 5440042 6137161 2878290.4
## 6 7972815 8574646 4126514.9
## 7 Inf Inf NaN
## 8 4830457 5356946 2536969.2
## 9 4127633 2814360 1664650.2
## 10 4809539 4082386 2204620.3
## 11 3735501 2978470 1646639.2
## 12 4091819 4863272 2220769.0
## 13 5910714 7710461 3334390.2
## 14 4968568 3615930 2082111.6
## 15 4423524 4634236 2260183.4
## 16 4766337 5151876 2472301.2
## 17 5868201 5766621 2905309.5
## 18 Inf Inf NaN
## 19 Inf Inf NaN
## 20 2138116 2492452 1149988.3
## 21 3345460 3350848 1669112.6
## 22 8252392 8334966 4146619.6
## 23 6486564 4517470 2643842.1
## 24 10944659 11727025 5640500.8
## 25 4856744 4582940 2354665.9
## 26 Inf Inf NaN
## 27 7481083 5478902 3158548.0
## 28 Inf Inf NaN
## 29 3046012 2840464 1462342.4
## 30 10878292 10943991 5445197.4
## 31 2203889 2507538 1171512.5
## 32 3675499 3237948 1715023.8
## 33 4001007 4579972 2129151.2
## 34 Inf Inf NaN
## 35 Inf Inf NaN
## 36 1822887 1929381 933550.4
## 37 5528810 6604961 3004304.1
## 38 10764413 10309504 5259880.6 i. f = 20%res2_2sc_0.2_0.2## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1433981 1568561 NaN Inf Inf NaN
## 2 10 2260716 2254727 NaN Inf Inf NaN
## 3 11 1492489 1307961 NaN Inf Inf NaN
## 4 12 1512705 1557990 NaN Inf Inf NaN
## 5 13 1852117 1696136 NaN Inf Inf NaN
## 6 14 2595039 2522857 NaN Inf Inf NaN
## 7 15 4196157 4093781 NaN Inf Inf NaN
## 8 16 2802059 2462385 NaN Inf Inf NaN
## 9 17 2044691 2604661 NaN Inf Inf NaN
## 10 18 1674359 1992998 NaN Inf Inf NaN
## 11 19 1872288 2334142 NaN Inf Inf NaN
## 12 2 1620708 1541808 NaN Inf Inf NaN
## 13 20 1968137 1569962 NaN Inf Inf NaN
## 14 21 1473338 1996094 NaN Inf Inf NaN
## 15 22 1924211 1975112 NaN Inf Inf NaN
## 16 23 2173965 1942761 NaN Inf Inf NaN
## 17 24 2663297 2825912 NaN Inf Inf NaN
## 18 25 3824371 3434742 NaN Inf Inf NaN
## 19 26 2157053 1542029 NaN Inf Inf NaN
## 20 27 1994931 1419642 NaN Inf Inf NaN
## 21 28 1451899 1708893 NaN Inf Inf NaN
## 22 29 1757749 1743382 NaN Inf Inf NaN
## 23 3 1731292 2347761 NaN Inf Inf NaN
## 24 4 2623709 2343380 NaN Inf Inf NaN
## 25 5 1934975 2044576 NaN Inf Inf NaN
## 26 6 2276835 2254195 NaN Inf Inf NaN
## 27 7 2533516 3329541 NaN Inf Inf NaN
## 28 71 2006234 2096386 NaN Inf Inf NaN
## 29 72 1493890 1885636 NaN Inf Inf NaN
## 30 73 3150922 2984196 NaN Inf Inf NaN
## 31 74 2138921 1858499 1996207 6.082540e+12 5.943340e+12 3.005781e+12
## 32 75 1740888 1929846 NaN Inf Inf NaN
## 33 76 1879651 1815536 1847692 6.690393e+12 6.614338e+12 3.325842e+12
## 34 77 1950094 2215518 NaN Inf Inf NaN
## 35 78 5988669 5482654 NaN Inf Inf NaN
## 36 79 2114429 1855137 NaN Inf Inf NaN
## 37 8 2039894 1700544 NaN Inf Inf NaN
## 38 9 2236321 2346558 NaN Inf Inf NaN
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 Inf Inf NaN
## 2 Inf Inf NaN
## 3 Inf Inf NaN
## 4 Inf Inf NaN
## 5 Inf Inf NaN
## 6 Inf Inf NaN
## 7 Inf Inf NaN
## 8 Inf Inf NaN
## 9 Inf Inf NaN
## 10 Inf Inf NaN
## 11 Inf Inf NaN
## 12 Inf Inf NaN
## 13 Inf Inf NaN
## 14 Inf Inf NaN
## 15 Inf Inf NaN
## 16 Inf Inf NaN
## 17 Inf Inf NaN
## 18 Inf Inf NaN
## 19 Inf Inf NaN
## 20 Inf Inf NaN
## 21 Inf Inf NaN
## 22 Inf Inf NaN
## 23 Inf Inf NaN
## 24 Inf Inf NaN
## 25 Inf Inf NaN
## 26 Inf Inf NaN
## 27 Inf Inf NaN
## 28 Inf Inf NaN
## 29 Inf Inf NaN
## 30 Inf Inf NaN
## 31 2752798 3174597 1471398
## 32 Inf Inf NaN
## 33 2567993 2818886 1337693
## 34 Inf Inf NaN
## 35 Inf Inf NaN
## 36 Inf Inf NaN
## 37 Inf Inf NaN
## 38 Inf Inf NaN 2. Metode PPSWR-SRSWOR res_2sc_1_1## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1612935 1560421 1583113 7845037888 6246258886 3400621228
## 2 10 2415381 2395026 2408539 7399776586 7100693339 3211221114
## 3 11 1431381 1199537 1395184 8245932030 51450562921 6875915199
## 4 12 1637684 1623131 1630113 8212889932 5248596301 3073317783
## 5 13 1665954 1542576 1633569 5690831695 23912890042 3815132250
## 6 14 2703353 2599541 2670319 12564929983 21881264350 6422952083
## 7 15 5014684 4559697 4861729 74159897116 161774665282 46847594703
## 8 16 2677256 2581864 2637922 7115205101 11210156075 3527281453
## 9 17 2068793 2754404 2079597 7137834371 469647155626 7003908201
## 10 18 1814688 2075180 1841400 5324417946 72010122784 4210197645
## 11 19 1933547 2456105 1947615 7093823887 271617910785 6880491460
## 12 2 1629336 1424421 1591631 7129361218 40864249772 5616990304
## 13 20 1929510 1560682 1899958 10605830742 135051953779 9708669813
## 14 21 1527249 2059804 1540420 6618947823 285491617059 6422717935
## 15 22 2067771 1930708 2028893 7435599617 19945580764 4737781677
## 16 23 2141332 2018621 2096147 10602832315 21564165489 5647913648
## 17 24 2717762 2840261 2763307 9741654428 18319221993 4840682254
## 18 25 3860999 3691128 3728776 48343981657 21421440486 13206060283
## 19 26 2013072 1459284 1991464 11814088696 306662673204 11334122082
## 20 27 1832752 1422739 1813855 7756597238 160690193355 7396468758
## 21 28 1640965 1709445 1679961 11455772041 8346341303 4296256896
## 22 29 1919034 1886470 1901626 12771124653 11526975983 6049940641
## 23 3 1711264 2555005 1724200 10768567268 704322355600 10590556356
## 24 4 2341407 2451042 2399525 13826138430 13869245317 5602794085
## 25 5 2037005 2099347 2055821 8812445207 13435820352 4529464454
## 26 6 2538739 2346667 2451205 16257549057 23307609312 7627173269
## 27 7 2758099 3520204 2773491 11551388927 577832723711 11298364976
## 28 71 2344845 2302714 2319280 22042402134 22769529831 8505163797
## 29 72 1960940 2069616 2009523 20428995070 24950592172 8414440143
## 30 73 3460877 3346887 3396430 38430312126 33884304401 16854134099
## 31 74 2201707 2032261 2162854 6776932651 27354876799 4646999026
## 32 75 1880316 2110536 1944336 17048243570 50493748831 10171739397
## 33 76 2035745 2044475 2044500 12381994590 8659531696 4943100542
## 34 77 2305963 2417556 2362677 27535106786 30886003406 12865155474
## 35 78 6705015 6209264 6425746 142625267898 120436435604 58669838178
## 36 79 2233947 2035182 2185448 9802575281 37238130601 6854867216
## 37 8 1980133 1696584 1958682 6140168989 77257761384 5670954623
## 38 9 2288451 2451507 2357965 14256484441 24716816108 5515475572
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 4853.123 3982.360 2138.899
## 2 3058.464 2980.850 1329.276
## 3 5654.902 42920.342 4808.911
## 4 4969.290 3206.885 1866.567
## 5 3367.558 15776.638 2290.299
## 6 4646.869 8547.289 2386.683
## 7 14543.652 35177.178 9425.290
## 8 2650.016 4325.527 1321.727
## 9 3423.255 170456.864 3342.147
## 10 2913.504 34495.134 2298.969
## 11 3665.351 110579.185 3530.989
## 12 4351.572 28604.728 3476.766
## 13 5467.735 86123.496 5072.905
## 14 4298.073 138311.424 4136.055
## 15 3579.096 10496.876 2317.231
## 16 4937.018 10583.853 2633.832
## 17 3573.239 6312.330 1745.107
## 18 12433.804 5797.644 3532.765
## 19 5854.803 210063.932 5675.789
## 20 4231.085 113100.208 4077.104
## 21 6961.448 4857.838 2566.647
## 22 6600.296 6067.164 3157.358
## 23 6256.953 275828.004 6107.040
## 24 5907.004 5663.321 2352.740
## 25 4294.701 6377.271 2205.964
## 26 6391.706 9916.699 3070.557
## 27 4180.319 163551.862 4065.359
## 28 9217.863 9870.319 3594.686
## 29 10232.454 12117.726 4174.558
## 30 11074.596 10152.668 4955.987
## 31 3069.758 13470.253 2119.298
## 32 9007.164 23978.600 5287.960
## 33 6071.746 4244.421 2416.221
## 34 11899.931 12732.266 5468.407
## 35 21252.793 19181.910 9028.978
## 36 4377.652 18470.963 3101.269
## 37 3080.686 45616.840 2876.323
## 38 6165.905 10109.823 2372.969 b. f = 90%res_2sc_0.9_0.9## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1626390 1551992 1581193 9693264202 8084784326 3555413958
## 2 10 2438008 2402067 2421305 9504756871 8279525699 3785294846
## 3 11 1445241 1129410 1416190 8799788621 99816421234 7862613247
## 4 12 1627870 1583333 1618532 9672503774 14810376681 3698555441
## 5 13 1671326 1565277 1642161 7644948444 21237897485 4292079236
## 6 14 2676250 2576179 2628186 12325100140 12989749626 5419947765
## 7 15 4849016 4533742 4686664 65408777824 72908418579 27410372117
## 8 16 2740643 2517961 2695991 9702601088 55524955051 6766332331
## 9 17 2133617 2751973 2153915 11929637740 379522918137 11437448499
## 10 18 1751679 2047034 1771727 5646295014 84311825337 5203802785
## 11 19 1914542 2444926 1931651 8824383932 278665476965 8512399490
## 12 2 1626493 1422611 1580976 9436364125 39463570586 6734634517
## 13 20 1933575 1527446 1900755 12865725791 163325884209 11740265951
## 14 21 1505518 2069853 1517607 6736525370 316126071294 6586498623
## 15 22 2058827 1876060 2007870 9274561977 33511750065 6014481169
## 16 23 2113374 1908075 2051763 11541798094 38643750683 6825056119
## 17 24 2652589 2726637 2688958 10624772799 18718500623 5321921403
## 18 25 3782932 3639673 3684043 48971260960 30531781701 17410566252
## 19 26 1939943 1453541 1914022 12445784034 234984667757 11759664380
## 20 27 1852038 1418125 1831759 8701005530 183158379033 8274782503
## 21 28 1642727 1747255 1681993 12248691133 18705151760 5682681416
## 22 29 1920792 1878593 1896831 13986110627 11839758714 6356882908
## 23 3 1814498 2625662 1834447 15864476906 645638858894 15406633959
## 24 4 2300902 2454046 2344366 16418309189 35556152300 7400483379
## 25 5 2012606 2137012 2050166 8185003200 22872836105 4005253562
## 26 6 2549791 2392911 2487680 17096913042 27211177758 8892485131
## 27 7 2674921 3545569 2687254 10457479158 761138458803 10283022757
## 28 71 2271595 2272801 2287224 20900350803 18499167536 9021176759
## 29 72 1903462 2040349 1984453 24067754065 19813306599 9316830674
## 30 73 3357280 3410888 3382421 44415206110 37221399075 18480317171
## 31 74 2195767 2001503 2148334 7903443651 37578881231 5124098195
## 32 75 1912463 2074170 1969691 22469334661 46437698143 9816705356
## 33 76 2050117 1989632 2008814 13133684172 8960921175 5090301788
## 34 77 2248827 2432414 2303401 27302387196 52958278326 15701046817
## 35 78 7181761 6497787 6839894 233666216773 259744872617 112906666347
## 36 79 2179322 2003532 2138407 11602119333 33088807199 7559026193
## 37 8 2001876 1681306 1972617 9011793429 101807265637 8103840097
## 38 9 2281095 2508325 2336995 13493990740 49428084954 8651406626
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 5929.495 5207.431 2216.246
## 2 3877.320 3426.418 1542.799
## 3 5989.550 90178.644 5455.720
## 4 5907.765 9448.227 2223.496
## 5 4553.822 14694.665 2614.370
## 6 4587.888 5062.649 2044.895
## 7 13288.966 16036.640 5712.011
## 8 3529.949 22129.563 2464.697
## 9 5551.114 138303.129 5278.766
## 10 3216.852 41172.534 2932.834
## 11 4606.794 113708.612 4405.128
## 12 5799.813 28091.834 4226.071
## 13 6623.230 107287.635 6140.172
## 14 4482.174 152846.291 4348.453
## 15 4457.748 18254.728 2937.780
## 16 5434.696 21213.358 3290.356
## 17 3974.668 7061.931 1995.110
## 18 12789.420 8400.411 4711.457
## 19 6380.646 162388.254 6107.215
## 20 4703.095 129560.729 4522.291
## 21 7333.772 10732.766 3408.000
## 22 7158.860 6192.413 3294.664
## 23 8674.527 245851.318 8327.512
## 24 7026.182 14618.388 3227.725
## 25 3993.914 10684.427 1958.707
## 26 6646.680 11501.205 3516.713
## 27 3890.381 213720.383 3806.846
## 28 9039.350 8089.771 3857.087
## 29 12593.408 9675.661 4722.900
## 30 13150.935 11021.017 5508.029
## 31 3599.474 18964.343 2349.257
## 32 11539.899 22421.940 5005.121
## 33 6394.532 4508.227 2538.479
## 34 11990.422 21944.621 6844.125
## 35 32197.454 39158.949 16144.638
## 36 5295.369 16528.165 3488.563
## 37 4480.942 61485.070 4086.515
## 38 5842.736 19589.820 3698.567 c. f = 80%res_2sc_0.8_0.8## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1613295 1592433 1598206 11107350123 8679139946 4663807768
## 2 10 2468849 2415596 2435858 11515759602 7478448811 4132475210
## 3 11 1405268 1144375 1365687 10097716352 60423440119 8416343851
## 4 12 1620193 1582832 1605858 12456311701 11072235603 4733853459
## 5 13 1668813 1584199 1616466 8937121944 6119417044 3353826851
## 6 14 2673893 2579573 2628412 14895135498 14831564884 5893180353
## 7 15 4881001 4519874 4716148 82092901976 107104381360 40391098199
## 8 16 2700579 2519860 2639789 11164272662 30620188821 6389074071
## 9 17 2026064 2750283 2042202 11380865574 518199949778 11089274252
## 10 18 1791337 2037516 1828603 7794524478 59264604101 6049138253
## 11 19 1929901 2453867 1949979 10279822544 266486521592 9831248241
## 12 2 1596118 1447619 1542332 11802503170 26843845716 6430800364
## 13 20 1954354 1533640 1919244 13494477156 190624883731 12119377550
## 14 21 1559104 2070385 1586135 10969158423 264887354686 9622791666
## 15 22 2032290 1904096 1969626 10197008683 12822212328 4441111548
## 16 23 2092671 1915266 2047874 13373967296 38083819652 7799022988
## 17 24 2746883 2800358 2767634 14445460222 22988022659 6972313628
## 18 25 3845604 3694678 3754430 78743828382 66824580261 32537314257
## 19 26 2058832 1464991 2031154 16158769646 347983166574 15346907239
## 20 27 1833237 1416024 1811104 9142025026 168896483073 8631031683
## 21 28 1637936 1764262 1713659 15085258258 12186927090 5744797820
## 22 29 1883531 1845320 1862295 16834038897 14973850871 7887785111
## 23 3 1830552 2620892 1860194 22986703161 606965446210 22015035761
## 24 4 2321493 2455610 2344542 20368015458 52671445686 8907345863
## 25 5 2088145 2123686 2087551 12461226158 15548537492 5291064073
## 26 6 2500373 2354786 2434691 21348690376 34558666956 9161053168
## 27 7 2690238 3520241 2705170 12136473112 693698752480 11894536946
## 28 71 2291240 2281929 2263349 29167051634 26106157095 11474036523
## 29 72 1901274 2050209 1932294 27116631953 60331057086 11464605694
## 30 73 3464057 3391080 3418055 51231965117 38225460926 21539664183
## 31 74 2217843 2012933 2163904 11372843873 39999728249 7358219360
## 32 75 1902639 2086858 1956142 20298421920 51562641925 10304421654
## 33 76 1976430 1980957 1987809 15627653154 10346545446 5609502687
## 34 77 2307293 2446674 2348535 39726044858 66188098825 20388868995
## 35 78 7007097 6417657 6660987 231473974917 172859805452 87671828918
## 36 79 2143495 2013226 2085396 12382751982 22613719291 6069444513
## 37 8 2050152 1710251 2012089 12390276059 112892961394 10804018038
## 38 9 2290529 2498570 2353657 19924567987 45446338487 11239799760
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 6880.149 5439.533 2912.516
## 2 4645.364 3081.614 1686.677
## 3 7164.359 53812.815 6148.385
## 4 7626.855 6978.069 2909.282
## 5 5249.913 3790.164 2021.117
## 6 5546.680 5676.561 2207.918
## 7 16695.025 23733.349 8447.071
## 8 4134.017 12033.473 2372.610
## 9 5581.908 188374.988 5397.074
## 10 4335.893 29030.250 3331.324
## 11 5303.914 108675.095 5021.086
## 12 7301.589 18816.764 4031.310
## 13 6842.332 125375.756 6226.054
## 14 6966.020 127979.980 6129.486
## 15 5008.658 6875.001 2248.359
## 16 6379.204 20805.989 3795.680
## 17 5238.941 8170.362 2534.951
## 18 20093.748 17850.425 8466.097
## 19 7817.649 236752.240 7525.209
## 20 4964.274 119706.293 4744.267
## 21 9204.267 6773.411 3365.451
## 22 8728.801 7882.238 4124.955
## 23 12424.484 231140.724 11704.544
## 24 8615.071 21473.499 3834.045
## 25 5938.018 7262.916 2525.580
## 26 8416.481 15400.881 3688.758
## 27 4499.002 196143.184 4385.019
## 28 12429.058 11385.728 5051.366
## 29 13756.880 29322.600 5902.740
## 30 14763.890 11279.967 6303.106
## 31 5127.787 19956.595 3357.043
## 32 10366.359 24628.853 5195.847
## 33 7882.562 5195.567 2798.369
## 34 16957.061 27223.838 8727.081
## 35 32724.228 26432.754 12825.580
## 36 5689.100 11178.374 2786.435
## 37 5991.990 66111.785 5313.105
## 38 8520.071 17860.603 4714.879 d. f = 70%res_2sc_0.7_0.7## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1650667 1576934 1609812 15175811976 13924521116 5465101545
## 2 10 2470594 2368112 2422062 16391186550 15865472506 7148652007
## 3 11 1402838 1142862 1355278 12755456565 65194782283 10150876444
## 4 12 1663858 1610783 1634355 12641405476 9712117450 5083554682
## 5 13 1672665 1551920 1625292 9932265013 17267424858 4737970424
## 6 14 2735300 2603447 2671134 18997751231 22230791204 9124452957
## 7 15 5055832 4614164 4838063 110747193821 144281056538 51547507164
## 8 16 2753737 2549555 2702174 14104875112 49107037454 8928650448
## 9 17 2073268 2764231 2095978 14839380956 474696564202 14232107128
## 10 18 1826460 2079805 1868669 9946880061 57075419345 7933834495
## 11 19 1931728 2479566 1956302 12753946259 295270827083 12081752544
## 12 2 1692038 1444265 1628233 14274626290 54291678165 9592545956
## 13 20 1945349 1543967 1900068 18144111822 149403865066 16016558236
## 14 21 1566762 2101602 1593541 12981356190 284622632937 12050382955
## 15 22 2003796 1905972 1962195 12153172022 15531712243 5400462243
## 16 23 2142764 1991084 2086952 17206709546 34543160393 8431561337
## 17 24 2698532 2732142 2703822 18679149047 20758306301 8188636944
## 18 25 3755091 3655901 3689803 68792963196 48045160438 27369049905
## 19 26 2043069 1460508 2011620 18364383294 334353613961 17347936293
## 20 27 1826372 1454128 1792601 12107134406 136462569054 10883731534
## 21 28 1582503 1727328 1655995 14065439946 15181890017 5910731087
## 22 29 1918438 1871745 1892040 20011597064 17674096656 9336812987
## 23 3 1796635 2676627 1822794 21118273374 769134214463 19992256885
## 24 4 2210458 2474114 2270977 20119513495 78694054189 11718605317
## 25 5 2169663 2154478 2158538 17588819417 18461006317 8111398386
## 26 6 2576101 2406993 2498720 26808747793 37164285671 13247120618
## 27 7 2656107 3584777 2671275 13901506452 855533782399 13657017598
## 28 71 2249232 2309992 2262167 26501745249 34697548095 13198260009
## 29 72 1903112 2074411 1984464 36886680719 40078734050 15619977041
## 30 73 3629204 3484440 3542261 79350280735 69525999412 31401444943
## 31 74 2275883 2029183 2222698 13930015008 70952837280 9425829994
## 32 75 1883769 2105734 1951251 26266705493 59983859661 13178391671
## 33 76 2101824 2073315 2095773 28488217416 29709491551 11467821352
## 34 77 2268690 2455772 2322493 49438536748 79738349972 25461176085
## 35 78 7150378 6545754 6726305 307670389845 149323599618 84888053747
## 36 79 2223103 2035821 2153463 18357841678 32689182371 10596473573
## 37 8 2045992 1675646 2005022 14439548966 137928524092 12659181881
## 38 9 2243524 2478256 2315628 19357483152 44888412335 12016936402
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 9161.777 9115.900 3378.635
## 2 6630.442 6733.932 2933.299
## 3 8922.941 58396.879 7313.139
## 4 7548.265 5985.642 3072.186
## 5 5887.241 11784.172 2915.224
## 6 6946.123 8562.665 3395.154
## 7 21762.390 30975.572 10378.224
## 8 5137.434 19567.063 3262.523
## 9 7087.280 171475.943 6725.731
## 10 5424.362 27456.477 4245.052
## 11 6565.482 118948.693 6145.030
## 12 8439.864 37706.305 5780.405
## 13 9240.482 96584.132 8346.860
## 14 8193.905 135296.623 7500.733
## 15 5991.719 8136.912 2684.511
## 16 8011.472 17505.609 3948.790
## 17 6870.426 7498.245 3027.370
## 18 18086.398 13135.939 7388.975
## 19 8891.680 226934.835 8527.462
## 20 6607.508 93604.152 6041.468
## 21 8810.497 8676.907 3569.424
## 22 10172.365 9127.942 4784.762
## 23 11359.708 287732.360 10616.034
## 24 8956.086 31526.455 5198.851
## 25 8009.820 8565.226 3740.976
## 26 10222.636 15433.109 5145.202
## 27 5224.861 238052.106 5103.160
## 28 11611.235 15091.473 5797.193
## 29 19102.155 19260.155 7878.186
## 30 21592.681 19989.092 8794.440
## 31 6122.890 35108.257 4174.769
## 32 13620.493 28686.928 6793.699
## 33 13352.260 14319.015 5278.243
## 34 21389.295 32370.122 10961.004
## 35 42706.781 22156.273 12180.627
## 36 8132.863 15939.900 4795.149
## 37 7015.729 83610.555 6258.972
## 38 8510.906 17927.150 5134.976 e. f = 60%res_2sc_0.6_0.6## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1644434 1605805 1618934 18118965801 11949649554 7023357798
## 2 10 2497052 2433711 2460268 19968921372 23339639130 8207712008
## 3 11 1425782 1167113 1361554 16945580889 54300054233 12359863223
## 4 12 1633397 1636190 1624188 18200330658 13059936022 6900357702
## 5 13 1684709 1501594 1647103 13449610535 49948551396 8257116505
## 6 14 2577647 2604326 2566234 20734481903 24804178303 8566635076
## 7 15 5162766 4609724 4914921 136687720705 195203629277 70644004977
## 8 16 2653264 2543531 2624575 15074242059 30926311319 7372253708
## 9 17 2026422 2745347 2043948 12361370402 510203644070 12026940202
## 10 18 1776763 2074343 1827156 11713563035 89077939120 8136289114
## 11 19 1924716 2441876 1957438 16098587612 258982760314 14874998952
## 12 2 1618041 1429545 1559941 16079287539 36543020764 10428436563
## 13 20 1941690 1545158 1877893 22134821219 170194974755 17375995797
## 14 21 1601335 2066416 1644278 17450859539 210925108310 14958084835
## 15 22 2051407 1935512 1980222 15662183533 15042714724 5312842522
## 16 23 2159135 2056218 2119482 21673799185 21784820905 10066933445
## 17 24 2666080 2815640 2726869 22847470405 36362056425 11374901102
## 18 25 3922016 3727320 3784415 102652652035 63848901131 36729658032
## 19 26 2079482 1442713 2040417 24342016144 400541052515 22752464849
## 20 27 1819708 1451631 1778261 14766082669 142722543235 12958598648
## 21 28 1681970 1704999 1700319 23663246581 18064493841 9793363477
## 22 29 1986561 1959552 1972701 35329435729 34477916904 17445071390
## 23 3 1756587 2622066 1790883 28370110712 739395732340 26989851465
## 24 4 2289956 2450584 2348317 31656441499 44200926803 14061646085
## 25 5 2041752 2101780 2058044 16239386080 18579405487 7951715442
## 26 6 2558940 2340476 2443869 34246427824 44003570868 13473887637
## 27 7 2724844 3523994 2752653 20953149688 633925300130 20041356949
## 28 71 2439866 2295940 2361983 55507761599 53222062123 23698271544
## 29 72 1881701 2062662 1916537 45115906920 91681934675 22082995874
## 30 73 3560345 3356046 3458741 85280956992 79112569267 35938618588
## 31 74 2165633 2026927 2111914 16153500600 29064198420 7713263740
## 32 75 1883045 2101507 1993528 35637659659 42310447835 11111629932
## 33 76 2067558 2058725 2075195 26510055884 21292614511 11196442799
## 34 77 2158966 2397653 2227691 43430483121 79989820379 20407497781
## 35 78 6715301 6247672 6343277 279328935965 79986207632 57022556543
## 36 79 2200540 2029856 2148497 22046545303 46557431464 11511400617
## 37 8 2031803 1714617 1963834 20106028785 90571967350 15034328217
## 38 9 2287244 2469690 2369439 27034897760 38152314106 10962484595
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 10987.510 7409.177 4321.343
## 2 7950.961 9495.406 3266.791
## 3 11584.265 46587.521 8849.544
## 4 11020.596 7941.098 4224.552
## 5 7928.522 35199.227 4892.796
## 6 7944.986 9374.868 3331.819
## 7 26107.961 42232.922 14040.785
## 8 5704.411 12152.259 2753.093
## 9 6076.537 185593.820 5859.976
## 10 6478.775 43527.668 4489.727
## 11 8346.358 105563.638 7583.469
## 12 9776.326 26171.384 6565.973
## 13 11274.411 109703.398 8962.944
## 14 10715.252 102401.656 9031.857
## 15 7633.088 8118.853 2652.167
## 16 10026.977 10757.133 4771.544
## 17 8532.262 12665.168 4218.175
## 18 25561.956 17211.804 9673.793
## 19 11620.016 276446.800 11057.192
## 20 7970.037 99929.565 7117.768
## 21 13978.789 10534.969 5727.363
## 22 17217.022 17001.029 8552.419
## 23 16035.060 282172.237 14984.547
## 24 13537.451 18072.965 5961.012
## 25 7867.850 8782.378 3844.025
## 26 13216.900 18286.318 5270.354
## 27 7612.013 179262.967 7208.633
## 28 21902.421 23241.216 9646.731
## 29 23003.766 44773.888 11451.464
## 30 23576.879 23703.013 10127.314
## 31 7415.678 14376.677 3553.028
## 32 18537.824 20126.705 5565.072
## 33 12968.660 10369.124 5407.060
## 34 19635.564 33813.136 9224.596
## 35 40773.467 12923.823 9007.746
## 36 9943.901 22873.267 5284.198
## 37 9841.760 53580.257 7555.512
## 38 11538.057 15575.373 4550.398 f. f = 50%res_2sc_0.5_0.5## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1627088 1602180 1609435 24638308146 14507014424 8804381078
## 2 10 2396439 2402405 2395748 20428579382 15136050029 8559609218
## 3 11 1427259 1130374 1361695 24133274715 94438930602 17234969511
## 4 12 1671947 1599374 1634432 25578778963 21525017440 10926713652
## 5 13 1719898 1614243 1665976 21231994533 20896136098 8413788503
## 6 14 2657524 2570184 2607055 36638947394 29547060240 15267892794
## 7 15 4683909 4494463 4603438 128540434459 137430535261 57269999599
## 8 16 2705955 2482806 2630994 27616283786 60744345267 15102987351
## 9 17 2035929 2731415 2072198 19620597435 511028086505 17458566517
## 10 18 1803288 2046903 1868488 16403932993 66197124842 9339715428
## 11 19 1842748 2447135 1881186 20299155434 370088900730 18398638943
## 12 2 1628224 1443629 1552644 21747148144 37142714931 11381039563
## 13 20 1924077 1527688 1827527 34414509622 160210802731 23983148870
## 14 21 1567167 2080915 1611520 19955064841 261518533280 17419926118
## 15 22 2023752 1916978 1973794 20734612315 22282056803 9048744550
## 16 23 2234872 1920026 2148888 31364432960 93542245421 21430335954
## 17 24 2641165 2720813 2669879 26174791664 31578804394 13560419587
## 18 25 3963680 3646357 3787167 118232550672 102525329831 50864763955
## 19 26 2090212 1465403 2045118 28218134857 383651524179 26059381513
## 20 27 1843741 1457139 1794378 18912379958 143207523537 16138432350
## 21 28 1615498 1789639 1677694 28066681986 47440081134 14139534422
## 22 29 2011279 1969307 1989396 41213935916 39052120720 20041408075
## 23 3 1877050 2715030 1929909 40883092253 688194523116 37425282188
## 24 4 2189046 2445423 2244859 31435367194 89021500165 17985575979
## 25 5 2115468 2138222 2108826 30328267446 27746327553 10254765417
## 26 6 2582670 2388805 2481174 45701624326 38185179771 18591833632
## 27 7 2600831 3501617 2630375 22995147323 827982276227 21695162150
## 28 71 2406074 2287452 2319698 61282292406 46482701100 21686985644
## 29 72 1897002 2065856 1951530 51803711263 68335984519 19879388591
## 30 73 3350928 3482212 3369799 85592044177 124838523775 38587839856
## 31 74 2220538 2020928 2128199 25636944324 34499328043 11079474226
## 32 75 1808320 2090615 1884344 48450556597 99080828241 25535510698
## 33 76 2156520 2037983 2076553 37902423067 32756279203 15412357582
## 34 77 2351942 2475444 2445714 82659687369 38211263810 22729998439
## 35 78 7251027 6506432 6807994 374456497721 271288423748 138990896286
## 36 79 2200908 2026501 2105509 32482045469 30911187531 13787752259
## 37 8 2019609 1690557 1955817 20802512192 107857372992 16381659190
## 38 9 2298544 2509675 2389346 28885719260 41847808991 11715982424
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 14457.782 8817.833 5273.225
## 2 8416.485 6271.273 3548.823
## 3 16310.666 87652.857 12093.610
## 4 15082.653 13625.302 6641.266
## 5 11869.205 12774.859 4724.378
## 6 13722.912 11395.807 5786.525
## 7 27403.815 30599.118 12316.632
## 8 10181.812 24574.086 5637.571
## 9 9311.775 186392.722 8199.927
## 10 9031.299 32233.167 5136.906
## 11 10867.528 150463.246 9687.997
## 12 13180.151 25403.170 7039.685
## 13 17533.459 106229.099 12504.998
## 14 12577.104 124331.627 10703.200
## 15 10137.809 12004.853 4531.523
## 16 13964.157 50098.023 9861.951
## 17 9935.803 11549.942 5103.756
## 18 29574.466 27816.177 13206.958
## 19 13348.111 266660.135 12604.572
## 20 10105.296 98654.375 8851.030
## 21 17016.015 26337.994 8403.667
## 22 19826.586 19170.568 9742.615
## 23 21433.947 253629.514 19129.818
## 24 13712.120 36619.341 7914.289
## 25 13958.431 12930.582 4811.340
## 26 17133.674 16003.722 7359.350
## 27 8711.219 235191.043 8156.458
## 28 24453.570 19929.174 8880.629
## 29 26456.201 33364.133 10292.207
## 30 25242.138 35698.696 11596.253
## 31 11565.156 17236.513 5133.015
## 32 25702.460 46891.403 13330.075
## 33 17444.269 16260.939 7446.956
## 34 34314.781 15060.064 9089.761
## 35 50260.417 40614.817 19722.037
## 36 14660.930 15160.099 6413.565
## 37 9971.189 63697.166 8027.529
## 38 12336.271 16542.885 4939.157 g. f = 40%res_2sc_0.4_0.4## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1624062 1610125 1611581 42335533016 38575187434 19755534607
## 2 10 2479069 2461803 2468638 34426746416 24803179498 12909685619
## 3 11 1499455 1161265 1408674 42402911093 115590232145 30118997835
## 4 12 1679679 1621295 1621835 33304424110 48486875152 15137950051
## 5 13 1784151 1643738 1715999 34501616632 41572718457 16885599211
## 6 14 2705845 2648689 2672490 42466700434 32020520623 16560914128
## 7 15 4923058 4579037 4742821 229403573446 180410612030 81288746326
## 8 16 2787417 2501716 2706385 36217997625 97026950553 21455303873
## 9 17 2136820 2751764 2221458 39959962745 375122245910 26220287228
## 10 18 1858921 2053757 1946242 30244348831 40316236114 14386382772
## 11 19 1930788 2451055 2014183 36139125533 260639741269 26731573484
## 12 2 1647160 1446888 1541730 33783322878 39489449200 10988296871
## 13 20 1916950 1524705 1807547 54573260932 200349869722 36739170081
## 14 21 1586760 2079890 1640271 26849682209 243635035915 21306380960
## 15 22 2080393 1960247 2016956 36035606567 37290083091 15316833579
## 16 23 2160809 1982158 2092511 34084617398 49809794180 18408992685
## 17 24 2635495 2836815 2732224 48129982434 48874687701 20356113145
## 18 25 4086344 3748657 3874671 201716689324 128934906803 72626477890
## 19 26 2085978 1456431 2021362 41339732043 370611594153 36861639716
## 20 27 1928809 1464905 1866983 28015910883 213487360237 23949832352
## 21 28 1663028 1785490 1697620 37586489296 102505536691 22461523953
## 22 29 1780064 1751559 1764546 35985423406 34147369849 17484572412
## 23 3 1877959 2710853 1957209 57989860716 661629043959 48060000933
## 24 4 2203916 2451725 2249148 48921943958 133650555150 27377613260
## 25 5 2042554 2137493 2041138 42670728339 55628243767 16931886520
## 26 6 2475236 2388919 2433375 54810855534 47364338723 24544092797
## 27 7 2652497 3575790 2693490 33000681427 865756924496 30525240800
## 28 71 2186213 2272527 2155727 70247098472 132223501976 31638785268
## 29 72 1902369 2042838 1859923 104916372028 150507981349 40111908501
## 30 73 3420300 3456732 3423730 122581130222 115353484328 55450415609
## 31 74 2242295 2002173 2147694 41213679501 57974230102 20864379770
## 32 75 1767592 2073931 1858178 40762320245 91245925385 20012559397
## 33 76 2129257 2017940 2098226 59848770554 68228066059 29992881881
## 34 77 2319782 2440900 2322747 108731049788 159131134040 52342924149
## 35 78 7156354 6497941 6690391 642665400534 266480064799 185826466150
## 36 79 2131160 2004951 2056905 50042038306 40398106731 21161751857
## 37 8 2016835 1720347 1939391 31080297793 95051074683 20729081487
## 38 9 2331335 2519657 2379058 42902457571 88268784074 21403064284
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 25611.25 23251.36 11918.847
## 2 13880.41 10239.23 5263.390
## 3 26632.60 104237.66 20318.432
## 4 18952.95 30573.01 9096.893
## 5 18817.85 25530.07 9614.057
## 6 15633.40 11813.90 6082.568
## 7 44987.05 40041.04 16823.824
## 8 12918.46 39053.09 7782.679
## 9 18005.19 136459.12 11866.304
## 10 15922.60 19427.47 7323.641
## 11 18532.96 105566.65 13398.008
## 12 20008.90 26959.48 6806.972
## 13 26923.32 134934.86 18345.240
## 14 16714.03 117094.81 13111.099
## 15 17076.81 19451.19 7393.814
## 16 15632.94 24860.19 8629.752
## 17 18075.89 16952.12 7407.987
## 18 47122.79 33828.87 18018.032
## 19 19685.80 257672.31 18143.448
## 20 14149.81 146171.24 12408.517
## 21 21295.24 54973.06 12336.470
## 22 19545.06 18790.37 9563.162
## 23 29521.86 245950.59 23587.581
## 24 21265.55 54514.17 12005.804
## 25 19466.77 26067.89 8128.306
## 26 21595.22 19739.02 9868.543
## 27 12302.04 239160.31 11241.075
## 28 29399.01 57133.84 13991.892
## 29 48432.77 73324.76 20689.065
## 30 35700.78 33568.32 16291.207
## 31 17906.39 28647.47 9286.591
## 32 22342.23 44756.36 10759.344
## 33 27495.20 34143.83 13974.629
## 34 43954.47 63731.01 21799.145
## 35 88278.69 40137.10 27206.718
## 36 22636.78 19751.57 9917.559
## 37 15207.32 55976.80 10437.168
## 38 17354.92 32964.62 8645.484 h. f = 30%res_2sc_0.3_0.3## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1624639 1579174 1595911 44530005096 22620075081 12878936450
## 2 10 2387273 2418662 2327284 61798258333 59329236866 20565765585
## 3 11 1403232 1152388 1305709 35891732608 57424925164 20448136707
## 4 12 1739856 1639604 1668772 58517903647 77394382751 26076208851
## 5 13 1611529 1491111 1554023 30236250427 28746417901 10778465723
## 6 14 2548011 2614995 2576464 66906337608 59876134199 28340174082
## 7 15 4948970 4618411 4791534 273771467164 267257629182 129535758338
## 8 16 2751068 2535160 2658500 54106458913 67769184936 24026905210
## 9 17 2039624 2736983 2115059 41675011607 498590919507 32111167667
## 10 18 1777182 2077274 1881535 37424206177 89813169383 15641393441
## 11 19 1830241 2423691 1891020 49658855560 421637488390 31997188756
## 12 2 1532007 1395422 1478495 41127863023 53431743363 19834093191
## 13 20 1981689 1544907 1850715 88611235021 204717933679 58941201237
## 14 21 1562061 2045099 1627769 33533220943 259215784662 22857081155
## 15 22 2040182 1934367 1973427 48732459235 40539013416 21090289402
## 16 23 2294484 2066911 2188374 77115449679 118971836949 37132905059
## 17 24 2688075 2816120 2725171 104334055221 104285023341 45218240664
## 18 25 4063702 3740971 3848599 294252991209 166630031680 88487692278
## 19 26 2025403 1415611 1936662 57340378222 389525933343 47613522161
## 20 27 1733731 1446834 1671786 23537765928 88736707388 17200097816
## 21 28 1622626 1692880 1667058 70943377460 57178861639 30955221542
## 22 29 1987797 1966309 1976558 69293448003 68056775444 34325094153
## 23 3 1874390 2627257 2029436 94538842659 527640195567 37551093952
## 24 4 2191937 2435005 2282276 68889164498 96532899439 32385669127
## 25 5 2183429 2084299 2098777 77498505452 47406506673 27433039237
## 26 6 2551732 2350061 2449114 121238635126 108618695290 53882005157
## 27 7 2740455 3492120 2816715 57508983492 577014048274 48272188890
## 28 71 2268367 2285092 2263603 120700404326 124455298833 57817580597
## 29 72 1880086 2054281 1906332 117976109176 189513850083 61576490796
## 30 73 3353908 3337212 3336897 210601114063 146771851474 84793461493
## 31 74 2210732 2018395 2101613 74348280508 48878695863 24749298663
## 32 75 2062440 2092900 2038330 180749935651 151967422533 74173823464
## 33 76 2136951 2064367 2065322 114493582725 65970087696 38186858644
## 34 77 2256702 2377777 2306555 185756627976 161085087756 78441009369
## 35 78 6734125 6181152 6375989 841366598898 514998726048 312269516759
## 36 79 2019687 2022049 1993514 51070234588 38246575176 17913763280
## 37 8 1950235 1691770 1838626 50301768298 75637719867 25011071223
## 38 9 2258358 2400294 2282672 67516190632 97130329143 33182962026
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 26643.60 14136.32 7941.412
## 2 24335.78 24094.35 8705.503
## 3 23352.02 48477.60 14412.000
## 4 31701.48 48472.95 14670.883
## 5 17386.87 21675.13 6853.371
## 6 25917.50 22584.28 10885.962
## 7 54507.51 56936.90 26425.360
## 8 19342.62 26670.11 8830.204
## 9 20181.47 179932.28 15357.522
## 10 20650.98 42882.58 8432.782
## 11 25212.06 173213.49 16668.548
## 12 26162.32 40491.11 13039.093
## 13 40889.32 133427.08 29171.203
## 14 20502.00 125322.42 13768.049
## 15 23431.12 20880.34 10554.197
## 16 32839.16 58402.28 16753.945
## 17 34700.73 35936.41 15692.278
## 18 68254.01 44373.10 22318.769
## 19 27620.59 274674.65 23616.387
## 20 13368.59 64546.27 10149.494
## 21 42209.72 32753.74 17891.935
## 22 33297.22 32935.92 16553.690
## 23 46150.30 200435.78 18866.753
## 24 29637.53 39660.81 13784.300
## 25 34034.69 22412.75 12831.580
## 26 46341.69 45692.00 21492.202
## 27 20921.99 162275.98 17201.976
## 28 51112.28 54205.33 25100.551
## 29 50445.52 91591.37 28116.845
## 30 58420.65 42564.20 24244.831
## 31 32883.58 24004.41 11411.052
## 32 80835.47 72223.55 35397.563
## 33 49188.94 30586.22 17642.209
## 34 80531.14 67292.40 33743.068
## 35 123838.68 82384.19 48454.348
## 36 24537.71 18432.58 8928.474
## 37 24958.20 46326.01 13054.802
## 38 25259.95 38503.63 12716.513 i. f = 20%res_2sc_0.2_0.2## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1691223 1574218 1603619 7.577977e+10 3.486898e+10 21516248873
## 2 10 2293102 2301545 2291751 1.084384e+11 1.059325e+11 52589996567
## 3 11 1346176 1158896 1247246 5.313367e+10 4.730405e+10 23197976175
## 4 12 1492624 1572130 1523598 7.399138e+10 5.662975e+10 27957631689
## 5 13 1734708 1636610 1664041 1.236993e+11 1.016917e+11 51766751547
## 6 14 2868064 2565653 2691508 2.055454e+11 1.461616e+11 83724059101
## 7 15 4848865 4427208 4518063 4.542614e+11 1.805892e+11 115200416222
## 8 16 2643377 2503426 2556438 1.862958e+11 1.287279e+11 64199083973
## 9 17 2107456 2698580 2181965 1.210220e+11 5.703685e+11 68605251756
## 10 18 1931667 2066745 1959811 1.255895e+11 1.631769e+11 60929789996
## 11 19 1866043 2419891 1979566 1.231638e+11 3.691553e+11 64153886158
## 12 2 1737344 1476247 1583907 1.299916e+11 8.400595e+10 48140500901
## 13 20 1841945 1531704 1655416 1.369415e+11 8.742261e+10 51481328363
## 14 21 1494139 2057788 1663289 8.490444e+10 2.463458e+11 40782116928
## 15 22 2068179 2028447 2046412 1.097180e+11 9.674786e+10 47596481197
## 16 23 2004740 2020789 2020338 8.650822e+10 1.310292e+11 44715671038
## 17 24 2724481 2785793 2762568 2.659311e+11 2.146186e+11 115421138765
## 18 25 3689118 3582622 3613083 3.147385e+11 2.408683e+11 132544532808
## 19 26 1978657 1488345 1814156 8.841907e+10 2.191686e+11 55365597526
## 20 27 1709439 1456432 1620009 5.150174e+10 1.055725e+11 26782473651
## 21 28 1559796 1733618 1622170 6.587145e+10 1.140149e+11 26826637558
## 22 29 1942713 1909214 1925588 1.531271e+11 1.514382e+11 76129177911
## 23 3 1875462 2641669 2041306 1.134837e+11 5.651923e+11 68089034738
## 24 4 2177980 2411545 2242176 7.379274e+10 1.284122e+11 37506716754
## 25 5 2170758 2095684 2085917 1.415214e+11 8.929052e+10 49673152917
## 26 6 2569195 2371119 2409040 2.539411e+11 1.858272e+11 99739580703
## 27 7 2673686 3462306 2820187 1.313935e+11 6.747393e+11 87833879471
## 28 71 2220466 2224050 2141850 2.530707e+11 1.908492e+11 95760210067
## 29 72 1920737 2005012 1898844 1.792261e+11 2.094025e+11 79320572591
## 30 73 3692267 3415619 3470288 6.130850e+11 3.904150e+11 221990946270
## 31 74 2212960 1954435 2081900 9.688368e+10 1.355539e+11 40700476732
## 32 75 1837394 2027898 1870922 3.067533e+11 2.987428e+11 117791938163
## 33 76 2113056 2037668 2002137 1.823575e+11 1.986296e+11 80021405206
## 34 77 2137946 2396144 2168164 3.191715e+11 4.317586e+11 161649404618
## 35 78 6633325 6224540 6393108 1.247178e+12 1.027193e+12 559594025162
## 36 79 2085323 1965400 1965441 1.625899e+11 8.972342e+10 48368992349
## 37 8 2075568 1644202 1886123 1.245285e+11 2.249349e+11 56789909979
## 38 9 2349878 2446664 2415553 1.173966e+11 9.722789e+10 52018154754
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 43092.06 22250.13 13354.51
## 2 46162.91 45991.34 22621.91
## 3 38532.87 42746.08 18260.58
## 4 45609.69 33163.67 16685.92
## 5 63959.34 57381.49 27850.37
## 6 71223.39 58016.13 31358.86
## 7 90619.58 40220.27 24836.97
## 8 63289.84 50492.88 23336.39
## 9 50421.51 208319.79 29324.25
## 10 61056.36 79932.04 30053.66
## 11 59179.86 154113.42 29834.45
## 12 68253.85 56672.19 28979.78
## 13 67830.51 59835.71 30562.07
## 14 53440.14 123833.95 24959.11
## 15 47731.03 46010.47 21157.45
## 16 40889.66 72418.05 21371.74
## 17 84096.60 71709.94 37457.18
## 18 83615.43 66801.94 36165.85
## 19 43275.07 147987.68 28634.91
## 20 28761.86 81829.44 15754.42
## 21 38344.24 64918.34 15684.36
## 22 69953.13 70020.78 34991.35
## 23 56806.68 206369.50 32288.76
## 24 33240.73 52030.03 16748.90
## 25 61760.76 41721.29 22905.14
## 26 96740.68 78552.46 41101.49
## 27 47513.13 190020.07 30927.86
## 28 96099.42 85360.36 41632.20
## 29 78586.44 102638.30 37465.96
## 30 143513.88 108543.83 58044.15
## 31 42246.58 73298.54 18266.08
## 32 130948.83 143198.51 56269.55
## 33 79198.12 92671.15 37610.34
## 34 131284.18 179668.33 70283.82
## 35 185418.86 161738.97 85990.31
## 36 71729.67 45490.86 23919.05
## 37 54403.50 137528.02 26232.65
## 38 47557.49 37911.88 20482.16MEnggunakan rumus MSE synthetic yang perbaikan
res2_2sc2_1_1## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1637409 1593637 1637409 0 955703889270 0
## 2 10 2452226 2411973 2452226 0 955693693914 0
## 3 11 1432925 1146047 1432925 0 958475729411 0
## 4 12 1628141 1604098 1628141 0 955657754859 0
## 5 13 1660654 1539506 1660654 0 956143923359 0
## 6 14 2677707 2613785 2677707 0 955778716546 0
## 7 15 5059790 4628255 5059790 0 962059273603 0
## 8 16 2702258 2538366 2702258 0 956564054318 0
## 9 17 2063216 2755581 2063216 0 972167779798 0
## 10 18 1802921 2055582 1802921 0 957839116898 0
## 11 19 1893080 2458794 1893080 0 966673452782 0
## 12 2 1631372 1434795 1631372 0 956970329890 0
## 13 20 1933963 1528077 1933963 0 961318643432 0
## 14 21 1543883 2075747 1543883 0 965392305402 0
## 15 22 2024290 1912770 2024290 0 956066666522 0
## 16 23 2145480 1996068 2145480 0 956407607501 0
## 17 24 2721467 2809820 2721467 0 955907001210 0
## 18 25 3923463 3735418 3923463 0 956857168758 0
## 19 26 2049321 1442285 2049321 0 968344465890 0
## 20 27 1850200 1431308 1850200 0 961688550294 0
## 21 28 1632444 1718061 1632444 0 955890588991 0
## 22 29 1945828 1910735 1945828 0 955680287983 0
## 23 3 1794553 2646763 1794553 0 980681363141 0
## 24 4 2297582 2461884 2297582 0 956568683966 0
## 25 5 2083579 2122123 2083579 0 955689051423 0
## 26 6 2526531 2357735 2526531 0 956620312116 0
## 27 7 2716794 3573111 2716794 0 980923248393 0
## 28 71 2297609 2288027 2297609 0 955640987735 0
## 29 72 1872873 2049375 1872873 0 956712061744 0
## 30 73 3484448 3408041 3484448 0 955839136750 0
## 31 74 2197198 2010784 2197198 0 956836105612 0
## 32 75 1867766 2087481 1867766 0 957302461609 0
## 33 76 2088732 2024320 2088732 0 955780884996 0
## 34 77 2275392 2420130 2275392 0 956360202086 0
## 35 78 7020112 6437879 7020112 0 967327314729 0
## 36 79 2178314 2013176 2178314 0 956578186333 0
## 37 8 2025630 1700059 2025630 0 959292868838 0
## 38 9 2279804 2477361 2279804 0 956983634820 0
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 0 599699.8 0
## 2 0 396229.0 0
## 3 0 836332.2 0
## 4 0 595760.3 0
## 5 0 621072.0 0
## 6 0 365668.4 0
## 7 0 207866.5 0
## 8 0 376842.5 0
## 9 0 352799.5 0
## 10 0 465969.9 0
## 11 0 393149.4 0
## 12 0 666973.7 0
## 13 0 629103.6 0
## 14 0 465081.8 0
## 15 0 499833.5 0
## 16 0 479145.8 0
## 17 0 340202.2 0
## 18 0 256158.0 0
## 19 0 671396.2 0
## 20 0 671895.1 0
## 21 0 556377.5 0
## 22 0 500163.7 0
## 23 0 370521.0 0
## 24 0 388551.5 0
## 25 0 450345.7 0
## 26 0 405737.1 0
## 27 0 274529.2 0
## 28 0 417670.4 0
## 29 0 466831.2 0
## 30 0 280465.9 0
## 31 0 475852.3 0
## 32 0 458592.1 0
## 33 0 472149.1 0
## 34 0 395169.0 0
## 35 0 150255.6 0
## 36 0 475158.7 0
## 37 0 564270.2 0
## 38 0 386291.6 0 b. f = 90%res2_2sc2_0.9_0.9## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1636986 1592285 1605817 2.095245e+12 939501537163 643510665115
## 2 10 2448930 2418047 2422009 6.140818e+12 939462547113 813948391981
## 3 11 1421352 1144085 1226305 2.252747e+12 942068621255 659452050664
## 4 12 1625446 1606477 1610865 2.489491e+12 939503056796 680656957020
## 5 13 1663225 1527483 1561391 2.848169e+12 940059147445 704702440040
## 6 14 2680657 2610325 2620477 5.579538e+12 939621730716 803176108336
## 7 15 5055180 4612941 4625779 3.156969e+13 946271924793 918664056694
## 8 16 2684384 2536543 2565974 3.788269e+12 940179118017 752517095146
## 9 17 2080308 2757703 2602674 3.231330e+12 955247032782 736848360721
## 10 18 1809275 2061249 1990655 2.464443e+12 941626994904 679583084160
## 11 19 1911671 2451197 2308891 2.682271e+12 949465974927 699424972517
## 12 2 1644132 1438254 1496551 2.387403e+12 940866763360 673250667903
## 13 20 1941041 1537992 1617432 3.865736e+12 945069892624 756605159480
## 14 21 1536454 2074466 1922837 2.425896e+12 949398817123 681850591435
## 15 22 2026726 1909133 1937094 2.969411e+12 939933292407 712451370127
## 16 23 2139697 1995460 2026329 3.467823e+12 940162285043 738705729621
## 17 24 2710160 2799417 2786535 5.668209e+12 939741428711 804760403369
## 18 25 3903980 3722677 3735154 1.262227e+13 940627384862 874708646552
## 19 26 2034148 1446449 1619942 2.283182e+12 951340872126 669128429781
## 20 27 1844046 1433882 1559125 2.159767e+12 945232808953 655844969874
## 21 28 1631040 1718602 1694491 2.557203e+12 939729767540 685922472615
## 22 29 1929624 1894598 1901683 3.802542e+12 939439618136 749785904963
## 23 3 1783932 2630155 2473691 4.280966e+12 964141097420 786129658556
## 24 4 2270448 2462961 2431266 5.282075e+12 940853345086 793503956869
## 25 5 2099107 2124879 2119922 4.026698e+12 939432221771 760542518144
## 26 6 2533881 2365748 2383069 8.171273e+12 940408913516 842628437618
## 27 7 2700917 3557805 3434911 5.786313e+12 964776266568 826505204196
## 28 71 2288659 2280277 2281355 4.940105e+12 939429768484 788369170039
## 29 72 1871972 2042432 1995729 2.781307e+12 940480845755 693393042171
## 30 73 3421763 3370910 3374841 1.011121e+13 939651417946 859270816826
## 31 74 2198368 2006751 2080333 1.503782e+12 940725816065 577162985111
## 32 75 1880351 2082021 2019916 2.233723e+12 940884249933 657722482385
## 33 76 2087826 2010660 2030568 2.522888e+12 939688747389 680715846932
## 34 77 2255208 2403001 2374909 4.372781e+12 940304530298 770237334056
## 35 78 6965927 6386184 6392925 8.087929e+13 951061150201 939981083270
## 36 79 2169969 2007866 2063647 1.747177e+12 940411818888 607495255119
## 37 8 2022379 1706131 1777624 3.245834e+12 942869432604 729327194838
## 38 9 2288095 2470224 2445293 6.118439e+12 940577538672 813331835413
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 1276953.8 590324.7 400609.5
## 2 2505777.4 388616.7 336116.0
## 3 1580943.4 824187.1 538039.5
## 4 1529183.7 585094.5 422615.1
## 5 1710044.0 615894.0 451353.8
## 6 2081490.0 359918.5 306467.6
## 7 6240168.4 205078.1 198540.5
## 8 1410527.4 370677.8 293252.2
## 9 1552709.6 346380.7 283081.8
## 10 1360170.6 456931.1 341386.5
## 11 1402125.9 387326.7 302886.9
## 12 1451074.9 654406.0 449887.6
## 13 1987388.7 614899.4 467957.1
## 14 1578426.3 457650.0 354588.4
## 15 1463386.2 492618.5 367887.2
## 16 1620064.9 471537.3 364698.9
## 17 2090763.7 335813.2 288879.2
## 18 3229057.7 252671.9 234167.0
## 19 1121090.2 658112.7 413094.3
## 20 1169786.0 659558.7 420620.5
## 21 1566176.4 547195.3 404960.8
## 22 1965360.0 496522.6 394442.2
## 23 2398289.3 366529.7 317770.1
## 24 2317903.7 381986.6 326362.4
## 25 1917879.5 442108.5 358747.5
## 26 3220566.8 397569.4 353600.4
## 27 2141335.3 271077.5 240534.6
## 28 2157065.1 412059.4 345619.6
## 29 1479264.4 460645.3 347335.7
## 30 2951257.3 278686.1 254548.2
## 31 683655.4 468980.5 277469.5
## 32 1183724.0 452074.7 325492.5
## 33 1204412.5 467483.2 335124.0
## 34 1932122.6 391418.7 324333.4
## 35 11605777.3 148861.8 146972.4
## 36 803206.2 468556.3 294365.1
## 37 1604556.7 552887.0 410331.4
## 38 2670595.0 380820.2 332643.3 c. f = 80%res2_2sc2_0.8_0.8## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1625221 1582494 1591626 3.324757e+12 9.900212e+11 747249861282
## 2 10 2448295 2419685 2422050 1.118283e+13 9.899656e+11 908923529155
## 3 11 1448698 1146046 1201177 4.488851e+12 9.930908e+11 810941372164
## 4 12 1632550 1606001 1610799 4.405111e+12 9.899696e+11 804884839182
## 5 13 1669021 1555428 1573418 5.245564e+12 9.903802e+11 828935384346
## 6 14 2722098 2644459 2650625 1.117575e+13 9.902662e+11 908833096804
## 7 15 5089721 4670848 4677488 6.186224e+13 9.961176e+11 980065543804
## 8 16 2720981 2529620 2553870 6.896086e+12 9.911813e+11 865203248392
## 9 17 2055816 2754683 2645222 5.697145e+12 1.006868e+12 850047667690
## 10 18 1808414 2048147 2004435 4.566822e+12 9.919135e+11 811813506708
## 11 19 1903053 2454138 2357241 4.810766e+12 1.000456e+12 825111266700
## 12 2 1604995 1420621 1456379 4.155904e+12 9.910886e+11 797004609678
## 13 20 1914344 1512889 1563729 6.955817e+12 9.955065e+11 866769584709
## 14 21 1520171 2072041 1959106 4.083676e+12 1.000459e+12 797633369770
## 15 22 2006261 1904671 1919928 5.282760e+12 9.903430e+11 829657009444
## 16 23 2152017 1998388 2018527 6.518963e+12 9.908096e+11 859182306868
## 17 24 2715240 2815421 2805620 1.029070e+13 9.902738e+11 897143873642
## 18 25 3966398 3761964 3769981 2.444537e+13 9.913540e+11 952189053582
## 19 26 2055269 1425197 1543220 4.373192e+12 1.003657e+12 814608202617
## 20 27 1835848 1431458 1521817 3.569666e+12 9.955668e+11 771095482303
## 21 28 1620923 1710788 1694024 4.453330e+12 9.902220e+11 807174970356
## 22 29 1943464 1911379 1915475 7.221622e+12 9.899241e+11 864000760059
## 23 3 1842051 2670374 2571705 8.109869e+12 1.013642e+12 894728884753
## 24 4 2297315 2459346 2443368 1.078545e+13 9.909759e+11 902079284795
## 25 5 2089675 2126603 2120018 6.712039e+12 9.900910e+11 856802615039
## 26 6 2528929 2368646 2379090 1.409092e+13 9.908406e+11 924244361545
## 27 7 2737518 3617489 3540056 1.070168e+13 1.016763e+12 927580312881
## 28 71 2300918 2272372 2274173 8.865780e+12 9.901215e+11 886199756724
## 29 72 1833508 2030200 1990911 4.686714e+12 9.914610e+11 805616638993
## 30 73 3484788 3420673 3423518 1.936067e+13 9.902466e+11 941267639221
## 31 74 2181289 1991418 2043753 2.602018e+12 9.912002e+11 713897765377
## 32 75 1821274 2068658 2008762 3.444609e+12 9.922193e+11 763371578575
## 33 76 2086997 2012329 2024476 4.360058e+12 9.902597e+11 800368937225
## 34 77 2258396 2403600 2385173 7.738031e+12 9.910241e+11 875191774659
## 35 78 7099338 6509233 6513244 1.463126e+14 1.002110e+12 995267493958
## 36 79 2150902 1993512 2028968 3.117607e+12 9.910122e+11 743341804649
## 37 8 2027749 1692762 1740687 6.022789e+12 9.938101e+11 850404035122
## 38 9 2288287 2474153 2459107 1.173991e+13 9.911885e+11 911923113453
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 2032685 626087.8 469353.4
## 2 4566724 409267.7 375364.8
## 3 3092752 867690.0 675467.5
## 4 2696601 616678.7 499781.4
## 5 3131391 637275.1 526994.8
## 6 4100788 374461.7 342840.1
## 7 12139687 213179.9 209445.7
## 8 2533812 391947.6 338843.3
## 9 2759983 365563.5 321406.5
## 10 2521602 484580.9 405228.6
## 11 2522982 407575.1 349937.1
## 12 2586338 698444.6 547540.3
## 13 3616762 659017.4 554729.8
## 14 2674211 482797.0 406927.8
## 15 2624289 520410.1 432364.7
## 16 3025185 496486.5 426083.0
## 17 3777222 351853.6 319793.9
## 18 6157205 263479.5 252523.7
## 19 2125012 705016.8 528177.2
## 20 1941429 695835.2 506866.6
## 21 2741535 579276.0 476639.1
## 22 3692732 519505.5 451711.7
## 23 4374149 379387.3 347536.3
## 24 4661471 402968.9 369104.7
## 25 3194154 465555.0 404053.5
## 26 5557894 418530.4 388637.3
## 27 3905090 280929.4 261934.4
## 28 3836700 435701.1 389626.0
## 29 2535705 488420.6 404382.4
## 30 5543322 289318.0 274770.2
## 31 1191545 497838.3 349363.5
## 32 1877458 479712.5 379777.5
## 33 2078207 492096.2 395178.1
## 34 3413150 412288.1 366910.3
## 35 20582448 153827.5 152682.7
## 36 1440611 497214.9 366312.9
## 37 2963341 587486.4 488644.6
## 38 5120505 400707.4 370888.7 d. f = 70%res2_2sc2_0.7_0.7## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1616307 1578886 1583843 4.456075e+12 915428980751 743958223930
## 2 10 2418382 2405519 2405951 1.495908e+13 915337581447 861684627967
## 3 11 1449058 1169190 1209808 5.753112e+12 918199859387 779087519389
## 4 12 1638281 1600037 1604724 6.218126e+12 915338075814 795711041545
## 5 13 1651059 1537062 1550116 6.845293e+12 915742872086 801008517423
## 6 14 2635718 2594061 2596059 1.436788e+13 915538563580 858220839704
## 7 15 4956678 4570370 4574847 7.702523e+13 920975030152 909648300877
## 8 16 2689935 2548547 2561433 8.829726e+12 916150912358 828099183793
## 9 17 2116274 2749459 2686778 8.831415e+12 929263087089 838098271974
## 10 18 1812652 2048644 2018374 6.352914e+12 917178103484 799627497004
## 11 19 1884646 2449682 2385580 7.330976e+12 926255206987 821515006920
## 12 2 1602153 1433267 1456297 5.594439e+12 916386860208 775709771733
## 13 20 1904186 1540388 1569605 1.039723e+13 920440479967 839798212712
## 14 21 1542188 2062208 1986324 5.825114e+12 924760986111 792849869298
## 15 22 2031159 1896481 1910457 7.387406e+12 915989967952 810753126239
## 16 23 2139800 1971372 1987674 8.400902e+12 916394611044 823833320812
## 17 24 2734362 2815637 2810414 1.378755e+13 915545702606 856571750712
## 18 25 3883256 3707509 3712766 2.826353e+13 916586371492 886392524354
## 19 26 2010666 1451640 1530751 5.658565e+12 926151346266 793467712454
## 20 27 1810387 1412662 1480694 4.600338e+12 920900999802 753531783879
## 21 28 1634796 1714692 1704484 6.505102e+12 915493837152 801516220524
## 22 29 1967574 1933639 1936348 1.059863e+13 915242383712 839781153446
## 23 3 1783744 2579527 2510546 1.064550e+13 937262181465 857529503483
## 24 4 2309992 2454956 2443414 1.585419e+13 916563979345 859545984487
## 25 5 2087482 2109798 2106771 1.085794e+13 915477578372 840359967976
## 26 6 2493951 2339984 2346828 1.882024e+13 916289905163 871745394533
## 27 7 2745183 3539126 3492524 1.521725e+13 937077461345 882112911354
## 28 71 2345345 2281799 2285016 1.245700e+13 915752228290 849267434779
## 29 72 1937661 2044848 2029065 7.676137e+12 916108083923 806942631494
## 30 73 3367689 3338240 3338419 2.451089e+13 916190073046 881969820716
## 31 74 2166577 2010224 2041592 3.593633e+12 916141913995 726857140387
## 32 75 1860847 2087107 2054943 5.748848e+12 917112309451 788161578167
## 33 76 2032231 1994677 1998412 5.312315e+12 915585901017 767105371870
## 34 77 2282816 2405977 2393619 1.120625e+13 916215022740 842287091538
## 35 78 6878159 6308461 6311451 1.787230e+14 926669693655 921804256123
## 36 79 2151334 2010293 2034663 4.157226e+12 916074345011 744737786718
## 37 8 2023873 1705086 1738444 7.942447e+12 918880593135 818471021287
## 38 9 2271166 2466509 2454641 1.488440e+13 916659104156 860356970088
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 2728019 580577.6 469832.4
## 2 6178160 380759.1 358370.3
## 3 3914544 787762.6 644888.4
## 4 3788135 573049.6 496503.6
## 5 4130301 596469.1 517378.2
## 6 5437616 352970.9 330598.8
## 7 15461588 201448.5 198772.8
## 8 3269661 359607.2 323333.9
## 9 4163760 337984.4 311993.5
## 10 3502501 448044.8 396433.2
## 11 3890461 378076.1 344332.1
## 12 3450043 639797.0 532761.7
## 13 5395920 598437.0 535537.9
## 14 3752715 448407.7 398934.8
## 15 3620623 483707.6 424851.8
## 16 3908042 465705.4 415033.7
## 17 5037709 325573.8 305106.0
## 18 7219994 247361.1 238829.5
## 19 2806124 638546.9 518633.0
## 20 2506955 652993.4 509315.6
## 21 3972698 534270.9 470432.1
## 22 5377194 474404.9 434513.9
## 23 5908922 363487.4 341599.0
## 24 6776416 373341.4 351788.7
## 25 5170901 433868.6 398818.5
## 26 7495904 391773.3 371573.3
## 27 5540779 264631.8 252452.0
## 28 5276906 401467.1 371721.7
## 29 3921581 448315.1 397888.0
## 30 7231743 274438.7 264163.4
## 31 1654849 456064.5 356152.7
## 32 3086749 439691.5 383749.1
## 33 2593562 459714.7 384056.8
## 34 4871898 380931.6 352025.7
## 35 25888108 146765.6 145923.7
## 36 1921725 456027.6 366179.6
## 37 3903917 539446.8 471099.2
## 38 6535503 371771.3 350554.7 e. f = 60%res2_2sc2_0.6_0.6## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1587576 1578601 1574877 5.432394e+12 842786483622 702201043815
## 2 10 2427684 2399478 2400469 1.797021e+13 842692034966 803217458031
## 3 11 1401387 1181179 1207242 6.607451e+12 844312092553 738752634672
## 4 12 1605746 1601218 1599878 6.851222e+12 842845108550 743358843326
## 5 13 1644666 1545833 1556270 8.084347e+12 842950403804 757535344817
## 6 14 2701051 2586828 2592115 1.859297e+13 843135283833 805962108710
## 7 15 4825377 4477796 4481493 7.614969e+13 847675746515 837704961703
## 8 16 2696272 2511261 2521896 1.309290e+13 843950298139 788863342508
## 9 17 2028061 2713785 2645650 8.428132e+12 858971802784 774623323280
## 10 18 1826595 2047310 2022304 7.699168e+12 844408736137 754401713416
## 11 19 1840162 2399841 2336071 7.552226e+12 853457373200 757651131604
## 12 2 1588838 1447062 1464085 5.969805e+12 843392361791 735693665022
## 13 20 1986656 1545138 1571414 1.401575e+13 849729141162 795430437022
## 14 21 1506534 2037784 1965883 6.340563e+12 852736712659 742848747250
## 15 22 2027951 1933005 1941386 8.547682e+12 842921308531 765287188898
## 16 23 2142823 2001454 2012220 9.828678e+12 843510924936 771242613473
## 17 24 2728370 2813293 2806926 1.590379e+13 843218323542 794612791553
## 18 25 3850988 3657485 3661183 3.611787e+13 844533559062 823073472291
## 19 26 2077600 1458565 1526776 6.997647e+12 855859227243 759759297081
## 20 27 1831280 1429147 1481939 5.853109e+12 848390870391 724934950705
## 21 28 1612772 1706487 1695854 7.609203e+12 843199007088 755527077878
## 22 29 1905717 1873519 1876515 9.832287e+12 842526404040 763202581869
## 23 3 1773547 2540345 2472482 1.020944e+13 863104203610 789070982081
## 24 4 2333637 2421249 2416017 1.869100e+13 843326137283 804695076372
## 25 5 2080418 2087722 2085841 1.207070e+13 842816851120 783718921190
## 26 6 2414774 2332135 2334981 1.926514e+13 842911299153 803415956210
## 27 7 2695261 3475228 3433925 1.621017e+13 863789507986 818637299022
## 28 71 2249927 2220809 2220262 1.388323e+13 843016143249 786862681700
## 29 72 1849993 1991523 1969539 6.394576e+12 843520487749 730543101315
## 30 73 3349716 3227808 3231286 2.696894e+13 843512683898 816658833181
## 31 74 2182707 1961266 1998995 3.990798e+12 844522305698 694366148102
## 32 75 1805328 2033623 1991820 5.827641e+12 845172773849 721094867360
## 33 76 1990009 1954471 1957073 4.930874e+12 842860103372 709402703541
## 34 77 2122599 2326346 2302666 1.017614e+13 844582944355 765208738134
## 35 78 6605093 6065277 6067509 2.048828e+14 852933353955 849348616465
## 36 79 2100811 1960406 1981150 4.788615e+12 843399799653 707739691201
## 37 8 2030356 1707254 1734864 9.174525e+12 846353501322 772493686111
## 38 9 2257793 2434545 2425184 1.779779e+13 843713314775 798734868997
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 3311768 533562.4 444925.2
## 2 7384362 352044.2 335365.6
## 3 4682943 719897.6 614812.4
## 4 4225859 527712.7 465388.5
## 5 4880990 548246.6 488438.9
## 6 6875284 325655.1 310616.3
## 7 15586780 188911.4 186518.6
## 8 4841793 336325.5 313050.3
## 9 4125485 316452.2 292816.0
## 10 4176776 412852.4 373352.9
## 11 4080624 355324.0 324009.9
## 12 3715849 582560.8 501915.6
## 13 7005101 551354.1 507037.9
## 14 4148414 418047.4 377585.5
## 15 4208514 436973.6 394896.0
## 16 4559335 421745.4 383414.5
## 17 5784689 299873.4 283238.2
## 18 9300757 230673.9 224570.9
## 19 3359730 588969.3 498879.5
## 20 3161816 593967.0 488952.8
## 21 4689484 495788.8 446714.9
## 22 5087226 451252.6 407291.6
## 23 5684633 339073.3 318387.7
## 24 7940934 348017.5 332737.9
## 25 5776207 403414.0 375544.9
## 26 7908411 361728.6 344302.7
## 27 5983669 248054.8 237943.4
## 28 6078456 379854.3 354553.1
## 29 3388443 424319.3 371374.5
## 30 8000599 261135.4 252525.5
## 31 1825532 431437.0 347865.3
## 32 3140194 416216.2 361857.3
## 33 2445368 431927.0 362644.7
## 34 4708199 363627.2 332607.6
## 35 30942714 140299.8 139659.1
## 36 2269119 431080.2 357723.2
## 37 4506254 496324.4 445645.5
## 38 7766151 348137.7 330544.5 f. f = 50%res2_2sc2_0.5_0.5## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1627076 1628308 1624899 6.813104e+12 928978726628 798002254499
## 2 10 2463634 2403272 2405793 2.060266e+13 929064026036 887906222564
## 3 11 1528357 1173926 1205490 1.050616e+13 933395199323 840330808541
## 4 12 1578191 1597892 1593998 6.840831e+12 928808188079 809556584680
## 5 13 1668720 1535740 1551487 9.661709e+12 929679246033 835478209484
## 6 14 2577312 2589395 2588452 1.758382e+13 929045887219 880049730740
## 7 15 5005324 4578876 4582501 1.091351e+14 936171704033 927436457610
## 8 16 2676179 2499590 2508146 1.367715e+13 930752689302 861574860511
## 9 17 2034597 2703588 2638054 9.974913e+12 944541387897 854480693384
## 10 18 1799060 1995348 1971227 9.088492e+12 930489294600 832097419406
## 11 19 1893053 2440797 2376364 8.655055e+12 939404967147 833408876374
## 12 2 1595747 1433322 1450838 7.713578e+12 929725753382 822432248744
## 13 20 1915551 1501319 1531642 1.236782e+13 935311267136 863434238619
## 14 21 1560879 2051373 2001553 8.992676e+12 937154024731 843795752983
## 15 22 1980581 1868042 1879220 8.104444e+12 929271904428 830462460268
## 16 23 2129981 1964478 1975894 1.106724e+13 930537919170 851995473783
## 17 24 2690581 2832940 2823055 1.276272e+13 929654944051 865022883405
## 18 25 4017313 3757820 3763603 3.993094e+13 931314840076 908662113805
## 19 26 2068217 1422399 1498149 7.404874e+12 943659241765 829656121780
## 20 27 1851414 1410402 1467613 7.038104e+12 935540465413 812191703176
## 21 28 1626755 1700916 1691823 8.253430e+12 929159949604 826415419107
## 22 29 1992458 1962151 1964644 1.360507e+13 928721905481 842557382848
## 23 3 1771144 2615418 2540866 1.396588e+13 953834970204 871255714112
## 24 4 2249290 2429616 2413681 1.962964e+13 930546797185 870200317989
## 25 5 2055809 2082080 2073835 1.164174e+13 929865078601 847177835457
## 26 6 2366850 2256948 2262279 1.723462e+13 929421911227 880177518491
## 27 7 2661085 3536624 3490839 1.801744e+13 955887079242 906527824435
## 28 71 2271207 2279727 2274382 1.597097e+13 929681476751 864225992488
## 29 72 1865281 2041868 2001694 8.992426e+12 931147708963 799012183781
## 30 73 3443270 3348332 3350054 3.472188e+13 929876886564 903865706508
## 31 74 2182305 2006444 2033326 5.092680e+12 930446033728 778659072142
## 32 75 1852538 2084857 2050674 7.901642e+12 931597825255 821003315535
## 33 76 2045557 1991975 1998415 7.323012e+12 929073454249 807338128766
## 34 77 2225991 2408124 2392836 1.164182e+13 931098262949 857653426885
## 35 78 6904907 6313135 6315434 2.437727e+14 941441570390 937767158614
## 36 79 2163777 2005724 2026227 6.107438e+12 929877851983 797068478302
## 37 8 2050571 1718862 1747487 9.938578e+12 933304715951 845488684998
## 38 9 2226086 2445019 2433270 1.715294e+13 930555392949 880842708510
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 4078432 572017.6 491768.8
## 2 8341444 386585.2 369028.1
## 3 6712699 797605.5 697816.7
## 4 4308286 582021.0 508376.7
## 5 5748686 612831.6 542719.5
## 6 6824182 358736.1 339941.2
## 7 21655610 204102.6 202027.5
## 8 5036517 372324.5 343489.5
## 9 4845044 349024.5 323643.7
## 10 4985654 466250.4 422273.9
## 11 4485846 384576.8 350437.2
## 12 4791645 649593.8 567228.4
## 13 6402746 624664.7 564807.0
## 14 5719407 456214.2 420938.3
## 15 4087356 499551.6 443526.7
## 16 5093724 474426.5 431363.8
## 17 4730269 329408.5 307458.5
## 18 9913736 248072.2 241653.3
## 19 3549338 663909.9 553860.6
## 20 3770060 664117.3 553546.9
## 21 5028355 549272.9 490552.1
## 22 6543959 476272.8 429009.4
## 23 7641053 364922.7 343081.5
## 24 8371104 382565.1 360044.0
## 25 5505147 445922.1 407962.8
## 26 7226514 411393.2 388613.0
## 27 6736195 269527.2 258965.8
## 28 6828918 407921.7 380021.7
## 29 4540154 456573.4 398710.4
## 30 10000325 277390.4 269501.9
## 31 2320888 464325.1 383055.4
## 32 4170085 447364.2 400580.3
## 33 3553118 466722.1 403626.5
## 34 5259236 387181.2 358963.2
## 35 35047495 148695.7 148059.4
## 36 2798946 464175.2 393506.4
## 37 4768211 543608.5 484231.8
## 38 7637790 380015.1 361273.3 g. f = 40%res2_2sc2_0.4_0.4## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1608938 1557611 1558547 5.190718e+12 967409313505 793743386389
## 2 10 2416227 2402915 2400820 2.147376e+13 967505716442 918358982359
## 3 11 1401256 1101907 1143363 8.043403e+12 970456750878 829231746429
## 4 12 1586585 1536892 1542301 8.344274e+12 966908404197 844226636462
## 5 13 1630510 1525620 1535718 9.510893e+12 967120625183 855248998386
## 6 14 2655271 2580683 2580901 2.276281e+13 967643347866 914249915954
## 7 15 4938873 4576759 4580224 9.295250e+13 974846239465 962687562172
## 8 16 2677275 2532628 2539626 1.568204e+13 967935532621 900098900006
## 9 17 2093857 2758605 2676693 1.287228e+13 983281668604 874779914195
## 10 18 1780489 2033975 1995930 7.856432e+12 969801309415 845277747189
## 11 19 1913709 2430644 2371134 9.315265e+12 976003401025 864691008646
## 12 2 1571062 1433270 1446140 7.905083e+12 967467446060 836089966587
## 13 20 1989441 1510772 1549123 1.153874e+13 975202542415 883651496723
## 14 21 1521635 2044751 1971050 7.548474e+12 976112106769 837331940986
## 15 22 2080768 1917682 1929628 1.190246e+13 968380238760 891791778101
## 16 23 2106348 1935569 1952039 9.360562e+12 967939269600 872907428595
## 17 24 2685221 2840781 2824863 1.692074e+13 968394878931 900862245340
## 18 25 4021220 3726820 3731211 5.396468e+13 972016489184 952428038584
## 19 26 2051055 1443206 1509826 8.205829e+12 981033039655 865228780783
## 20 27 1857142 1371470 1432948 7.689773e+12 975428763729 827937022923
## 21 28 1615073 1707875 1687394 7.094971e+12 967575832417 828580897317
## 22 29 1917013 1879589 1883232 1.127190e+13 966469894229 859146964814
## 23 3 1628841 2564287 2421307 7.308381e+12 997014263926 855189546878
## 24 4 2293359 2452534 2430914 1.973957e+13 970192984060 895244747402
## 25 5 2063273 2098413 2091147 1.520771e+13 967474597656 896042683191
## 26 6 2459345 2346421 2349466 2.558885e+13 967610948408 922072088837
## 27 7 2682244 3542710 3499430 2.048295e+13 992801655827 943200268860
## 28 71 2147234 2230905 2212756 1.683686e+13 969008762883 871857549879
## 29 72 1882179 1988176 1961953 1.017501e+13 968064671517 816622445035
## 30 73 3382618 3316332 3316473 3.702574e+13 969142318820 941565346717
## 31 74 2237472 1954930 2001194 5.111396e+12 969788346901 804956748229
## 32 75 1741644 2031275 1975459 4.721770e+12 970355604248 785953793154
## 33 76 1991465 1918001 1921147 7.371362e+12 968035724098 819448628980
## 34 77 2231703 2372365 2351149 1.277247e+13 969048895000 883462727317
## 35 78 7006792 6381056 6383398 2.580246e+14 981511722968 977577492341
## 36 79 2040377 1953384 1970451 3.763868e+12 967205174100 756628825038
## 37 8 1978746 1685276 1708194 1.123883e+13 969913682281 882517442749
## 38 9 2287518 2487688 2476785 2.331137e+13 968307051653 917536399218
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 3153396 619713.9 508152.1
## 2 8669037 402552.0 382426.4
## 3 5602197 907527.9 736804.6
## 4 5226260 635385.4 551443.7
## 5 5708137 640786.5 560981.4
## 6 8469026 373310.8 353066.3
## 7 17897437 211212.4 208395.7
## 8 5767054 382712.0 354775.8
## 9 5905521 355497.7 326432.8
## 10 4250082 478474.2 423960.5
## 11 4816394 400612.9 364451.6
## 12 4926134 674766.6 577351.1
## 13 5672017 653807.6 576524.8
## 14 4854094 476362.6 424567.6
## 15 5687518 505571.4 462335.6
## 16 4437752 502117.9 448539.2
## 17 6078669 339622.1 317574.2
## 18 13014893 258933.5 253361.2
## 19 3953713 685582.0 575604.1
## 20 4014553 715263.0 578852.2
## 21 4286485 572514.5 494553.1
## 22 5753529 518053.7 458004.5
## 23 4262099 386576.2 349516.2
## 24 7968161 394369.1 367265.4
## 25 7195689 460464.4 427711.6
## 26 9976121 413397.6 393184.9
## 27 7560726 278260.7 267726.1
## 28 7432775 435848.6 395276.2
## 29 5051538 490870.6 417380.3
## 30 10737960 291135.3 282760.9
## 31 2267684 500293.7 404564.7
## 32 2694388 481062.9 400114.8
## 33 3544299 509863.8 429724.4
## 34 5599295 410137.5 377365.1
## 35 36287010 152420.5 151751.4
## 36 1834830 499461.5 385917.4
## 37 5587364 576656.9 517138.3
## 38 9996805 388921.9 369721.1 h. f = 30%res2_2sc2_0.3_0.3## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1557814 1551710 1549607 7.663532e+12 862447872741 720586586430
## 2 10 2418841 2380705 2379866 2.720478e+13 863018614202 828325176836
## 3 11 1434265 1162108 1195130 6.855437e+12 864946480780 744862700038
## 4 12 1684398 1614045 1615423 9.015394e+12 863678185729 777594452293
## 5 13 1518427 1409623 1413589 9.472012e+12 863861356172 717959483148
## 6 14 2663220 2542298 2545103 2.182192e+13 863799320858 819951954125
## 7 15 4968500 4518861 4520724 Inf 872963898178 867350024817
## 8 16 2733113 2494468 2510205 1.332008e+13 865182525867 804831648328
## 9 17 1921628 2674501 2575846 8.169789e+12 882450240794 771698187576
## 10 18 1751501 2035811 2005929 8.618005e+12 865462576687 770175382544
## 11 19 1876156 2366768 2293985 7.418348e+12 871622256074 757958409190
## 12 2 1685573 1393469 1432704 7.182960e+12 865368166928 744663180202
## 13 20 1902736 1511355 1541664 1.218339e+13 869480028277 766657561152
## 14 21 1470119 2016231 1937660 7.598645e+12 873193034580 756278973200
## 15 22 1945898 1866243 1872252 8.841754e+12 862976999063 767775626238
## 16 23 2183698 2051402 2060014 1.075137e+13 863592441695 783598932527
## 17 24 2572417 2705935 2693831 1.631171e+13 863443435114 792426251872
## 18 25 3771490 3589609 3593979 Inf 863899691451 843664332127
## 19 26 1943330 1406841 1449398 Inf 873282543316 803778578892
## 20 27 1708308 1425428 1480623 3.663351e+12 865589334556 689480776799
## 21 28 1586629 1615170 1605486 5.586135e+12 863102212695 720698588828
## 22 29 1973889 1950889 1953186 1.736921e+13 862126860662 787579191633
## 23 3 1758122 2505354 2419328 1.222021e+13 882403644894 785985118275
## 24 4 2347631 2390382 2374994 2.845884e+13 865792902569 810716576413
## 25 5 1940240 2054353 2042372 9.457369e+12 863232810594 781676053306
## 26 6 2453888 2284372 2287929 Inf 865528211247 842072336754
## 27 7 2630302 3474058 3431353 1.984480e+13 887617178641 843132811825
## 28 71 2042854 2221886 2212271 Inf 865394557120 829396420807
## 29 72 1795388 1989355 1941028 5.961947e+12 865028504516 715398430980
## 30 73 3128778 3204950 3199962 3.552441e+13 864483464059 836881813986
## 31 74 2209428 1959381 1992305 5.018326e+12 865734716774 709682713525
## 32 75 1767609 2031637 1979040 6.719013e+12 865585711557 736163834980
## 33 76 2249655 2003144 2022163 9.591414e+12 866416353159 771022342745
## 34 77 2338948 2276602 2270834 Inf 864329038067 811761082019
## 35 78 6572536 6018285 6018431 Inf 874013373966 873580172890
## 36 79 1992904 1961236 1942123 3.811316e+12 864404102476 641303418335
## 37 8 1968301 1675271 1701705 1.136806e+13 866676050676 765367731946
## 38 9 2296794 2431599 2423823 2.499399e+13 863540823659 825807439717
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 4534153 566217.2 468474.2
## 2 10819093 365254.6 350486.6
## 3 4648710 748037.9 624935.3
## 4 5174193 537225.3 482259.9
## 5 5440042 622356.6 510616.2
## 6 7972815 339359.3 321762.0
## 7 Inf 192701.1 191349.8
## 8 4830457 347632.5 321184.1
## 9 4127633 330816.4 300291.3
## 10 4809539 428549.1 385779.1
## 11 3735501 368064.4 329463.8
## 12 4091819 631576.5 522945.3
## 13 5910714 581288.1 499692.5
## 14 4968568 432930.0 390539.6
## 15 4423524 467022.4 412685.3
## 16 4766337 424551.3 383012.8
## 17 5868201 321256.3 295442.7
## 18 Inf 241153.5 235040.5
## 19 Inf 625888.6 560924.2
## 20 2138116 608328.3 465735.2
## 21 3345460 543199.4 454095.1
## 22 8252392 448185.9 406136.6
## 23 6486564 352376.5 324443.8
## 24 10944659 362762.8 341713.2
## 25 4856744 421042.1 383462.7
## 26 Inf 385327.5 372898.1
## 27 7481083 254694.0 244975.7
## 28 Inf 389519.9 374860.1
## 29 3046012 435403.4 367833.6
## 30 10878292 269491.1 261123.3
## 31 2203889 442956.3 356276.4
## 32 3675499 426822.7 370757.9
## 33 4001007 434383.7 381604.9
## 34 Inf 380267.7 358383.9
## 35 Inf 144712.3 144639.1
## 36 1822887 441448.4 330668.3
## 37 5528810 522888.4 451792.3
## 38 10764413 355791.1 341297.3 i. f = 20%res2_2sc2_0.2_0.2## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1433981 1568561 1568041 Inf 699297309759 628255006000
## 2 10 2260716 2254727 2254727 Inf 700191429368 700191429368
## 3 11 1492489 1307961 1314142 Inf 701879901671 677947206218
## 4 12 1512705 1557990 1557930 Inf 700618842586 646525145030
## 5 13 1852117 1696136 1696136 Inf 702470109218 702470109218
## 6 14 2595039 2522857 2522857 Inf 701674743585 701674743585
## 7 15 4196157 4093781 4093781 Inf 704549221215 704549221215
## 8 16 2802059 2462385 2464740 Inf 707274324580 692165517563
## 9 17 2044691 2604661 2599024 Inf 713570070461 701883911982
## 10 18 1674359 1992998 1948882 Inf 707927991344 639068692583
## 11 19 1872288 2334142 2323064 Inf 707169371439 690516352319
## 12 2 1620708 1541808 1548523 Inf 699230347662 623409568518
## 13 20 1968137 1569962 1587546 Inf 708734204366 647256529434
## 14 21 1473338 1996094 1996094 Inf 709919153558 709919153558
## 15 22 1924211 1975112 1974623 Inf 698614640538 694475987140
## 16 23 2173965 1942761 1960420 Inf 703122201857 659455719109
## 17 24 2663297 2825912 2825625 Inf 699446686476 697534727788
## 18 25 3824371 3434742 3443650 Inf 706862664182 691163838796
## 19 26 2157053 1542029 1555484 Inf 714152578747 684541005548
## 20 27 1994931 1419642 1445179 Inf 711563663492 672494408479
## 21 28 1451899 1708893 1677066 Inf 701270284784 604000398323
## 22 29 1757749 1743382 1742365 Inf 697954336372 636092976163
## 23 3 1731292 2347761 2327069 Inf 715133394691 690333737152
## 24 4 2623709 2343380 2330951 Inf 711502248924 667931414773
## 25 5 1934975 2044576 2044576 Inf 702133299566 702133299566
## 26 6 2276835 2254195 2254195 Inf 701402509190 701402509190
## 27 7 2533516 3329541 3329541 Inf 723274631363 723274631363
## 28 71 2006234 2096386 2096386 Inf 710909560940 710909560940
## 29 72 1493890 1885636 1885636 Inf 705164980820 705164980820
## 30 73 3150922 2984196 2984196 Inf 704678019113 704678019113
## 31 74 2138921 1858499 1880132 6.082540e+12 702700453263 604839233636
## 32 75 1740888 1929846 1929012 Inf 706489074098 701395962354
## 33 76 1879651 1815536 1824253 6.690393e+12 700523026124 501759876968
## 34 77 1950094 2215518 2215518 Inf 710548225610 710548225610
## 35 78 5988669 5482654 5482654 Inf 710396013546 710396013546
## 36 79 2114429 1855137 1855266 Inf 711209716098 657529940189
## 37 8 2039894 1700544 1706439 Inf 704853475588 687051895621
## 38 9 2236321 2346558 2346691 Inf 701604690216 700813135422
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 Inf 455046.9 404266.2
## 2 Inf 310660.5 310660.5
## 3 Inf 556338.5 532070.4
## 4 Inf 452680.9 415400.0
## 5 Inf 439052.1 439052.1
## 6 Inf 276495.3 276495.3
## 7 Inf 170843.8 170843.8
## 8 Inf 290211.3 283527.0
## 9 Inf 274963.2 270945.3
## 10 Inf 361366.1 331461.8
## 11 Inf 303800.9 297954.8
## 12 Inf 474974.2 418781.7
## 13 Inf 464707.6 416470.5
## 14 Inf 358313.9 358313.9
## 15 Inf 371683.2 369922.5
## 16 Inf 375153.8 347058.6
## 17 Inf 249349.3 248718.4
## 18 Inf 205117.5 199925.8
## 19 Inf 488664.1 456539.5
## 20 Inf 512244.6 477377.5
## 21 Inf 421494.4 366619.1
## 22 Inf 415969.5 379775.8
## 23 Inf 307461.8 298134.7
## 24 Inf 304429.4 287709.8
## 25 Inf 345347.7 345347.7
## 26 Inf 318028.9 318028.9
## 27 Inf 215713.6 215713.6
## 28 Inf 340354.5 340354.5
## 29 Inf 378603.4 378603.4
## 30 Inf 241067.3 241067.3
## 31 2752798 381333.2 323897.4
## 32 Inf 368679.5 365864.8
## 33 2567993 400642.5 278297.4
## 34 Inf 327788.6 327788.6
## 35 Inf 128197.1 128197.1
## 36 Inf 387853.5 359762.2
## 37 Inf 417036.8 405317.0
## 38 Inf 301351.2 301034.5 2. Metode PPSWR-SRSWOR res_2sc2_1_1## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1612935 1560421 1612492 7845037888 886390683405 7770821543
## 2 10 2415381 2395026 2415192 7399776586 886492986512 7336749142
## 3 11 1431381 1199537 1429027 8245932030 888298285791 8158489725
## 4 12 1637684 1623131 1637399 8212889932 886463683960 8132660852
## 5 13 1665954 1542576 1665011 5690831695 887238449554 5650504039
## 6 14 2703353 2599541 2701691 12564929983 887161771311 12382595960
## 7 15 5014684 4559697 4977639 74159897116 894259889673 67912089963
## 8 16 2677256 2581864 2676464 7115205101 886708060103 7057605914
## 9 17 2068793 2754404 2074053 7137834371 902680647689 7078235249
## 10 18 1814688 2075180 1816129 5324417946 888906494266 5290491988
## 11 19 1933547 2456105 1937630 7093823887 895850535437 7036205945
## 12 2 1629336 1424421 1627616 7129361218 887894738068 7069932667
## 13 20 1929510 1560682 1925047 10605830742 891262468190 10477414664
## 14 21 1527249 2059804 1531038 6618947823 896312564065 6567189028
## 15 22 2067771 1930708 2066569 7435599617 887090822863 7368367928
## 16 23 2141332 2018621 2139784 10602832315 887118411120 10476065477
## 17 24 2717762 2840261 2719021 9741654428 887016204855 9629988331
## 18 25 3860999 3691128 3850785 48343981657 887653011035 45502831374
## 19 26 2013072 1459284 2005677 11814088696 897221743272 11654708888
## 20 27 1832752 1422739 1829195 7756597238 892048295986 7688538472
## 21 28 1640965 1709445 1641798 11455772041 886650155914 11305286114
## 22 29 1919034 1886470 1918549 12771124653 886282687644 12566344894
## 23 3 1711264 2555005 1721090 10768567268 910898093650 10636073546
## 24 4 2341407 2451042 2343053 13826138430 886831648408 13603803896
## 25 5 2037005 2099347 2037426 8812445207 886610972886 8718452079
## 26 6 2538739 2346667 2535248 16257549057 887590967383 15953040737
## 27 7 2758099 3520204 2767674 11551388927 906563376056 11396350027
## 28 71 2344845 2302714 2342360 22042402134 887266431101 21424985689
## 29 72 1960940 2069616 1962577 20428995070 887168654312 19905509969
## 30 73 3460877 3346887 3455279 38430312126 887387836520 36733178489
## 31 74 2201707 2032261 2200390 6776932651 887385416865 6723008361
## 32 75 1880316 2110536 1884376 17048243570 888423339507 16702339110
## 33 76 2035745 2044475 2035926 12381994590 886424702428 12207519677
## 34 77 2305963 2417556 2308824 27535106786 887542723918 26664124637
## 35 78 6705015 6209264 6636513 142625267898 895077168157 122736665150
## 36 79 2233947 2035182 2231721 9802575281 887823360481 9692357766
## 37 8 1980133 1696584 1978106 6140168989 889115576668 6095189514
## 38 9 2288451 2451507 2290737 14256484441 887365223104 14006967594
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 4853.123 568436.3 4808.694
## 2 3058.464 370652.3 3032.713
## 3 5654.902 742396.1 5606.738
## 4 4969.290 546690.6 4922.297
## 5 3367.558 577584.0 3346.242
## 6 4646.869 341186.8 4582.424
## 7 14543.652 195970.8 13443.053
## 8 2650.016 343363.7 2629.417
## 9 3423.255 327608.9 3386.141
## 10 2913.504 428215.6 2892.775
## 11 3665.351 364609.5 3627.818
## 12 4351.572 623419.5 4320.315
## 13 5467.735 571160.1 5414.861
## 14 4298.073 434869.9 4253.922
## 15 3579.096 460441.4 3549.125
## 16 4937.018 439525.3 4881.790
## 17 3573.239 312750.9 3530.866
## 18 12433.804 240578.5 11747.024
## 19 5854.803 615251.9 5798.067
## 20 4231.085 627241.1 4202.258
## 21 6961.448 520177.7 6867.011
## 22 6600.296 469940.0 6496.927
## 23 6256.953 356688.2 6144.351
## 24 5907.004 361615.9 5807.709
## 25 4294.701 422048.6 4248.558
## 26 6391.706 378155.9 6281.131
## 27 4180.319 257334.3 4109.307
## 28 9217.863 385414.9 8978.263
## 29 10232.454 429048.3 9968.514
## 30 11074.596 265335.7 10607.336
## 31 3069.758 437003.1 3047.311
## 32 9007.164 421217.4 8805.899
## 33 6071.746 434329.9 5985.623
## 34 11899.931 367618.2 11511.238
## 35 21252.793 144035.1 18480.951
## 36 4377.652 436624.4 4333.042
## 37 3080.686 524354.1 3061.640
## 38 6165.905 362174.3 6052.795 b. f = 90%res_2sc2_0.9_0.9## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1626390 1551992 1625550 9693264202 961673504665 9581648555
## 2 10 2438008 2402067 2437511 9504756871 961670789381 9406676321
## 3 11 1445241 1129410 1442244 8799788621 965090417143 8711458658
## 4 12 1627870 1583333 1627300 9672503774 962068912236 9570883569
## 5 13 1671326 1565277 1670382 7644948444 962187266485 7577304462
## 6 14 2676250 2576179 2674864 12325100140 962000064587 12164623125
## 7 15 4849016 4533742 4826333 65408777824 965687290214 60710663314
## 8 16 2740643 2517961 2738275 9702601088 963539613707 9601061717
## 9 17 2133617 2751973 2140924 11929637740 974843394592 11776245549
## 10 18 1751679 2047034 1753399 5646295014 964447034746 5610680678
## 11 19 1914542 2444926 1919285 8824383932 971258474075 8741459130
## 12 2 1626493 1422611 1624495 9436364125 962995067204 9341861157
## 13 20 1933575 1527446 1927962 12865725791 967420603355 12684364665
## 14 21 1505518 2069853 1509426 6736525370 972478223585 6689686865
## 15 22 2058827 1876060 2056986 9274561977 962781835700 9180308681
## 16 23 2113374 1908075 2110956 11541798094 963054310942 11396878019
## 17 24 2652589 2726637 2653288 10624772799 962140357933 10504482652
## 18 25 3782932 3639673 3774509 48971260960 962645660307 46391546023
## 19 26 1939943 1453541 1933541 12445784034 969877115141 12279896295
## 20 27 1852038 1418125 1848094 8701005530 967960871443 8619482119
## 21 28 1642727 1747255 1643724 12248691133 962092639585 12079014735
## 22 29 1920792 1878593 1920179 13986110627 961419042731 13755147570
## 23 3 1814498 2625662 1827190 15864476906 984155490451 15588261689
## 24 4 2300902 2454046 2302894 16418309189 962872115124 16110290963
## 25 5 2012606 2137012 2013386 8185003200 962225712366 8106689966
## 26 6 2549791 2392911 2546841 17096913042 962585097609 16777044676
## 27 7 2674921 3545569 2683969 10457479158 987951787078 10335068105
## 28 71 2271595 2272801 2270697 20900350803 961921393580 20382135515
## 29 72 1903462 2040349 1906569 24067754065 962166842153 23445254550
## 30 73 3357280 3410888 3357851 44415206110 962543482767 42313063925
## 31 74 2195767 2001503 2194208 7903443651 962905471979 7837211004
## 32 75 1912463 2074170 1915330 22469334661 962903092567 21882041091
## 33 76 2050117 1989632 2049288 13133684172 961591732624 12950648057
## 34 77 2248827 2432414 2252935 27302387196 963320423982 26467528404
## 35 78 7181761 6497787 7048424 233666216773 978359206092 187225909367
## 36 79 2179322 2003532 2177075 11602119333 962822459543 11457461351
## 37 8 2001876 1681306 1998734 9011793429 965166377461 8921875302
## 38 9 2281095 2508325 2283892 13493990740 963395409305 13275527024
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 5929.495 619561.1 5865.003
## 2 3877.320 400711.6 3838.549
## 3 5989.550 863345.5 5943.875
## 4 5907.765 609676.0 5848.426
## 5 4553.822 617554.6 4516.564
## 6 4587.888 373542.8 4530.851
## 7 13288.966 212877.9 12416.309
## 8 3529.949 382815.9 3496.336
## 9 5551.114 353899.4 5461.348
## 10 3216.852 471399.1 3193.351
## 11 4606.794 396798.8 4552.060
## 12 5799.813 676866.0 5748.877
## 13 6623.230 633772.9 6550.680
## 14 4482.174 469388.8 4439.146
## 15 4457.748 515410.9 4417.156
## 16 5434.696 508624.7 5373.265
## 17 3974.668 353609.7 3928.997
## 18 12789.420 264528.8 12153.914
## 19 6380.646 667789.6 6317.757
## 20 4703.095 683760.3 4669.456
## 21 7333.772 550579.5 7230.584
## 22 7158.860 514485.6 7044.554
## 23 8674.527 374456.0 8460.770
## 24 7026.182 391902.2 6891.484
## 25 3993.914 449711.9 3955.070
## 26 6646.680 402094.7 6531.446
## 27 3890.381 278481.4 3831.198
## 28 9039.350 422960.0 8828.206
## 29 12593.408 471341.8 12249.671
## 30 13150.935 281997.8 12532.991
## 31 3599.474 480934.9 3571.838
## 32 11539.899 464063.5 11225.090
## 33 6394.532 484453.3 6308.319
## 34 11990.422 395527.7 11610.579
## 35 32197.454 150228.3 26344.161
## 36 5295.369 480410.0 5235.266
## 37 4480.942 574579.5 4443.875
## 38 5842.736 384530.3 5742.571 c. f = 80%res_2sc2_0.8_0.8## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1613295 1592433 1612984 11107350123 938558016051 10962130744
## 2 10 2468849 2415596 2468114 11515759602 938785852843 11371488821
## 3 11 1405268 1144375 1402462 10097716352 940898666591 9983573977
## 4 12 1620193 1582832 1619505 12456311701 938877104201 12283080064
## 5 13 1668813 1584199 1667974 8937121944 938780668638 8831104601
## 6 14 2673893 2579573 2672245 14895135498 939173720597 14654420537
## 7 15 4881001 4519874 4850243 82092901976 944022479184 75122229128
## 8 16 2700579 2519860 2698327 11164272662 939849274200 11030733720
## 9 17 2026064 2750283 2034488 11380865574 956728310001 11240710660
## 10 18 1791337 2037516 1793347 7794524478 940779286319 7726849770
## 11 19 1929901 2453867 1935484 10279822544 948010569614 10162556603
## 12 2 1596118 1447619 1593827 11802503170 939665375544 11642145485
## 13 20 1954354 1533640 1948253 13494477156 945505501226 13291304044
## 14 21 1559104 2070385 1564705 10969158423 947979196130 10837612910
## 15 22 2032290 1904096 2030915 10197008683 939158609264 10083081686
## 16 23 2092671 1915266 2090176 13373967296 940081975130 13181790746
## 17 24 2746883 2800358 2747398 14445460222 939249453056 14216353685
## 18 25 3845604 3694678 3829489 78743828382 940718988069 72002146654
## 19 26 2058832 1464991 2048791 16158769646 951023521065 15872638712
## 20 27 1833237 1416024 1829154 9142025026 944606161475 9049949910
## 21 28 1637936 1764262 1639968 15085258258 939162561403 14841432211
## 22 29 1883531 1845320 1882897 16834038897 938531046989 16440937790
## 23 3 1830552 2620892 1848734 22986703161 960189390484 22407535595
## 24 4 2321493 2455610 2323227 20368015458 940279689130 19891703529
## 25 5 2088145 2123686 2088414 12461226158 938934521244 12287622701
## 26 6 2500373 2354786 2496324 21348690376 939876951892 20801243032
## 27 7 2690238 3520241 2700501 12136473112 962806048285 11977549362
## 28 71 2291240 2281929 2289102 29167051634 939360329498 28072824786
## 29 72 1901274 2050209 1903284 27116631953 940703642130 26177161008
## 30 73 3464057 3391080 3459958 51231965117 939283532747 48509316205
## 31 74 2217843 2012933 2215410 11372843873 940180480492 11236126791
## 32 75 1902639 2086858 1905395 20298421920 940259354415 19712929137
## 33 76 1976430 1980957 1976396 15627653154 938954765163 15344083445
## 34 77 2307293 2446674 2311051 39726044858 940173246058 37927611076
## 35 78 7007097 6417657 6891262 231473974917 951346362258 185457009614
## 36 79 2143495 2013226 2141321 12382751982 939458001742 12196806972
## 37 8 2050152 1710251 2045583 12390276059 942787014211 12213064485
## 38 9 2290529 2498570 2294228 19924567987 940311926681 19408815082
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 6880.149 589758.6 6791.895
## 2 4645.364 389084.1 4588.909
## 3 7164.359 827872.7 7098.603
## 4 7626.855 594837.3 7525.478
## 5 5249.913 593216.8 5191.922
## 6 5546.680 363983.3 5460.929
## 7 16695.025 208744.2 15392.702
## 8 4134.017 373475.7 4088.196
## 9 5581.908 347832.0 5490.385
## 10 4335.893 462216.2 4293.367
## 11 5303.914 386138.4 5228.043
## 12 7301.589 649175.0 7215.196
## 13 6842.332 618292.3 6763.283
## 14 6966.020 457548.5 6858.188
## 15 5008.658 494291.0 4956.320
## 16 6379.204 493910.7 6295.458
## 17 5238.941 336553.5 5155.618
## 18 20093.748 254895.2 18490.649
## 19 7817.649 649574.0 7719.312
## 20 4964.274 667946.1 4925.820
## 21 9204.267 532759.6 9043.909
## 22 8728.801 511992.2 8536.494
## 23 12424.484 366021.9 11990.543
## 24 8615.071 382727.6 8411.522
## 25 5938.018 441785.5 5854.624
## 26 8416.481 399190.8 8219.011
## 27 4499.002 273098.1 4422.345
## 28 12429.058 412006.9 11990.603
## 29 13756.880 459472.2 13284.070
## 30 14763.890 276670.7 13997.146
## 31 5127.787 467808.3 5071.733
## 32 10366.359 451268.8 10060.867
## 33 7882.562 476081.8 7740.827
## 34 16957.061 384084.2 16181.254
## 35 32724.228 147908.5 26707.698
## 36 5689.100 467392.4 5612.794
## 37 5991.990 551774.9 5921.568
## 38 8520.071 376406.2 8291.865 d. f = 70%res_2sc2_0.7_0.7## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1650667 1576934 1649380 15175811976 9.791660e+11 14908162102
## 2 10 2470594 2368112 2468887 16391186550 9.795049e+11 16113246720
## 3 11 1402838 1142862 1399211 12755456565 9.813624e+11 12571786960
## 4 12 1663858 1610783 1663014 12641405476 9.789735e+11 12470341887
## 5 13 1672665 1551920 1671325 9932265013 9.794714e+11 9816262016
## 6 14 2735300 2603447 2732714 18997751231 9.796481e+11 18625801175
## 7 15 5055832 4614164 5008389 110747193821 9.867601e+11 98894047446
## 8 16 2753737 2549555 2750770 14104875112 9.806745e+11 13900447743
## 9 17 2073268 2764231 2083135 14839380956 9.955700e+11 14602189390
## 10 18 1826460 2079805 1828953 9946880061 9.810006e+11 9843350166
## 11 19 1931728 2479566 1938547 12753946259 9.893110e+11 12581032501
## 12 2 1692038 1444265 1688520 14274626290 9.810539e+11 14061522886
## 13 20 1945349 1543967 1937734 18144111822 9.844670e+11 17790409555
## 14 21 1566762 2101602 1573370 12981356190 9.889517e+11 12796227945
## 15 22 2003796 1905972 2002460 12153172022 9.794383e+11 11987932363
## 16 23 2142764 1991084 2139861 17206709546 9.801919e+11 16900358363
## 17 24 2698532 2732142 2698816 18679149047 9.792439e+11 18309266514
## 18 25 3755091 3655901 3746345 68792963196 9.795707e+11 63724178240
## 19 26 2043069 1460508 2032159 18364383294 9.908522e+11 18013381103
## 20 27 1826372 1454128 1821599 12107134406 9.838126e+11 11943708386
## 21 28 1582503 1727328 1584468 14065439946 9.794960e+11 13853717717
## 22 29 1918438 1871745 1917562 20011597064 9.787701e+11 19530953690
## 23 3 1796635 2676627 1813548 21118273374 1.005940e+12 20544761080
## 24 4 2210458 2474114 2214863 20119513495 9.817864e+11 19631117643
## 25 5 2169663 2154478 2168879 17588819417 9.792189e+11 17235661468
## 26 6 2576101 2406993 2570345 26808747793 9.804237e+11 25978599289
## 27 7 2656107 3584777 2668687 13901506452 1.008670e+12 13703524240
## 28 71 2249232 2309992 2250080 26501745249 9.794967e+11 25695327841
## 29 72 1903112 2074411 1908211 36886680719 9.803565e+11 35398030327
## 30 73 3629204 3484440 3614874 79350280735 9.820434e+11 73132392440
## 31 74 2275883 2029183 2272355 13930015008 9.815374e+11 13729028819
## 32 75 1883769 2105734 1888295 26266705493 9.809350e+11 25277295413
## 33 76 2101824 2073315 2099633 28488217416 9.801771e+11 27597481123
## 34 77 2268690 2455772 2275437 49438536748 9.812054e+11 46863465850
## 35 78 7150378 6545754 7004739 307670389845 9.924581e+11 233684030891
## 36 79 2223103 2035821 2218843 18357841678 9.804056e+11 17970024871
## 37 8 2045992 1675646 2040272 14439548966 9.839436e+11 14216588154
## 38 9 2243524 2478256 2248271 19357483152 9.808250e+11 18871401274
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 9161.777 621101.4 9009.030
## 2 6630.442 413411.3 6522.413
## 3 8922.941 864851.6 8822.800
## 4 7548.265 609006.5 7451.286
## 5 5887.241 632647.4 5824.698
## 6 6946.123 376373.4 6816.797
## 7 21762.390 213704.0 19648.112
## 8 5137.434 385176.1 5068.432
## 9 7087.280 360041.7 6940.812
## 10 5424.362 471826.1 5360.733
## 11 6565.482 398854.7 6454.226
## 12 8439.864 679493.8 8331.019
## 13 9240.482 639587.8 9099.250
## 14 8193.905 470256.2 8044.246
## 15 5991.719 514474.1 5915.392
## 16 8011.472 493485.2 7880.487
## 17 6870.426 359183.4 6735.447
## 18 18086.398 268036.9 16812.579
## 19 8891.680 678849.8 8772.417
## 20 6607.508 676722.6 6537.273
## 21 8810.497 567390.5 8668.291
## 22 10172.365 525160.6 9936.412
## 23 11359.708 375519.6 10948.010
## 24 8956.086 396643.8 8725.851
## 25 8009.820 454215.2 7855.350
## 26 10222.636 407070.9 9940.549
## 27 5224.861 280996.6 5125.515
## 28 11611.235 424140.8 11257.409
## 29 19102.155 472843.4 18292.742
## 30 21592.681 281544.5 20008.008
## 31 6122.890 484094.3 6043.970
## 32 13620.493 466244.7 13086.890
## 33 13352.260 473212.0 12963.566
## 34 21389.295 399374.2 20237.155
## 35 42706.781 151298.5 33168.963
## 36 8132.863 481892.6 7981.594
## 37 7015.729 587681.1 6929.111
## 38 8510.906 395404.2 8274.357 e. f = 60%res_2sc2_0.6_0.6## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1644434 1605805 1643332 18118965801 908200785137 17716933384
## 2 10 2497052 2433711 2495378 19968921372 908854904609 19522671037
## 3 11 1425782 1167113 1420655 16945580889 910312525908 16580106058
## 4 12 1633397 1636190 1633063 18200330658 908234214045 17817567965
## 5 13 1684709 1501594 1682348 13449610535 909735160916 13224223948
## 6 14 2577647 2604326 2577429 20734481903 908771697285 20208813865
## 7 15 5162766 4609724 5087633 136687720705 919285513255 117652324728
## 8 16 2653264 2543531 2651538 15074242059 909129059809 14816752487
## 9 17 2026422 2745347 2035811 12361370402 925875263127 12190486506
## 10 18 1776763 2074343 1780209 11713563035 911331348909 11540955315
## 11 19 1924716 2441876 1933489 16098587612 917341343591 15801456276
## 12 2 1618041 1429545 1614247 16079287539 909535866008 15768498274
## 13 20 1941690 1545158 1931409 22134821219 914487841799 21563622467
## 14 21 1601335 2066416 1609561 17450859539 915730813243 17081265944
## 15 22 2051407 1935512 2049459 15662183533 908603798844 15382151163
## 16 23 2159135 2056218 2156842 21673799185 908793944495 21158668870
## 17 24 2666080 2815640 2669453 22847470405 909190012109 22233358919
## 18 25 3922016 3727320 3897536 102652652035 910415858155 90967073402
## 19 26 2079482 1442713 2062629 24342016144 922506920168 23679813529
## 20 27 1819708 1451631 1813343 14766082669 913286422142 14509942577
## 21 28 1681970 1704999 1682421 23663246581 908057355025 23022531493
## 22 29 1986561 1959552 1985577 35329435729 907885142587 33701810504
## 23 3 1756587 2622066 1781706 28370110712 934330464112 27445734649
## 24 4 2289956 2450584 2293771 31656441499 909483334639 30273360324
## 25 5 2041752 2101780 2042411 16239386080 908342076897 15926110912
## 26 6 2558940 2340476 2549958 34246427824 910144070620 32877120303
## 27 7 2724844 3523994 2741912 20953149688 930437795380 20451704419
## 28 71 2439866 2295940 2425826 55507761599 910186956254 51515629838
## 29 72 1881701 2062662 1885853 45115906920 910754896330 42466802116
## 30 73 3560345 3356046 3535053 85280956992 912266300799 77020748073
## 31 74 2165633 2026927 2162953 16153500600 909223496832 15862825031
## 32 75 1883045 2101507 1890094 35637659659 909855499854 33893974307
## 33 76 2067558 2058725 2067385 26510055884 908715021800 25728025677
## 34 77 2158966 2397653 2167141 43430483121 911058823080 41054332073
## 35 78 6715301 6247672 6603031 279328935965 916194821165 207949416125
## 36 79 2200540 2029856 2196157 22046545303 909647430980 21319140673
## 37 8 2031803 1714617 2024611 20106028785 911672262494 19651108330
## 38 9 2287244 2469690 2292037 27034897760 909522791525 26065844439
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 10987.510 566738.5 10754.419
## 2 7950.961 372970.6 7779.901
## 3 11584.265 784496.7 11385.704
## 4 11020.596 556128.8 10794.888
## 5 7928.522 610884.8 7806.363
## 6 7944.986 348780.8 7748.943
## 7 26107.961 199130.6 22862.662
## 8 5704.411 358367.4 5610.061
## 9 6076.537 337424.3 5963.994
## 10 6478.775 441279.6 6373.380
## 11 8346.358 375753.3 8154.033
## 12 9776.326 639651.2 9615.367
## 13 11274.411 595380.5 11048.747
## 14 10715.252 443403.3 10439.265
## 15 7633.088 472028.3 7504.542
## 16 10026.977 443928.2 9799.067
## 17 8532.262 323410.0 8294.022
## 18 25561.956 244055.0 22877.524
## 19 11620.016 642797.3 11402.216
## 20 7970.037 631156.8 7864.903
## 21 13978.789 533986.1 13600.721
## 22 17217.022 468378.4 16449.508
## 23 16035.060 356061.1 15288.799
## 24 13537.451 371138.1 12930.666
## 25 7867.850 432386.1 7715.985
## 26 13216.900 390200.0 12747.119
## 27 7612.013 263520.6 7383.962
## 28 21902.421 397210.6 20542.322
## 29 23003.766 443031.5 21659.438
## 30 23576.879 271824.3 21518.207
## 31 7415.678 449953.9 7292.600
## 32 18537.824 434180.9 17567.193
## 33 12968.660 444285.1 12583.123
## 34 19635.564 380689.5 18515.615
## 35 40773.467 146313.1 31054.181
## 36 9943.901 449699.0 9643.022
## 37 9841.760 531630.9 9656.502
## 38 11538.057 368597.2 11101.552 f. f = 50%res_2sc2_0.5_0.5## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1627088 1602180 1625084 24638308146 950240582731 23565601328
## 2 10 2396439 2402405 2396052 20428579382 950198959850 19925194015
## 3 11 1427259 1130374 1418735 24133274715 953762012648 23357999588
## 4 12 1671947 1599374 1669092 25578778963 950908668409 24826842918
## 5 13 1719898 1614243 1715921 21231994533 950895205664 20525362847
## 6 14 2657524 2570184 2654105 36638947394 950712258989 35171731642
## 7 15 4683909 4494463 4658476 128540434459 953813919585 111001250825
## 8 16 2705955 2482806 2699452 27616283786 952319845413 26783421570
## 9 17 2035929 2731415 2048240 19620597435 968121995991 19089810981
## 10 18 1803288 2046903 1807165 16403932993 952538526071 16091711742
## 11 19 1842748 2447135 1854610 20299155434 963285422620 19813706631
## 12 2 1628224 1443629 1623393 21747148144 951550105666 21214580669
## 13 20 1924077 1527688 1908948 34414509622 956534983178 33035823283
## 14 21 1567167 2080915 1577216 19955064841 959529751653 19464738874
## 15 22 2023752 1916978 2021091 20734612315 950671405696 20198658273
## 16 23 2234872 1920026 2224789 31364432960 953948209570 30320581349
## 17 24 2641165 2720813 2643676 26174791664 951018569819 25437609674
## 18 25 3963680 3646357 3924355 118232550672 955183544326 103605269025
## 19 26 2090212 1465403 2071855 28218134857 964026167546 27343194492
## 20 27 1843741 1457139 1835598 18912379958 955414107010 18477822506
## 21 28 1615498 1789639 1619387 28066681986 951406187877 27066134869
## 22 29 2011279 1969307 2009477 41213935916 949898310862 39023656915
## 23 3 1877050 2715030 1909222 40883092253 974964372937 39061120716
## 24 4 2189046 2445423 2194085 31435367194 953002882747 30041946417
## 25 5 2115468 2138222 2114369 30328267446 950770102388 29137946637
## 26 6 2582670 2388805 2571506 45701624326 951957163337 42968382739
## 27 7 2600831 3501617 2620459 22995147323 979167814598 22395049508
## 28 71 2406074 2287452 2393703 61282292406 952011993313 56756985986
## 29 72 1897002 2065856 1902611 51803711263 952496357532 48582471525
## 30 73 3350928 3482212 3356252 85592044177 953154684842 77491884975
## 31 74 2220538 2020928 2215001 25636944324 951597543985 24929052604
## 32 75 1808320 2090615 1816212 48450556597 953837090274 45085614515
## 33 76 2156520 2037983 2151094 37902423067 951238602992 36257492543
## 34 77 2351942 2475444 2358761 82659687369 951552845800 75338309452
## 35 78 7251027 6506432 7044850 374456497721 970842815337 261597427918
## 36 79 2200908 2026501 2194946 32482045469 951168719623 31246442770
## 37 8 2019609 1690557 2010638 20802512192 954260313275 20157129193
## 38 9 2298544 2509675 2304854 28885719260 951833140499 27788805303
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 14457.782 594785.3 13906.588
## 2 8416.485 397001.9 8214.584
## 3 16310.666 862239.8 15925.256
## 4 15082.653 601013.5 14676.735
## 5 11869.205 596945.6 11538.636
## 6 13722.912 371236.5 13194.140
## 7 27403.815 211731.1 23768.939
## 8 10181.812 383666.4 9899.544
## 9 9311.775 354037.4 9012.319
## 10 9031.299 466535.2 8840.898
## 11 10867.528 393091.4 10542.689
## 12 13180.151 661878.5 12904.549
## 13 17533.459 630377.1 16992.354
## 14 12577.104 460851.4 12190.012
## 15 10137.809 499210.3 9895.007
## 16 13964.157 500777.7 13564.413
## 17 9935.803 350389.2 9642.837
## 18 29574.466 261838.8 26242.353
## 19 13348.111 660337.8 13062.464
## 20 10105.296 660882.8 9927.567
## 21 17016.015 533482.4 16389.311
## 22 19826.586 483371.1 18818.439
## 23 21433.947 358888.9 20122.071
## 24 13712.120 389267.1 13109.283
## 25 13958.431 444552.4 13434.923
## 26 17133.674 399421.8 16230.390
## 27 8711.219 279045.6 8424.119
## 28 24453.570 416227.8 22897.068
## 29 26456.201 461782.9 24770.601
## 30 25242.138 273063.7 22861.583
## 31 11565.156 471891.7 11274.609
## 32 25702.460 456944.6 23953.616
## 33 17444.269 468526.5 16738.388
## 34 34314.781 383859.1 31262.690
## 35 50260.417 148436.8 36502.442
## 36 14660.930 470304.3 14145.668
## 37 9971.189 566950.4 9736.195
## 38 12336.271 379316.9 11833.854 g. f = 40%res_2sc2_0.4_0.4## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1624062 1610125 1623985 42335533016 9.836363e+11 39879417906
## 2 10 2479069 2461803 2477983 34426746416 9.843533e+11 33157596585
## 3 11 1499455 1161265 1481299 42402911093 9.887748e+11 39848666624
## 4 12 1679679 1621295 1674917 33304424110 9.855595e+11 31966989866
## 5 13 1784151 1643738 1777435 34501616632 9.852589e+11 32816054340
## 6 14 2705845 2648689 2703184 42466700434 9.845007e+11 40526312715
## 7 15 4923058 4579037 4836938 229403573446 9.928755e+11 175926159122
## 8 16 2787417 2501716 2776641 36217997625 9.875094e+11 34661863925
## 9 17 2136820 2751764 2158202 39959962745 9.978198e+11 37630657826
## 10 18 1858921 2053757 1864428 30244348831 9.852152e+11 29153007362
## 11 19 1930788 2451055 1948422 36139125533 9.937404e+11 34642355499
## 12 2 1647160 1446888 1639250 33783322878 9.856569e+11 32381365824
## 13 20 1916950 1524705 1889532 54573260932 9.919048e+11 50884923307
## 14 21 1586760 2079890 1599018 26849682209 9.926638e+11 26084578460
## 15 22 2080393 1960247 2075459 36035606567 9.849253e+11 34423824810
## 16 23 2160809 1982158 2153334 34084617398 9.854719e+11 32790067700
## 17 24 2635495 2836815 2645730 48129982434 9.854321e+11 44782347759
## 18 25 4086344 3748657 4020934 201716689324 9.901776e+11 162608681751
## 19 26 2085978 1456431 2059805 41339732043 9.977119e+11 39465632263
## 20 27 1928809 1464905 1912923 28015910883 9.918343e+11 27004608915
## 21 28 1663028 1785490 1659873 37586489296 9.874300e+11 35410787902
## 22 29 1780064 1751559 1778970 35985423406 9.835700e+11 34313758940
## 23 3 1877959 2710853 1919551 57989860716 1.008321e+12 53757811082
## 24 4 2203916 2451725 2210206 48921943958 9.875830e+11 45285721447
## 25 5 2042554 2137493 2040694 42670728339 9.861225e+11 40303670622
## 26 6 2475236 2388919 2466967 54810855534 9.847339e+11 51203662840
## 27 7 2652497 3575790 2680534 33000681427 1.014498e+12 31779005623
## 28 71 2186213 2272527 2175217 70247098472 9.882858e+11 62760135127
## 29 72 1902369 2042838 1888999 104916372028 9.891450e+11 88221234798
## 30 73 3420300 3456732 3417536 122581130222 9.863301e+11 106975723871
## 31 74 2242295 2002173 2228977 41213679501 9.868512e+11 39157096787
## 32 75 1767592 2073931 1778347 40762320245 9.871782e+11 38278552423
## 33 76 2129257 2017940 2121342 59848770554 9.857352e+11 55505942009
## 34 77 2319782 2440900 2313409 108731049788 9.873413e+11 91666179605
## 35 78 7156354 6497941 6898771 642665400534 9.998544e+11 383979669711
## 36 79 2131160 2004951 2120419 50042038306 9.848179e+11 46339851121
## 37 8 2016835 1720347 2006115 31080297793 9.875237e+11 30012021061
## 38 9 2331335 2519657 2335460 42902457571 9.862735e+11 40252974319
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 25611.25 609221.1 24152.72
## 2 13880.41 398912.3 13376.04
## 3 26632.60 872330.2 25547.47
## 4 18952.95 613033.8 18289.46
## 5 18817.85 601990.0 18001.73
## 6 15633.40 371301.5 14932.69
## 7 44987.05 216057.6 35530.06
## 8 12918.46 398294.9 12415.81
## 9 18005.19 362762.3 16802.87
## 10 15922.60 484495.9 15317.67
## 11 18532.96 405462.9 17598.90
## 12 20008.90 681201.6 19304.82
## 13 26923.32 661311.8 25652.73
## 14 16714.03 476334.8 16116.91
## 15 17076.81 509012.5 16360.05
## 16 15632.94 506420.9 15111.20
## 17 18075.89 350149.6 16752.27
## 18 47122.79 264153.9 39001.78
## 19 19685.80 691237.1 19069.50
## 20 14149.81 679421.7 13785.56
## 21 21295.24 554141.2 20238.75
## 22 19545.06 568595.0 18671.17
## 23 29521.86 370004.0 26671.83
## 24 21265.55 402296.0 19699.69
## 25 19466.77 459652.0 18496.25
## 26 21595.22 413029.7 20302.62
## 27 12302.04 281764.8 11728.03
## 28 29399.01 437003.9 26783.40
## 29 48432.77 488161.8 41841.98
## 30 35700.78 284761.8 31192.58
## 31 17906.39 496791.8 17155.03
## 32 22342.23 479746.7 20856.84
## 33 27495.20 496984.7 25662.11
## 34 43954.47 405745.3 37664.07
## 35 88278.69 152421.9 54941.09
## 36 22636.78 495490.2 21195.17
## 37 15207.32 573791.9 14781.00
## 38 17354.92 390304.0 16316.77 h. f = 30%res_2sc2_0.3_0.3## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1624639 1579174 1620746 44530005096 916186184673 42197959363
## 2 10 2387273 2418662 2379291 61798258333 918092429610 56193022544
## 3 11 1403232 1152388 1390029 35891732608 918118650267 33407072192
## 4 12 1739856 1639604 1727747 58517903647 918526088843 53472244701
## 5 13 1611529 1491111 1605575 30236250427 916310798780 28460129713
## 6 14 2548011 2614995 2550696 66906337608 916542135437 61692225804
## 7 15 4948970 4618411 4858169 273771467164 925580437181 206382544193
## 8 16 2751068 2535160 2736913 54106458913 918083171994 50502146937
## 9 17 2039624 2736983 2067405 41675011607 933733961635 39604654734
## 10 18 1777182 2077274 1788488 37424206177 919044310776 35316690263
## 11 19 1830241 2423691 1850744 49658855560 929981916150 45884989045
## 12 2 1532007 1395422 1526181 41127863023 916520532527 38658086511
## 13 20 1981689 1544907 1929995 88611235021 924816841475 77301913781
## 14 21 1562061 2045099 1575773 33533220943 924848819445 31942216002
## 15 22 2040182 1934367 2033367 48732459235 916278592294 45793815473
## 16 23 2294484 2066911 2273351 77115449679 919911829459 69667890916
## 17 24 2688075 2816120 2683897 104334055221 918666697218 85198012040
## 18 25 4063702 3740971 3954898 294252991209 924249629913 206938813806
## 19 26 2025403 1415611 1985935 57340378222 930412940334 53552242481
## 20 27 1733731 1446834 1725227 23537765928 918857053791 22685264852
## 21 28 1622626 1692880 1625856 70943377460 915796136436 64075096857
## 22 29 1987797 1966309 1986421 69293448003 915146745823 61176184105
## 23 3 1874390 2627257 1928033 94538842659 936308715446 80061311880
## 24 4 2191937 2435005 2200750 68889164498 918395405402 60484300444
## 25 5 2183429 2084299 2169605 77498505452 917067408516 70320150282
## 26 6 2551732 2350061 2524533 121238635126 918292960176 103063674424
## 27 7 2740455 3492120 2783539 57508983492 936777006571 53815738030
## 28 71 2268367 2285092 2258576 120700404326 917306491794 99040498310
## 29 72 1880086 2054281 1849811 117976109176 923903681282 88584891198
## 30 73 3353908 3337212 3331585 210601114063 918385542264 163284679436
## 31 74 2210732 2018395 2193384 74348280508 918093260044 67974745395
## 32 75 2062440 2092900 2045331 180749935651 917997706004 139251324186
## 33 76 2136951 2064367 2111737 114493582725 917671927777 92516460952
## 34 77 2256702 2377777 2262769 185756627976 919077440958 143995027567
## 35 78 6734125 6181152 6469802 841366598898 926358166522 431531780989
## 36 79 2019687 2022049 2017160 51070234588 916391638605 47879441579
## 37 8 1950235 1691770 1930723 50301768298 918805769785 46242608076
## 38 9 2258358 2400294 2257980 67516190632 917246728582 55784225550
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 26643.60 581369.2 25365.50
## 2 24335.78 380808.8 22389.38
## 3 23352.02 803318.9 22251.42
## 4 31701.48 567473.7 29339.60
## 5 17386.87 617748.3 16562.46
## 6 25917.50 350900.4 23893.32
## 7 54507.51 199829.2 42009.57
## 8 19342.62 361843.0 18174.92
## 9 20181.47 341120.3 18933.80
## 10 20650.98 444510.3 19340.25
## 11 25212.06 382430.6 23187.08
## 12 26162.32 665388.4 24771.41
## 13 40889.32 606392.5 37331.53
## 14 20502.00 452272.9 19375.52
## 15 23431.12 474268.1 22130.20
## 16 32839.16 449691.9 30099.75
## 17 34700.73 324799.4 29201.84
## 18 68254.01 246010.0 50312.74
## 19 27620.59 664363.4 26424.95
## 20 13368.59 643930.4 13011.03
## 21 42209.72 544132.3 38171.56
## 22 33297.22 476470.0 29567.81
## 23 46150.30 356253.5 38138.91
## 24 29637.53 377254.6 26163.52
## 25 34034.69 441910.8 31244.22
## 26 46341.69 394730.0 40007.64
## 27 20921.99 267428.4 19249.17
## 28 51112.28 401129.2 42637.08
## 29 50445.52 450069.8 41801.22
## 30 58420.65 274074.7 46573.26
## 31 32883.58 455709.9 30363.69
## 32 80835.47 438982.0 64107.16
## 33 49188.94 447399.7 41194.03
## 34 80531.14 385945.1 62317.48
## 35 123838.68 149036.6 66234.12
## 36 24537.71 453920.8 23086.49
## 37 24958.20 545305.3 23369.74
## 38 25259.95 383969.9 21502.51 i. f = 20%res_2sc2_0.2_0.2## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 1691223 1574218 1673958 7.577977e+10 916386345572 67575792177
## 2 10 2293102 2301545 2287710 1.084384e+11 915489919015 87622801863
## 3 11 1346176 1158896 1335428 5.313367e+10 916949333998 49567807606
## 4 12 1492624 1572130 1492546 7.399138e+10 916342203303 61820129713
## 5 13 1734708 1636610 1716846 1.236993e+11 917271329927 100284688872
## 6 14 2868064 2565653 2805841 2.055454e+11 919208516947 150145412328
## 7 15 4848865 4427208 4679411 4.542614e+11 927120434307 291185156195
## 8 16 2643377 2503426 2578429 1.862958e+11 921454022447 134802748750
## 9 17 2107456 2698580 2135559 1.210220e+11 934295863562 94165450401
## 10 18 1931667 2066745 1933382 1.255895e+11 917968588968 104149473250
## 11 19 1866043 2419891 1903801 1.231638e+11 927811963391 101518144704
## 12 2 1737344 1476247 1695437 1.299916e+11 918475805566 102980725503
## 13 20 1841945 1531704 1773866 1.369415e+11 921186004143 101893231257
## 14 21 1494139 2057788 1538261 8.490444e+10 926272923295 74970492795
## 15 22 2068179 2028447 2048522 1.097180e+11 919219675746 91471715614
## 16 23 2004740 2020789 1998590 8.650822e+10 918948246764 75146470181
## 17 24 2724481 2785793 2698966 2.659311e+11 917763157243 148942223614
## 18 25 3689118 3582622 3658210 3.147385e+11 917397751844 229003973434
## 19 26 1978657 1488345 1929275 8.841907e+10 925152043602 79381559442
## 20 27 1709439 1456432 1692824 5.150174e+10 919398092629 47826815411
## 21 28 1559796 1733618 1558209 6.587145e+10 919838713964 58811638208
## 22 29 1942713 1909214 1938862 1.531271e+11 914908738390 110920705182
## 23 3 1875462 2641669 1951207 1.134837e+11 937972843151 96246853715
## 24 4 2177980 2411545 2192406 7.379274e+10 919459469365 66423266007
## 25 5 2170758 2095684 2144348 1.415214e+11 917977163838 115422985467
## 26 6 2569195 2371119 2507841 2.539411e+11 920046009103 186243345456
## 27 7 2673686 3462306 2755825 1.313935e+11 940689391398 109417766527
## 28 71 2220466 2224050 2154568 2.530707e+11 921890907659 158640320818
## 29 72 1920737 2005012 1879214 1.792261e+11 924076092495 132572212302
## 30 73 3692267 3415619 3477402 6.130850e+11 936680993094 308510662210
## 31 74 2212960 1954435 2178248 9.688368e+10 921294326858 84920672980
## 32 75 1837394 2027898 1792507 3.067533e+11 923565471561 166017445174
## 33 76 2113056 2037668 2078184 1.823575e+11 920386059854 137071162902
## 34 77 2137946 2396144 2133245 3.191715e+11 924827961029 191044266144
## 35 78 6633325 6224540 6419633 1.247178e+12 927285344056 497150313527
## 36 79 2085323 1965400 2037173 1.625899e+11 919024224178 125493992606
## 37 8 2075568 1644202 2009639 1.245285e+11 926334113515 103595432859
## 38 9 2349878 2446664 2356802 1.173966e+11 915634004783 99798552968
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 43092.06 586407.5 39147.09
## 2 46162.91 403551.6 37690.27
## 3 38532.87 826040.6 36318.42
## 4 45609.69 597541.7 39098.07
## 5 63959.34 577537.0 53111.00
## 6 71223.39 360325.9 53651.22
## 7 90619.58 209249.7 61306.46
## 8 63289.84 372587.1 49373.75
## 9 50421.51 348088.0 40113.09
## 10 61056.36 450688.2 51189.62
## 11 59179.86 383155.2 48353.08
## 12 68253.85 640518.7 57172.64
## 13 67830.51 622905.1 55600.68
## 14 53440.14 450853.0 45763.27
## 15 47731.03 463274.3 41003.54
## 16 40889.66 482778.5 35961.94
## 17 84096.60 333800.3 52369.89
## 18 83615.43 256424.5 61776.61
## 19 43275.07 635838.1 40131.82
## 20 28761.86 643400.8 27177.80
## 21 38344.24 534541.9 34792.95
## 22 69953.13 500737.2 52123.95
## 23 56806.68 354606.3 46269.98
## 24 33240.73 383064.3 29780.71
## 25 61760.76 440590.1 51804.94
## 26 96740.68 393217.5 73389.94
## 27 47513.13 270360.4 38701.93
## 28 96099.42 413884.0 66653.96
## 29 78586.44 460751.8 62166.72
## 30 143513.88 272411.8 83914.40
## 31 42246.58 472625.5 37910.04
## 32 130948.83 455942.6 80874.25
## 33 79198.12 461667.6 61174.96
## 34 131284.18 383952.0 84685.16
## 35 185418.86 146985.1 76631.31
## 36 71729.67 468321.9 58030.99
## 37 54403.50 571444.5 47705.61
## 38 47557.49 373388.0 40597.15