title: “skripsweet3” author: “Nadra” date: “March 26, 2018” output: word_document
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"))#Agregasi variabel X kedua
x1 <- aggregate(Xnew_ya$X1, list(Xnew_ya$R102),sum,na.rm=TRUE)
x2 <- aggregate(Xnew_ya$X2, list(Xnew_ya$R102),sum,na.rm=TRUE)
x3 <- aggregate(Xnew_ya$X3, list(Xnew_ya$R102),sum,na.rm=TRUE)
x4 <- aggregate(Xnew_ya$X4, list(Xnew_ya$R102),sum,na.rm=TRUE)
x5 <- aggregate(Xnew_ya$X5, list(Xnew_ya$R102),sum,na.rm=TRUE)
x6 <- aggregate(Xnew_ya$X6, list(Xnew_ya$R102),sum,na.rm=TRUE)
x7 <- aggregate(Xnew_ya$X7, list(Xnew_ya$R102),sum,na.rm=TRUE)
x8 <- aggregate(Xnew_ya$X8, list(Xnew_ya$R102),sum,na.rm=TRUE)
x9 <- aggregate(Xnew_ya$X9, list(Xnew_ya$R102),sum,na.rm=TRUE)
x10 <- aggregate(Xnew_ya$X10, list(Xnew_ya$R102),sum,na.rm=TRUE)
x11 <- aggregate(Xnew_ya$X11, list(Xnew_ya$R102),sum,na.rm=TRUE)
x12 <- aggregate(Xnew_ya$X12, 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)
x15 <- aggregate(Xnew_ya$X15, list(Xnew_ya$R102),sum,na.rm=TRUE)
x16 <- aggregate(Xnew_ya$X16, list(Xnew_ya$R102),sum,na.rm=TRUE)
x17 <- aggregate(Xnew_ya$X17, 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)
x22 <- aggregate(Xnew_ya$X22, list(Xnew_ya$R102),sum,na.rm=TRUE)
x23 <- aggregate(Xnew_ya$X23, 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)
x26 <- aggregate(Xnew_ya$X26, list(Xnew_ya$R102),sum,na.rm=TRUE)
x27 <- aggregate(Xnew_ya$X27, list(Xnew_ya$R102),sum,na.rm=TRUE)
x28 <- aggregate(Xnew_ya$X28, list(Xnew_ya$R102),sum,na.rm=TRUE)
x29 <- aggregate(Xnew_ya$X29, list(Xnew_ya$R102),sum,na.rm=TRUE)
x30 <- aggregate(Xnew_ya$X30, list(Xnew_ya$R102),sum,na.rm=TRUE)
x31 <- aggregate(Xnew_ya$X31, list(Xnew_ya$R102),sum,na.rm=TRUE)
x32 <- aggregate(Xnew_ya$X32, list(Xnew_ya$R102),sum,na.rm=TRUE)
x33 <- aggregate(Xnew_ya$X33, list(Xnew_ya$R102),sum,na.rm=TRUE)
x34 <- aggregate(Xnew_ya$X34, list(Xnew_ya$R102),sum,na.rm=TRUE)
x35 <- aggregate(Xnew_ya$X35, list(Xnew_ya$R102),sum,na.rm=TRUE)
# variabel terpilih : x7,x13,x14,x18,x19,x20,x21,x24,x25,x34
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)
dt1 <- data.frame(cbind("bs_o"=Jatim_pengeluaran$`No Kode Sampel`,"pop"=Jatim_pengeluaran$Exp_krb,"kab"=Jatim_pengeluaran$R102))#1. Without weight
#Pakai rumus awal
resampling = function(data, dataX, proporsi,iterasi ){
myDataY = data
prop = proporsi
propinsi <- function(data, prop){
cb_sampling <- function(data, prop){
mydata= data
bs_u <-unique(mydata$bs_o)
bs_u <- data.frame(bs_u)
samp_bs <- sample(1:nrow(bs_u), size = prop*nrow(bs_u), replace = FALSE)
sampi= bs_u[samp_bs,]
sampi2 = as.matrix(sampi)
dat = as.matrix(mydata)
hasil = c()
for(i in 1: length(sampi2)){
hasil = rbind(hasil,dat[dat[,1]==sampi[i],])
}
hasil <- data.frame(hasil)
hasil$pop = as.numeric(hasil$pop)
est_ybar_1cs <- function(data,f){
hs_sam <- data
Mi <- data.frame(table(hs_sam$bs_o))
y_bar_i <- aggregate(hs_sam$pop, by=list(hs_sam$bs_o), FUN=mean)
y_bar_n <- (sum(y_bar_i$x*Mi$Freq))/sum(Mi$Freq)
M_bar <- sum(Mi$Freq)/nrow(Mi)
pnyb <- sum((Mi$Freq^2)*((y_bar_i$x-rep(c(y_bar_n),nrow(Mi)))^2))
pmbg <- nrow(Mi)*(M_bar)^2*(nrow(Mi)-1)
var_y_bar_n <- (1-f)*pnyb/pmbg
Y_cap_est <- (nrow(Mi)/f)*M_bar*y_bar_n
Var_cap_est <- ((nrow(Mi)/f)*M_bar)^2*var_y_bar_n
rel_var_y <- Var_cap_est/Y_cap_est
hasil <- cbind(Y_cap_est,Var_cap_est,rel_var_y)
return(hasil)
}
result = est_ybar_1cs(hasil, prop)
}
r102 = unique(dt1$kab)
datakoeh = list()
for(b in r102){
datakoeh[[b]] <- subset(dt1, kab %in% b)
}
try_ye <- list()
for(c in r102){
try_ye[[c]] <- cb_sampling(data = datakoeh[[c]], prop = prop)
}
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]] = propinsi(myDataY, prop)
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 kedua #Pakai rumus kedua
resampling_2 <- function(data, dataX, Ni, proporsi,iterasi ){
myDataY <- data
prop <- proporsi
propinsi <- function(data, prop){
cb_sampling <- function(data, prop){
mydata<- data
bs_u <-unique(mydata$bs_o)
bs_u <- data.frame(bs_u)
samp_bs <- sample(1:nrow(bs_u), size = prop*nrow(bs_u), replace = FALSE)
sampi=bs_u[samp_bs,]
sampi2 = as.matrix(sampi)
dat <- as.matrix(mydata)
hasil <- c()
for(i in 1: length(sampi2)){
hasil = rbind(hasil,dat[dat[,1]==sampi[i],])
}
hasil <- data.frame(hasil)
hasil$pop <- as.numeric(hasil$pop)
est_ybar_1cs <- function(data,f){
hs_sam <- data
Mi <- data.frame(table(hs_sam$bs_o))
y_bar_i <- aggregate(hs_sam$pop, by=list(hs_sam$bs_o), FUN=mean)
y_bar_n <- (sum(y_bar_i$x*Mi$Freq))/sum(Mi$Freq)
M_bar <- sum(Mi$Freq)/nrow(Mi)
pnyb <- sum((Mi$Freq^2)*((y_bar_i$x-rep(c(y_bar_n),nrow(Mi)))^2))
pmbg <- nrow(Mi)*(M_bar)^2*(nrow(Mi)-1)
var_y_bar_n <- (1-f)*pnyb/pmbg
Y_cap_est <- (nrow(Mi)/f)*M_bar*y_bar_n
Var_cap_est <- ((nrow(Mi)/f)*M_bar)^2*var_y_bar_n
rel_var_y <- Var_cap_est/Y_cap_est
hasil <- cbind(Y_cap_est,Var_cap_est,rel_var_y)
return(hasil)
}
result <- est_ybar_1cs(hasil, prop)
}
r102 = unique(dt1$kab)
datakoeh = list()
for(b in r102){
datakoeh[[b]] <- subset(dt1, kab %in% b)
}
try_ye <- list()
for(c in r102){
try_ye[[c]] <- cb_sampling(data = datakoeh[[c]], prop = prop)
}
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]] = propinsi(myDataY, prop)
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)
}#2. With weight
#Pakai rumus awal
resampling_weight = function(data, dataX, proporsi,iterasi ){
myDataY = data
prop = proporsi
propinsi <- function(data, prop){
cb_sampling <- function(data, prop){
mydata= data
bs_u <-unique(mydata$bs_o)
bs_u <- data.frame(bs_u)
samp_bs <- sample(1:nrow(bs_u), size = prop*nrow(bs_u), replace = FALSE)
sampi=bs_u[samp_bs,]
sampi2 = as.matrix(sampi)
dat = as.matrix(mydata)
hasil = c()
for(i in 1: length(sampi2)){
hasil = rbind(hasil,dat[dat[,1]==sampi[i],])
}
hasil <- data.frame(hasil)
hasil$pop = as.numeric(hasil$pop)
est_ybar_1cs <- function(data,f){
hs_sam <- data
Mi <- data.frame(table(hs_sam$bs_o))
y_bar_i <- aggregate(hs_sam$pop, by=list(hs_sam$bs_o), FUN=mean)
y_bar_n <- (sum(y_bar_i$x*Mi$Freq))/sum(Mi$Freq)
M_bar <- sum(Mi$Freq)/nrow(Mi)
pnyb <- sum((Mi$Freq^2)*((y_bar_i$x-rep(c(y_bar_n),nrow(Mi)))^2))
pmbg <- nrow(Mi)*(M_bar)^2*(nrow(Mi)-1)
var_y_bar_n <- (1-f)*pnyb/pmbg
Y_cap_est <- (nrow(Mi)/f)*M_bar*y_bar_n
Y_cap_est_weight <- 1/f*Y_cap_est #yang pakai weight
Var_cap_est <- ((nrow(Mi)/f)*M_bar)^2*var_y_bar_n
rel_var_y <- Var_cap_est/Y_cap_est
hasil <- cbind(Y_cap_est,Y_cap_est_weight,Var_cap_est,rel_var_y)
return(hasil)
}
result = est_ybar_1cs(hasil, prop)
}
r102 = unique(dt1$kab)
datakoeh = list()
for(b in r102){
datakoeh[[b]] <- subset(dt1, kab %in% b)
}
try_ye <- list()
for(c in r102){
try_ye[[c]] <- cb_sampling(data = datakoeh[[c]], prop = prop)
}
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]] = propinsi(myDataY, prop)
est_dataset [[j]] = Y_synth_comp(dataX,dataset[[j]]$Y_cap_est_weight,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 kedua#Pakai rumus kedua
resampling_weight_2 <- function(data, dataX, Ni, proporsi,iterasi ){
myDataY <- data
prop <- proporsi
propinsi <- function(data, prop){
cb_sampling <- function(data, prop){
mydata<- data
bs_u <-unique(mydata$bs_o)
bs_u <- data.frame(bs_u)
samp_bs <- sample(1:nrow(bs_u), size = prop*nrow(bs_u), replace = FALSE)
sampi=bs_u[samp_bs,]
sampi2 = as.matrix(sampi)
dat <- as.matrix(mydata)
hasil <- c()
for(i in 1: length(sampi2)){
hasil = rbind(hasil,dat[dat[,1]==sampi[i],])
}
hasil <- data.frame(hasil)
hasil$pop <- as.numeric(hasil$pop)
est_ybar_1cs <- function(data,f){
hs_sam <- data
Mi <- data.frame(table(hs_sam$bs_o))
y_bar_i <- aggregate(hs_sam$pop, by=list(hs_sam$bs_o), FUN=mean)
y_bar_n <- (sum(y_bar_i$x*Mi$Freq))/sum(Mi$Freq)
M_bar <- sum(Mi$Freq)/nrow(Mi)
pnyb <- sum((Mi$Freq^2)*((y_bar_i$x-rep(c(y_bar_n),nrow(Mi)))^2))
pmbg <- nrow(Mi)*(M_bar)^2*(nrow(Mi)-1)
var_y_bar_n <- (1-f)*pnyb/pmbg
Y_cap_est <- (nrow(Mi)/f)*M_bar*y_bar_n
Y_cap_est_weight <- 1/f*Y_cap_est #yang pakai weight
Var_cap_est <- ((nrow(Mi)/f)*M_bar)^2*var_y_bar_n
rel_var_y <- Var_cap_est/Y_cap_est
hasil <- cbind(Y_cap_est,Y_cap_est_weight,Var_cap_est,rel_var_y)
return(hasil)
}
result <- est_ybar_1cs(hasil, prop)
}
r102 = unique(dt1$kab)
datakoeh = list()
for(b in r102){
datakoeh[[b]] <- subset(dt1, kab %in% b)
}
try_ye <- list()
for(c in r102){
try_ye[[c]] <- cb_sampling(data = datakoeh[[c]], prop = prop)
}
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]] = propinsi(myDataY, prop)
est_dataset [[j]] = Y_synth_comp(dataX,dataset[[j]]$Y_cap_est_weight,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)
}Running the data
res_1sc_1 <- resampling(data = dt1, dataX = mat_x, proporsi = 1, iterasi = 10)
res_1sc2_1 <- resampling_2(data = dt1, dataX = mat_x, Ni = Ni_desa, proporsi = 1, iterasi = 10)
res_1sc_0.9 <- resampling(data = dt1, dataX = mat_x, proporsi = 0.9, iterasi = 10)
res_1sc2_0.9 <- resampling_2(data = dt1, dataX = mat_x, Ni = Ni_desa , proporsi = 0.9, iterasi = 10)
res_1sc_0.8 <- resampling(data = dt1, dataX = mat_x, proporsi = 0.8, iterasi = 10)
res_1sc2_0.8 <- resampling_2(data = dt1, dataX = mat_x, Ni = Ni_desa, proporsi = 0.8, iterasi = 10)
res_1sc_0.7 <- resampling(data = dt1, dataX = mat_x, proporsi = 0.7, iterasi = 10)
res_1sc2_0.7 <- resampling_2(data = dt1, dataX = mat_x, Ni = Ni_desa, proporsi = 0.7, iterasi = 10)
res_1sc_0.6 <- resampling(data = dt1, dataX = mat_x, proporsi = 0.6, iterasi = 10)
res_1sc2_0.6 <- resampling_2(data = dt1, dataX = mat_x, Ni = Ni_desa, proporsi = 0.6, iterasi = 10)
res_1sc_0.5 <- resampling(data = dt1, dataX = mat_x, proporsi = 0.5, iterasi = 10)
res_1sc2_0.5 <- resampling_2(data = dt1, dataX = mat_x, Ni = Ni_desa, proporsi = 0.5, iterasi = 10)
ress_1sc_1 <- resampling_weight(data = dt1, dataX = mat_x, proporsi = 1, iterasi = 10)
ress_1sc2_1 <- resampling_weight_2(data = dt1, dataX = mat_x, Ni = Ni_desa, proporsi = 1, iterasi = 10)
ress_1sc_0.9 <- resampling_weight(data = dt1, dataX = mat_x, proporsi = 0.9, iterasi = 10)
ress_1sc2_0.9 <- resampling_weight_2(data = dt1, dataX = mat_x, Ni = Ni_desa , proporsi = 0.9, iterasi = 10)
ress_1sc_0.8 <- resampling_weight(data = dt1, dataX = mat_x, proporsi = 0.8, iterasi = 10)
ress_1sc2_0.8 <- resampling_weight_2(data = dt1, dataX = mat_x, Ni = Ni_desa, proporsi = 0.8, iterasi = 10)
ress_1sc_0.7 <- resampling_weight(data = dt1, dataX = mat_x, proporsi = 0.7, iterasi = 10)
ress_1sc2_0.7 <- resampling_weight_2(data = dt1, dataX = mat_x, Ni = Ni_desa, proporsi = 0.7, iterasi = 10)
ress_1sc_0.6 <- resampling_weight(data = dt1, dataX = mat_x, proporsi = 0.6, iterasi = 10)
ress_1sc2_0.6 <- resampling_weight_2(data = dt1, dataX = mat_x, Ni = Ni_desa, proporsi = 0.6, iterasi = 10)
ress_1sc_0.5 <- resampling_weight(data = dt1, dataX = mat_x, proporsi = 0.5, iterasi = 10)
ress_1sc2_0.5 <- resampling_weight_2(data = dt1, dataX = mat_x, Ni = Ni_desa, proporsi = 0.5, iterasi = 10)Hasil Resampling a. without weight 1. Resampling 100%
res_1sc_1## kab Y dir Y syn Y com MSE dir MSE syn MSE com Relatif MSE dir
## 1 1 226801 222986.8 226801 0 1.454819e+07 0 0
## 2 10 481671 491371.2 481671 0 9.409479e+07 0 0
## 3 11 284635 300358.1 284635 0 2.472163e+08 0 0
## 4 12 316410 289505.1 316410 0 7.238717e+08 0 0
## 5 13 343206 349175.9 343206 0 3.563970e+07 0 0
## 6 14 434778 420501.9 434778 0 2.038059e+08 0 0
## 7 15 461280 450573.8 461280 0 1.146229e+08 0 0
## 8 16 314028 309741.7 314028 0 1.837214e+07 0 0
## 9 17 351541 365175.3 351541 0 1.858928e+08 0 0
## 10 18 317206 333449.3 317206 0 2.638447e+08 0 0
## 11 19 285390 241795.8 285390 0 1.900454e+09 0 0
## 12 2 286903 287012.9 286903 0 1.208006e+04 0 0
## 13 20 258840 245472.0 258840 0 1.787042e+08 0 0
## 14 21 308461 290790.1 308461 0 3.122619e+08 0 0
## 15 22 379756 380949.8 379756 0 1.425268e+06 0 0
## 16 23 350703 360940.2 350703 0 1.048012e+08 0 0
## 17 24 349030 354571.2 349030 0 3.070515e+07 0 0
## 18 25 340725 346014.6 340725 0 2.797959e+07 0 0
## 19 26 278631 305628.1 278631 0 7.288427e+08 0 0
## 20 27 287661 302653.5 287661 0 2.247761e+08 0 0
## 21 28 286903 283530.1 286903 0 1.137645e+07 0 0
## 22 29 349733 350486.5 349733 0 5.678016e+05 0 0
## 23 3 285390 281784.1 285390 0 1.300264e+07 0 0
## 24 4 319600 321717.8 319600 0 4.485242e+06 0 0
## 25 5 341551 347863.8 341551 0 3.985086e+07 0 0
## 26 6 441330 428562.5 441330 0 1.630078e+08 0 0
## 27 7 561270 564932.5 561270 0 1.341375e+07 0 0
## 28 71 163306 152747.8 163306 0 1.114753e+08 0 0
## 29 72 113050 132366.7 113050 0 3.731344e+08 0 0
## 30 73 266815 231109.6 266815 0 1.274876e+09 0 0
## 31 74 157080 145013.1 157080 0 1.456100e+08 0 0
## 32 75 123753 136477.7 123753 0 1.619178e+08 0 0
## 33 76 115440 144147.8 115440 0 8.241364e+08 0 0
## 34 77 121771 125802.1 121771 0 1.625002e+07 0 0
## 35 78 542361 567311.8 542361 0 6.225448e+08 0 0
## 36 79 132355 139290.8 132355 0 4.810515e+07 0 0
## 37 8 318003 317330.1 318003 0 4.527983e+05 0 0
## 38 9 617716 595940.7 617716 0 4.741638e+08 0 0
## Relatif MSE syn Relatif MSE com
## 1 6.524240e+01 0
## 2 1.914943e+02 0
## 3 8.230718e+02 0
## 4 2.500376e+03 0
## 5 1.020681e+02 0
## 6 4.846729e+02 0
## 7 2.543931e+02 0
## 8 5.931437e+01 0
## 9 5.090509e+02 0
## 10 7.912589e+02 0
## 11 7.859749e+03 0
## 12 4.208892e-02 0
## 13 7.280026e+02 0
## 14 1.073840e+03 0
## 15 3.741353e+00 0
## 16 2.903562e+02 0
## 17 8.659796e+01 0
## 18 8.086247e+01 0
## 19 2.384737e+03 0
## 20 7.426846e+02 0
## 21 4.012433e+01 0
## 22 1.620038e+00 0
## 23 4.614397e+01 0
## 24 1.394154e+01 0
## 25 1.145588e+02 0
## 26 3.803595e+02 0
## 27 2.374400e+01 0
## 28 7.298000e+02 0
## 29 2.818945e+03 0
## 30 5.516329e+03 0
## 31 1.004116e+03 0
## 32 1.186405e+03 0
## 33 5.717302e+03 0
## 34 1.291712e+02 0
## 35 1.097359e+03 0
## 36 3.453577e+02 0
## 37 1.426900e+00 0
## 38 7.956560e+02 0 res_1sc2_1## kab Y dir Y syn Y com MSE dir MSE syn MSE com Relatif MSE dir
## 1 1 226801 222986.8 226801 0 14599353633 0 0
## 2 10 481671 491371.2 481671 0 14602096619 0 0
## 3 11 284635 300358.1 284635 0 14607376670 0 0
## 4 12 316410 289505.1 316410 0 14623813064 0 0
## 5 13 343206 349175.9 343206 0 14600080926 0 0
## 6 14 434778 420501.9 434778 0 14605879761 0 0
## 7 15 461280 450573.8 461280 0 14602804483 0 0
## 8 16 314028 309741.7 314028 0 14599485493 0 0
## 9 17 351541 365175.3 351541 0 14605262067 0 0
## 10 18 317206 333449.3 317206 0 14607950065 0 0
## 11 19 285390 241795.8 285390 0 14664384875 0 0
## 12 2 286903 287012.9 286903 0 14598852388 0 0
## 13 20 258840 245472.0 258840 0 14605014186 0 0
## 14 21 308461 290790.1 308461 0 14609619622 0 0
## 15 22 379756 380949.8 379756 0 14598901118 0 0
## 16 23 350703 360940.2 350703 0 14602465807 0 0
## 17 24 349030 354571.2 349030 0 14599910769 0 0
## 18 25 340725 346014.6 340725 0 14599816785 0 0
## 19 26 278631 305628.1 278631 0 14623984478 0 0
## 20 27 287661 302653.5 287661 0 14606602872 0 0
## 21 28 286903 283530.1 286903 0 14599244263 0 0
## 22 29 349733 350486.5 349733 0 14598871550 0 0
## 23 3 285390 281784.1 285390 0 14599300338 0 0
## 24 4 319600 321717.8 319600 0 14599006635 0 0
## 25 5 341551 347863.8 341551 0 14600226139 0 0
## 26 6 441330 428562.5 441330 0 14604472931 0 0
## 27 7 561270 564932.5 561270 0 14599314514 0 0
## 28 71 163306 152747.8 163306 0 14602695949 0 0
## 29 72 113050 132366.7 113050 0 14611718676 0 0
## 30 73 266815 231109.6 266815 0 14642813228 0 0
## 31 74 157080 145013.1 157080 0 14603873006 0 0
## 32 75 123753 136477.7 123753 0 14604435342 0 0
## 33 76 115440 144147.8 115440 0 14627270468 0 0
## 34 77 121771 125802.1 121771 0 14599412316 0 0
## 35 78 542361 567311.8 542361 0 14620319033 0 0
## 36 79 132355 139290.8 132355 0 14600510769 0 0
## 37 8 318003 317330.1 318003 0 14598867585 0 0
## 38 9 617716 595940.7 617716 0 14615202446 0 0
## Relatif MSE syn Relatif MSE com
## 1 65471.83 0
## 2 29717.04 0
## 3 48633.20 0
## 4 50513.14 0
## 5 41812.97 0
## 6 34734.39 0
## 7 32409.35 0
## 8 47134.38 0
## 9 39995.21 0
## 10 43808.61 0
## 11 60647.81 0
## 12 50864.79 0
## 13 59497.69 0
## 14 50241.12 0
## 15 38322.37 0
## 16 40456.74 0
## 17 41176.24 0
## 18 42194.23 0
## 19 47848.95 0
## 20 48261.79 0
## 21 51490.99 0
## 22 41653.16 0
## 23 51810.24 0
## 24 45378.29 0
## 25 41971.10 0
## 26 34077.81 0
## 27 25842.58 0
## 28 95600.03 0
## 29 110388.19 0
## 30 63358.74 0
## 31 100707.27 0
## 32 107009.69 0
## 33 101474.13 0
## 34 116050.60 0
## 35 25771.22 0
## 36 104820.36 0
## 37 46005.30 0
## 38 24524.59 0 2. Resampling 90% res_1sc_0.9## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 202742.9 198505.2 202045.3 2896549 15336488 2396373.5
## 2 10 433848.4 439377.9 435571.5 9293068 21991026 6245777.2
## 3 11 253275.2 267866.1 253558.5 4104925 210415038 4024098.4
## 4 12 285047.7 259335.3 284830.9 5567810 656044366 5520633.2
## 5 13 304140.3 309283.3 305320.5 5982514 20788405 4560119.6
## 6 14 387627.0 375849.0 386550.8 12550986 127574986 11380806.5
## 7 15 414338.9 403440.8 413259.4 11538425 109580748 10350275.3
## 8 16 282903.3 276068.1 282207.4 4724879 42380497 4236910.9
## 9 17 311251.1 325755.2 311662.3 5959340 204581551 5789282.2
## 10 18 285445.8 297955.8 285752.2 3819672 153163165 3725230.9
## 11 19 253650.2 214880.0 253518.5 5107239 1498146478 5089788.5
## 12 2 255229.4 255232.0 255231.4 3546116 3361836 1724513.3
## 13 20 227413.3 219266.5 227050.1 2944371 63817302 2811027.3
## 14 21 277258.3 259161.5 276947.1 5633829 322210758 5536389.7
## 15 22 339909.6 338815.9 339272.8 9179344 7497880 4091005.1
## 16 23 310669.7 323508.0 311233.1 7177815 158943959 6857440.5
## 17 24 309092.8 315051.4 310581.2 8776320 27252029 6541147.9
## 18 25 301844.0 308566.3 303338.6 9648618 37259969 7312202.2
## 19 26 247235.1 271791.6 247515.4 6866347 596894708 6787427.2
## 20 27 255983.0 270591.9 256336.5 5160873 208528587 5035535.5
## 21 28 255078.7 251805.4 253730.8 4571648 6772017 2583616.2
## 22 29 309802.2 310475.7 310149.7 8287884 7814824 4021549.4
## 23 3 253875.7 250747.1 252263.9 4878761 5570792 2094237.2
## 24 4 287680.0 286861.7 287244.5 7053811 6326847 3332278.3
## 25 5 302826.8 310393.9 303493.3 4937735 53560556 4483563.1
## 26 6 393964.7 382842.8 393186.5 8564663 116563923 7952852.9
## 27 7 502340.0 504239.2 503686.5 14386731 9453904 5379741.2
## 28 71 147635.4 135791.7 147477.1 1860017 139484451 1834691.1
## 29 72 100725.3 117525.0 100790.2 1089293 281234154 1085061.2
## 30 73 237766.2 203639.3 237610.0 5328707 1160833492 5304257.9
## 31 74 139166.7 129046.9 138880.5 2889981 99632606 2807522.8
## 32 75 107557.4 121347.9 107633.3 1042327 189783975 1036527.9
## 33 76 102961.1 128547.1 102997.7 935502 653946476 934156.7
## 34 77 106436.2 110096.4 106756.7 1071555 13635809 959914.3
## 35 78 479480.1 503209.5 480170.5 16265513 549829990 15783290.0
## 36 79 115064.4 123866.9 115228.0 1438460 76169424 1411671.6
## 37 8 286162.9 283356.3 283894.6 7311009 1804414 1365077.5
## 38 9 547253.8 530589.3 546240.7 16844535 261336303 15814610.2
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 14.288052 77.277801 11.860528
## 2 21.419807 50.068009 14.346492
## 3 16.207530 785.141808 15.870766
## 4 19.533320 2529.988255 19.382584
## 5 19.668846 67.194473 14.936791
## 6 32.381089 339.371025 29.441780
## 7 27.847582 271.492808 25.040961
## 8 16.701246 153.486575 15.012528
## 9 19.146407 627.975922 18.575384
## 10 13.381463 514.002192 13.036657
## 11 20.134755 6972.104662 20.076436
## 12 13.893845 13.171931 6.756722
## 13 12.946960 291.129743 12.380415
## 14 20.318131 1243.322115 19.989277
## 15 27.004996 22.128425 12.058664
## 16 23.104742 491.097821 22.033705
## 17 28.393459 86.446709 21.062025
## 18 31.962946 120.697655 24.119066
## 19 27.770355 2196.065055 27.420137
## 20 20.161583 770.592380 19.644769
## 21 17.922280 26.919733 10.181548
## 22 26.751680 25.170570 12.966349
## 23 19.217092 22.253727 8.291393
## 24 24.519859 22.054786 11.600842
## 25 16.306158 172.531157 14.775386
## 26 21.739759 304.335636 20.225433
## 27 28.639688 18.753526 10.683100
## 28 12.600390 1027.154746 12.441997
## 29 10.814333 2392.807169 10.765355
## 30 22.410720 5702.469422 22.322572
## 31 20.766329 772.208680 20.215581
## 32 9.691321 1564.054136 9.630649
## 33 9.085974 5086.547132 9.069671
## 34 10.068042 123.612052 8.997215
## 35 33.918428 1093.716098 32.867239
## 36 12.501863 614.833837 12.251634
## 37 25.548795 6.368328 4.805900
## 38 30.779853 492.580583 28.951481 res_1sc2_0.9## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 202742.9 198505.2 202741.8 2896549 11643481916 2895825.6
## 2 10 433848.4 439377.9 433852.9 9293068 11643931953 9285651.9
## 3 11 253275.2 267866.1 253280.4 4104925 11650250431 4103476.6
## 4 12 285047.7 259335.3 285035.4 5567810 11665667404 5565147.0
## 5 13 304140.3 309283.3 304143.0 5982514 11643776326 5979433.3
## 6 14 387627.0 375849.0 387614.3 12550986 11647685121 12537463.4
## 7 15 414338.9 403440.8 414328.1 11538425 11647029714 11526994.4
## 8 16 282903.3 276068.1 282900.6 4724879 11644477514 4722958.7
## 9 17 311251.1 325755.2 311258.5 5959340 11650113222 5956279.1
## 10 18 285445.8 297955.8 285449.9 3819672 11648266392 3818417.9
## 11 19 253650.2 214880.0 253633.3 5107239 11694689526 5104999.6
## 12 2 255229.4 255232.0 255229.4 3546116 11642859545 3545033.4
## 13 20 227413.3 219266.5 227411.3 2944371 11645155317 2943622.5
## 14 21 277258.3 259161.5 277249.6 5633829 11654158177 5631087.6
## 15 22 339909.6 338815.9 339908.7 9179344 11642911172 9172099.6
## 16 23 310669.7 323508.0 310677.6 7177815 11648581528 7173386.3
## 17 24 309092.8 315051.4 309097.3 8776320 11644095548 8769692.8
## 18 25 301844.0 308566.3 301849.6 9648618 11644470728 9640615.6
## 19 26 247235.1 271791.6 247249.5 6866347 11663672538 6862301.9
## 20 27 255983.0 270591.9 255989.5 5160873 11650221793 5158580.5
## 21 28 255078.7 251805.4 255077.4 4571648 11643235454 4569849.1
## 22 29 309802.2 310475.7 309802.7 8287884 11642869503 8281975.4
## 23 3 253875.7 250747.1 253874.4 4878761 11643210425 4876712.0
## 24 4 287680.0 286861.7 287679.5 7053811 11642878259 7049536.6
## 25 5 302826.8 310393.9 302830.0 4937735 11644870373 4935638.2
## 26 6 393964.7 382842.8 393956.5 8564663 11647167970 8558366.9
## 27 7 502340.0 504239.2 502342.3 14386731 11643049423 14368963.1
## 28 71 147635.4 135791.7 147633.6 1860017 11647727138 1859719.8
## 29 72 100725.3 117525.0 100726.9 1089293 11652588482 1089191.0
## 30 73 237766.2 203639.3 237750.6 5328707 11683065680 5326273.4
## 31 74 139166.7 129046.9 139164.2 2889981 11646388452 2889260.8
## 32 75 107557.4 121347.9 107558.7 1042327 11649433408 1042233.6
## 33 76 102961.1 128547.1 102963.2 935502 11665435328 935426.7
## 34 77 106436.2 110096.4 106436.6 1071555 11643360341 1071455.9
## 35 78 479480.1 503209.5 479513.1 16265513 11662373725 16242800.2
## 36 79 115064.4 123866.9 115065.5 1438460 11645529325 1438281.9
## 37 8 286162.9 283356.3 286161.1 7311009 11643129051 7306418.6
## 38 9 547253.8 530589.3 547229.7 16844535 11652445634 16820206.6
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 14.288052 58656.53 14.284558
## 2 21.419807 26500.96 21.402497
## 3 16.207530 43493.69 16.201480
## 4 19.533320 44983.65 19.524818
## 5 19.668846 37647.56 19.658548
## 6 32.381089 30990.44 32.347256
## 7 27.847582 28869.39 27.820722
## 8 16.701246 42179.99 16.694622
## 9 19.146407 35763.32 19.136114
## 10 13.381463 39094.04 13.376876
## 11 20.134755 54424.32 20.127275
## 12 13.893845 45616.91 13.889605
## 13 12.946960 53110.16 12.943787
## 14 20.318131 44968.64 20.308890
## 15 27.004996 34363.61 26.983752
## 16 23.104742 36007.96 23.089895
## 17 28.393459 36959.94 28.371602
## 18 31.962946 37737.92 31.935852
## 19 27.770355 42914.02 27.752377
## 20 20.161583 43054.76 20.152114
## 21 17.922280 46238.90 17.915319
## 22 26.751680 37500.33 26.732565
## 23 19.217092 46433.87 19.209121
## 24 24.519859 40587.01 24.505040
## 25 16.306158 37516.27 16.299059
## 26 21.739759 30422.72 21.724229
## 27 28.639688 23090.24 28.604185
## 28 12.600390 85777.67 12.598536
## 29 10.814333 99151.73 10.813149
## 30 22.410720 57372.35 22.401956
## 31 20.766329 90250.97 20.761531
## 32 9.691321 96002.53 9.690341
## 33 9.085974 90751.32 9.085061
## 34 10.068042 105758.04 10.067081
## 35 33.918428 23176.27 33.868738
## 36 12.501863 94018.32 12.500200
## 37 25.548795 41090.30 25.532912
## 38 30.779853 21961.21 30.736748 3. Resampling 80% res_1sc_0.8## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 179158.62 175290.21 177713.91 5322159 10324741 3112830
## 2 10 384749.25 391126.27 387367.87 16477656 24760754 9512476
## 3 11 221861.50 236618.39 222334.04 6957809 211716581 6732867
## 4 12 252889.62 228837.53 252450.93 10550147 568543571 10357073
## 5 13 273489.25 273380.59 273394.18 11208944 9976202 5244359
## 6 14 341604.75 333608.72 338719.60 21881568 44908242 13011675
## 7 15 369128.88 360454.71 366544.55 20243316 63047755 12882922
## 8 16 251619.62 246121.72 249950.90 8921334 22289750 6058551
## 9 17 279725.62 290101.20 280857.19 11713645 96202703 10423896
## 10 18 253765.00 263998.66 254421.78 6711901 98282371 6279478
## 11 19 222606.50 190783.31 222337.49 8555933 1004373643 8482590
## 12 2 223652.12 225376.25 224907.29 6161111 2964889 1879663
## 13 20 202920.00 193816.08 202339.13 5285092 77712702 4944435
## 14 21 246297.75 229135.30 245674.99 10665525 284285011 10274813
## 15 22 300069.88 300866.18 300508.03 16023132 14562587 7607427
## 16 23 279225.12 286372.79 281055.14 12857730 39227283 9450668
## 17 24 277639.00 282777.59 280734.95 15693841 11159256 5991133
## 18 25 272666.75 277214.07 275776.94 18190288 9868956 5586942
## 19 26 216606.38 240415.17 217115.56 12089974 556104434 11828550
## 20 27 224175.75 237981.66 224856.01 9384909 181597158 8916900
## 21 28 223502.62 222904.63 223122.05 8103627 7014985 3690196
## 22 29 277892.25 278437.11 278168.08 14900557 14580814 7369204
## 23 3 222457.62 221021.34 221552.41 8190352 5767696 3307735
## 24 4 255920.00 254896.13 255375.74 12672186 11462378 6012454
## 25 5 272332.75 274453.05 273831.57 8573372 4383357 2782593
## 26 6 346241.12 338899.02 344187.47 14810349 39894107 10480629
## 27 7 444156.62 448228.17 447230.23 25925417 11011622 7199252
## 28 71 130077.38 120436.02 129756.07 3077105 91262765 2972863
## 29 72 88595.50 103882.62 88723.20 1952580 232062519 1936232
## 30 73 208225.50 180012.02 207912.05 8845466 788014949 8747131
## 31 74 121275.00 114227.67 120560.40 5018921 44761789 4502780
## 32 75 96580.75 107518.99 96749.45 1837925 118096170 1809103
## 33 76 90487.50 112478.78 90563.61 1670692 482308110 1664809
## 34 77 94967.62 97698.63 95777.53 1956878 6316260 1213066
## 35 78 428702.88 448015.24 430189.93 28614720 345857047 26391944
## 36 79 102870.50 109526.14 103271.82 2660601 41795595 2498448
## 37 8 254482.38 251345.12 251905.28 13135237 3725806 2639860
## 38 9 487526.50 471888.81 485633.26 29539233 215493460 25940495
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 29.70201 58.98426 17.490196
## 2 42.82857 63.32299 24.572035
## 3 31.36497 894.59312 30.286401
## 4 41.71476 2484.71609 41.022705
## 5 40.98667 36.47934 19.180720
## 6 64.06263 134.63456 38.325989
## 7 54.84136 174.62107 35.017455
## 8 35.45519 90.60215 24.222743
## 9 41.87533 331.58110 37.114297
## 10 26.44928 372.24574 24.681483
## 11 38.43832 5264.53401 38.155494
## 12 27.54715 13.16055 8.361348
## 13 26.04441 401.01532 24.436257
## 14 43.30539 1240.62765 41.824870
## 15 53.39785 48.40155 25.313389
## 16 46.05130 136.92975 33.640599
## 17 56.52478 39.41421 21.353990
## 18 66.71290 35.58177 20.289149
## 19 55.81344 2313.14608 54.478615
## 20 41.86774 762.97546 39.658467
## 21 36.25790 31.44821 16.531796
## 22 53.61898 52.36596 26.491465
## 23 36.81867 26.08053 14.924864
## 24 49.51687 44.96771 23.543927
## 25 31.48164 15.97766 10.164087
## 26 42.78020 117.64642 30.436962
## 27 58.37256 24.57987 16.106344
## 28 23.63938 757.54935 22.893952
## 29 22.04345 2233.55897 21.827358
## 30 42.47486 4378.43404 42.066045
## 31 41.38463 391.99915 37.350128
## 32 19.02638 1098.40519 18.695364
## 33 18.46323 4286.84793 18.382552
## 34 20.61227 64.63241 12.729578
## 35 66.75133 772.36638 61.357558
## 36 25.86461 381.53617 24.194050
## 37 51.61875 14.81681 10.468421
## 38 60.58973 456.69802 53.415319 res_1sc2_0.8## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 179158.62 175290.21 179156.40 5322159 9277548285 5319086
## 2 10 384749.25 391126.27 384760.57 16477656 9278430751 16448328
## 3 11 221861.50 236618.39 221872.57 6957809 9284549233 6952564
## 4 12 252889.62 228837.53 252862.34 10550147 9296977486 10538145
## 5 13 273489.25 273380.59 273489.16 11208944 9277051245 11195386
## 6 14 341604.75 333608.72 341585.97 21881568 9279311834 21829915
## 7 15 369128.88 360454.71 369109.97 20243316 9279880843 20199095
## 8 16 251619.62 246121.72 251614.32 8921334 9278084981 8912721
## 9 17 279725.62 290101.20 279738.70 11713645 9280729991 11698812
## 10 18 253765.00 263998.66 253772.40 6711901 9280629229 6707040
## 11 19 222606.50 190783.31 222577.31 8555933 9311937343 8547985
## 12 2 223652.12 225376.25 223653.27 6161111 9277118952 6156970
## 13 20 202920.00 193816.08 202914.81 5285092 9279870730 5282054
## 14 21 246297.75 229135.30 246278.08 10665525 9287179445 10653237
## 15 22 300069.88 300866.18 300071.25 16023132 9277059101 15995481
## 16 23 279225.12 286372.79 279235.05 12857730 9278804772 12839884
## 17 24 277639.00 282777.59 277647.70 15693841 9277934706 15667280
## 18 25 272666.75 277214.07 272675.64 18190288 9277788785 18154634
## 19 26 216606.38 240415.17 216637.28 12089974 9296601648 12074195
## 20 27 224175.75 237981.66 224189.71 9384909 9283594326 9375328
## 21 28 223502.62 222904.63 223502.11 8103627 9277046277 8096506
## 22 29 277892.25 278437.11 277893.13 14900557 9277019763 14876578
## 23 3 222457.62 221021.34 222456.37 8190352 9277092277 8183082
## 24 4 255920.00 254896.13 255918.61 12672186 9277050455 12654861
## 25 5 272332.75 274453.05 272334.71 8573372 9277183426 8565413
## 26 6 346241.12 338899.02 346229.45 14810349 9278895098 14786615
## 27 7 444156.62 448228.17 444167.98 25925417 9277637640 25852974
## 28 71 130077.38 120436.02 130074.15 3077105 9280261836 3076073
## 29 72 88595.50 103882.62 88598.72 1952580 9285078223 1952167
## 30 73 208225.50 180012.02 208198.66 8845466 9304486683 8837043
## 31 74 121275.00 114227.67 121271.20 5018921 9278725313 5016189
## 32 75 96580.75 107518.99 96582.91 1837925 9281144396 1837558
## 33 76 90487.50 112478.78 90491.45 1670692 9293697661 1670388
## 34 77 94967.62 97698.63 94968.20 1956878 9277294018 1956462
## 35 78 428702.88 448015.24 428762.21 28614720 9289921557 28526695
## 36 79 102870.50 109526.14 102872.41 2660601 9278541709 2659837
## 37 8 254482.38 251345.12 254477.94 13135237 9277358375 13116576
## 38 9 487526.50 471888.81 487476.92 29539233 9285458140 29445411
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 29.70201 52928.04 29.68524
## 2 42.82857 23722.29 42.75108
## 3 31.36497 39238.99 31.33975
## 4 41.71476 40627.18 41.67182
## 5 40.98667 33934.53 40.93710
## 6 64.06263 27814.80 63.91491
## 7 54.84136 25744.85 54.72437
## 8 35.45519 37697.45 35.42171
## 9 41.87533 31991.22 41.82034
## 10 26.44928 35154.16 26.42935
## 11 38.43832 48809.14 38.40771
## 12 27.54715 41163.18 27.52849
## 13 26.04441 47880.24 26.03011
## 14 43.30539 40531.43 43.25897
## 15 53.39785 30834.60 53.30546
## 16 46.05130 32401.75 45.98573
## 17 56.52478 32810.36 56.42734
## 18 66.71290 33468.14 66.57996
## 19 55.81344 38669.09 55.73261
## 20 41.86774 39009.83 41.82236
## 21 36.25790 41619.14 36.22612
## 22 53.61898 33318.38 53.53252
## 23 36.81867 41973.69 36.78619
## 24 49.51687 36395.33 49.44944
## 25 31.48164 33802.27 31.45218
## 26 42.78020 27379.45 42.71307
## 27 58.37256 20698.35 58.20794
## 28 23.63938 77057.14 23.63205
## 29 22.04345 89383.05 22.03799
## 30 42.47486 51688.66 42.43991
## 31 41.38463 81232.16 41.36341
## 32 19.02638 86323.63 19.02216
## 33 18.46323 82630.51 18.45906
## 34 20.61227 94961.82 20.60775
## 35 66.75133 20735.65 66.53674
## 36 25.86461 84717.63 25.85670
## 37 51.61875 36911.03 51.54631
## 38 60.58973 19677.14 60.40344 4. Resampling 70% res_1sc_0.7## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 155045.29 153800.31 154218.48 6954616 4743247 2719589
## 2 10 337514.43 342379.06 341239.04 22431652 7647797 4997572
## 3 11 198529.57 207676.58 199596.20 9672909 74819683 8517231
## 4 12 221606.71 201704.72 220932.58 13432871 383312135 12974030
## 5 13 233785.29 240163.36 235940.47 13542047 27739121 8766923
## 6 14 303089.57 290951.77 300545.29 30265416 119205691 23571663
## 7 15 320262.29 311345.45 317098.19 27144723 55733840 16713041
## 8 16 220494.43 214135.29 218680.67 11524474 29105896 8215928
## 9 17 239706.14 253489.17 240725.25 14021460 176175762 12964854
## 10 18 222483.29 231595.23 223436.01 8569218 75339847 7633392
## 11 19 198378.00 166962.71 198004.70 11728571 975500756 11587168
## 12 2 199207.14 197989.19 198514.87 8360667 6732383 3710623
## 13 20 178285.71 170882.17 177285.79 7339708 47804797 6320118
## 14 21 216030.29 201592.95 215139.54 12775069 196870252 11958032
## 15 22 261187.86 261710.21 261487.31 21294772 19595511 10166549
## 16 23 238714.86 249015.19 240366.49 16942707 89541625 14197936
## 17 24 237724.29 242292.14 241108.14 19421438 8112222 4914402
## 18 25 232975.00 237711.59 236027.03 22266252 13955098 7267296
## 19 26 194414.00 211548.89 195383.79 16303262 281835329 15338469
## 20 27 200264.71 210298.94 201443.41 11816053 88978226 10407700
## 21 28 199584.57 196982.34 197666.45 11264669 4883945 3211512
## 22 29 238098.00 238722.58 238413.78 18773759 18376816 9286498
## 23 3 198077.29 194470.42 195403.16 11548273 4334317 2782356
## 24 4 224000.00 223113.85 223513.32 16352727 14935694 7782305
## 25 5 232480.29 241911.22 233727.71 11568061 78703827 9974919
## 26 6 307291.29 298454.68 305062.61 19537992 59351954 14496020
## 27 7 385881.14 390661.15 389050.79 33126843 20578271 11931396
## 28 71 114532.71 106174.39 113992.41 4228771 72250369 3921947
## 29 72 76519.00 92180.95 76682.17 2552894 242928308 2525908
## 30 73 185937.86 158334.50 185480.34 12607804 751427473 12396840
## 31 74 108921.43 101225.39 108017.67 6925858 52514747 6100584
## 32 75 85369.00 95084.76 85642.04 2584289 93734393 2506565
## 33 76 78021.43 99634.38 78127.24 2285941 465140263 2274627
## 34 77 83797.86 86670.96 84556.77 2628376 9341921 1626871
## 35 78 374831.00 392670.23 377008.48 38714974 280283234 33936189
## 36 79 91138.14 97140.49 91719.62 3434729 32994900 3085256
## 37 8 222722.00 220554.64 221378.21 17100765 11994859 6913855
## 38 9 427802.43 413472.42 425153.79 37775363 168422635 30672238
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 44.85426 30.82338 17.63115
## 2 66.46439 22.33727 14.67318
## 3 48.73628 360.14831 42.68502
## 4 60.60581 1900.56919 58.71606
## 5 57.93889 115.39419 37.17048
## 6 99.88496 409.59391 78.41954
## 7 84.75476 179.11548 52.62602
## 8 52.25767 135.91285 37.56265
## 9 58.49315 694.95364 53.85359
## 10 38.51145 325.24847 34.16084
## 11 59.11861 5843.08581 58.51757
## 12 41.97277 34.00429 18.69340
## 13 41.17064 279.90118 35.65205
## 14 59.13256 976.51828 55.58242
## 15 81.52177 74.89565 38.88162
## 16 70.97654 359.59093 59.07334
## 17 81.70222 33.44473 20.41394
## 18 95.58263 58.78820 30.91861
## 19 83.84718 1332.38354 78.50876
## 20 59.00525 423.10936 51.66598
## 21 56.43621 24.78269 16.23612
## 22 78.85695 76.98783 38.95519
## 23 58.29237 22.34778 14.23437
## 24 73.00530 66.92816 34.81717
## 25 49.76427 325.34335 42.69287
## 26 63.57202 198.75234 47.49430
## 27 85.84627 52.65929 30.67442
## 28 36.85918 679.22058 34.31060
## 29 33.35265 2635.33909 32.92935
## 30 67.76691 4745.44968 66.79792
## 31 63.58582 519.12509 56.47941
## 32 30.25628 985.69532 29.25711
## 33 29.29889 4667.27677 29.11399
## 34 31.38108 107.41172 19.38817
## 35 103.28147 714.09031 90.02112
## 36 37.68895 339.60161 33.64392
## 37 76.78445 54.36782 31.22454
## 38 88.29142 407.37521 72.13551 res_1sc2_0.7## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 155045.29 153800.31 155044.03 6954616 6992781355 6947663
## 2 10 337514.43 342379.06 337530.06 22431652 6993556415 22359249
## 3 11 198529.57 207676.58 198542.22 9672909 6995618639 9659434
## 4 12 221606.71 201704.72 221568.52 13432871 7006385963 13406797
## 5 13 233785.29 240163.36 233797.60 13542047 6994128590 13515587
## 6 14 303089.57 290951.77 303037.61 30265416 6997859277 30134099
## 7 15 320262.29 311345.45 320227.77 27144723 6995562983 27039094
## 8 16 220494.43 214135.29 220483.93 11524474 6994106148 11505354
## 9 17 239706.14 253489.17 239733.68 14021460 6999263626 13993015
## 10 18 222483.29 231595.23 222494.41 8569218 6995598517 8558575
## 11 19 198378.00 166962.71 198325.63 11728571 7026747492 11708733
## 12 2 199207.14 197989.19 199205.70 8360667 6992761249 8350571
## 13 20 178285.71 170882.17 178277.97 7339708 6994606636 7331948
## 14 21 216030.29 201592.95 216004.01 12775069 6999934250 12751257
## 15 22 261187.86 261710.21 261189.39 21294772 6992763696 21229790
## 16 23 238714.86 249015.19 238739.73 16942707 6996376975 16901704
## 17 24 237724.29 242292.14 237736.88 19421438 6993469673 19367511
## 18 25 232975.00 237711.59 232990.04 22266252 6993582073 22194685
## 19 26 194414.00 211548.89 194453.68 16303262 7002985742 16265268
## 20 27 200264.71 210298.94 200281.64 11816053 6996180766 11795878
## 21 28 199584.57 196982.34 199580.36 11264669 6992956127 11246467
## 22 29 238098.00 238722.58 238099.67 18773759 6992718789 18723098
## 23 3 198077.29 194470.42 198071.32 11548273 6993171257 11529030
## 24 4 224000.00 223113.85 223997.95 16352727 6992753964 16314271
## 25 5 232480.29 241911.22 232495.86 11568061 6995817925 11548882
## 26 6 307291.29 298454.68 307266.60 19537992 6995425444 19483004
## 27 7 385881.14 390661.15 385903.59 33126843 6993668580 32970453
## 28 71 114532.71 106174.39 114527.51 4228771 6995342313 4226180
## 29 72 76519.00 92180.95 76524.71 2552894 7001169970 2551952
## 30 73 185937.86 158334.50 185888.32 12607804 7019051145 12585048
## 31 74 108921.43 101225.39 108913.82 6925858 6994754777 6918973
## 32 75 85369.00 95084.76 85372.57 2584289 6996026435 2583324
## 33 76 78021.43 99634.38 78028.48 2285941 7008823246 2285188
## 34 77 83797.86 86670.96 83798.93 2628376 6993088339 2627385
## 35 78 374831.00 392670.23 374929.02 38714974 7003705039 38501499
## 36 79 91138.14 97140.49 91141.08 3434729 6993961295 3433028
## 37 8 222722.00 220554.64 222716.75 17100765 6992881167 17058828
## 38 9 427802.43 413472.42 427725.55 37775363 6999815377 37571285
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 44.85426 45467.53 44.80979
## 2 66.46439 20426.47 66.24673
## 3 48.73628 33685.64 48.66523
## 4 60.60581 34736.07 60.49879
## 5 57.93889 29122.66 57.82258
## 6 99.88496 24051.70 99.46854
## 7 84.75476 22468.98 84.43412
## 8 52.25767 32662.28 52.17350
## 9 58.49315 27611.78 58.36770
## 10 38.51145 30206.34 38.46167
## 11 59.11861 42085.94 59.03445
## 12 41.97277 35319.04 41.92238
## 13 41.17064 40932.94 41.12891
## 14 59.13256 34723.09 59.02968
## 15 81.52177 26719.94 81.27257
## 16 70.97654 28096.52 70.79738
## 17 81.70222 28864.16 81.47103
## 18 95.58263 29421.23 95.26919
## 19 83.84718 33103.39 83.63471
## 20 59.00525 33267.97 58.89946
## 21 56.43621 35500.59 56.34622
## 22 78.85695 29292.62 78.64357
## 23 58.29237 35960.79 58.19703
## 24 73.00530 31341.74 72.83427
## 25 49.76427 28919.00 49.67841
## 26 63.57202 23439.02 63.39829
## 27 85.84627 17902.18 85.43603
## 28 36.85918 65886.60 36.83836
## 29 33.35265 75951.65 33.33784
## 30 67.76691 44333.18 67.66280
## 31 63.58582 69102.09 63.52706
## 32 30.25628 73578.60 30.24371
## 33 29.29889 70347.13 29.28656
## 34 31.38108 80690.71 31.36883
## 35 103.28147 17836.44 102.68510
## 36 37.68895 71999.58 37.66906
## 37 76.78445 31706.08 76.59794
## 38 88.29142 16929.43 87.83036 5. Resampling 60% res_1sc_0.6## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 131209.17 130546.80 130804.05 7790370 6703861 3547409
## 2 10 288712.00 294558.06 292489.46 24564277 14702548 7772658
## 3 11 166135.33 176101.76 167171.91 10335000 89796552 9244925
## 4 12 190165.00 171300.81 189336.53 15627389 340578932 14936737
## 5 13 203035.33 203043.92 203042.71 16407864 15380413 7930051
## 6 14 255308.67 248302.15 251352.13 33996993 30068945 12718652
## 7 15 276685.50 268579.56 272944.48 32397899 45887081 15633211
## 8 16 187863.17 183890.71 184861.73 13215098 4769015 3024282
## 9 17 207584.17 217640.85 209124.08 15413602 86073129 13010925
## 10 18 190483.83 197706.38 191933.19 9997570 43763623 7727681
## 11 19 166433.50 142425.65 165859.33 13785927 562904660 13452378
## 12 2 168001.67 167979.56 168008.40 8648007 7727332 4062807
## 13 20 153940.00 145142.91 153036.08 7933214 69716246 7108093
## 14 21 185474.00 171509.13 184227.53 17090833 180276893 15483147
## 15 22 221108.00 223526.84 222594.50 23178324 16685085 9514102
## 16 23 207752.00 212452.35 211012.52 19480639 9858849 5204307
## 17 24 206421.67 208912.90 207938.89 23409424 16519734 9542173
## 18 25 201641.83 205545.34 204159.13 27431043 18598126 10097936
## 19 26 162738.67 179983.89 163798.61 18112178 284114446 16972954
## 20 27 168300.83 177878.36 169772.00 13878449 78608667 11628901
## 21 28 167253.33 167484.43 167371.61 12571606 12248676 6203155
## 22 29 206618.17 207089.76 206855.23 22795775 22568102 11340647
## 23 3 166656.67 164502.48 165269.63 12732799 7819024 4756415
## 24 4 191680.00 190989.74 191323.99 18081133 17466611 8880555
## 25 5 202051.17 205899.12 204110.26 12210829 11174491 5044392
## 26 6 260126.83 254413.62 256522.36 22536140 15520755 7128552
## 27 7 327930.00 334830.51 331804.49 37742686 33196379 14885296
## 28 71 98569.00 90526.73 97894.66 5251516 63378387 4785661
## 29 72 64312.00 78221.17 64520.76 2883401 191435461 2838504
## 30 73 155092.83 134259.12 154432.65 13679542 423755151 13238236
## 31 74 91025.00 86134.69 89397.88 7724005 16560693 4908097
## 32 75 73756.67 80923.46 74185.83 3036826 48801293 2847857
## 33 76 65566.67 83662.37 65691.27 2250220 325755151 2234194
## 34 77 71549.67 73516.06 72236.58 2881784 5255019 1557109
## 35 78 323869.33 337361.81 327142.90 43855631 139823844 32942315
## 36 79 78337.50 82536.22 79224.68 3655034 14270838 2844592
## 37 8 190961.50 188536.72 189425.18 21349069 14574356 8450279
## 38 9 369298.17 355732.86 366123.82 42703155 142712899 32441775
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 59.38240 51.32543 27.11705
## 2 85.06906 50.01441 26.69132
## 3 62.21859 509.54227 55.30682
## 4 82.15989 1988.55796 78.87493
## 5 80.81001 75.77245 39.06098
## 6 133.20864 120.91727 50.32864
## 7 117.08111 170.24831 56.94585
## 8 70.31617 25.95448 16.31795
## 9 74.25086 395.45264 62.21708
## 10 52.45862 221.10782 40.26539
## 11 82.83908 3952.37375 81.11766
## 12 51.47694 46.02128 24.18606
## 13 51.53723 480.53026 46.45115
## 14 92.20024 1051.10623 84.08143
## 15 104.82710 74.66995 42.75041
## 16 93.76400 46.30186 24.74057
## 17 113.42121 79.08526 45.88975
## 18 136.04056 90.51993 49.57343
## 19 111.41376 1579.00515 103.74991
## 20 82.44847 441.77108 68.49124
## 21 75.17471 73.13942 37.06580
## 22 110.29899 108.94896 54.80969
## 23 76.40021 47.49312 28.76597
## 24 94.33011 91.43304 46.41165
## 25 60.42776 54.27746 24.76888
## 26 86.66289 60.81575 27.61559
## 27 115.11700 99.09822 45.02055
## 28 53.16688 699.00573 48.74097
## 29 44.82219 2445.76695 43.97944
## 30 88.17425 3157.99367 85.69411
## 31 84.85586 192.83937 54.87768
## 32 41.16572 602.76755 38.38251
## 33 34.31957 3891.11647 34.00859
## 34 40.28553 71.34514 21.65621
## 35 135.43797 415.25190 100.77764
## 36 46.64461 172.79699 35.91495
## 37 111.79449 77.23727 44.59134
## 38 115.61941 401.42999 88.59829 res_1sc2_0.6## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 131209.17 130546.80 131208.19 7790370 5213493816 7778649
## 2 10 288712.00 294558.06 288739.12 24564277 5214684070 24448487
## 3 11 166135.33 176101.76 166155.21 10335000 5216909162 10314074
## 4 12 190165.00 171300.81 190108.72 15627389 5225739327 15580462
## 5 13 203035.33 203043.92 203035.22 16407864 5213491779 16355488
## 6 14 255308.67 248302.15 255263.46 33996993 5215344906 33773952
## 7 15 276685.50 268579.56 276635.14 32397899 5215996973 32197408
## 8 16 187863.17 183890.71 187852.91 13215098 5214025571 13181214
## 9 17 207584.17 217640.85 207613.75 15413602 5216955893 15367922
## 10 18 190483.83 197706.38 190497.47 9997570 5215310184 9978113
## 11 19 166433.50 142425.65 166370.35 13785927 5233342232 13749208
## 12 2 168001.67 167979.56 168001.61 8648007 5213488098 8633482
## 13 20 153940.00 145142.91 153926.63 7933214 5216133918 7921061
## 14 21 185474.00 171509.13 185428.69 17090833 5220262134 17034408
## 15 22 221108.00 223526.84 221118.65 23178324 5213680255 23075039
## 16 23 207752.00 212452.35 207769.42 19480639 5214275961 19407649
## 17 24 206421.67 208912.90 206433.02 23409424 5213693926 23303689
## 18 25 201641.83 205545.34 201661.79 27431043 5214130827 27286142
## 19 26 162738.67 179983.89 162798.73 18112178 5223877958 18048842
## 20 27 168300.83 177878.36 168326.09 13878449 5216645561 13841088
## 21 28 167253.33 167484.43 167253.93 12571606 5213467486 12540773
## 22 29 206618.17 207089.76 206620.22 22795775 5213464201 22694511
## 23 3 166656.67 164502.48 166651.43 12732799 5213625791 12701520
## 24 4 191680.00 190989.74 191677.72 18081133 5213477541 18017161
## 25 5 202051.17 205899.12 202060.05 12210829 5214052574 12181730
## 26 6 260126.83 254413.62 260102.37 22536140 5214656653 22438596
## 27 7 327930.00 334830.51 327979.97 37742686 5215372539 37469373
## 28 71 98569.00 90526.73 98560.67 5251516 5215822899 5246126
## 29 72 64312.00 78221.17 64319.66 2883401 5220157001 2881776
## 30 73 155092.83 134259.12 155038.32 13679542 5228540305 13643502
## 31 74 91025.00 86134.69 91017.87 7724005 5214293753 7712363
## 32 75 73756.67 80923.46 73760.82 3036826 5215243871 3035042
## 33 76 65566.67 83662.37 65574.46 2250220 5224766880 2249223
## 34 77 71549.67 73516.06 71550.72 2881784 5213669760 2880172
## 35 78 323869.33 337361.81 323981.93 43855631 5219790125 43485524
## 36 79 78337.50 82536.22 78340.44 3655034 5214074484 3652456
## 37 8 190961.50 188536.72 190951.50 21349069 5213689961 21261213
## 38 9 369298.17 355732.86 369188.17 42703155 5219850007 42353627
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 59.38240 39937.30 59.29347
## 2 85.06906 17704.01 84.66022
## 3 62.21859 29625.11 62.08490
## 4 82.15989 30507.84 81.93767
## 5 80.81001 25677.13 80.55218
## 6 133.20864 21004.06 132.35803
## 7 117.08111 19420.72 116.37793
## 8 70.31617 28353.96 70.13989
## 9 74.25086 23970.51 74.02021
## 10 52.45862 26379.89 52.35282
## 11 82.83908 36744.00 82.65022
## 12 51.47694 31036.42 51.39049
## 13 51.53723 35939.56 51.46279
## 14 92.20024 30436.96 91.91828
## 15 104.82710 23325.36 104.35497
## 16 93.76400 24544.04 93.40483
## 17 113.42121 24956.36 112.90243
## 18 136.04056 25368.17 135.30875
## 19 111.41376 29024.34 110.98189
## 20 82.44847 29327.71 82.21411
## 21 75.17471 31129.27 74.99003
## 22 110.29899 25176.08 109.80804
## 23 76.40021 31693.85 76.21494
## 24 94.33011 27297.02 93.99746
## 25 60.42776 25323.22 60.28111
## 26 86.66289 20497.33 86.29579
## 27 115.11700 15575.88 114.26552
## 28 53.16688 57618.60 53.11710
## 29 44.82219 66741.65 44.79150
## 30 88.17425 38943.92 87.97324
## 31 84.85586 60542.16 84.73463
## 32 41.16572 64451.87 41.13922
## 33 34.31957 62458.37 34.30017
## 34 40.28553 70922.64 40.26238
## 35 135.43797 15472.52 134.24738
## 36 46.64461 63179.13 46.60997
## 37 111.79449 27654.02 111.34039
## 38 115.61941 14673.56 114.70762 6. Resampling 50% res_1sc_0.5## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 113907.8 111726.27 112283.23 8717251 3918291 2524733
## 2 10 241477.0 246170.72 244890.54 24388779 10232366 6622293
## 3 11 142582.2 150303.29 143954.67 10476229 50831606 8551929
## 4 12 158563.4 144982.34 157402.94 15710842 169215037 14347238
## 5 13 172060.4 175018.24 173978.07 16573695 9497335 5802966
## 6 14 217533.2 210989.52 213513.58 33418170 22952251 12629905
## 7 15 231665.8 226783.12 228750.73 33342537 25667225 12872834
## 8 16 157057.0 155080.11 155598.14 13801797 8127544 4760299
## 9 17 176064.0 183021.57 178605.38 17152244 32226928 10340083
## 10 18 159042.2 166837.14 160452.18 10488447 52019646 8333879
## 11 19 142507.2 121383.93 141789.50 15116645 432285029 14584940
## 12 2 143641.8 143667.54 143666.87 9125392 8552724 4405351
## 13 20 129744.0 122626.68 128540.55 8412391 43611978 6883555
## 14 21 155138.0 145874.91 153031.62 18845098 68871697 14086788
## 15 22 190358.6 190365.83 190354.74 24896365 24192939 12259804
## 16 23 175814.4 180587.81 178552.16 20980634 16461286 8255394
## 17 24 174640.6 177365.25 176442.54 25182792 16550363 9550161
## 18 25 171080.8 173737.87 173165.08 29766390 18469183 10237406
## 19 26 138397.8 152856.99 139801.87 19366298 195389747 17228630
## 20 27 143716.8 151601.49 145505.32 13971628 48591350 10694578
## 21 28 143641.6 141947.78 142540.55 12007602 8419319 4757781
## 22 29 175103.4 175423.52 175263.99 22923651 22805470 11432159
## 23 3 143413.8 141878.38 142368.89 13405935 9772510 5407828
## 24 4 159920.0 161319.48 160717.38 20014139 17311229 9218141
## 25 5 171483.4 174541.82 173532.40 14959313 9271200 5023870
## 26 6 221289.8 214785.79 218396.75 22268352 29657683 11090950
## 27 7 282816.4 284105.99 283583.21 42041524 37208536 19559474
## 28 71 82320.2 76645.92 81424.92 4502617 30591222 3556847
## 29 72 56502.4 66333.56 56826.17 3150978 93927889 3041622
## 30 73 134256.8 116809.43 133320.03 15943107 294659689 15031560
## 31 74 78680.0 72626.40 77272.56 8481580 28351376 6475760
## 32 75 62560.0 68380.91 63192.49 3344154 33327567 2890520
## 33 76 57840.0 72298.17 58021.77 2627943 207088344 2594028
## 34 77 61608.4 63294.80 62338.13 2852508 3361140 1376904
## 35 78 275134.2 287216.49 279440.77 49526635 101596532 29532918
## 36 79 65947.0 69766.46 67032.71 4044496 10946822 2814765
## 37 8 159281.6 158868.55 159060.79 21726878 20574436 10561390
## 38 9 309248.6 298816.52 304798.81 46369056 62847861 26411352
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 76.54833 35.04214 22.47097
## 2 101.00767 41.57272 27.07811
## 3 73.49470 337.36305 59.42808
## 4 99.09476 1167.64790 91.16507
## 5 96.28595 54.22969 33.34086
## 6 153.52464 108.82914 59.03068
## 7 144.06951 113.33139 56.25816
## 8 87.89356 52.45150 30.61265
## 9 97.41515 175.77079 57.96048
## 10 65.94059 311.27442 51.97281
## 11 106.09375 3562.33477 102.89354
## 12 63.52780 59.52688 30.65936
## 13 64.83134 356.58571 53.54768
## 14 121.51768 472.25510 91.95640
## 15 130.77582 127.05134 64.39499
## 16 119.28821 90.97370 46.29926
## 17 144.20558 93.43765 54.18268
## 18 174.24196 106.15015 59.07931
## 19 139.76481 1279.34497 123.22992
## 20 97.20702 320.35104 73.49379
## 21 83.57939 59.24046 33.34691
## 22 130.91045 129.99564 65.22536
## 23 93.47642 68.70245 37.92491
## 24 125.13317 107.39655 57.37322
## 25 87.16677 53.10235 28.98781
## 26 100.55761 137.75698 50.49326
## 27 148.65778 130.97325 68.96378
## 28 54.72042 399.20616 43.48728
## 29 55.75816 1415.38949 53.50968
## 30 118.80400 2526.40894 112.77655
## 31 107.79842 390.83519 83.80854
## 32 53.42789 486.92483 45.81432
## 33 45.43470 2861.02393 44.70413
## 34 46.33829 53.04248 22.15785
## 35 179.95812 355.97083 106.58544
## 36 61.33572 156.70676 42.04308
## 37 136.42256 129.47607 66.39614
## 38 149.94792 210.32436 86.64447 res_1sc2_0.5## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 113907.8 111726.27 113902.58 8717251 3695179685 8696473
## 2 10 241477.0 246170.72 241507.70 24388779 3695808509 24227091
## 3 11 142582.2 150303.29 142604.39 10476229 3697112215 10445799
## 4 12 158563.4 144982.34 158506.04 15710842 3701374906 15643844
## 5 13 172060.4 175018.24 172073.54 16573695 3695332694 16497206
## 6 14 217533.2 210989.52 217473.16 33418170 3696637934 33111410
## 7 15 231665.8 226783.12 231622.73 33342537 3696132360 33038548
## 8 16 157057.0 155080.11 157049.48 13801797 3695194238 13749770
## 9 17 176064.0 183021.57 176095.99 17152244 3696700882 17072179
## 10 18 159042.2 166837.14 159064.11 10488447 3697153603 10458438
## 11 19 142507.2 121383.93 142421.41 15116645 3710425795 15053358
## 12 2 143641.8 143667.54 143641.92 9125392 3695017899 9102569
## 13 20 129744.0 122626.68 129727.70 8412391 3696792095 8392455
## 14 21 155138.0 145874.91 155091.38 18845098 3698022868 18747793
## 15 22 190358.6 190365.83 190358.85 24896365 3695022408 24728624
## 16 23 175814.4 180587.81 175840.88 20980634 3695895768 20859797
## 17 24 174640.6 177365.25 174659.23 25182792 3695295821 25009771
## 18 25 171080.8 173737.87 171103.15 29766390 3695422865 29525092
## 19 26 138397.8 152856.99 138471.93 19366298 3702403532 19264291
## 20 27 143716.8 151601.49 143746.37 13971628 3697155495 13918383
## 21 28 143641.6 141947.78 143636.05 12007602 3695121885 11967727
## 22 29 175103.4 175423.52 175105.39 22923651 3695002227 22780468
## 23 3 143413.8 141878.38 143408.32 13405935 3695123442 13356511
## 24 4 159920.0 161319.48 159927.14 20014139 3695091355 19903332
## 25 5 171483.4 174541.82 171495.33 14959313 3695396042 14898027
## 26 6 221289.8 214785.79 221249.79 22268352 3696632203 22130601
## 27 7 282816.4 284105.99 282830.53 42041524 3695164806 41565148
## 28 71 82320.2 76645.92 82313.33 4502617 3696208284 4497007
## 29 72 56502.4 66333.56 56510.72 3150978 3698345699 3148228
## 30 73 134256.8 116809.43 134181.95 15943107 3705708593 15873908
## 31 74 78680.0 72626.40 78666.14 8481580 3696268253 8461986
## 32 75 62560.0 68380.91 62565.19 3344154 3696262694 3341060
## 33 76 57840.0 72298.17 57850.30 2627943 3702229748 2626019
## 34 77 61608.4 63294.80 61609.67 2852508 3695157584 2850138
## 35 78 275134.2 287216.49 275292.76 49526635 3700209295 48866970
## 36 79 65947.0 69766.46 65951.21 4044496 3695515093 4039963
## 37 8 159281.6 158868.55 159279.49 21726878 3695037891 21597426
## 38 9 309248.6 298816.52 309119.13 46369056 3698764252 45789624
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 76.54833 33076.13 76.36942
## 2 101.00767 15013.59 100.32508
## 3 73.49470 24600.18 73.26930
## 4 99.09476 25531.66 98.70814
## 5 96.28595 21114.39 95.83422
## 6 153.52464 17520.38 152.15959
## 7 144.06951 16297.85 142.78174
## 8 87.89356 23828.11 87.56651
## 9 97.41515 20198.13 96.94275
## 10 65.94059 22161.40 65.74285
## 11 106.09375 30566.98 105.71509
## 12 63.52780 25719.52 63.36883
## 13 64.83134 30149.56 64.68626
## 14 121.51768 25350.21 120.92670
## 15 130.77582 19411.32 129.89450
## 16 119.28821 20467.93 118.58344
## 17 144.20558 20835.21 143.19907
## 18 174.24196 21270.84 172.80333
## 19 139.76481 24222.08 138.95533
## 20 97.20702 24388.08 96.81650
## 21 83.57939 26032.52 83.30522
## 22 130.91045 21064.04 130.09131
## 23 93.47642 26044.13 93.13537
## 24 125.13317 22905.10 124.43503
## 25 87.16677 21171.56 86.80404
## 26 100.55761 17211.32 99.95511
## 27 148.65778 13005.82 146.96587
## 28 54.72042 48226.73 54.65686
## 29 55.75816 55759.61 55.70115
## 30 118.80400 31724.62 118.35502
## 31 107.79842 50899.74 107.56849
## 32 53.42789 54059.50 53.37401
## 33 45.43470 51219.28 45.39307
## 34 46.33829 58382.05 46.29871
## 35 179.95812 12883.38 177.45931
## 36 61.33572 52975.88 61.26294
## 37 136.42256 23259.26 135.61133
## 38 149.94792 12378.18 148.13694 ress_1sc_1## kab Y dir Y syn Y com MSE dir MSE syn MSE com Relatif MSE dir
## 1 1 226801 222986.8 226801 0 1.454819e+07 0 0
## 2 10 481671 491371.2 481671 0 9.409479e+07 0 0
## 3 11 284635 300358.1 284635 0 2.472163e+08 0 0
## 4 12 316410 289505.1 316410 0 7.238717e+08 0 0
## 5 13 343206 349175.9 343206 0 3.563970e+07 0 0
## 6 14 434778 420501.9 434778 0 2.038059e+08 0 0
## 7 15 461280 450573.8 461280 0 1.146229e+08 0 0
## 8 16 314028 309741.7 314028 0 1.837214e+07 0 0
## 9 17 351541 365175.3 351541 0 1.858928e+08 0 0
## 10 18 317206 333449.3 317206 0 2.638447e+08 0 0
## 11 19 285390 241795.8 285390 0 1.900454e+09 0 0
## 12 2 286903 287012.9 286903 0 1.208006e+04 0 0
## 13 20 258840 245472.0 258840 0 1.787042e+08 0 0
## 14 21 308461 290790.1 308461 0 3.122619e+08 0 0
## 15 22 379756 380949.8 379756 0 1.425268e+06 0 0
## 16 23 350703 360940.2 350703 0 1.048012e+08 0 0
## 17 24 349030 354571.2 349030 0 3.070515e+07 0 0
## 18 25 340725 346014.6 340725 0 2.797959e+07 0 0
## 19 26 278631 305628.1 278631 0 7.288427e+08 0 0
## 20 27 287661 302653.5 287661 0 2.247761e+08 0 0
## 21 28 286903 283530.1 286903 0 1.137645e+07 0 0
## 22 29 349733 350486.5 349733 0 5.678016e+05 0 0
## 23 3 285390 281784.1 285390 0 1.300264e+07 0 0
## 24 4 319600 321717.8 319600 0 4.485242e+06 0 0
## 25 5 341551 347863.8 341551 0 3.985086e+07 0 0
## 26 6 441330 428562.5 441330 0 1.630078e+08 0 0
## 27 7 561270 564932.5 561270 0 1.341375e+07 0 0
## 28 71 163306 152747.8 163306 0 1.114753e+08 0 0
## 29 72 113050 132366.7 113050 0 3.731344e+08 0 0
## 30 73 266815 231109.6 266815 0 1.274876e+09 0 0
## 31 74 157080 145013.1 157080 0 1.456100e+08 0 0
## 32 75 123753 136477.7 123753 0 1.619178e+08 0 0
## 33 76 115440 144147.8 115440 0 8.241364e+08 0 0
## 34 77 121771 125802.1 121771 0 1.625002e+07 0 0
## 35 78 542361 567311.8 542361 0 6.225448e+08 0 0
## 36 79 132355 139290.8 132355 0 4.810515e+07 0 0
## 37 8 318003 317330.1 318003 0 4.527983e+05 0 0
## 38 9 617716 595940.7 617716 0 4.741638e+08 0 0
## Relatif MSE syn Relatif MSE com
## 1 6.524240e+01 0
## 2 1.914943e+02 0
## 3 8.230718e+02 0
## 4 2.500376e+03 0
## 5 1.020681e+02 0
## 6 4.846729e+02 0
## 7 2.543931e+02 0
## 8 5.931437e+01 0
## 9 5.090509e+02 0
## 10 7.912589e+02 0
## 11 7.859749e+03 0
## 12 4.208892e-02 0
## 13 7.280026e+02 0
## 14 1.073840e+03 0
## 15 3.741353e+00 0
## 16 2.903562e+02 0
## 17 8.659796e+01 0
## 18 8.086247e+01 0
## 19 2.384737e+03 0
## 20 7.426846e+02 0
## 21 4.012433e+01 0
## 22 1.620038e+00 0
## 23 4.614397e+01 0
## 24 1.394154e+01 0
## 25 1.145588e+02 0
## 26 3.803595e+02 0
## 27 2.374400e+01 0
## 28 7.298000e+02 0
## 29 2.818945e+03 0
## 30 5.516329e+03 0
## 31 1.004116e+03 0
## 32 1.186405e+03 0
## 33 5.717302e+03 0
## 34 1.291712e+02 0
## 35 1.097359e+03 0
## 36 3.453577e+02 0
## 37 1.426900e+00 0
## 38 7.956560e+02 0 ress_1sc2_1## kab Y dir Y syn Y com MSE dir MSE syn MSE com Relatif MSE dir
## 1 1 226801 222986.8 226801 0 14599353633 0 0
## 2 10 481671 491371.2 481671 0 14602096619 0 0
## 3 11 284635 300358.1 284635 0 14607376670 0 0
## 4 12 316410 289505.1 316410 0 14623813064 0 0
## 5 13 343206 349175.9 343206 0 14600080926 0 0
## 6 14 434778 420501.9 434778 0 14605879761 0 0
## 7 15 461280 450573.8 461280 0 14602804483 0 0
## 8 16 314028 309741.7 314028 0 14599485493 0 0
## 9 17 351541 365175.3 351541 0 14605262067 0 0
## 10 18 317206 333449.3 317206 0 14607950065 0 0
## 11 19 285390 241795.8 285390 0 14664384875 0 0
## 12 2 286903 287012.9 286903 0 14598852388 0 0
## 13 20 258840 245472.0 258840 0 14605014186 0 0
## 14 21 308461 290790.1 308461 0 14609619622 0 0
## 15 22 379756 380949.8 379756 0 14598901118 0 0
## 16 23 350703 360940.2 350703 0 14602465807 0 0
## 17 24 349030 354571.2 349030 0 14599910769 0 0
## 18 25 340725 346014.6 340725 0 14599816785 0 0
## 19 26 278631 305628.1 278631 0 14623984478 0 0
## 20 27 287661 302653.5 287661 0 14606602872 0 0
## 21 28 286903 283530.1 286903 0 14599244263 0 0
## 22 29 349733 350486.5 349733 0 14598871550 0 0
## 23 3 285390 281784.1 285390 0 14599300338 0 0
## 24 4 319600 321717.8 319600 0 14599006635 0 0
## 25 5 341551 347863.8 341551 0 14600226139 0 0
## 26 6 441330 428562.5 441330 0 14604472931 0 0
## 27 7 561270 564932.5 561270 0 14599314514 0 0
## 28 71 163306 152747.8 163306 0 14602695949 0 0
## 29 72 113050 132366.7 113050 0 14611718676 0 0
## 30 73 266815 231109.6 266815 0 14642813228 0 0
## 31 74 157080 145013.1 157080 0 14603873006 0 0
## 32 75 123753 136477.7 123753 0 14604435342 0 0
## 33 76 115440 144147.8 115440 0 14627270468 0 0
## 34 77 121771 125802.1 121771 0 14599412316 0 0
## 35 78 542361 567311.8 542361 0 14620319033 0 0
## 36 79 132355 139290.8 132355 0 14600510769 0 0
## 37 8 318003 317330.1 318003 0 14598867585 0 0
## 38 9 617716 595940.7 617716 0 14615202446 0 0
## Relatif MSE syn Relatif MSE com
## 1 65471.83 0
## 2 29717.04 0
## 3 48633.20 0
## 4 50513.14 0
## 5 41812.97 0
## 6 34734.39 0
## 7 32409.35 0
## 8 47134.38 0
## 9 39995.21 0
## 10 43808.61 0
## 11 60647.81 0
## 12 50864.79 0
## 13 59497.69 0
## 14 50241.12 0
## 15 38322.37 0
## 16 40456.74 0
## 17 41176.24 0
## 18 42194.23 0
## 19 47848.95 0
## 20 48261.79 0
## 21 51490.99 0
## 22 41653.16 0
## 23 51810.24 0
## 24 45378.29 0
## 25 41971.10 0
## 26 34077.81 0
## 27 25842.58 0
## 28 95600.03 0
## 29 110388.19 0
## 30 63358.74 0
## 31 100707.27 0
## 32 107009.69 0
## 33 101474.13 0
## 34 116050.60 0
## 35 25771.22 0
## 36 104820.36 0
## 37 46005.30 0
## 38 24524.59 0 2. Resampling 90% ress_1sc_0.9## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 202742.9 220561.3 224642.1 2896549 19613373 2491406.8
## 2 10 433848.4 488197.7 483604.6 9293068 29329270 6824762.4
## 3 11 253275.2 297629.0 281671.8 4104925 260734536 4039455.5
## 4 12 285047.7 288150.4 316524.5 5567810 811237346 5529596.8
## 5 13 304140.3 343648.1 338995.8 5982514 27068003 4830374.5
## 6 14 387627.0 417610.0 429728.1 12550986 160444041 11603140.6
## 7 15 414338.9 448267.5 459405.0 11538425 137991418 10576023.8
## 8 16 282903.3 306742.3 313710.7 4724879 53429906 4329624.9
## 9 17 311251.1 361950.3 346204.7 5959340 253967686 5821593.3
## 10 18 285445.8 331062.0 317437.8 3819672 189986299 3743174.7
## 11 19 253650.2 238755.6 281715.0 5107239 1850761547 5093104.1
## 12 2 255229.4 283591.2 283590.7 3546116 3318609 1712347.7
## 13 20 227413.3 243629.4 252354.5 2944371 79477447 2836362.6
## 14 21 277258.3 287957.2 307784.7 5633829 399112575 5554903.2
## 15 22 339909.6 376462.1 376939.9 9179344 7103462 3944719.0
## 16 23 310669.7 359453.4 345695.6 7177815 197910795 6918311.7
## 17 24 309092.8 350057.1 344776.0 8776320 35703123 6965830.6
## 18 25 301844.0 342851.4 336727.3 9648618 48263217 7756121.2
## 19 26 247235.1 301990.6 274957.9 6866347 738517671 6802421.9
## 20 27 255983.0 300657.7 284743.7 5160873 258653275 5059349.7
## 21 28 255078.7 279783.8 282157.8 4571648 9230490 2850124.1
## 22 29 309802.2 344973.0 344613.8 8287884 7703859 3991583.7
## 23 3 253875.7 278607.9 280599.8 4878761 7911113 2550076.5
## 24 4 287680.0 318735.3 319153.1 7053811 6156325 3282511.6
## 25 5 302826.8 344882.1 337074.0 4937735 67282378 4569855.9
## 26 6 393964.7 425380.9 437038.1 8564663 145915072 8069096.7
## 27 7 502340.0 560265.7 559528.4 14386731 9646669 5666384.7
## 28 71 147635.4 150879.7 163896.9 1860017 172639326 1839503.1
## 29 72 100725.3 130583.3 111975.4 1089293 347458172 1085865.3
## 30 73 237766.2 226265.8 264044.1 5328707 1434377711 5308903.2
## 31 74 139166.7 143385.4 154372.1 2889981 123681115 2823189.8
## 32 75 107557.4 134831.0 119576.6 1042327 234545700 1037629.7
## 33 76 102961.1 142830.1 114434.2 935502 807560767 934412.3
## 34 77 106436.2 122329.3 118550.9 1071555 17085685 981126.0
## 35 78 479480.1 559121.7 533377.1 16265513 682617823 15874912.4
## 36 79 115064.4 137629.9 127996.6 1438460 94373742 1416761.3
## 37 8 286162.9 314840.3 315701.3 7311009 2930032 1894844.1
## 38 9 547253.8 589543.7 607148.0 16844535 326588599 16010296.0
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 14.288052 88.944640 11.090807
## 2 21.419807 60.096177 14.117308
## 3 16.207530 875.614857 14.341223
## 4 19.533320 2815.630508 17.470118
## 5 19.668846 78.743833 14.249644
## 6 32.381089 384.128931 27.001381
## 7 27.847582 307.696455 23.017856
## 8 16.701246 174.153811 13.800632
## 9 19.146407 701.613086 16.815385
## 10 13.381463 573.819934 11.791901
## 11 20.134755 7751.800358 18.078769
## 12 13.893845 11.702425 6.038178
## 13 12.946960 326.312324 11.239381
## 14 20.318131 1386.058263 18.046635
## 15 27.004996 18.867757 10.465650
## 16 23.104742 550.348267 20.013209
## 17 28.393459 101.932343 20.204544
## 18 31.962946 140.709546 23.042633
## 19 27.770355 2445.405665 24.737955
## 20 20.161583 860.240130 17.768568
## 21 17.922280 33.024559 10.098965
## 22 26.751680 22.332010 11.582699
## 23 19.217092 28.437195 9.079668
## 24 24.519859 19.314207 10.285059
## 25 16.306158 195.059606 13.559124
## 26 21.739759 342.873454 18.462264
## 27 28.639688 17.220035 10.126880
## 28 12.600390 1144.174807 11.224907
## 29 10.814333 2660.631441 9.697190
## 30 22.410720 6341.601803 20.105396
## 31 20.766329 862.737168 18.288360
## 32 9.691321 1739.651339 8.677953
## 33 9.085974 5653.255465 8.165491
## 34 10.068042 139.401231 8.280090
## 35 33.918428 1222.063189 29.760055
## 36 12.501863 685.600330 11.069210
## 37 25.548795 9.309839 5.998490
## 38 30.779853 554.013700 26.369457 ress_1sc2_0.9## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 202742.9 220561.3 225268.9 2896549 14374669032 2895963
## 2 10 433848.4 488197.7 482057.8 9293068 14375224633 9287060
## 3 11 253275.2 297629.0 281421.5 4104925 14383025223 4103752
## 4 12 285047.7 288150.4 316708.6 5567810 14402058523 5565653
## 5 13 304140.3 343648.1 337936.1 5982514 14375032501 5980018
## 6 14 387627.0 417610.0 430685.3 12550986 14379858174 12540031
## 7 15 414338.9 448267.5 460366.8 11538425 14379049030 11529164
## 8 16 282903.3 306742.3 314334.5 4724879 14375898166 4723323
## 9 17 311251.1 361950.3 345841.2 5959340 14382855829 5956861
## 10 18 285445.8 331062.0 317165.7 3819672 14380575793 3818656
## 11 19 253650.2 238755.6 281818.4 5107239 14437888303 5105425
## 12 2 255229.4 283591.2 283588.3 3546116 14373900673 3545239
## 13 20 227413.3 243629.4 252679.6 2944371 14376734960 2943765
## 14 21 277258.3 287957.2 308056.9 5633829 14387849601 5631608
## 15 22 339909.6 376462.1 377676.5 9179344 14373964410 9173475
## 16 23 310669.7 359453.4 345195.6 7177815 14380964849 7174227
## 17 24 309092.8 350057.1 343440.5 8776320 14375426602 8770951
## 18 25 301844.0 342851.4 335387.2 9648618 14375889788 9642135
## 19 26 247235.1 301990.6 274718.7 6866347 14399595725 6863070
## 20 27 255983.0 300657.7 284431.4 5160873 14382989868 5159016
## 21 28 255078.7 279783.8 283419.6 4571648 14374364758 4570191
## 22 29 309802.2 344973.0 344225.1 8287884 14373912967 8283097
## 23 3 253875.7 278607.9 282082.9 4878761 14374333858 4877101
## 24 4 287680.0 318735.3 319644.0 7053811 14373923776 7050348
## 25 5 302826.8 344882.1 336477.1 4937735 14376383177 4936037
## 26 6 393964.7 425380.9 437731.2 8564663 14379219716 8559562
## 27 7 502340.0 560265.7 558157.7 14386731 14374135090 14372336
## 28 71 147635.4 150879.7 164037.7 1860017 14379910047 1859776
## 29 72 100725.3 130583.3 111918.4 1089293 14385911706 1089210
## 30 73 237766.2 226265.8 264170.7 5328707 14423537877 5326736
## 31 74 139166.7 143385.4 154627.4 2889981 14378257349 2889398
## 32 75 107557.4 134831.0 119509.4 1042327 14382016553 1042251
## 33 76 102961.1 142830.1 114403.1 935502 14401772010 935441
## 34 77 106436.2 122329.3 118262.8 1071555 14374518940 1071475
## 35 78 479480.1 559121.7 532785.4 16265513 14397992253 16247111
## 36 79 115064.4 137629.9 127850.4 1438460 14377196697 1438316
## 37 8 286162.9 314840.3 317957.2 7311009 14374233396 7307290
## 38 9 547253.8 589543.7 608038.1 16844535 14385735350 16824824
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 14.288052 65173.93 12.856700
## 2 21.419807 29445.51 19.265205
## 3 16.207530 48326.33 14.582366
## 4 19.533320 49981.83 17.573790
## 5 19.668846 41830.62 17.694453
## 6 32.381089 34433.83 29.118311
## 7 27.847582 32077.10 25.043239
## 8 16.701246 46866.66 15.026293
## 9 19.146407 39737.03 17.224262
## 10 13.381463 43437.82 12.039973
## 11 20.134755 60471.46 18.115826
## 12 13.893845 50685.45 12.501369
## 13 12.946960 59011.29 11.649951
## 14 20.318131 49965.16 18.279581
## 15 27.004996 38181.79 24.289007
## 16 23.104742 40008.85 20.783443
## 17 28.393459 41066.60 25.538177
## 18 31.962946 41931.02 28.746896
## 19 27.770355 47682.24 24.980212
## 20 20.161583 47838.62 18.138521
## 21 17.922280 51376.55 16.124977
## 22 26.751680 41667.03 24.062576
## 23 19.217092 51593.19 17.289572
## 24 24.519859 45096.68 22.057069
## 25 16.306158 41684.74 14.670366
## 26 21.739759 33803.02 19.554460
## 27 28.639688 25655.83 25.749831
## 28 12.600390 95308.53 11.339000
## 29 10.814333 110168.58 9.732037
## 30 22.410720 63747.05 20.163259
## 31 20.766329 100278.86 18.686198
## 32 9.691321 106669.48 8.721474
## 33 9.085974 100834.80 8.176711
## 34 10.068042 117508.93 9.060537
## 35 33.918428 25751.41 30.490351
## 36 12.501863 104464.80 11.250464
## 37 25.548795 45655.89 22.982335
## 38 30.779853 24401.34 27.670436 3. Resampling 80% ress_1sc_0.8## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 179158.62 219112.8 222792.5 5322159 19126122 3908188
## 2 10 384749.25 488907.8 483031.5 16477656 47957360 12019941
## 3 11 221861.50 295773.0 277704.9 6957809 334720925 6813846
## 4 12 252889.62 286046.9 315761.1 10550147 894283788 10426580
## 5 13 273489.25 341725.7 341700.7 11208944 9282785 4978351
## 6 14 341604.75 417010.9 424697.8 21881568 82477510 16204837
## 7 15 369128.88 450568.4 459343.6 20243316 109898983 15532664
## 8 16 251619.62 307652.1 313189.5 8921334 39845985 7089153
## 9 17 279725.62 362626.5 350562.3 11713645 156905648 10888206
## 10 18 253765.00 329998.3 317731.7 6711901 157341648 6435150
## 11 19 222606.50 238479.1 278042.9 8555933 1574146529 8508993
## 12 2 223652.12 281720.3 281240.0 6161111 2240020 1514246
## 13 20 202920.00 242270.1 253185.3 5285092 124398961 5067072
## 14 21 246297.75 286419.1 307374.0 10665525 450194688 10415469
## 15 22 300069.88 376082.7 375669.4 16023132 13741030 7337910
## 16 23 279225.12 357966.0 350495.4 12857730 68525102 10677210
## 17 24 277639.00 353472.0 349525.5 15693841 26264123 9484108
## 18 25 272666.75 346517.6 343540.4 18190288 21130602 8894269
## 19 26 216606.38 300519.0 271165.3 12089974 875713788 11922663
## 20 27 224175.75 297477.1 280763.9 9384909 289024570 9085383
## 21 28 223502.62 278630.8 278852.8 8103627 6580727 3500767
## 22 29 277892.25 348046.4 347712.6 14900557 14400959 7322570
## 23 3 222457.62 276276.7 276753.9 8190352 4509702 2686799
## 24 4 255920.00 318620.2 319194.1 12672186 10781861 5809315
## 25 5 272332.75 343066.3 342325.8 8573372 4065564 2390683
## 26 6 346241.12 423623.8 431158.5 14810349 70665364 12039328
## 27 7 444156.62 560285.2 558739.2 25925417 13179687 7723208
## 28 71 130077.38 150545.0 162339.7 3077105 144328942 3010390
## 29 72 88595.50 129853.3 110846.5 1952580 363696012 1942117
## 30 73 208225.50 225015.0 260031.1 8845466 1236248932 8782532
## 31 74 121275.00 142784.6 151022.1 5018921 72763438 4688590
## 32 75 96580.75 134398.7 120860.9 1837925 185559098 1819479
## 33 76 90487.50 140598.5 113170.3 1670692 754546186 1666927
## 34 77 94967.62 122123.3 119357.5 1956878 10969900 1480838
## 35 78 428702.88 560019.1 537068.2 28614720 556497416 27192143
## 36 79 102870.50 136907.7 128909.2 2660601 66802206 2556823
## 37 8 254482.38 314181.4 315252.2 13135237 5132261 3375556
## 38 9 487526.50 589861.0 607893.5 29539233 353324350 27236040
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 29.70201 87.391788 17.529680
## 2 42.82857 98.111837 24.892037
## 3 31.36497 1131.473695 24.539321
## 4 41.71476 3126.641454 33.017698
## 5 40.98667 27.150298 14.566610
## 6 64.06263 197.810672 38.119186
## 7 54.84136 243.547431 33.754418
## 8 35.45519 129.565458 22.627740
## 9 41.87533 432.646262 31.059025
## 10 26.44928 476.747469 20.253466
## 11 38.43832 6600.849009 30.605915
## 12 27.54715 7.953235 5.386075
## 13 26.04441 513.540173 20.013020
## 14 43.30539 1571.730283 33.886989
## 15 53.39785 36.536687 19.530829
## 16 46.05130 191.365701 30.471328
## 17 56.52478 74.240547 27.138930
## 18 66.71290 61.015767 25.936062
## 19 55.81344 2914.061759 43.966735
## 20 41.86774 971.465366 32.361783
## 21 36.25790 23.596550 12.545027
## 22 53.61898 41.376251 21.058953
## 23 36.81867 16.308103 9.699578
## 24 49.51687 33.837824 18.200286
## 25 31.48164 11.853463 6.991899
## 26 42.78020 166.724987 27.918753
## 27 58.37256 23.532894 13.842666
## 28 23.63938 958.431479 18.530191
## 29 22.04345 2800.406309 17.524064
## 30 42.47486 5495.156153 33.770690
## 31 41.38463 509.769672 31.046530
## 32 19.02638 1380.698457 15.051515
## 33 18.46323 5365.245509 14.729265
## 34 20.61227 89.807044 12.438718
## 35 66.75133 994.199795 50.635982
## 36 25.86461 487.852209 19.835150
## 37 51.61875 16.340664 10.699007
## 38 60.58973 599.041167 44.803678 ress_1sc2_0.8## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 179158.62 219112.8 223946.5 5322159 14496169196 5320192
## 2 10 384749.25 488907.8 480945.6 16477656 14497548049 16458874
## 3 11 221861.50 295773.0 277335.7 6957809 14507108177 6954451
## 4 12 252889.62 286046.9 316090.2 10550147 14526527322 10542462
## 5 13 273489.25 341725.7 341861.5 11208944 14495392571 11200263
## 6 14 341604.75 417010.9 426990.9 21881568 14498924741 21848482
## 7 15 369128.88 450568.4 461396.0 20243316 14499813817 20214992
## 8 16 251619.62 307652.1 314520.3 8921334 14497007783 8915820
## 9 17 279725.62 362626.5 349667.5 11713645 14501140610 11704147
## 10 18 253765.00 329998.3 317212.2 6711901 14500983171 6708789
## 11 19 222606.50 238479.1 278234.8 8555933 14549902099 8550845
## 12 2 223652.12 281720.3 279566.1 6161111 14495498362 6158460
## 13 20 202920.00 242270.1 253645.8 5285092 14499798016 5283147
## 14 21 246297.75 286419.1 307856.4 10665525 14511217883 10657658
## 15 22 300069.88 376082.7 375088.4 16023132 14495404845 16005424
## 16 23 279225.12 357966.0 349039.4 12857730 14498132456 12846303
## 17 24 277639.00 353472.0 347055.7 15693841 14496772978 15676831
## 18 25 272666.75 346517.6 340840.6 18190288 14496544977 18167453
## 19 26 216606.38 300519.0 270782.7 12089974 14525940075 12079870
## 20 27 224175.75 297477.1 280230.9 9384909 14505616134 9378775
## 21 28 223502.62 278630.8 279377.9 8103627 14495384807 8099068
## 22 29 277892.25 348046.4 347366.0 14900557 14495343379 14885202
## 23 3 222457.62 276276.7 278071.0 8190352 14495456683 8185698
## 24 4 255920.00 318620.2 319898.9 12672186 14495391335 12661093
## 25 5 272332.75 343066.3 340417.5 8573372 14495599103 8568276
## 26 6 346241.12 423623.8 432792.1 14810349 14498273590 14795150
## 27 7 444156.62 560285.2 555204.9 25925417 14496308813 25879006
## 28 71 130077.38 150545.0 162594.1 3077105 14500409119 3076444
## 29 72 88595.50 129853.3 110747.0 1952580 14507934724 1952316
## 30 73 208225.50 225015.0 260260.4 8845466 14538260441 8840074
## 31 74 121275.00 142784.6 151590.7 5018921 14498008302 5017172
## 32 75 96580.75 134398.7 120727.7 1837925 14501788118 1837690
## 33 76 90487.50 140598.5 113112.5 1670692 14521402596 1670497
## 34 77 94967.62 122123.3 118710.0 1956878 14495771903 1956612
## 35 78 428702.88 560019.1 535926.1 28614720 14515502432 28558321
## 36 79 102870.50 136907.7 128589.7 2660601 14497721421 2660112
## 37 8 254482.38 314181.4 318099.4 13135237 14495872461 13123287
## 38 9 487526.50 589861.0 609368.4 29539233 14508528344 29479118
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 29.70201 66160.05 23.75302
## 2 42.82857 29652.87 34.22315
## 3 31.36497 49048.74 25.07906
## 4 41.71476 50783.98 33.34981
## 5 40.98667 42418.16 32.76395
## 6 64.06263 34768.50 51.17441
## 7 54.84136 32181.07 43.81314
## 8 35.45519 47121.81 28.34701
## 9 41.87533 39989.03 33.47210
## 10 26.44928 43942.71 21.14922
## 11 38.43832 61011.42 30.73498
## 12 27.54715 51453.98 22.02817
## 13 26.04441 59850.30 20.82821
## 14 43.30539 50664.28 34.62054
## 15 53.39785 38543.25 42.67095
## 16 46.05130 40502.19 36.80745
## 17 56.52478 41012.95 45.16990
## 18 66.71290 41835.18 53.30220
## 19 55.81344 48336.36 44.60934
## 20 41.86774 48762.28 33.47095
## 21 36.25790 52023.92 28.99004
## 22 53.61898 41647.97 42.85089
## 23 36.81867 52467.11 29.43830
## 24 49.51687 45494.17 39.57896
## 25 31.48164 42252.84 25.17022
## 26 42.78020 34224.31 34.18977
## 27 58.37256 25872.93 46.61368
## 28 23.63938 96321.42 18.90775
## 29 22.04345 111728.82 17.63197
## 30 42.47486 64610.83 33.96199
## 31 41.38463 101540.20 33.09684
## 32 19.02638 107904.53 15.21894
## 33 18.46323 103288.14 14.76845
## 34 20.61227 118702.27 16.48750
## 35 66.75133 25919.56 53.29107
## 36 25.86461 105897.03 20.68764
## 37 51.61875 46138.79 41.25789
## 38 60.58973 24596.42 48.37630 4. Resampling 70% ress_1sc_0.7## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 155045.29 219714.7 220135.7 6954616 3617897 2145215
## 2 10 337514.43 489112.9 485431.0 22431652 28443020 11011174
## 3 11 198529.57 296680.8 284360.3 9672909 162760953 9106627
## 4 12 221606.71 288149.6 316109.1 13432871 796250814 13208039
## 5 13 233785.29 343090.5 335487.6 13542047 70705234 11202236
## 6 14 303089.57 415645.4 431204.1 30265416 274777660 26985477
## 7 15 320262.29 444779.2 455302.7 27144723 141995200 22033199
## 8 16 220494.43 305907.6 313722.4 11524474 71394649 9903286
## 9 17 239706.14 362127.4 343150.7 14021460 374136136 13503723
## 10 18 222483.29 330850.3 318500.2 8569218 162673772 8110664
## 11 19 198378.00 238518.2 283135.8 11728571 2003025158 11659283
## 12 2 199207.14 282841.7 283415.2 8360667 5037639 3028066
## 13 20 178285.71 244117.4 253993.9 7339708 105200099 6840109
## 14 21 216030.29 287989.9 307991.2 12775069 415072524 12374721
## 15 22 261187.86 373871.7 373642.6 21294772 17826892 9482485
## 16 23 238714.86 355736.0 342177.4 16942707 200372256 15597769
## 17 24 237724.29 346131.6 342339.3 19421438 28296359 10399299
## 18 25 232975.00 339588.0 335821.7 22266252 33663509 10773931
## 19 26 194414.00 302212.7 278413.1 16303262 592142843 15830513
## 20 27 200264.71 300427.1 286917.5 11816053 193886558 11125960
## 21 28 199584.57 281403.3 282592.0 11264669 6074806 3268360
## 22 29 238098.00 341032.3 340596.4 18773759 17963671 9179525
## 23 3 198077.29 277814.9 280647.9 11548273 16040567 5939164
## 24 4 224000.00 318734.1 319220.8 16352727 13460822 7261484
## 25 5 232480.29 345587.5 332987.9 11568061 172660282 10787421
## 26 6 307291.29 426363.8 437427.5 19537992 141461897 17067426
## 27 7 385881.14 558087.4 554295.1 33126843 41281316 16046351
## 28 71 114532.71 151677.7 163239.9 4228771 151851107 4078427
## 29 72 76519.00 131687.1 109427.1 2552894 498429150 2539671
## 30 73 185937.86 226192.1 265305.2 12607804 1546647862 12504431
## 31 74 108921.43 144607.7 154969.4 6925858 114381500 6521474
## 32 75 85369.00 135835.4 122146.8 2584289 193984450 2546204
## 33 76 78021.43 142334.8 111533.2 2285941 951645087 2280397
## 34 77 83797.86 123815.7 120293.3 2628376 20766070 1970417
## 35 78 374831.00 560957.5 536997.1 38714974 612301776 36373369
## 36 79 91138.14 138772.1 130604.4 3434729 70911453 3263487
## 37 8 222722.00 315078.1 316000.8 17100765 8590512 5497269
## 38 9 427802.43 590674.9 609292.3 37775363 383036878 34294832
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 44.85426 16.44563 9.732906
## 2 66.46439 58.15313 22.721235
## 3 48.73628 548.43874 32.033967
## 4 60.60581 2763.62430 41.776909
## 5 57.93889 205.94090 33.398363
## 6 99.88496 660.92301 62.590759
## 7 84.75476 319.39770 48.368805
## 8 52.25767 233.36905 31.561444
## 9 58.49315 1033.09152 39.350276
## 10 38.51145 491.59785 25.462306
## 11 59.11861 8398.44572 41.177076
## 12 41.97277 17.81070 10.684432
## 13 41.17064 431.15179 26.932225
## 14 59.13256 1441.18921 40.177808
## 15 81.52177 47.71332 25.388595
## 16 70.97654 563.27457 45.586262
## 17 81.70222 81.59690 30.397872
## 18 95.58263 99.40195 32.246814
## 19 83.84718 1959.55345 56.857102
## 20 59.00525 645.37966 38.778566
## 21 56.43621 21.61207 11.545720
## 22 78.85695 52.68001 26.954121
## 23 58.29237 57.86715 21.143030
## 24 73.00530 42.21273 22.744493
## 25 49.76427 499.61763 32.402572
## 26 63.57202 331.62323 39.006883
## 27 85.84627 74.02856 29.009147
## 28 36.85918 999.32887 24.929644
## 29 33.35265 3784.94244 23.201502
## 30 67.76691 6837.19440 47.104898
## 31 63.58582 791.45511 42.083160
## 32 30.25628 1427.93302 20.836085
## 33 29.29889 6684.25355 20.445758
## 34 31.38108 167.16362 16.468884
## 35 103.28147 1091.95848 67.735104
## 36 37.68895 510.91127 24.990011
## 37 76.78445 27.25148 17.389309
## 38 88.29142 648.51890 56.280523 ress_1sc2_0.7## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 155045.29 219714.7 221492.4 6954616 14270982358 6951207
## 2 10 337514.43 489112.9 482174.4 22431652 14272564112 22396115
## 3 11 198529.57 296680.8 283622.5 9672909 14276772732 9666302
## 4 12 221606.71 288149.6 316554.3 13432871 14298746864 13420082
## 5 13 233785.29 343090.5 333987.6 13542047 14273731816 13529068
## 6 14 303089.57 415645.4 432948.7 30265416 14281345464 30200927
## 7 15 320262.29 444779.2 457493.3 27144723 14276659148 27092861
## 8 16 220494.43 305907.6 314984.7 11524474 14273686017 11515097
## 9 17 239706.14 362127.4 342456.6 14021460 14284211482 14007507
## 10 18 222483.29 330850.3 317841.1 8569218 14276731668 8564000
## 11 19 198378.00 238518.2 283360.5 11728571 14340301004 11718842
## 12 2 199207.14 282841.7 284580.6 8360667 14270941324 8355717
## 13 20 178285.71 244117.4 254688.5 7339708 14274707420 7335903
## 14 21 216030.29 287989.9 308596.3 12775069 14285580102 12763390
## 15 22 261187.86 373871.7 373126.6 21294772 14270946319 21262881
## 16 23 238714.86 355736.0 341038.7 16942707 14278320356 16922590
## 17 24 237724.29 346131.6 339614.9 19421438 14272387087 19394976
## 18 25 232975.00 339588.0 332832.0 22266252 14272616475 22231126
## 19 26 194414.00 302212.7 277762.1 16303262 14291807637 16284623
## 20 27 200264.71 300427.1 286104.3 11816053 14277919931 11806159
## 21 28 199584.57 281403.3 285117.9 11264669 14271339035 11255743
## 22 29 238098.00 341032.3 340141.2 18773759 14270854671 18748901
## 23 3 198077.29 277814.9 282963.4 11548273 14271778076 11538836
## 24 4 224000.00 318734.1 319998.6 16352727 14270926458 16333861
## 25 5 232480.29 345587.5 332125.6 11568061 14277179439 11558655
## 26 6 307291.29 426363.8 438970.2 19537992 14276378457 19511009
## 27 7 385881.14 558087.4 551274.5 33126843 14272793021 33050027
## 28 71 114532.71 151677.7 163614.5 4228771 14276208802 4227501
## 29 72 76519.00 131687.1 109316.9 2552894 14288101980 2552432
## 30 73 185937.86 226192.1 265590.8 12607804 14324594174 12596643
## 31 74 108921.43 144607.7 155596.7 6925858 14275009749 6922482
## 32 75 85369.00 135835.4 121958.2 2584289 14277604969 2583816
## 33 76 78021.43 142334.8 111464.1 2285941 14303720910 2285572
## 34 77 83797.86 123815.7 119712.0 2628376 14271608856 2627890
## 35 78 374831.00 560957.5 535541.7 38714974 14293275590 38610075
## 36 79 91138.14 138772.1 130199.4 3434729 14273390398 3433896
## 37 8 222722.00 315078.1 318170.6 17100765 14271186055 17080190
## 38 9 427802.43 590674.9 611092.4 37775363 14285337504 37675087
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 44.85426 64953.62 31.38272
## 2 66.46439 29180.68 46.45029
## 3 48.73628 48122.34 34.09101
## 4 60.60581 49622.95 42.38733
## 5 57.93889 41603.80 40.51729
## 6 99.88496 34359.57 69.77633
## 7 84.75476 32098.54 59.21814
## 8 52.25767 46660.40 36.55148
## 9 58.49315 39445.40 40.90213
## 10 38.51145 43151.91 26.94093
## 11 59.11861 60122.78 41.35414
## 12 41.97277 50455.77 29.36364
## 13 41.17064 58475.63 28.80513
## 14 59.13256 49604.42 41.35747
## 15 81.52177 38171.35 56.97963
## 16 70.97654 40137.88 49.62205
## 17 81.70222 41234.51 57.11214
## 18 95.58263 42030.33 66.80015
## 19 83.84718 47290.56 58.62006
## 20 59.00525 47525.67 41.26736
## 21 56.43621 50715.13 39.47445
## 22 78.85695 41846.61 55.12657
## 23 58.29237 51372.55 40.77193
## 24 73.00530 44773.91 51.04498
## 25 49.76427 41312.86 34.80551
## 26 63.57202 33484.32 44.44074
## 27 85.84627 25574.54 59.95133
## 28 36.85918 94123.71 25.79428
## 29 33.35265 108502.36 23.34177
## 30 67.76691 63333.12 47.40110
## 31 63.58582 98717.28 44.48991
## 32 30.25628 105112.29 21.17508
## 33 29.29889 100495.90 20.50499
## 34 31.38108 115272.44 21.96255
## 35 103.28147 25480.63 72.09187
## 36 37.68895 102856.55 26.37544
## 37 76.78445 45294.41 53.68506
## 38 88.29142 24184.91 61.64542 5. Resampling 60% ress_1sc_0.6## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 131209.17 217578.0 217974.8 7790370 5997205 3257929
## 2 10 288712.00 490930.1 483856.1 24564277 74335376 16858614
## 3 11 166135.33 293502.9 277514.2 10335000 267808200 9942573
## 4 12 190165.00 285501.4 316444.6 15627389 973834612 15378754
## 5 13 203035.33 338406.5 338423.1 16407864 13553836 7333530
## 6 14 255308.67 413836.9 422230.5 33996993 118271574 21578502
## 7 15 276685.50 447632.6 458619.5 32397899 173937059 24543040
## 8 16 187863.17 306484.5 311054.9 13215098 32638774 8933298
## 9 17 207584.17 362734.7 346897.6 15413602 266493985 14548639
## 10 18 190483.83 329510.6 318342.7 9997570 139339077 9180410
## 11 19 166433.50 237376.1 277044.7 13785927 1588132370 13665849
## 12 2 168001.67 279965.9 280110.2 8648007 6160098 3396893
## 13 20 153940.00 241904.9 256024.3 7933214 207759729 7636170
## 14 21 185474.00 285848.5 308375.5 17090833 531152849 16512066
## 15 22 221108.00 372544.7 371398.0 23178324 11293705 6927133
## 16 23 207752.00 354087.2 348914.4 19480639 46543566 11943560
## 17 24 206421.67 348188.2 347214.0 23409424 9394424 5942973
## 18 25 201641.83 342575.6 339159.1 27431043 38861295 12147664
## 19 26 162738.67 299973.2 271867.1 18112178 821406221 17702057
## 20 27 168300.83 296463.9 281384.1 13878449 243030208 13068612
## 21 28 167253.33 279140.7 278961.5 12571606 11674579 6045985
## 22 29 206618.17 345149.6 344762.5 22795775 22163350 11237264
## 23 3 166656.67 274170.8 275186.3 12732799 5439858 3582999
## 24 4 191680.00 318316.2 318836.0 18081133 16374129 8559357
## 25 5 202051.17 343165.2 338773.1 12210829 37014170 6862501
## 26 6 260126.83 424022.7 431004.3 22536140 74154953 15530117
## 27 7 327930.00 558050.9 549838.3 37742686 121500356 24063881
## 28 71 98569.00 150877.9 163877.1 5251516 185387105 5083808
## 29 72 64312.00 130368.6 107311.9 2883401 536891215 2867238
## 30 73 155092.83 223765.2 258091.9 13679542 1201416827 13520672
## 31 74 91025.00 143557.8 150732.1 7724005 59733488 6710278
## 32 75 73756.67 134872.4 123185.3 3036826 140957948 2968797
## 33 76 65566.67 139437.3 109352.5 2250220 908875810 2244450
## 34 77 71549.67 122526.8 119841.6 2881784 15447471 1784842
## 35 78 323869.33 562269.7 541746.4 43855631 466365131 39926837
## 36 79 78337.50 137560.4 131094.8 3655034 46139055 3363275
## 37 8 190961.50 314227.9 315434.7 21349069 11001943 6594215
## 38 9 369298.17 592888.1 613592.3 42703155 472341441 39009058
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 59.38240 27.53285 14.92975
## 2 85.06906 151.59150 34.89052
## 3 62.21859 911.81591 35.83210
## 4 82.15989 3411.55953 48.58839
## 5 80.81001 40.07413 21.67672
## 6 133.20864 285.49602 50.88901
## 7 117.08111 387.44538 53.34425
## 8 70.31617 106.55272 28.68676
## 9 74.25086 734.63213 41.93859
## 10 52.45862 422.45328 28.82668
## 11 82.83908 6690.55747 49.33246
## 12 51.47694 22.03071 12.13840
## 13 51.53723 859.19343 29.82796
## 14 92.20024 1858.14294 53.57461
## 15 104.82710 30.28391 18.64885
## 16 93.76400 131.23424 34.28749
## 17 113.42121 26.97053 17.12804
## 18 136.04056 113.19721 35.98125
## 19 111.41376 2739.02221 65.18582
## 20 82.44847 819.52350 46.43621
## 21 75.17471 41.82927 21.67575
## 22 110.29899 64.19672 32.58558
## 23 76.40021 19.84180 13.01471
## 24 94.33011 51.42321 26.84158
## 25 60.42776 107.85662 20.35095
## 26 86.66289 174.46171 35.95918
## 27 115.11700 217.79727 43.93937
## 28 53.16688 1226.86569 30.94860
## 29 44.82219 4115.60551 26.71088
## 30 88.17425 5371.98657 52.37034
## 31 84.85586 416.98059 44.52289
## 32 41.16572 1044.64713 24.09592
## 33 34.31957 6513.88824 20.52448
## 34 40.28553 125.80026 15.00327
## 35 135.43797 830.75747 73.72635
## 36 46.64461 335.22432 25.65150
## 37 111.79449 35.00503 20.88304
## 38 115.61941 797.07649 63.56848 ress_1sc2_0.6## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 131209.17 217578.0 218681.4 7790370 14481927266 7786146
## 2 10 288712.00 490930.1 481203.0 24564277 14485233528 24522466
## 3 11 166135.33 293502.9 276904.2 10335000 14491414339 10327457
## 4 12 190165.00 285501.4 316907.8 15627389 14515942574 15610462
## 5 13 203035.33 338406.5 338392.2 16407864 14481921609 16388969
## 6 14 255308.67 413836.9 425487.2 33996993 14487069185 33916355
## 7 15 276685.50 447632.6 461112.2 32397899 14488880481 32325435
## 8 16 187863.17 306484.5 313099.1 13215098 14483404363 13202879
## 9 17 207584.17 362734.7 345991.4 15413602 14491544146 15397126
## 10 18 190483.83 329510.6 317481.2 9997570 14486972734 9990557
## 11 19 166433.50 237376.1 277351.2 13785927 14537061757 13772686
## 12 2 168001.67 279965.9 280002.7 8648007 14481911382 8642772
## 13 20 153940.00 241904.9 256558.6 7933214 14489260882 7928834
## 14 21 185474.00 285848.5 309096.1 17090833 14500728149 17070476
## 15 22 221108.00 372544.7 368519.7 23178324 14482445152 23141034
## 16 23 207752.00 354087.2 346263.8 19480639 14484099891 19454299
## 17 24 206421.67 348188.2 344042.9 23409424 14482483127 23371248
## 18 25 201641.83 342575.6 336081.7 27431043 14483696741 27378700
## 19 26 162738.67 299973.2 271267.2 18112178 14510772105 18089325
## 20 27 168300.83 296463.9 280516.6 13878449 14490682114 13864976
## 21 28 167253.33 279140.7 278755.9 12571606 14481854127 12560488
## 22 29 206618.17 345149.6 344364.8 22795775 14481845002 22759214
## 23 3 166656.67 274170.8 277758.0 12732799 14482293863 12721521
## 24 4 191680.00 318316.2 319465.3 18081133 14481882057 18058050
## 25 5 202051.17 343165.2 336757.3 12210829 14483479372 12200337
## 26 6 260126.83 424022.7 433530.0 22536140 14485157370 22500926
## 27 7 327930.00 558050.9 546580.1 37742686 14487145942 37643831
## 28 71 98569.00 150877.9 164276.7 5251516 14488396940 5249574
## 29 72 64312.00 130368.6 107191.3 2883401 14500436112 2882816
## 30 73 155092.83 223765.2 258455.3 13679542 14523723069 13666545
## 31 74 91025.00 143557.8 151704.0 7724005 14484149315 7719809
## 32 75 73756.67 134872.4 122930.3 3036826 14486788532 3036183
## 33 76 65566.67 139437.3 109282.5 2250220 14513241334 2249861
## 34 77 71549.67 122526.8 119250.1 2881784 14482416001 2881203
## 35 78 323869.33 562269.7 539850.2 43855631 14499417014 43721660
## 36 79 78337.50 137560.4 130564.3 3655034 14483540232 3654105
## 37 8 190961.50 314227.9 318263.1 21349069 14482472113 21317356
## 38 9 369298.17 592888.1 615430.6 42703155 14499583353 42576656
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 59.38240 66562.16 35.61021
## 2 85.06906 29506.68 50.95285
## 3 62.21859 49375.18 37.30224
## 4 82.15989 50846.41 49.24785
## 5 80.81001 42795.21 48.43020
## 6 133.20864 35006.76 79.74069
## 7 117.08111 32367.86 70.09619
## 8 70.31617 47256.60 42.15156
## 9 74.25086 39950.85 44.50060
## 10 52.45862 43966.48 31.45229
## 11 82.83908 61240.01 49.66259
## 12 51.47694 51727.37 30.86747
## 13 51.53723 59899.26 30.90624
## 14 92.20024 50728.27 55.25912
## 15 104.82710 38875.60 62.79399
## 16 93.76400 40906.73 56.18063
## 17 113.42121 41593.94 67.94034
## 18 136.04056 42280.28 81.46572
## 19 111.41376 48373.91 66.75474
## 20 82.44847 48879.52 49.41837
## 21 75.17471 51882.12 45.06487
## 22 110.29899 41960.14 66.07304
## 23 76.40021 52823.08 45.80004
## 24 94.33011 45495.04 56.52605
## 25 60.42776 42205.37 36.22493
## 26 86.66289 34162.21 51.91822
## 27 115.11700 25959.80 68.88540
## 28 53.16688 96031.00 31.88937
## 29 44.82219 111236.09 26.88668
## 30 88.17425 64906.54 52.86107
## 31 84.85586 100903.59 50.88730
## 32 41.16572 107419.79 24.69371
## 33 34.31957 104097.28 20.58755
## 34 40.28553 118204.40 24.16632
## 35 135.43797 25787.54 81.00414
## 36 46.64461 105298.54 27.97928
## 37 111.79449 46090.03 66.97835
## 38 115.61941 24455.94 69.17369 6. Resampling 50% ress_1sc_0.5## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 113907.8 223452.5 225939.7 8717251 13417966 4534225
## 2 10 241477.0 492341.4 485736.1 24388779 69612650 16024492
## 3 11 142582.2 300606.6 285850.6 10476229 234755111 9995154
## 4 12 158563.4 289964.7 316546.6 15710842 723992674 15369941
## 5 13 172060.4 350036.5 346749.1 16573695 24378226 8182938
## 6 14 217533.2 421979.0 432613.9 33418170 160819356 25546472
## 7 15 231665.8 453566.2 460561.5 33342537 108196747 19361050
## 8 16 157057.0 310160.2 312112.5 13801797 14709754 6044108
## 9 17 176064.0 366043.1 353398.7 17152244 180364446 15449204
## 10 18 159042.2 333674.3 318789.4 10488447 239543927 9949805
## 11 19 142507.2 242767.9 284655.5 15116645 1774490049 14983719
## 12 2 143641.8 287335.1 287418.6 9125392 7134182 3867150
## 13 20 129744.0 245253.4 258886.3 8412391 199685086 8030182
## 14 21 155138.0 291749.8 309222.8 18845098 332022081 17655520
## 15 22 190358.6 380731.7 380633.3 24896365 22082658 11469399
## 16 23 175814.4 361175.6 354040.6 20980634 84422649 13427586
## 17 24 174640.6 354730.5 352535.8 25182792 19913598 9099559
## 18 25 171080.8 347475.7 344347.6 29766390 42205377 15083132
## 19 26 138397.8 305714.0 277497.6 19366298 839657882 18831881
## 20 27 143716.8 303203.0 288327.9 13971628 236280282 13152365
## 21 28 143641.6 283895.6 285254.7 12007602 9299534 4507520
## 22 29 175103.4 350847.0 350531.2 22923651 22450928 11341677
## 23 3 143413.8 283756.8 285191.5 13405935 12772696 5977135
## 24 4 159920.0 322639.0 321783.3 20014139 12781638 7238155
## 25 5 171483.4 349083.6 345090.1 14959313 35185267 8306718
## 26 6 221289.8 429571.6 440655.0 22268352 167281631 18080130
## 27 7 282816.4 568212.0 567257.4 42041524 32567509 17658463
## 28 71 82320.2 153291.8 164192.8 4502617 135872740 4266174
## 29 72 56502.4 132667.1 113166.7 3150978 385164489 3123639
## 30 73 134256.8 233618.9 268045.2 15943107 1226468077 15715220
## 31 74 78680.0 145252.8 156656.3 8481580 138850244 7980125
## 32 75 62560.0 136761.8 125436.2 3344154 143342732 3230746
## 33 76 57840.0 144596.3 115770.9 2627943 836237208 2619464
## 34 77 61608.4 126589.6 123548.8 2852508 15991623 1924609
## 35 78 275134.2 574433.0 552421.7 49526635 554966034 44528206
## 36 79 65947.0 139532.9 132436.9 4044496 55920776 3737063
## 37 8 159281.6 317737.1 318005.1 21726878 17117109 9437413
## 38 9 309248.6 597633.0 616272.3 46369056 390498614 41379630
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 76.54833 60.16135 20.03730
## 2 101.00767 141.56504 33.06432
## 3 73.49470 779.22303 34.97584
## 4 99.09476 2497.85659 48.56177
## 5 96.28595 69.56145 23.63074
## 6 153.52464 381.46102 58.89866
## 7 144.06951 238.54967 41.85054
## 8 87.89356 47.37535 19.31765
## 9 97.41515 492.13661 43.71634
## 10 65.94059 716.86878 31.21070
## 11 106.09375 7311.48905 52.64846
## 12 63.52780 24.83832 13.45641
## 13 64.83134 816.12645 31.01618
## 14 121.51768 1138.25859 57.10807
## 15 130.77582 57.96872 30.12256
## 16 119.28821 232.77264 38.01100
## 17 144.20558 56.00244 25.86188
## 18 174.24196 121.60850 43.95976
## 19 139.76481 2748.71804 67.79761
## 20 97.20702 778.95745 45.61056
## 21 83.57939 32.81937 15.77618
## 22 130.91045 63.98642 32.35394
## 23 93.47642 45.05633 20.94153
## 24 125.13317 39.67578 22.51904
## 25 87.16677 100.94394 24.15022
## 26 100.55761 388.41636 40.88473
## 27 148.65778 57.32157 31.13328
## 28 54.72042 886.50288 25.97270
## 29 55.75816 2902.08160 27.59671
## 30 118.80400 5257.61700 58.65288
## 31 107.79842 956.85403 50.94198
## 32 53.42789 1047.22461 25.75342
## 33 45.43470 5776.54243 22.62580
## 34 46.33829 126.01854 15.66211
## 35 179.95812 970.54248 80.68590
## 36 61.33572 400.34893 28.22429
## 37 136.42256 53.85421 29.67587
## 38 149.94792 653.44410 67.14807 ress_1sc2_0.5## kab Y dir Y syn Y com MSE dir MSE syn MSE com
## 1 1 113907.8 223452.5 227813.0 8717251 14780718740 8712047
## 2 10 241477.0 492341.4 482969.4 24388779 14783234034 24348153
## 3 11 142582.2 300606.6 285175.5 10476229 14788448859 10468605
## 4 12 158563.4 289964.7 317098.0 15710842 14805499624 15694038
## 5 13 172060.4 350036.5 344127.4 16573695 14781330776 16554504
## 6 14 217533.2 421979.0 435036.2 33418170 14786551735 33340938
## 7 15 231665.8 453566.2 463309.9 33342537 14784529442 33266007
## 8 16 157057.0 310160.2 314110.2 13801797 14780776953 13788753
## 9 17 176064.0 366043.1 352144.1 17152244 14786803526 17132157
## 10 18 159042.2 333674.3 318095.4 10488447 14788614412 10480929
## 11 19 142507.2 242767.9 284971.4 15116645 14841703181 15100772
## 12 2 143641.8 287335.1 287283.7 9125392 14780071594 9119676
## 13 20 129744.0 245253.4 259479.8 8412391 14787168381 8407398
## 14 21 155138.0 291749.8 310252.6 18845098 14792091474 18820675
## 15 22 190358.6 380731.7 380717.3 24896365 14780089630 24854216
## 16 23 175814.4 361175.6 351642.1 20980634 14783583070 20950291
## 17 24 174640.6 354730.5 349290.6 25182792 14781183286 25139309
## 18 25 171080.8 347475.7 342172.8 29766390 14781691460 29705692
## 19 26 138397.8 305714.0 276832.8 19366298 14809614129 19340694
## 20 27 143716.8 303203.0 287448.4 13971628 14788621982 13958278
## 21 28 143641.6 283895.6 287280.4 12007602 14780487541 11997607
## 22 29 175103.4 350847.0 350207.8 22923651 14780008907 22887685
## 23 3 143413.8 283756.8 286824.9 13405935 14780493768 13393544
## 24 4 159920.0 322639.0 319843.6 20014139 14780365421 19986319
## 25 5 171483.4 349083.6 342972.8 14959313 14781584168 14943944
## 26 6 221289.8 429571.6 442559.5 22268352 14786528812 22233750
## 27 7 282816.4 568212.0 565639.9 42041524 14780659225 41921402
## 28 71 82320.2 153291.8 164637.0 4502617 14784833135 4501213
## 29 72 56502.4 132667.1 113009.0 3150978 14793382794 3150290
## 30 73 134256.8 233618.9 268476.1 15943107 14822834371 15925750
## 31 74 78680.0 145252.8 157353.1 8481580 14785073013 8476673
## 32 75 62560.0 136761.8 125122.6 3344154 14785050774 3343380
## 33 76 57840.0 144596.3 115685.2 2627943 14808918990 2627462
## 34 77 61608.4 126589.6 123217.4 2852508 14780630337 2851915
## 35 78 275134.2 574433.0 550348.5 49526635 14800837181 49360041
## 36 79 65947.0 139532.9 131896.1 4044496 14782060374 4043361
## 37 8 159281.6 317737.1 318562.1 21726878 14780151563 21694367
## 38 9 309248.6 597633.0 618431.8 46369056 14795057008 46222814
## Relatif MSE dir Relatif MSE syn Relatif MSE com
## 1 76.54833 66152.27 38.25176
## 2 101.00767 30027.17 50.41807
## 3 73.49470 49200.37 36.71911
## 4 99.09476 51063.33 49.49891
## 5 96.28595 42228.78 48.08630
## 6 153.52464 35040.77 76.59052
## 7 144.06951 32595.69 71.87268
## 8 87.89356 47656.22 43.90578
## 9 97.41515 40396.27 48.64831
## 10 65.94059 44322.79 32.94552
## 11 106.09375 61133.95 52.99941
## 12 63.52780 51439.04 31.74399
## 13 64.83134 60299.12 32.39750
## 14 121.51768 50700.43 60.68469
## 15 130.77582 38822.64 65.27718
## 16 119.28821 40935.86 59.55561
## 17 144.20558 41670.42 71.97630
## 18 174.24196 42541.69 86.94001
## 19 139.76481 48444.16 69.78077
## 20 97.20702 48776.16 48.55455
## 21 83.57939 52065.03 41.75534
## 22 130.91045 42128.08 65.35234
## 23 93.47642 52088.26 46.69546
## 24 125.13317 45810.20 62.47894
## 25 87.16677 42343.13 43.53790
## 26 100.55761 34422.64 50.20315
## 27 148.65778 26011.64 74.11557
## 28 54.72042 96453.46 27.35226
## 29 55.75816 111519.22 27.87195
## 30 118.80400 63449.25 59.34571
## 31 107.79842 101799.47 53.87042
## 32 53.42789 108119.01 26.70720
## 33 45.43470 102438.55 22.71214
## 34 46.33829 116764.11 23.16419
## 35 179.95812 25766.76 89.66339
## 36 61.33572 105951.77 30.65875
## 37 136.42256 46518.52 68.10941
## 38 149.94792 24756.37 74.74550Penentuan X signifikan
#Cek X signifikan
urut_kab <- as.numeric(c(1,10,11,12,13,14,15,16,17,18,19,2,20,21,22,23,24,25,26,27,28,29,3,4,5,6,7,71,72,73,74,75,76,77,78,79,8,9))
y_kab_x <- res_1sc_1$`Y dir`
gab_9 <- cbind(urut_kab,y_kab_x)
gab_10 <- gab_9[order(gab_9[,1]),]
y_x <- gab_10[,2]
# variabel terpilih : x7,x13,x14,x18,x19,x20,x21,x24,x25,x34
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)
#Cari yang paling signifikan
# 1. Cek semua variabel
hsl_br_x <- lm(y_x~x1$x+x2$x+x3$x+x4$x+x5$x+x6$x+x7$x+x8$x+x9$x+x10$x
+x11$x+x12$x+x13$x+x14$x+x15$x+x16$x+x17$x+x18$x+x19$x+x20$x
+x21$x+x22$x+x23$x+x24$x+x25$x+x26$x+x27$x+x28$x+x29$x+x30$x
+x31$x+x32$x+x33$x+x34$x+x35$x)
summary(hsl_br_x)##
## Call:
## lm(formula = y_x ~ x1$x + x2$x + x3$x + x4$x + x5$x + x6$x +
## x7$x + x8$x + x9$x + x10$x + x11$x + x12$x + x13$x + x14$x +
## x15$x + x16$x + x17$x + x18$x + x19$x + x20$x + x21$x + x22$x +
## x23$x + x24$x + x25$x + x26$x + x27$x + x28$x + x29$x + x30$x +
## x31$x + x32$x + x33$x + x34$x + x35$x)
##
## Residuals:
## Min 1Q Median 3Q Max
## -26447 -6066 -1168 5921 22270
##
## Coefficients: (5 not defined because of singularities)
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.163e+05 1.429e+04 8.133 8.2e-05 ***
## x1$x 1.616e+01 1.963e+01 0.823 0.438
## x2$x -4.456e+01 4.866e+01 -0.916 0.390
## x3$x 7.560e+01 9.714e+01 0.778 0.462
## x4$x 5.263e+03 7.708e+03 0.683 0.517
## x5$x NA NA NA NA
## x6$x NA NA NA NA
## x7$x 1.850e+04 1.543e+04 1.199 0.270
## x8$x -5.485e+02 1.621e+03 -0.338 0.745
## x9$x 1.766e-01 3.679e+00 0.048 0.963
## x10$x 3.132e+00 3.624e+00 0.864 0.416
## x11$x NA NA NA NA
## x12$x 2.327e+02 3.479e+02 0.669 0.525
## x13$x 9.796e-01 7.885e-01 1.242 0.254
## x14$x -1.049e+01 1.256e+01 -0.836 0.431
## x15$x NA NA NA NA
## x16$x -1.076e+01 1.959e+01 -0.549 0.600
## x17$x 1.928e+03 1.409e+03 1.368 0.213
## x18$x -1.929e+04 1.727e+04 -1.117 0.301
## x19$x -1.990e+04 1.752e+04 -1.135 0.294
## x20$x -1.918e+02 2.425e+02 -0.791 0.455
## x21$x 5.152e+02 6.421e+02 0.802 0.449
## x22$x 7.323e+01 3.956e+02 0.185 0.858
## x23$x 1.625e+02 3.994e+02 0.407 0.696
## x24$x 6.364e+01 3.419e+01 1.861 0.105
## x25$x -5.354e+01 3.050e+01 -1.756 0.123
## x26$x 2.991e+00 5.256e+01 0.057 0.956
## x27$x 3.679e+01 5.674e+01 0.648 0.537
## x28$x -3.040e+01 3.918e+01 -0.776 0.463
## x29$x 2.819e+01 1.989e+02 0.142 0.891
## x30$x NA NA NA NA
## x31$x 8.888e+00 7.195e+01 0.124 0.905
## x32$x 1.379e+02 3.073e+02 0.449 0.667
## x33$x -1.314e+02 1.273e+02 -1.033 0.336
## x34$x 2.438e+00 4.653e+00 0.524 0.616
## x35$x 4.386e-01 2.124e+00 0.206 0.842
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 24600 on 7 degrees of freedom
## Multiple R-squared: 0.9923, Adjusted R-squared: 0.9593
## F-statistic: 30.08 on 30 and 7 DF, p-value: 5.436e-05# 2. Seleksi variabel awal
#Hasil dengan backward
step(lm(y_x~x1$x+x2$x+x3$x+x4$x+x5$x+x6$x+x7$x+x8$x+x9$x+x10$x
+x11$x+x12$x+x13$x+x14$x+x15$x+x16$x+x17$x+x18$x+x19$x+x20$x
+x21$x+x22$x+x23$x+x24$x+x25$x+x26$x+x27$x+x28$x+x29$x+x30$x
+x31$x+x32$x+x33$x+x34$x+x35$x),direction="backward")## Start: AIC=766.1
## y_x ~ x1$x + x2$x + x3$x + x4$x + x5$x + x6$x + x7$x + x8$x +
## x9$x + x10$x + x11$x + x12$x + x13$x + x14$x + x15$x + x16$x +
## x17$x + x18$x + x19$x + x20$x + x21$x + x22$x + x23$x + x24$x +
## x25$x + x26$x + x27$x + x28$x + x29$x + x30$x + x31$x + x32$x +
## x33$x + x34$x + x35$x
##
##
## Step: AIC=766.1
## y_x ~ x1$x + x2$x + x3$x + x4$x + x5$x + x6$x + x7$x + x8$x +
## x9$x + x10$x + x11$x + x12$x + x13$x + x14$x + x15$x + x16$x +
## x17$x + x18$x + x19$x + x20$x + x21$x + x22$x + x23$x + x24$x +
## x25$x + x26$x + x27$x + x28$x + x29$x + x31$x + x32$x + x33$x +
## x34$x + x35$x
##
##
## Step: AIC=766.1
## y_x ~ x1$x + x2$x + x3$x + x4$x + x5$x + x6$x + x7$x + x8$x +
## x9$x + x10$x + x11$x + x12$x + x13$x + x14$x + x16$x + x17$x +
## x18$x + x19$x + x20$x + x21$x + x22$x + x23$x + x24$x + x25$x +
## x26$x + x27$x + x28$x + x29$x + x31$x + x32$x + x33$x + x34$x +
## x35$x
##
##
## Step: AIC=766.1
## y_x ~ x1$x + x2$x + x3$x + x4$x + x5$x + x6$x + x7$x + x8$x +
## x9$x + x10$x + x12$x + x13$x + x14$x + x16$x + x17$x + x18$x +
## x19$x + x20$x + x21$x + x22$x + x23$x + x24$x + x25$x + x26$x +
## x27$x + x28$x + x29$x + x31$x + x32$x + x33$x + x34$x + x35$x
##
##
## Step: AIC=766.1
## y_x ~ x1$x + x2$x + x3$x + x4$x + x5$x + x7$x + x8$x + x9$x +
## x10$x + x12$x + x13$x + x14$x + x16$x + x17$x + x18$x + x19$x +
## x20$x + x21$x + x22$x + x23$x + x24$x + x25$x + x26$x + x27$x +
## x28$x + x29$x + x31$x + x32$x + x33$x + x34$x + x35$x
##
##
## Step: AIC=766.1
## y_x ~ x1$x + x2$x + x3$x + x4$x + x7$x + x8$x + x9$x + x10$x +
## x12$x + x13$x + x14$x + x16$x + x17$x + x18$x + x19$x + x20$x +
## x21$x + x22$x + x23$x + x24$x + x25$x + x26$x + x27$x + x28$x +
## x29$x + x31$x + x32$x + x33$x + x34$x + x35$x
##
## Df Sum of Sq RSS AIC
## - x9$x 1 1394069 4236087902 764.11
## - x26$x 1 1959467 4236653299 764.12
## - x31$x 1 9232767 4243926599 764.18
## - x29$x 1 12158194 4246852026 764.21
## - x22$x 1 20729763 4255423595 764.29
## - x35$x 1 25785512 4260479345 764.33
## - x8$x 1 69309711 4304003543 764.72
## - x23$x 1 100209110 4334902943 764.99
## - x32$x 1 121833222 4356527054 765.18
## - x34$x 1 166142605 4400836437 765.56
## - x16$x 1 182361517 4417055350 765.70
## <none> 4234693833 766.10
## - x27$x 1 254408621 4489102454 766.32
## - x12$x 1 270782407 4505476239 766.46
## - x4$x 1 282047791 4516741624 766.55
## - x28$x 1 364364428 4599058261 767.24
## - x3$x 1 366431585 4601125418 767.26
## - x20$x 1 378542172 4613236005 767.36
## - x21$x 1 389513187 4624207019 767.45
## - x1$x 1 409777539 4644471371 767.61
## - x14$x 1 422594204 4657288036 767.72
## - x10$x 1 451971179 4686665011 767.96
## - x2$x 1 507401761 4742095593 768.40
## - x33$x 1 645040032 4879733865 769.49
## - x18$x 1 754988836 4989682669 770.34
## - x19$x 1 779726454 5014420286 770.52
## - x7$x 1 868997538 5103691371 771.19
## - x13$x 1 933697774 5168391606 771.67
## - x17$x 1 1132796745 5367490578 773.11
## - x25$x 1 1864682363 6099376196 777.97
## - x24$x 1 2096104599 6330798432 779.38
##
## Step: AIC=764.11
## y_x ~ x1$x + x2$x + x3$x + x4$x + x7$x + x8$x + x10$x + x12$x +
## x13$x + x14$x + x16$x + x17$x + x18$x + x19$x + x20$x + x21$x +
## x22$x + x23$x + x24$x + x25$x + x26$x + x27$x + x28$x + x29$x +
## x31$x + x32$x + x33$x + x34$x + x35$x
##
## Df Sum of Sq RSS AIC
## - x26$x 1 1988159 4238076061 762.13
## - x31$x 1 8138874 4244226776 762.19
## - x29$x 1 19929604 4256017506 762.29
## - x22$x 1 23199854 4259287756 762.32
## - x35$x 1 38055896 4274143797 762.45
## - x8$x 1 101148758 4337236660 763.01
## - x32$x 1 129486289 4365574190 763.26
## - x23$x 1 133714530 4369802431 763.30
## - x16$x 1 186551170 4422639072 763.75
## <none> 4236087902 764.11
## - x27$x 1 269773831 4505861733 764.46
## - x12$x 1 295341259 4531429160 764.68
## - x34$x 1 302783663 4538871565 764.74
## - x20$x 1 399669633 4635757535 765.54
## - x4$x 1 409641276 4645729178 765.62
## - x3$x 1 426488113 4662576014 765.76
## - x28$x 1 432717139 4668805041 765.81
## - x21$x 1 592052023 4828139925 767.09
## - x14$x 1 601073012 4837160913 767.16
## - x1$x 1 729889558 4965977460 768.16
## - x2$x 1 881813479 5117901381 769.30
## - x10$x 1 911363511 5147451412 769.52
## - x33$x 1 1020417606 5256505507 770.32
## - x18$x 1 1186014255 5422102156 771.49
## - x19$x 1 1229782179 5465870080 771.80
## - x17$x 1 1289890464 5525978366 772.22
## - x7$x 1 1397097774 5633185676 772.95
## - x13$x 1 1775729050 6011816952 775.42
## - x25$x 1 2588693214 6824781116 780.24
## - x24$x 1 2738599692 6974687594 781.06
##
## Step: AIC=762.13
## y_x ~ x1$x + x2$x + x3$x + x4$x + x7$x + x8$x + x10$x + x12$x +
## x13$x + x14$x + x16$x + x17$x + x18$x + x19$x + x20$x + x21$x +
## x22$x + x23$x + x24$x + x25$x + x27$x + x28$x + x29$x + x31$x +
## x32$x + x33$x + x34$x + x35$x
##
## Df Sum of Sq RSS AIC
## - x22$x 1 21649885 4259725946 760.33
## - x29$x 1 25096588 4263172649 760.36
## - x35$x 1 36091115 4274167176 760.45
## - x8$x 1 111797952 4349874013 761.12
## - x23$x 1 158690865 4396766926 761.53
## - x32$x 1 159261381 4397337442 761.53
## - x16$x 1 188646717 4426722778 761.79
## <none> 4238076061 762.13
## - x31$x 1 231269601 4469345662 762.15
## - x27$x 1 273042208 4511118269 762.50
## - x12$x 1 293448423 4531524485 762.68
## - x34$x 1 348481100 4586557162 763.13
## - x20$x 1 399855202 4637931263 763.56
## - x28$x 1 431448762 4669524823 763.82
## - x4$x 1 446780372 4684856433 763.94
## - x3$x 1 450667689 4688743750 763.97
## - x21$x 1 688917736 4926993797 765.86
## - x1$x 1 734753020 4972829081 766.21
## - x2$x 1 906154334 5144230395 767.50
## - x10$x 1 909836945 5147913006 767.52
## - x14$x 1 931826879 5169902940 767.68
## - x17$x 1 1308556098 5546632159 770.36
## - x18$x 1 1440643576 5678719637 771.25
## - x19$x 1 1495643702 5733719763 771.62
## - x7$x 1 1664730353 5902806414 772.72
## - x33$x 1 1861052047 6099128108 773.97
## - x25$x 1 3059640151 7297716212 780.78
## - x24$x 1 3255330296 7493406357 781.79
## - x13$x 1 3268400563 7506476624 781.85
##
## Step: AIC=760.33
## y_x ~ x1$x + x2$x + x3$x + x4$x + x7$x + x8$x + x10$x + x12$x +
## x13$x + x14$x + x16$x + x17$x + x18$x + x19$x + x20$x + x21$x +
## x23$x + x24$x + x25$x + x27$x + x28$x + x29$x + x31$x + x32$x +
## x33$x + x34$x + x35$x
##
## Df Sum of Sq RSS AIC
## - x35$x 1 31716636 4291442582 758.61
## - x29$x 1 37485186 4297211131 758.66
## - x8$x 1 95599915 4355325861 759.17
## - x23$x 1 146294866 4406020812 759.61
## - x32$x 1 149675937 4409401883 759.64
## <none> 4259725946 760.33
## - x27$x 1 251445730 4511171676 760.50
## - x31$x 1 297048647 4556774593 760.89
## - x16$x 1 320481955 4580207901 761.08
## - x34$x 1 328777523 4588503469 761.15
## - x12$x 1 492922815 4752648761 762.49
## - x4$x 1 506797414 4766523359 762.60
## - x3$x 1 523508002 4783233947 762.73
## - x28$x 1 665701163 4925427109 763.84
## - x20$x 1 696490868 4956216814 764.08
## - x1$x 1 824810413 5084536359 765.05
## - x2$x 1 888634073 5148360019 765.53
## - x14$x 1 940173231 5199899177 765.90
## - x10$x 1 1019684885 5279410831 766.48
## - x21$x 1 1057620668 5317346613 766.75
## - x17$x 1 1838197759 6097923705 771.96
## - x33$x 1 2185429919 6445155865 774.06
## - x18$x 1 2354986896 6614712842 775.05
## - x19$x 1 2450546894 6710272840 775.59
## - x7$x 1 2801300987 7061026932 777.53
## - x25$x 1 3039096997 7298822942 778.79
## - x24$x 1 3245242447 7504968393 779.85
## - x13$x 1 3856269520 8115995466 782.82
##
## Step: AIC=758.61
## y_x ~ x1$x + x2$x + x3$x + x4$x + x7$x + x8$x + x10$x + x12$x +
## x13$x + x14$x + x16$x + x17$x + x18$x + x19$x + x20$x + x21$x +
## x23$x + x24$x + x25$x + x27$x + x28$x + x29$x + x31$x + x32$x +
## x33$x + x34$x
##
## Df Sum of Sq RSS AIC
## - x29$x 1 37382969 4328825551 756.94
## - x8$x 1 64851172 4356293753 757.18
## - x23$x 1 114660139 4406102721 757.61
## - x32$x 1 200931180 4492373762 758.35
## <none> 4291442582 758.61
## - x31$x 1 266527179 4557969761 758.90
## - x27$x 1 272196357 4563638939 758.94
## - x16$x 1 302241796 4593684378 759.19
## - x12$x 1 477023368 4768465950 760.61
## - x34$x 1 504375945 4795818527 760.83
## - x3$x 1 521462635 4812905217 760.97
## - x4$x 1 560034584 4851477166 761.27
## - x20$x 1 674213408 4965655990 762.15
## - x28$x 1 767188753 5058631335 762.86
## - x2$x 1 900901346 5192343928 763.85
## - x1$x 1 1045916317 5337358899 764.90
## - x14$x 1 1060877791 5352320373 765.00
## - x21$x 1 1082164348 5373606930 765.15
## - x10$x 1 1086644127 5378086709 765.18
## - x17$x 1 1867647847 6159090429 770.34
## - x33$x 1 2354167260 6645609842 773.23
## - x18$x 1 2561767567 6853210148 774.40
## - x19$x 1 2642636980 6934079562 774.84
## - x7$x 1 2932508533 7223951115 776.40
## - x25$x 1 3014779957 7306222539 776.83
## - x24$x 1 3288469958 7579912540 778.22
## - x13$x 1 5339672227 9631114809 787.33
##
## Step: AIC=756.94
## y_x ~ x1$x + x2$x + x3$x + x4$x + x7$x + x8$x + x10$x + x12$x +
## x13$x + x14$x + x16$x + x17$x + x18$x + x19$x + x20$x + x21$x +
## x23$x + x24$x + x25$x + x27$x + x28$x + x31$x + x32$x + x33$x +
## x34$x
##
## Df Sum of Sq RSS AIC
## - x8$x 1 152595558 4481421109 756.25
## <none> 4328825551 756.94
## - x23$x 1 238020410 4566845961 756.97
## - x32$x 1 263628990 4592454541 757.18
## - x31$x 1 283436346 4612261897 757.35
## - x27$x 1 329480106 4658305657 757.72
## - x12$x 1 444678672 4773504223 758.65
## - x34$x 1 489092237 4817917788 759.00
## - x16$x 1 516249368 4845074919 759.22
## - x20$x 1 642415208 4971240759 760.20
## - x3$x 1 671119301 4999944852 760.41
## - x28$x 1 736340297 5065165848 760.91
## - x2$x 1 882342600 5211168151 761.99
## - x1$x 1 1009744469 5338570020 762.90
## - x10$x 1 1049715523 5378541074 763.19
## - x4$x 1 1097132447 5425957998 763.52
## - x14$x 1 1260051507 5588877058 764.65
## - x21$x 1 1262949432 5591774983 764.67
## - x17$x 1 1850048140 6178873691 768.46
## - x33$x 1 2339395613 6668221164 771.36
## - x25$x 1 3107660601 7436486152 775.50
## - x18$x 1 3137753242 7466578793 775.65
## - x19$x 1 3286215741 7615041292 776.40
## - x24$x 1 3324088663 7652914214 776.59
## - x7$x 1 3631333415 7960158966 778.08
## - x13$x 1 5489003823 9817829374 786.06
##
## Step: AIC=756.25
## y_x ~ x1$x + x2$x + x3$x + x4$x + x7$x + x10$x + x12$x + x13$x +
## x14$x + x16$x + x17$x + x18$x + x19$x + x20$x + x21$x + x23$x +
## x24$x + x25$x + x27$x + x28$x + x31$x + x32$x + x33$x + x34$x
##
## Df Sum of Sq RSS AIC
## - x23$x 1 126844344 4608265453 755.31
## - x32$x 1 207621122 4689042231 755.97
## <none> 4481421109 756.25
## - x31$x 1 268976258 4750397367 756.47
## - x12$x 1 305242094 4786663203 756.76
## - x27$x 1 368140831 4849561940 757.25
## - x16$x 1 386692188 4868113297 757.40
## - x3$x 1 554786128 5036207237 758.69
## - x28$x 1 587145905 5068567014 758.93
## - x20$x 1 735694265 5217115374 760.03
## - x2$x 1 764227482 5245648591 760.24
## - x1$x 1 863452954 5344874063 760.95
## - x34$x 1 932995336 5414416445 761.44
## - x10$x 1 960449023 5441870132 761.63
## - x4$x 1 967882247 5449303356 761.68
## - x17$x 1 1784959873 6266380981 766.99
## - x14$x 1 1814022446 6295443555 767.17
## - x21$x 1 1876482036 6357903145 767.54
## - x33$x 1 2195427173 6676848282 769.40
## - x25$x 1 3029418843 7510839952 773.88
## - x24$x 1 3277258538 7758679647 775.11
## - x18$x 1 4410441862 8891862971 780.29
## - x7$x 1 4444362250 8925783359 780.44
## - x19$x 1 4598392265 9079813374 781.09
## - x13$x 1 5385078125 9866499234 784.24
##
## Step: AIC=755.31
## y_x ~ x1$x + x2$x + x3$x + x4$x + x7$x + x10$x + x12$x + x13$x +
## x14$x + x16$x + x17$x + x18$x + x19$x + x20$x + x21$x + x24$x +
## x25$x + x27$x + x28$x + x31$x + x32$x + x33$x + x34$x
##
## Df Sum of Sq RSS AIC
## - x32$x 1 89913114 4.6982e+09 754.05
## - x31$x 1 145497228 4.7538e+09 754.50
## <none> 4.6083e+09 755.31
## - x16$x 1 294042349 4.9023e+09 755.66
## - x27$x 1 374571297 4.9828e+09 756.28
## - x3$x 1 428279819 5.0365e+09 756.69
## - x12$x 1 447593055 5.0559e+09 756.84
## - x28$x 1 540104128 5.1484e+09 757.53
## - x2$x 1 638708939 5.2470e+09 758.25
## - x1$x 1 754261206 5.3625e+09 759.07
## - x20$x 1 766047138 5.3743e+09 759.16
## - x34$x 1 818404052 5.4267e+09 759.53
## - x10$x 1 843781482 5.4520e+09 759.70
## - x4$x 1 844769028 5.4530e+09 759.71
## - x17$x 1 1736064572 6.3443e+09 765.46
## - x21$x 1 2097685378 6.7060e+09 767.57
## - x14$x 1 2114650054 6.7229e+09 767.67
## - x33$x 1 2168438060 6.7767e+09 767.97
## - x24$x 1 3333018255 7.9413e+09 773.99
## - x25$x 1 3333322769 7.9416e+09 774.00
## - x18$x 1 4663815664 9.2721e+09 779.88
## - x7$x 1 4774673705 9.3829e+09 780.33
## - x19$x 1 4836051609 9.4443e+09 780.58
## - x13$x 1 5826679650 1.0435e+10 784.37
##
## Step: AIC=754.05
## y_x ~ x1$x + x2$x + x3$x + x4$x + x7$x + x10$x + x12$x + x13$x +
## x14$x + x16$x + x17$x + x18$x + x19$x + x20$x + x21$x + x24$x +
## x25$x + x27$x + x28$x + x31$x + x33$x + x34$x
##
## Df Sum of Sq RSS AIC
## - x31$x 1 151567025 4.8497e+09 753.26
## - x16$x 1 250749907 4.9489e+09 754.02
## <none> 4.6982e+09 754.05
## - x27$x 1 287874355 4.9861e+09 754.31
## - x3$x 1 367029545 5.0652e+09 754.91
## - x12$x 1 387656715 5.0858e+09 755.06
## - x28$x 1 480515320 5.1787e+09 755.75
## - x2$x 1 559249005 5.2574e+09 756.32
## - x1$x 1 664914458 5.3631e+09 757.08
## - x20$x 1 733738117 5.4319e+09 757.56
## - x4$x 1 754892879 5.4531e+09 757.71
## - x10$x 1 810351669 5.5085e+09 758.10
## - x34$x 1 886377316 5.5846e+09 758.62
## - x17$x 1 1674352704 6.3725e+09 763.63
## - x21$x 1 2053267277 6.7514e+09 765.83
## - x14$x 1 2159146163 6.8573e+09 766.42
## - x33$x 1 2319380481 7.0176e+09 767.30
## - x25$x 1 3446691889 8.1449e+09 772.96
## - x24$x 1 3557695347 8.2559e+09 773.47
## - x18$x 1 4574814944 9.2730e+09 777.89
## - x7$x 1 4689800911 9.3880e+09 778.35
## - x19$x 1 4746388133 9.4446e+09 778.58
## - x13$x 1 7060241079 1.1758e+10 786.91
##
## Step: AIC=753.26
## y_x ~ x1$x + x2$x + x3$x + x4$x + x7$x + x10$x + x12$x + x13$x +
## x14$x + x16$x + x17$x + x18$x + x19$x + x20$x + x21$x + x24$x +
## x25$x + x27$x + x28$x + x33$x + x34$x
##
## Df Sum of Sq RSS AIC
## - x16$x 1 255096224 5.1048e+09 753.20
## <none> 4.8497e+09 753.26
## - x3$x 1 320483217 5.1702e+09 753.69
## - x27$x 1 398007742 5.2478e+09 754.25
## - x28$x 1 424977505 5.2747e+09 754.45
## - x12$x 1 488371583 5.3381e+09 754.90
## - x1$x 1 574683489 5.4244e+09 755.51
## - x2$x 1 599760629 5.4495e+09 755.69
## - x20$x 1 715960899 5.5657e+09 756.49
## - x10$x 1 717820453 5.5676e+09 756.50
## - x34$x 1 745228016 5.5950e+09 756.69
## - x4$x 1 757122240 5.6069e+09 756.77
## - x17$x 1 1945476689 6.7952e+09 764.07
## - x21$x 1 1974787719 6.8245e+09 764.24
## - x14$x 1 2032629069 6.8824e+09 764.56
## - x33$x 1 2169128101 7.0189e+09 765.30
## - x25$x 1 3394182933 8.2439e+09 771.42
## - x24$x 1 3464596817 8.3143e+09 771.74
## - x18$x 1 4831892750 9.6816e+09 777.52
## - x7$x 1 4870671098 9.7204e+09 777.68
## - x19$x 1 4980266323 9.8300e+09 778.10
## - x13$x 1 6912068507 1.1762e+10 784.92
##
## Step: AIC=753.2
## y_x ~ x1$x + x2$x + x3$x + x4$x + x7$x + x10$x + x12$x + x13$x +
## x14$x + x17$x + x18$x + x19$x + x20$x + x21$x + x24$x + x25$x +
## x27$x + x28$x + x33$x + x34$x
##
## Df Sum of Sq RSS AIC
## - x3$x 1 172146767 5.2770e+09 752.46
## <none> 5.1048e+09 753.20
## - x27$x 1 285212000 5.3901e+09 753.27
## - x12$x 1 463182889 5.5680e+09 754.50
## - x28$x 1 482528992 5.5874e+09 754.64
## - x2$x 1 550091527 5.6549e+09 755.09
## - x10$x 1 618691146 5.7235e+09 755.55
## - x34$x 1 646212113 5.7511e+09 755.73
## - x20$x 1 658387723 5.7632e+09 755.81
## - x4$x 1 917918545 6.0228e+09 757.49
## - x1$x 1 970166345 6.0750e+09 757.81
## - x17$x 1 1692303364 6.7971e+09 762.08
## - x21$x 1 2120408338 7.2253e+09 764.40
## - x33$x 1 2223160543 7.3280e+09 764.94
## - x24$x 1 3386340312 8.4912e+09 770.54
## - x25$x 1 3512674607 8.6175e+09 771.10
## - x14$x 1 4667504960 9.7723e+09 775.88
## - x18$x 1 4690297442 9.7951e+09 775.97
## - x7$x 1 4810609503 9.9155e+09 776.43
## - x19$x 1 4836492902 9.9413e+09 776.53
## - x13$x 1 6815510546 1.1920e+10 783.43
##
## Step: AIC=752.46
## y_x ~ x1$x + x2$x + x4$x + x7$x + x10$x + x12$x + x13$x + x14$x +
## x17$x + x18$x + x19$x + x20$x + x21$x + x24$x + x25$x + x27$x +
## x28$x + x33$x + x34$x
##
## Df Sum of Sq RSS AIC
## <none> 5.2770e+09 752.46
## - x2$x 1 379556768 5.6565e+09 753.10
## - x27$x 1 393047716 5.6700e+09 753.19
## - x28$x 1 484293413 5.7613e+09 753.80
## - x34$x 1 496530734 5.7735e+09 753.88
## - x10$x 1 549351183 5.8263e+09 754.23
## - x12$x 1 752326579 6.0293e+09 755.53
## - x4$x 1 767831983 6.0448e+09 755.63
## - x20$x 1 785370790 6.0624e+09 755.74
## - x1$x 1 1167207409 6.4442e+09 758.06
## - x17$x 1 1859734558 7.1367e+09 761.94
## - x33$x 1 2052242768 7.3292e+09 762.95
## - x24$x 1 3380845127 8.6578e+09 769.28
## - x25$x 1 3616547073 8.8935e+09 770.30
## - x21$x 1 3854815767 9.1318e+09 771.30
## - x14$x 1 4740562894 1.0018e+10 774.82
## - x18$x 1 4811377888 1.0088e+10 775.09
## - x19$x 1 4903647745 1.0181e+10 775.43
## - x7$x 1 4929257866 1.0206e+10 775.53
## - x13$x 1 7129247097 1.2406e+10 782.95##
## Call:
## lm(formula = y_x ~ x1$x + x2$x + x4$x + x7$x + x10$x + x12$x +
## x13$x + x14$x + x17$x + x18$x + x19$x + x20$x + x21$x + x24$x +
## x25$x + x27$x + x28$x + x33$x + x34$x)
##
## Coefficients:
## (Intercept) x1$x x2$x x4$x x7$x
## 114692.419 10.702 -16.665 1461.718 21251.091
## x10$x x12$x x13$x x14$x x17$x
## 1.158 223.326 1.033 -14.512 1657.706
## x18$x x19$x x20$x x21$x x24$x
## -22601.643 -23024.615 -193.912 718.521 48.719
## x25$x x27$x x28$x x33$x x34$x
## -40.576 34.673 -19.963 -120.294 1.668#hasil:
lm(formula = y_x ~ x1$x + x2$x + x4$x + x7$x + x10$x + x12$x +
x13$x + x14$x + x17$x + x18$x + x19$x + x20$x + x21$x + x24$x +
x25$x + x27$x + x28$x + x33$x + x34$x)##
## Call:
## lm(formula = y_x ~ x1$x + x2$x + x4$x + x7$x + x10$x + x12$x +
## x13$x + x14$x + x17$x + x18$x + x19$x + x20$x + x21$x + x24$x +
## x25$x + x27$x + x28$x + x33$x + x34$x)
##
## Coefficients:
## (Intercept) x1$x x2$x x4$x x7$x
## 114692.419 10.702 -16.665 1461.718 21251.091
## x10$x x12$x x13$x x14$x x17$x
## 1.158 223.326 1.033 -14.512 1657.706
## x18$x x19$x x20$x x21$x x24$x
## -22601.643 -23024.615 -193.912 718.521 48.719
## x25$x x27$x x28$x x33$x x34$x
## -40.576 34.673 -19.963 -120.294 1.668#Hasil dengan both
cek1<-step(lm(y_x~x1$x+x2$x+x3$x+x4$x+x5$x+x6$x+x7$x+x8$x+x9$x+x10$x
+x11$x+x12$x+x13$x+x14$x+x15$x+x16$x+x17$x+x18$x+x19$x+x20$x
+x21$x+x22$x+x23$x+x24$x+x25$x+x26$x+x27$x+x28$x+x29$x+x30$x
+x31$x+x32$x+x33$x+x34$x+x35$x),direction="both")## Start: AIC=766.1
## y_x ~ x1$x + x2$x + x3$x + x4$x + x5$x + x6$x + x7$x + x8$x +
## x9$x + x10$x + x11$x + x12$x + x13$x + x14$x + x15$x + x16$x +
## x17$x + x18$x + x19$x + x20$x + x21$x + x22$x + x23$x + x24$x +
## x25$x + x26$x + x27$x + x28$x + x29$x + x30$x + x31$x + x32$x +
## x33$x + x34$x + x35$x
##
##
## Step: AIC=766.1
## y_x ~ x1$x + x2$x + x3$x + x4$x + x5$x + x6$x + x7$x + x8$x +
## x9$x + x10$x + x11$x + x12$x + x13$x + x14$x + x15$x + x16$x +
## x17$x + x18$x + x19$x + x20$x + x21$x + x22$x + x23$x + x24$x +
## x25$x + x26$x + x27$x + x28$x + x29$x + x31$x + x32$x + x33$x +
## x34$x + x35$x
##
##
## Step: AIC=766.1
## y_x ~ x1$x + x2$x + x3$x + x4$x + x5$x + x6$x + x7$x + x8$x +
## x9$x + x10$x + x11$x + x12$x + x13$x + x14$x + x16$x + x17$x +
## x18$x + x19$x + x20$x + x21$x + x22$x + x23$x + x24$x + x25$x +
## x26$x + x27$x + x28$x + x29$x + x31$x + x32$x + x33$x + x34$x +
## x35$x
##
##
## Step: AIC=766.1
## y_x ~ x1$x + x2$x + x3$x + x4$x + x5$x + x6$x + x7$x + x8$x +
## x9$x + x10$x + x12$x + x13$x + x14$x + x16$x + x17$x + x18$x +
## x19$x + x20$x + x21$x + x22$x + x23$x + x24$x + x25$x + x26$x +
## x27$x + x28$x + x29$x + x31$x + x32$x + x33$x + x34$x + x35$x
##
##
## Step: AIC=766.1
## y_x ~ x1$x + x2$x + x3$x + x4$x + x5$x + x7$x + x8$x + x9$x +
## x10$x + x12$x + x13$x + x14$x + x16$x + x17$x + x18$x + x19$x +
## x20$x + x21$x + x22$x + x23$x + x24$x + x25$x + x26$x + x27$x +
## x28$x + x29$x + x31$x + x32$x + x33$x + x34$x + x35$x
##
##
## Step: AIC=766.1
## y_x ~ x1$x + x2$x + x3$x + x4$x + x7$x + x8$x + x9$x + x10$x +
## x12$x + x13$x + x14$x + x16$x + x17$x + x18$x + x19$x + x20$x +
## x21$x + x22$x + x23$x + x24$x + x25$x + x26$x + x27$x + x28$x +
## x29$x + x31$x + x32$x + x33$x + x34$x + x35$x
##
## Df Sum of Sq RSS AIC
## - x9$x 1 1394069 4236087902 764.11
## - x26$x 1 1959467 4236653299 764.12
## - x31$x 1 9232767 4243926599 764.18
## - x29$x 1 12158194 4246852026 764.21
## - x22$x 1 20729763 4255423595 764.29
## - x35$x 1 25785512 4260479345 764.33
## - x8$x 1 69309711 4304003543 764.72
## - x23$x 1 100209110 4334902943 764.99
## - x32$x 1 121833222 4356527054 765.18
## - x34$x 1 166142605 4400836437 765.56
## - x16$x 1 182361517 4417055350 765.70
## <none> 4234693833 766.10
## - x27$x 1 254408621 4489102454 766.32
## - x12$x 1 270782407 4505476239 766.46
## - x4$x 1 282047791 4516741624 766.55
## - x28$x 1 364364428 4599058261 767.24
## - x3$x 1 366431585 4601125418 767.26
## - x20$x 1 378542172 4613236005 767.36
## - x21$x 1 389513187 4624207019 767.45
## - x1$x 1 409777539 4644471371 767.61
## - x14$x 1 422594204 4657288036 767.72
## - x10$x 1 451971179 4686665011 767.96
## - x2$x 1 507401761 4742095593 768.40
## - x33$x 1 645040032 4879733865 769.49
## - x18$x 1 754988836 4989682669 770.34
## - x19$x 1 779726454 5014420286 770.52
## - x7$x 1 868997538 5103691371 771.19
## - x13$x 1 933697774 5168391606 771.67
## - x17$x 1 1132796745 5367490578 773.11
## - x25$x 1 1864682363 6099376196 777.97
## - x24$x 1 2096104599 6330798432 779.38
##
## Step: AIC=764.11
## y_x ~ x1$x + x2$x + x3$x + x4$x + x7$x + x8$x + x10$x + x12$x +
## x13$x + x14$x + x16$x + x17$x + x18$x + x19$x + x20$x + x21$x +
## x22$x + x23$x + x24$x + x25$x + x26$x + x27$x + x28$x + x29$x +
## x31$x + x32$x + x33$x + x34$x + x35$x
##
## Df Sum of Sq RSS AIC
## - x26$x 1 1988159 4238076061 762.13
## - x31$x 1 8138874 4244226776 762.19
## - x29$x 1 19929604 4256017506 762.29
## - x22$x 1 23199854 4259287756 762.32
## - x35$x 1 38055896 4274143797 762.45
## - x8$x 1 101148758 4337236660 763.01
## - x32$x 1 129486289 4365574190 763.26
## - x23$x 1 133714530 4369802431 763.30
## - x16$x 1 186551170 4422639072 763.75
## <none> 4236087902 764.11
## - x27$x 1 269773831 4505861733 764.46
## - x12$x 1 295341259 4531429160 764.68
## - x34$x 1 302783663 4538871565 764.74
## - x20$x 1 399669633 4635757535 765.54
## - x4$x 1 409641276 4645729178 765.62
## - x3$x 1 426488113 4662576014 765.76
## - x28$x 1 432717139 4668805041 765.81
## + x9$x 1 1394069 4234693833 766.10
## + x11$x 1 1394069 4234693833 766.10
## - x21$x 1 592052023 4828139925 767.09
## - x14$x 1 601073012 4837160913 767.16
## - x1$x 1 729889558 4965977460 768.16
## - x2$x 1 881813479 5117901381 769.30
## - x10$x 1 911363511 5147451412 769.52
## - x33$x 1 1020417606 5256505507 770.32
## - x18$x 1 1186014255 5422102156 771.49
## - x19$x 1 1229782179 5465870080 771.80
## - x17$x 1 1289890464 5525978366 772.22
## - x7$x 1 1397097774 5633185676 772.95
## - x13$x 1 1775729050 6011816952 775.42
## - x25$x 1 2588693214 6824781116 780.24
## - x24$x 1 2738599692 6974687594 781.06
##
## Step: AIC=762.13
## y_x ~ x1$x + x2$x + x3$x + x4$x + x7$x + x8$x + x10$x + x12$x +
## x13$x + x14$x + x16$x + x17$x + x18$x + x19$x + x20$x + x21$x +
## x22$x + x23$x + x24$x + x25$x + x27$x + x28$x + x29$x + x31$x +
## x32$x + x33$x + x34$x + x35$x
##
## Df Sum of Sq RSS AIC
## - x22$x 1 21649885 4259725946 760.33
## - x29$x 1 25096588 4263172649 760.36
## - x35$x 1 36091115 4274167176 760.45
## - x8$x 1 111797952 4349874013 761.12
## - x23$x 1 158690865 4396766926 761.53
## - x32$x 1 159261381 4397337442 761.53
## - x16$x 1 188646717 4426722778 761.79
## <none> 4238076061 762.13
## - x31$x 1 231269601 4469345662 762.15
## - x27$x 1 273042208 4511118269 762.50
## - x12$x 1 293448423 4531524485 762.68
## - x34$x 1 348481100 4586557162 763.13
## - x20$x 1 399855202 4637931263 763.56
## - x28$x 1 431448762 4669524823 763.82
## - x4$x 1 446780372 4684856433 763.94
## - x3$x 1 450667689 4688743750 763.97
## + x26$x 1 1988159 4236087902 764.11
## + x9$x 1 1422762 4236653299 764.12
## + x11$x 1 1422762 4236653299 764.12
## - x21$x 1 688917736 4926993797 765.86
## - x1$x 1 734753020 4972829081 766.21
## - x2$x 1 906154334 5144230395 767.50
## - x10$x 1 909836945 5147913006 767.52
## - x14$x 1 931826879 5169902940 767.68
## - x17$x 1 1308556098 5546632159 770.36
## - x18$x 1 1440643576 5678719637 771.25
## - x19$x 1 1495643702 5733719763 771.62
## - x7$x 1 1664730353 5902806414 772.72
## - x33$x 1 1861052047 6099128108 773.97
## - x25$x 1 3059640151 7297716212 780.78
## - x24$x 1 3255330296 7493406357 781.79
## - x13$x 1 3268400563 7506476624 781.85
##
## Step: AIC=760.33
## y_x ~ x1$x + x2$x + x3$x + x4$x + x7$x + x8$x + x10$x + x12$x +
## x13$x + x14$x + x16$x + x17$x + x18$x + x19$x + x20$x + x21$x +
## x23$x + x24$x + x25$x + x27$x + x28$x + x29$x + x31$x + x32$x +
## x33$x + x34$x + x35$x
##
## Df Sum of Sq RSS AIC
## - x35$x 1 31716636 4291442582 758.61
## - x29$x 1 37485186 4297211131 758.66
## - x8$x 1 95599915 4355325861 759.17
## - x23$x 1 146294866 4406020812 759.61
## - x32$x 1 149675937 4409401883 759.64
## <none> 4259725946 760.33
## - x27$x 1 251445730 4511171676 760.50
## - x31$x 1 297048647 4556774593 760.89
## - x16$x 1 320481955 4580207901 761.08
## - x34$x 1 328777523 4588503469 761.15
## + x22$x 1 21649885 4238076061 762.13
## + x9$x 1 3563260 4256162685 762.29
## + x11$x 1 3563260 4256162685 762.29
## + x26$x 1 438190 4259287756 762.32
## - x12$x 1 492922815 4752648761 762.49
## - x4$x 1 506797414 4766523359 762.60
## - x3$x 1 523508002 4783233947 762.73
## - x28$x 1 665701163 4925427109 763.84
## - x20$x 1 696490868 4956216814 764.08
## - x1$x 1 824810413 5084536359 765.05
## - x2$x 1 888634073 5148360019 765.53
## - x14$x 1 940173231 5199899177 765.90
## - x10$x 1 1019684885 5279410831 766.48
## - x21$x 1 1057620668 5317346613 766.75
## - x17$x 1 1838197759 6097923705 771.96
## - x33$x 1 2185429919 6445155865 774.06
## - x18$x 1 2354986896 6614712842 775.05
## - x19$x 1 2450546894 6710272840 775.59
## - x7$x 1 2801300987 7061026932 777.53
## - x25$x 1 3039096997 7298822942 778.79
## - x24$x 1 3245242447 7504968393 779.85
## - x13$x 1 3856269520 8115995466 782.82
##
## Step: AIC=758.61
## y_x ~ x1$x + x2$x + x3$x + x4$x + x7$x + x8$x + x10$x + x12$x +
## x13$x + x14$x + x16$x + x17$x + x18$x + x19$x + x20$x + x21$x +
## x23$x + x24$x + x25$x + x27$x + x28$x + x29$x + x31$x + x32$x +
## x33$x + x34$x
##
## Df Sum of Sq RSS AIC
## - x29$x 1 37382969 4328825551 756.94
## - x8$x 1 64851172 4356293753 757.18
## - x23$x 1 114660139 4406102721 757.61
## - x32$x 1 200931180 4492373762 758.35
## <none> 4291442582 758.61
## - x31$x 1 266527179 4557969761 758.90
## - x27$x 1 272196357 4563638939 758.94
## - x16$x 1 302241796 4593684378 759.19
## + x35$x 1 31716636 4259725946 760.33
## + x9$x 1 25269066 4266173515 760.38
## + x11$x 1 25269066 4266173515 760.38
## + x22$x 1 17275406 4274167176 760.45
## + x26$x 1 176157 4291266425 760.61
## - x12$x 1 477023368 4768465950 760.61
## - x34$x 1 504375945 4795818527 760.83
## - x3$x 1 521462635 4812905217 760.97
## - x4$x 1 560034584 4851477166 761.27
## - x20$x 1 674213408 4965655990 762.15
## - x28$x 1 767188753 5058631335 762.86
## - x2$x 1 900901346 5192343928 763.85
## - x1$x 1 1045916317 5337358899 764.90
## - x14$x 1 1060877791 5352320373 765.00
## - x21$x 1 1082164348 5373606930 765.15
## - x10$x 1 1086644127 5378086709 765.18
## - x17$x 1 1867647847 6159090429 770.34
## - x33$x 1 2354167260 6645609842 773.23
## - x18$x 1 2561767567 6853210148 774.40
## - x19$x 1 2642636980 6934079562 774.84
## - x7$x 1 2932508533 7223951115 776.40
## - x25$x 1 3014779957 7306222539 776.83
## - x24$x 1 3288469958 7579912540 778.22
## - x13$x 1 5339672227 9631114809 787.33
##
## Step: AIC=756.94
## y_x ~ x1$x + x2$x + x3$x + x4$x + x7$x + x8$x + x10$x + x12$x +
## x13$x + x14$x + x16$x + x17$x + x18$x + x19$x + x20$x + x21$x +
## x23$x + x24$x + x25$x + x27$x + x28$x + x31$x + x32$x + x33$x +
## x34$x
##
## Df Sum of Sq RSS AIC
## - x8$x 1 152595558 4481421109 756.25
## <none> 4328825551 756.94
## - x23$x 1 238020410 4566845961 756.97
## - x32$x 1 263628990 4592454541 757.18
## - x31$x 1 283436346 4612261897 757.35
## - x27$x 1 329480106 4658305657 757.72
## + x29$x 1 37382969 4291442582 758.61
## - x12$x 1 444678672 4773504223 758.65
## + x35$x 1 31614420 4297211131 758.66
## + x22$x 1 28556115 4300269436 758.69
## + x9$x 1 12894711 4315930840 758.82
## + x11$x 1 12894711 4315930840 758.82
## + x26$x 1 1001026 4327824525 758.93
## - x34$x 1 489092237 4817917788 759.00
## - x16$x 1 516249368 4845074919 759.22
## - x20$x 1 642415208 4971240759 760.20
## - x3$x 1 671119301 4999944852 760.41
## - x28$x 1 736340297 5065165848 760.91
## - x2$x 1 882342600 5211168151 761.99
## - x1$x 1 1009744469 5338570020 762.90
## - x10$x 1 1049715523 5378541074 763.19
## - x4$x 1 1097132447 5425957998 763.52
## - x14$x 1 1260051507 5588877058 764.65
## - x21$x 1 1262949432 5591774983 764.67
## - x17$x 1 1850048140 6178873691 768.46
## - x33$x 1 2339395613 6668221164 771.36
## - x25$x 1 3107660601 7436486152 775.50
## - x18$x 1 3137753242 7466578793 775.65
## - x19$x 1 3286215741 7615041292 776.40
## - x24$x 1 3324088663 7652914214 776.59
## - x7$x 1 3631333415 7960158966 778.08
## - x13$x 1 5489003823 9817829374 786.06
##
## Step: AIC=756.25
## y_x ~ x1$x + x2$x + x3$x + x4$x + x7$x + x10$x + x12$x + x13$x +
## x14$x + x16$x + x17$x + x18$x + x19$x + x20$x + x21$x + x23$x +
## x24$x + x25$x + x27$x + x28$x + x31$x + x32$x + x33$x + x34$x
##
## Df Sum of Sq RSS AIC
## - x23$x 1 126844344 4608265453 755.31
## - x32$x 1 207621122 4689042231 755.97
## <none> 4481421109 756.25
## - x31$x 1 268976258 4750397367 756.47
## - x12$x 1 305242094 4786663203 756.76
## + x8$x 1 152595558 4328825551 756.94
## + x29$x 1 125127356 4356293753 757.18
## - x27$x 1 368140831 4849561940 757.25
## - x16$x 1 386692188 4868113297 757.40
## + x35$x 1 21115925 4460305184 758.07
## + x9$x 1 16502820 4464918289 758.11
## + x11$x 1 16502820 4464918289 758.11
## + x22$x 1 12473051 4468948058 758.15
## + x26$x 1 12272255 4469148854 758.15
## - x3$x 1 554786128 5036207237 758.69
## - x28$x 1 587145905 5068567014 758.93
## - x20$x 1 735694265 5217115374 760.03
## - x2$x 1 764227482 5245648591 760.24
## - x1$x 1 863452954 5344874063 760.95
## - x34$x 1 932995336 5414416445 761.44
## - x10$x 1 960449023 5441870132 761.63
## - x4$x 1 967882247 5449303356 761.68
## - x17$x 1 1784959873 6266380981 766.99
## - x14$x 1 1814022446 6295443555 767.17
## - x21$x 1 1876482036 6357903145 767.54
## - x33$x 1 2195427173 6676848282 769.40
## - x25$x 1 3029418843 7510839952 773.88
## - x24$x 1 3277258538 7758679647 775.11
## - x18$x 1 4410441862 8891862971 780.29
## - x7$x 1 4444362250 8925783359 780.44
## - x19$x 1 4598392265 9079813374 781.09
## - x13$x 1 5385078125 9866499234 784.24
##
## Step: AIC=755.31
## y_x ~ x1$x + x2$x + x3$x + x4$x + x7$x + x10$x + x12$x + x13$x +
## x14$x + x16$x + x17$x + x18$x + x19$x + x20$x + x21$x + x24$x +
## x25$x + x27$x + x28$x + x31$x + x32$x + x33$x + x34$x
##
## Df Sum of Sq RSS AIC
## - x32$x 1 89913114 4.6982e+09 754.05
## - x31$x 1 145497228 4.7538e+09 754.50
## <none> 4.6083e+09 755.31
## - x16$x 1 294042349 4.9023e+09 755.66
## + x29$x 1 195470195 4.4128e+09 755.67
## + x23$x 1 126844344 4.4814e+09 756.25
## - x27$x 1 374571297 4.9828e+09 756.28
## - x3$x 1 428279819 5.0365e+09 756.69
## - x12$x 1 447593055 5.0559e+09 756.84
## + x26$x 1 44686128 4.5636e+09 756.94
## + x35$x 1 41721354 4.5665e+09 756.97
## + x8$x 1 41419491 4.5668e+09 756.97
## + x9$x 1 24785395 4.5835e+09 757.11
## + x11$x 1 24785395 4.5835e+09 757.11
## + x22$x 1 18127660 4.5901e+09 757.16
## - x28$x 1 540104128 5.1484e+09 757.53
## - x2$x 1 638708939 5.2470e+09 758.25
## - x1$x 1 754261206 5.3625e+09 759.07
## - x20$x 1 766047138 5.3743e+09 759.16
## - x34$x 1 818404052 5.4267e+09 759.53
## - x10$x 1 843781482 5.4520e+09 759.70
## - x4$x 1 844769028 5.4530e+09 759.71
## - x17$x 1 1736064572 6.3443e+09 765.46
## - x21$x 1 2097685378 6.7060e+09 767.57
## - x14$x 1 2114650054 6.7229e+09 767.67
## - x33$x 1 2168438060 6.7767e+09 767.97
## - x24$x 1 3333018255 7.9413e+09 773.99
## - x25$x 1 3333322769 7.9416e+09 774.00
## - x18$x 1 4663815664 9.2721e+09 779.88
## - x7$x 1 4774673705 9.3829e+09 780.33
## - x19$x 1 4836051609 9.4443e+09 780.58
## - x13$x 1 5826679650 1.0435e+10 784.37
##
## Step: AIC=754.05
## y_x ~ x1$x + x2$x + x3$x + x4$x + x7$x + x10$x + x12$x + x13$x +
## x14$x + x16$x + x17$x + x18$x + x19$x + x20$x + x21$x + x24$x +
## x25$x + x27$x + x28$x + x31$x + x33$x + x34$x
##
## Df Sum of Sq RSS AIC
## - x31$x 1 151567025 4.8497e+09 753.26
## - x16$x 1 250749907 4.9489e+09 754.02
## <none> 4.6982e+09 754.05
## - x27$x 1 287874355 4.9861e+09 754.31
## + x29$x 1 184889753 4.5133e+09 754.52
## - x3$x 1 367029545 5.0652e+09 754.91
## - x12$x 1 387656715 5.0858e+09 755.06
## + x26$x 1 96272088 4.6019e+09 755.26
## + x32$x 1 89913114 4.6083e+09 755.31
## + x8$x 1 59097457 4.6391e+09 755.57
## - x28$x 1 480515320 5.1787e+09 755.75
## + x9$x 1 12953615 4.6852e+09 755.94
## + x11$x 1 12953615 4.6852e+09 755.94
## + x22$x 1 10381713 4.6878e+09 755.96
## + x23$x 1 9136335 4.6890e+09 755.97
## + x35$x 1 1596879 4.6966e+09 756.04
## - x2$x 1 559249005 5.2574e+09 756.32
## - x1$x 1 664914458 5.3631e+09 757.08
## - x20$x 1 733738117 5.4319e+09 757.56
## - x4$x 1 754892879 5.4531e+09 757.71
## - x10$x 1 810351669 5.5085e+09 758.10
## - x34$x 1 886377316 5.5846e+09 758.62
## - x17$x 1 1674352704 6.3725e+09 763.63
## - x21$x 1 2053267277 6.7514e+09 765.83
## - x14$x 1 2159146163 6.8573e+09 766.42
## - x33$x 1 2319380481 7.0176e+09 767.30
## - x25$x 1 3446691889 8.1449e+09 772.96
## - x24$x 1 3557695347 8.2559e+09 773.47
## - x18$x 1 4574814944 9.2730e+09 777.89
## - x7$x 1 4689800911 9.3880e+09 778.35
## - x19$x 1 4746388133 9.4446e+09 778.58
## - x13$x 1 7060241079 1.1758e+10 786.91
##
## Step: AIC=753.26
## y_x ~ x1$x + x2$x + x3$x + x4$x + x7$x + x10$x + x12$x + x13$x +
## x14$x + x16$x + x17$x + x18$x + x19$x + x20$x + x21$x + x24$x +
## x25$x + x27$x + x28$x + x33$x + x34$x
##
## Df Sum of Sq RSS AIC
## - x16$x 1 255096224 5.1048e+09 753.20
## <none> 4.8497e+09 753.26
## - x3$x 1 320483217 5.1702e+09 753.69
## + x31$x 1 151567025 4.6982e+09 754.05
## + x29$x 1 132485839 4.7173e+09 754.20
## - x27$x 1 398007742 5.2478e+09 754.25
## + x8$x 1 113938581 4.7358e+09 754.35
## + x26$x 1 109626027 4.7401e+09 754.39
## - x28$x 1 424977505 5.2747e+09 754.45
## + x32$x 1 95982910 4.7538e+09 754.50
## + x22$x 1 64100120 4.7856e+09 754.75
## - x12$x 1 488371583 5.3381e+09 754.90
## + x23$x 1 14140559 4.8356e+09 755.14
## + x35$x 1 11216529 4.8385e+09 755.17
## + x9$x 1 2032731 4.8477e+09 755.24
## + x11$x 1 2032731 4.8477e+09 755.24
## - x1$x 1 574683489 5.4244e+09 755.51
## - x2$x 1 599760629 5.4495e+09 755.69
## - x20$x 1 715960899 5.5657e+09 756.49
## - x10$x 1 717820453 5.5676e+09 756.50
## - x34$x 1 745228016 5.5950e+09 756.69
## - x4$x 1 757122240 5.6069e+09 756.77
## - x17$x 1 1945476689 6.7952e+09 764.07
## - x21$x 1 1974787719 6.8245e+09 764.24
## - x14$x 1 2032629069 6.8824e+09 764.56
## - x33$x 1 2169128101 7.0189e+09 765.30
## - x25$x 1 3394182933 8.2439e+09 771.42
## - x24$x 1 3464596817 8.3143e+09 771.74
## - x18$x 1 4831892750 9.6816e+09 777.52
## - x7$x 1 4870671098 9.7204e+09 777.68
## - x19$x 1 4980266323 9.8300e+09 778.10
## - x13$x 1 6912068507 1.1762e+10 784.92
##
## Step: AIC=753.2
## y_x ~ x1$x + x2$x + x3$x + x4$x + x7$x + x10$x + x12$x + x13$x +
## x14$x + x17$x + x18$x + x19$x + x20$x + x21$x + x24$x + x25$x +
## x27$x + x28$x + x33$x + x34$x
##
## Df Sum of Sq RSS AIC
## - x3$x 1 172146767 5.2770e+09 752.46
## <none> 5.1048e+09 753.20
## + x16$x 1 255096224 4.8497e+09 753.26
## - x27$x 1 285212000 5.3901e+09 753.27
## + x29$x 1 219110159 4.8857e+09 753.54
## + x31$x 1 155913341 4.9489e+09 754.02
## + x26$x 1 117308627 4.9875e+09 754.32
## - x12$x 1 463182889 5.5680e+09 754.50
## - x28$x 1 482528992 5.5874e+09 754.64
## + x32$x 1 50728844 5.0541e+09 754.82
## + x35$x 1 48799129 5.0560e+09 754.84
## + x8$x 1 45896815 5.0589e+09 754.86
## + x23$x 1 32739206 5.0721e+09 754.96
## - x2$x 1 550091527 5.6549e+09 755.09
## + x9$x 1 3701586 5.1011e+09 755.18
## + x11$x 1 3701586 5.1011e+09 755.18
## + x22$x 1 628492 5.1042e+09 755.20
## - x10$x 1 618691146 5.7235e+09 755.55
## - x34$x 1 646212113 5.7511e+09 755.73
## - x20$x 1 658387723 5.7632e+09 755.81
## - x4$x 1 917918545 6.0228e+09 757.49
## - x1$x 1 970166345 6.0750e+09 757.81
## - x17$x 1 1692303364 6.7971e+09 762.08
## - x21$x 1 2120408338 7.2253e+09 764.40
## - x33$x 1 2223160543 7.3280e+09 764.94
## - x24$x 1 3386340312 8.4912e+09 770.54
## - x25$x 1 3512674607 8.6175e+09 771.10
## - x14$x 1 4667504960 9.7723e+09 775.88
## - x18$x 1 4690297442 9.7951e+09 775.97
## - x7$x 1 4810609503 9.9155e+09 776.43
## - x19$x 1 4836492902 9.9413e+09 776.53
## - x13$x 1 6815510546 1.1920e+10 783.43
##
## Step: AIC=752.46
## y_x ~ x1$x + x2$x + x4$x + x7$x + x10$x + x12$x + x13$x + x14$x +
## x17$x + x18$x + x19$x + x20$x + x21$x + x24$x + x25$x + x27$x +
## x28$x + x33$x + x34$x
##
## Df Sum of Sq RSS AIC
## <none> 5.2770e+09 752.46
## + x29$x 1 209703745 5.0673e+09 752.92
## - x2$x 1 379556768 5.6565e+09 753.10
## - x27$x 1 393047716 5.6700e+09 753.19
## + x3$x 1 172146767 5.1048e+09 753.20
## + x5$x 1 172146767 5.1048e+09 753.20
## + x6$x 1 172146767 5.1048e+09 753.20
## + x31$x 1 116310107 5.1607e+09 753.62
## + x16$x 1 106759774 5.1702e+09 753.69
## - x28$x 1 484293413 5.7613e+09 753.80
## + x26$x 1 87175097 5.1898e+09 753.83
## - x34$x 1 496530734 5.7735e+09 753.88
## + x35$x 1 75637798 5.2014e+09 753.91
## + x23$x 1 58004037 5.2190e+09 754.04
## + x8$x 1 41653935 5.2353e+09 754.16
## - x10$x 1 549351183 5.8263e+09 754.23
## + x32$x 1 23897015 5.2531e+09 754.29
## + x9$x 1 9917175 5.2671e+09 754.39
## + x11$x 1 9917175 5.2671e+09 754.39
## + x22$x 1 8413238 5.2686e+09 754.40
## - x12$x 1 752326579 6.0293e+09 755.53
## - x4$x 1 767831983 6.0448e+09 755.63
## - x20$x 1 785370790 6.0624e+09 755.74
## - x1$x 1 1167207409 6.4442e+09 758.06
## - x17$x 1 1859734558 7.1367e+09 761.94
## - x33$x 1 2052242768 7.3292e+09 762.95
## - x24$x 1 3380845127 8.6578e+09 769.28
## - x25$x 1 3616547073 8.8935e+09 770.30
## - x21$x 1 3854815767 9.1318e+09 771.30
## - x14$x 1 4740562894 1.0018e+10 774.82
## - x18$x 1 4811377888 1.0088e+10 775.09
## - x19$x 1 4903647745 1.0181e+10 775.43
## - x7$x 1 4929257866 1.0206e+10 775.53
## - x13$x 1 7129247097 1.2406e+10 782.95#hasil
lm(formula = y_x ~ x1$x + x2$x + x4$x + x7$x + x10$x + x12$x +
x13$x + x14$x + x17$x + x18$x + x19$x + x20$x + x21$x + x24$x +
x25$x + x27$x + x28$x + x33$x + x34$x)##
## Call:
## lm(formula = y_x ~ x1$x + x2$x + x4$x + x7$x + x10$x + x12$x +
## x13$x + x14$x + x17$x + x18$x + x19$x + x20$x + x21$x + x24$x +
## x25$x + x27$x + x28$x + x33$x + x34$x)
##
## Coefficients:
## (Intercept) x1$x x2$x x4$x x7$x
## 114692.419 10.702 -16.665 1461.718 21251.091
## x10$x x12$x x13$x x14$x x17$x
## 1.158 223.326 1.033 -14.512 1657.706
## x18$x x19$x x20$x x21$x x24$x
## -22601.643 -23024.615 -193.912 718.521 48.719
## x25$x x27$x x28$x x33$x x34$x
## -40.576 34.673 -19.963 -120.294 1.668# 3. Seleksi variabel hasil seleksi
var_sig_1 <- lm(y_x~x1$x+x2$x+x4$x+x7$x+x10$x+x12$x+x13$x+
x14$x+x17$x+x18$x+x19$x+x20$x+x21$x+x24$x+
x25$x+x27$x+x28$x+x33$x+x34$x+x35$x)
summary(var_sig_1) #signifikan x7,x13,x14,x18,x19,x21,x24,x25##
## Call:
## lm(formula = y_x ~ x1$x + x2$x + x4$x + x7$x + x10$x + x12$x +
## x13$x + x14$x + x17$x + x18$x + x19$x + x20$x + x21$x + x24$x +
## x25$x + x27$x + x28$x + x33$x + x34$x + x35$x)
##
## Residuals:
## Min 1Q Median 3Q Max
## -25207.7 -8027.3 87.7 8455.4 24511.9
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.155e+05 7.676e+03 15.051 2.93e-11 ***
## x1$x 9.337e+00 6.128e+00 1.524 0.14595
## x2$x -1.435e+01 1.567e+01 -0.916 0.37269
## x4$x 1.113e+03 1.159e+03 0.961 0.35014
## x7$x 2.083e+04 5.363e+03 3.884 0.00119 **
## x10$x 1.122e+00 8.676e-01 1.293 0.21330
## x12$x 2.063e+02 1.465e+02 1.408 0.17702
## x13$x 1.103e+00 2.569e-01 4.295 0.00049 ***
## x14$x -1.410e+01 3.779e+00 -3.731 0.00166 **
## x17$x 1.553e+03 7.045e+02 2.204 0.04156 *
## x18$x -2.204e+04 5.809e+03 -3.795 0.00145 **
## x19$x -2.247e+04 5.859e+03 -3.835 0.00133 **
## x20$x -1.740e+02 1.275e+02 -1.364 0.19023
## x21$x 6.805e+02 2.164e+02 3.145 0.00590 **
## x24$x 4.858e+01 1.466e+01 3.314 0.00410 **
## x25$x -4.031e+01 1.181e+01 -3.412 0.00332 **
## x27$x 3.153e+01 3.124e+01 1.009 0.32696
## x28$x -1.768e+01 1.652e+01 -1.070 0.29944
## x33$x -1.293e+02 4.984e+01 -2.594 0.01892 *
## x34$x 1.915e+00 1.400e+00 1.367 0.18929
## x35$x -2.962e-01 5.958e-01 -0.497 0.62541
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 17490 on 17 degrees of freedom
## Multiple R-squared: 0.9905, Adjusted R-squared: 0.9794
## F-statistic: 89.06 on 20 and 17 DF, p-value: 1.989e-13var_sig_2 <- lm(y_x~x1$x+x2$x+x4$x+x7$x+x10$x+x12$x+x13$x+
x14$x+x17$x+x18$x+x19$x+x20$x+x21$x+x24$x+
x25$x+x27$x+x28$x+x33$x+x34$x)
summary(var_sig_2) #signifikan x7,x13,x14,x18,x19,x21,x24,x25##
## Call:
## lm(formula = y_x ~ x1$x + x2$x + x4$x + x7$x + x10$x + x12$x +
## x13$x + x14$x + x17$x + x18$x + x19$x + x20$x + x21$x + x24$x +
## x25$x + x27$x + x28$x + x33$x + x34$x)
##
## Residuals:
## Min 1Q Median 3Q Max
## -25292.6 -8144.7 -391.5 7785.3 26311.2
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.147e+05 7.329e+03 15.650 6.33e-12 ***
## x1$x 1.070e+01 5.363e+00 1.995 0.061368 .
## x2$x -1.667e+01 1.465e+01 -1.138 0.270106
## x4$x 1.462e+03 9.032e+02 1.618 0.122973
## x7$x 2.125e+04 5.183e+03 4.100 0.000672 ***
## x10$x 1.158e+00 8.463e-01 1.369 0.187876
## x12$x 2.233e+02 1.394e+02 1.602 0.126572
## x13$x 1.032e+00 2.094e-01 4.931 0.000108 ***
## x14$x -1.451e+01 3.609e+00 -4.021 0.000801 ***
## x17$x 1.658e+03 6.582e+02 2.519 0.021459 *
## x18$x -2.260e+04 5.579e+03 -4.051 0.000749 ***
## x19$x -2.302e+04 5.630e+03 -4.090 0.000688 ***
## x20$x -1.939e+02 1.185e+02 -1.637 0.119048
## x21$x 7.185e+02 1.982e+02 3.626 0.001932 **
## x24$x 4.872e+01 1.435e+01 3.396 0.003221 **
## x25$x -4.058e+01 1.155e+01 -3.512 0.002488 **
## x27$x 3.467e+01 2.995e+01 1.158 0.262044
## x28$x -1.996e+01 1.553e+01 -1.285 0.214987
## x33$x -1.203e+02 4.547e+01 -2.646 0.016434 *
## x34$x 1.668e+00 1.281e+00 1.301 0.209528
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 17120 on 18 degrees of freedom
## Multiple R-squared: 0.9904, Adjusted R-squared: 0.9803
## F-statistic: 97.83 on 19 and 18 DF, p-value: 2.208e-14var_sig_3 <- lm(y_x~x1$x+x4$x+x7$x+x10$x+x12$x+x13$x+
x14$x+x17$x+x18$x+x19$x+x20$x+x21$x+x24$x+
x25$x+x27$x+x28$x+x33$x+x34$x)
summary(var_sig_3) #signifikan x7,x13,x14,x18,x19,x21,x24,x25##
## Call:
## lm(formula = y_x ~ x1$x + x4$x + x7$x + x10$x + x12$x + x13$x +
## x14$x + x17$x + x18$x + x19$x + x20$x + x21$x + x24$x + x25$x +
## x27$x + x28$x + x33$x + x34$x)
##
## Residuals:
## Min 1Q Median 3Q Max
## -26743.9 -7578.6 238.3 8728.9 28424.2
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.134e+05 7.291e+03 15.549 2.92e-12 ***
## x1$x 5.869e+00 3.300e+00 1.778 0.091346 .
## x4$x 1.179e+03 8.750e+02 1.347 0.193766
## x7$x 1.860e+04 4.664e+03 3.988 0.000788 ***
## x10$x 1.132e+00 8.525e-01 1.328 0.199939
## x12$x 1.711e+02 1.326e+02 1.290 0.212640
## x13$x 1.001e+00 2.091e-01 4.786 0.000128 ***
## x14$x -1.284e+01 3.321e+00 -3.866 0.001041 **
## x17$x 1.343e+03 6.018e+02 2.231 0.037895 *
## x18$x -1.968e+04 4.991e+03 -3.943 0.000873 ***
## x19$x -2.008e+04 5.039e+03 -3.985 0.000793 ***
## x20$x -1.601e+02 1.156e+02 -1.385 0.181977
## x21$x 6.231e+02 1.809e+02 3.444 0.002719 **
## x24$x 4.503e+01 1.408e+01 3.198 0.004740 **
## x25$x -3.679e+01 1.115e+01 -3.300 0.003765 **
## x27$x 3.321e+01 3.015e+01 1.101 0.284481
## x28$x -1.458e+01 1.491e+01 -0.978 0.340479
## x33$x -1.098e+02 4.487e+01 -2.448 0.024268 *
## x34$x 1.475e+00 1.280e+00 1.152 0.263540
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 17250 on 19 degrees of freedom
## Multiple R-squared: 0.9897, Adjusted R-squared: 0.98
## F-statistic: 101.6 on 18 and 19 DF, p-value: 3.96e-15var_sig_4 <- lm(y_x~x1$x+x4$x+x7$x+x10$x+x12$x+x13$x+
x14$x+x17$x+x18$x+x19$x+x20$x+x21$x+x24$x+
x25$x+x27$x+x33$x+x34$x)
summary(var_sig_4) #signifikan x7,x13,x14,x18,x19,x21,x24,x25##
## Call:
## lm(formula = y_x ~ x1$x + x4$x + x7$x + x10$x + x12$x + x13$x +
## x14$x + x17$x + x18$x + x19$x + x20$x + x21$x + x24$x + x25$x +
## x27$x + x33$x + x34$x)
##
## Residuals:
## Min 1Q Median 3Q Max
## -24044.9 -8264.8 -41.4 8474.1 29903.4
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.139e+05 7.264e+03 15.677 1.06e-12 ***
## x1$x 4.832e+00 3.122e+00 1.548 0.137321
## x4$x 8.068e+02 7.871e+02 1.025 0.317599
## x7$x 1.720e+04 4.434e+03 3.879 0.000934 ***
## x10$x 1.052e+00 8.476e-01 1.241 0.229065
## x12$x 1.233e+02 1.232e+02 1.001 0.328762
## x13$x 9.853e-01 2.083e-01 4.731 0.000128 ***
## x14$x -1.205e+01 3.218e+00 -3.744 0.001278 **
## x17$x 1.193e+03 5.814e+02 2.053 0.053433 .
## x18$x -1.820e+04 4.751e+03 -3.831 0.001044 **
## x19$x -1.863e+04 4.811e+03 -3.873 0.000947 ***
## x20$x -1.292e+02 1.110e+02 -1.163 0.258417
## x21$x 5.941e+02 1.783e+02 3.333 0.003318 **
## x24$x 3.925e+01 1.277e+01 3.074 0.005986 **
## x25$x -3.245e+01 1.022e+01 -3.177 0.004743 **
## x27$x 3.507e+01 3.005e+01 1.167 0.257018
## x33$x -1.052e+02 4.457e+01 -2.360 0.028538 *
## x34$x 1.822e+00 1.228e+00 1.483 0.153563
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 17240 on 20 degrees of freedom
## Multiple R-squared: 0.9892, Adjusted R-squared: 0.98
## F-statistic: 107.8 on 17 and 20 DF, p-value: 5.89e-16var_sig_5 <- lm(y_x~x1$x+x4$x+x7$x+x10$x+x13$x+
x14$x+x17$x+x18$x+x19$x+x20$x+x21$x+x24$x+
x25$x+x27$x+x33$x+x34$x)
summary(var_sig_5) #signifikan x7,x13,x14,x18,x19,x21,x24,x25,x34##
## Call:
## lm(formula = y_x ~ x1$x + x4$x + x7$x + x10$x + x13$x + x14$x +
## x17$x + x18$x + x19$x + x20$x + x21$x + x24$x + x25$x + x27$x +
## x33$x + x34$x)
##
## Residuals:
## Min 1Q Median 3Q Max
## -22815.5 -9634.5 -758.5 7766.0 30639.9
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.142e+05 7.258e+03 15.732 4.28e-13 ***
## x1$x 4.559e+00 3.110e+00 1.466 0.157493
## x4$x 4.329e+02 6.929e+02 0.625 0.538825
## x7$x 1.496e+04 3.829e+03 3.907 0.000811 ***
## x10$x 8.902e-01 8.321e-01 1.070 0.296857
## x13$x 9.852e-01 2.083e-01 4.730 0.000114 ***
## x14$x -1.084e+01 2.983e+00 -3.634 0.001553 **
## x17$x 1.023e+03 5.561e+02 1.840 0.079881 .
## x18$x -1.567e+04 4.024e+03 -3.894 0.000836 ***
## x19$x -1.610e+04 4.093e+03 -3.934 0.000761 ***
## x20$x -1.470e+02 1.096e+02 -1.341 0.194194
## x21$x 5.290e+02 1.660e+02 3.187 0.004436 **
## x24$x 3.500e+01 1.204e+01 2.907 0.008439 **
## x25$x -2.912e+01 9.660e+00 -3.015 0.006595 **
## x27$x 2.493e+01 2.830e+01 0.881 0.388256
## x33$x -1.064e+02 4.456e+01 -2.387 0.026473 *
## x34$x 2.636e+00 9.205e-01 2.864 0.009298 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 17240 on 21 degrees of freedom
## Multiple R-squared: 0.9887, Adjusted R-squared: 0.98
## F-statistic: 114.4 on 16 and 21 DF, p-value: < 2.2e-16var_sig_6 <- lm(y_x~x1$x+x4$x+x7$x+x10$x+x13$x+
x14$x+x17$x+x18$x+x19$x+x20$x+x21$x+x24$x+
x25$x+x33$x+x34$x)
summary(var_sig_6) #signifikan x7,x13,x14,x18,x19,x21,x25,x34##
## Call:
## lm(formula = y_x ~ x1$x + x4$x + x7$x + x10$x + x13$x + x14$x +
## x17$x + x18$x + x19$x + x20$x + x21$x + x24$x + x25$x + x33$x +
## x34$x)
##
## Residuals:
## Min 1Q Median 3Q Max
## -22589 -10456 155 6506 31826
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.146e+05 7.202e+03 15.917 1.48e-13 ***
## x1$x 4.030e+00 3.036e+00 1.328 0.197935
## x4$x 1.853e+02 6.301e+02 0.294 0.771455
## x7$x 1.469e+04 3.798e+03 3.869 0.000829 ***
## x10$x 1.012e+00 8.163e-01 1.240 0.227917
## x13$x 1.004e+00 2.061e-01 4.871 7.21e-05 ***
## x14$x -1.073e+01 2.965e+00 -3.618 0.001525 **
## x17$x 7.719e+02 4.748e+02 1.626 0.118225
## x18$x -1.512e+04 3.956e+03 -3.823 0.000927 ***
## x19$x -1.558e+04 4.029e+03 -3.866 0.000835 ***
## x20$x -1.784e+02 1.031e+02 -1.731 0.097531 .
## x21$x 5.970e+02 1.462e+02 4.083 0.000493 ***
## x24$x 3.331e+01 1.183e+01 2.816 0.010056 *
## x25$x -2.761e+01 9.458e+00 -2.919 0.007945 **
## x33$x -1.018e+02 4.403e+01 -2.312 0.030493 *
## x34$x 3.039e+00 7.947e-01 3.824 0.000925 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 17150 on 22 degrees of freedom
## Multiple R-squared: 0.9882, Adjusted R-squared: 0.9802
## F-statistic: 123.3 on 15 and 22 DF, p-value: < 2.2e-16var_sig_7 <- lm(y_x~x1$x+x4$x+x7$x+x13$x+x14$x+x17$x+x18$x+x19$x
+x20$x+x21$x+x24$x+x25$x+x33$x+x34$x)
summary(var_sig_7) #signifikan x7,x13,x14,x18,x19,x21,x25,x34##
## Call:
## lm(formula = y_x ~ x1$x + x4$x + x7$x + x13$x + x14$x + x17$x +
## x18$x + x19$x + x20$x + x21$x + x24$x + x25$x + x33$x + x34$x)
##
## Residuals:
## Min 1Q Median 3Q Max
## -21092 -10013 272 6927 34331
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.132e+05 7.184e+03 15.750 8.2e-14 ***
## x1$x 4.019e+00 3.071e+00 1.309 0.203616
## x4$x -1.762e+02 5.651e+02 -0.312 0.757950
## x7$x 1.304e+04 3.598e+03 3.625 0.001421 **
## x13$x 9.244e-01 1.981e-01 4.665 0.000107 ***
## x14$x -1.011e+01 2.957e+00 -3.419 0.002349 **
## x17$x 4.492e+02 4.017e+02 1.118 0.275102
## x18$x -1.317e+04 3.671e+03 -3.588 0.001556 **
## x19$x -1.357e+04 3.733e+03 -3.636 0.001384 **
## x20$x -2.046e+02 1.021e+02 -2.005 0.056921 .
## x21$x 5.422e+02 1.410e+02 3.845 0.000826 ***
## x24$x 3.233e+01 1.194e+01 2.708 0.012547 *
## x25$x -2.717e+01 9.561e+00 -2.842 0.009238 **
## x33$x -8.397e+01 4.210e+01 -1.995 0.058059 .
## x34$x 2.923e+00 7.983e-01 3.662 0.001299 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 17350 on 23 degrees of freedom
## Multiple R-squared: 0.9874, Adjusted R-squared: 0.9798
## F-statistic: 128.9 on 14 and 23 DF, p-value: < 2.2e-16var_sig_8 <- lm(y_x~x1$x+x7$x+x13$x+x14$x+x17$x+x18$x+x19$x
+x20$x+x21$x+x24$x+x25$x+x33$x+x34$x)
summary(var_sig_8) #signifikan x7,x13,x14,x18,x19,x21,x25,x34##
## Call:
## lm(formula = y_x ~ x1$x + x7$x + x13$x + x14$x + x17$x + x18$x +
## x19$x + x20$x + x21$x + x24$x + x25$x + x33$x + x34$x)
##
## Residuals:
## Min 1Q Median 3Q Max
## -21449 -10307 -1380 7324 33754
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.131e+05 7.047e+03 16.052 2.46e-14 ***
## x1$x 4.165e+00 2.978e+00 1.399 0.174731
## x7$x 1.329e+04 3.438e+03 3.867 0.000737 ***
## x13$x 9.197e-01 1.938e-01 4.745 7.94e-05 ***
## x14$x -1.026e+01 2.861e+00 -3.587 0.001484 **
## x17$x 5.008e+02 3.591e+02 1.394 0.175976
## x18$x -1.348e+04 3.469e+03 -3.886 0.000703 ***
## x19$x -1.388e+04 3.531e+03 -3.930 0.000629 ***
## x20$x -2.069e+02 9.990e+01 -2.071 0.049259 *
## x21$x 5.519e+02 1.349e+02 4.092 0.000417 ***
## x24$x 3.251e+01 1.170e+01 2.780 0.010409 *
## x25$x -2.724e+01 9.378e+00 -2.904 0.007781 **
## x33$x -8.327e+01 4.124e+01 -2.019 0.054772 .
## x34$x 2.909e+00 7.819e-01 3.720 0.001065 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 17020 on 24 degrees of freedom
## Multiple R-squared: 0.9874, Adjusted R-squared: 0.9805
## F-statistic: 144.3 on 13 and 24 DF, p-value: < 2.2e-16var_sig_9 <- lm(y_x~x1$x+x7$x+x13$x+x14$x+x18$x+x19$x
+x20$x+x21$x+x24$x+x25$x+x33$x+x34$x)
summary(var_sig_9) #signifikan x7,x13,x14,x18,x19,x21,x25,x34##
## Call:
## lm(formula = y_x ~ x1$x + x7$x + x13$x + x14$x + x18$x + x19$x +
## x20$x + x21$x + x24$x + x25$x + x33$x + x34$x)
##
## Residuals:
## Min 1Q Median 3Q Max
## -23027 -11137 -219 8748 35613
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.115e+05 7.080e+03 15.748 1.74e-14 ***
## x1$x 6.585e+00 2.465e+00 2.672 0.013092 *
## x7$x 1.233e+04 3.431e+03 3.594 0.001393 **
## x13$x 8.172e-01 1.827e-01 4.473 0.000146 ***
## x14$x -1.006e+01 2.911e+00 -3.458 0.001964 **
## x18$x -1.207e+04 3.379e+03 -3.570 0.001479 **
## x19$x -1.250e+04 3.453e+03 -3.619 0.001308 **
## x20$x -2.299e+02 1.004e+02 -2.291 0.030673 *
## x21$x 5.093e+02 1.338e+02 3.806 0.000814 ***
## x24$x 3.063e+01 1.184e+01 2.588 0.015866 *
## x25$x -2.492e+01 9.402e+00 -2.650 0.013741 *
## x33$x -6.767e+01 4.044e+01 -1.674 0.106693
## x34$x 2.619e+00 7.679e-01 3.411 0.002206 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 17340 on 25 degrees of freedom
## Multiple R-squared: 0.9863, Adjusted R-squared: 0.9798
## F-statistic: 150.5 on 12 and 25 DF, p-value: < 2.2e-16var_sig_10 <- lm(y_x~x1$x+x7$x+x13$x+x14$x+x18$x+x19$x
+x20$x+x21$x+x24$x+x25$x+x34$x)
summary(var_sig_10) #signifikan x7,x13,x14,x18,x19,x20,x21,x34##
## Call:
## lm(formula = y_x ~ x1$x + x7$x + x13$x + x14$x + x18$x + x19$x +
## x20$x + x21$x + x24$x + x25$x + x34$x)
##
## Residuals:
## Min 1Q Median 3Q Max
## -28451 -8927 -3035 10740 41140
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.075e+05 6.897e+03 15.589 1.04e-14 ***
## x1$x 4.553e+00 2.218e+00 2.053 0.05030 .
## x7$x 9.040e+03 2.907e+03 3.110 0.00450 **
## x13$x 5.170e-01 3.590e-02 14.400 6.67e-14 ***
## x14$x -7.808e+00 2.668e+00 -2.927 0.00702 **
## x18$x -8.901e+03 2.896e+03 -3.073 0.00492 **
## x19$x -9.262e+03 2.958e+03 -3.130 0.00428 **
## x20$x -3.122e+02 9.049e+01 -3.450 0.00193 **
## x21$x 4.474e+02 1.330e+02 3.364 0.00239 **
## x24$x 3.231e+01 1.219e+01 2.650 0.01353 *
## x25$x -2.489e+01 9.722e+00 -2.560 0.01663 *
## x34$x 2.494e+00 7.903e-01 3.156 0.00402 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 17930 on 26 degrees of freedom
## Multiple R-squared: 0.9848, Adjusted R-squared: 0.9784
## F-statistic: 153.3 on 11 and 26 DF, p-value: < 2.2e-16var_sig_11 <- lm(y_x~x7$x+x13$x+x14$x+x18$x+x19$x
+x20$x+x21$x+x24$x+x25$x+x34$x)
summary(var_sig_11) #signifikan x7,x13,x14,x18,x19,x20,x21,x24,x25,x34##
## Call:
## lm(formula = y_x ~ x7$x + x13$x + x14$x + x18$x + x19$x + x20$x +
## x21$x + x24$x + x25$x + x34$x)
##
## Residuals:
## Min 1Q Median 3Q Max
## -28708 -10103 -1656 10669 43594
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.078e+05 7.295e+03 14.772 1.86e-14 ***
## x7$x 8.252e+03 3.048e+03 2.707 0.01163 *
## x13$x 5.733e-01 2.453e-02 23.365 < 2e-16 ***
## x14$x -7.224e+00 2.806e+00 -2.575 0.01584 *
## x18$x -8.090e+03 3.035e+03 -2.666 0.01281 *
## x19$x -8.465e+03 3.103e+03 -2.728 0.01106 *
## x20$x -3.019e+02 9.558e+01 -3.159 0.00388 **
## x21$x 4.493e+02 1.407e+02 3.194 0.00355 **
## x24$x 3.818e+01 1.254e+01 3.044 0.00515 **
## x25$x -3.011e+01 9.926e+00 -3.033 0.00529 **
## x34$x 2.889e+00 8.108e-01 3.563 0.00139 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 18960 on 27 degrees of freedom
## Multiple R-squared: 0.9824, Adjusted R-squared: 0.9758
## F-statistic: 150.3 on 10 and 27 DF, p-value: < 2.2e-16cek2<-step(lm(y_x~x7$x+x13$x+x14$x+x18$x+x19$x
+x20$x+x21$x+x24$x+x25$x+x34$x),direction="backward")## Start: AIC=757.64
## y_x ~ x7$x + x13$x + x14$x + x18$x + x19$x + x20$x + x21$x +
## x24$x + x25$x + x34$x
##
## Df Sum of Sq RSS AIC
## <none> 9.7102e+09 757.64
## - x14$x 1 2.3840e+09 1.2094e+10 763.98
## - x18$x 1 2.5555e+09 1.2266e+10 764.51
## - x7$x 1 2.6352e+09 1.2345e+10 764.76
## - x19$x 1 2.6772e+09 1.2387e+10 764.89
## - x25$x 1 3.3093e+09 1.3020e+10 766.78
## - x24$x 1 3.3333e+09 1.3044e+10 766.85
## - x20$x 1 3.5879e+09 1.3298e+10 767.59
## - x21$x 1 3.6687e+09 1.3379e+10 767.82
## - x34$x 1 4.5650e+09 1.4275e+10 770.28
## - x13$x 1 1.9634e+11 2.0605e+11 871.72