data <- c(12, 45, 52, 58, 61, 63, 67, 70, 72, 75, 78, 82, 88, 95, 310)
# Fungsi untuk sorting
sort_function <- function(data){
sorted_data <- sort(data)
return(sorted_data)
}
# Fungsi untuk menghitung mean
mean_function <- function(data){
mean_value <- mean(data)
return(mean_value)
}
# Mengurutkan data
data <- sort_function(data)
# Fungsi Winsorization
winsor_function <- function(data, alpha){
data_sort <- sort_function(data)
n <- length(data_sort)
k <- floor(n * alpha)
winsor_data <- data_sort
# Winsorize bagian bawah
for(i in 1:k){
winsor_data[i] <- data_sort[k+1]
}
# Winsorize bagian atas
for(i in (n-k+1):n){
winsor_data[i] <- data_sort[n-k]
}
return(winsor_data)
}
# Fungsi untuk Winsorized mean
winsor_mean <- function(data, alpha){
data_winsor <- winsor_function(data, alpha)
mean_w <- mean_function(data_winsor)
return(list(
winsor_data = data_winsor,
winsor_mean = mean_w
))
}
# Hasil Winsorized dengan alpha 0.2 dan 0
result_02 <- winsor_mean(data, 0.2)
result_00 <- winsor_mean(data, 0)
# Menampilkan hasil
cat("Winsorized data alpha 0.2 :", result_02$winsor_data, "\n")
Winsorized data alpha 0.2 : 58 58 58 58 61 63 67 70 72 75 78 82 82 82 82
cat("Winsorized mean alpha 0.2 :", result_02$winsor_mean, "\n\n")
Winsorized mean alpha 0.2 : 69.73333
cat("Winsorized data alpha 0 :", result_00$winsor_data, "\n")
Winsorized data alpha 0 : 12 45 52 58 61 63 67 70 72 75 78 82 88 95 310 310
cat("Winsorized mean alpha 0 :", result_00$winsor_mean, "\n")
Winsorized mean alpha 0 : 96.125
Add a new chunk by clicking the Insert Chunk button on the
toolbar or by pressing Ctrl+Alt+I.
When you save the notebook, an HTML file containing the code and
output will be saved alongside it (click the Preview button or
press Ctrl+Shift+K to preview the HTML file).
The preview shows you a rendered HTML copy of the contents of the
editor. Consequently, unlike Knit, Preview does not
run any R code chunks. Instead, the output of the chunk when it was last
run in the editor is displayed.
LS0tDQp0aXRsZTogIlIgTm90ZWJvb2siDQpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sNCi0tLQ0KDQoNCmBgYHtyfQ0KZGF0YSA8LSBjKDEyLCA0NSwgNTIsIDU4LCA2MSwgNjMsIDY3LCA3MCwgNzIsIDc1LCA3OCwgODIsIDg4LCA5NSwgMzEwKQ0KDQojIEZ1bmdzaSB1bnR1ayBzb3J0aW5nDQpzb3J0X2Z1bmN0aW9uIDwtIGZ1bmN0aW9uKGRhdGEpew0KICBzb3J0ZWRfZGF0YSA8LSBzb3J0KGRhdGEpDQogIHJldHVybihzb3J0ZWRfZGF0YSkNCn0NCg0KIyBGdW5nc2kgdW50dWsgbWVuZ2hpdHVuZyBtZWFuDQptZWFuX2Z1bmN0aW9uIDwtIGZ1bmN0aW9uKGRhdGEpew0KICBtZWFuX3ZhbHVlIDwtIG1lYW4oZGF0YSkNCiAgcmV0dXJuKG1lYW5fdmFsdWUpDQp9DQoNCiMgTWVuZ3VydXRrYW4gZGF0YQ0KZGF0YSA8LSBzb3J0X2Z1bmN0aW9uKGRhdGEpDQoNCiMgRnVuZ3NpIFdpbnNvcml6YXRpb24NCndpbnNvcl9mdW5jdGlvbiA8LSBmdW5jdGlvbihkYXRhLCBhbHBoYSl7DQogIA0KICBkYXRhX3NvcnQgPC0gc29ydF9mdW5jdGlvbihkYXRhKQ0KICANCiAgbiA8LSBsZW5ndGgoZGF0YV9zb3J0KQ0KICBrIDwtIGZsb29yKG4gKiBhbHBoYSkNCiAgDQogIHdpbnNvcl9kYXRhIDwtIGRhdGFfc29ydA0KICANCiAgIyBXaW5zb3JpemUgYmFnaWFuIGJhd2FoDQogIGZvcihpIGluIDE6ayl7DQogICAgd2luc29yX2RhdGFbaV0gPC0gZGF0YV9zb3J0W2srMV0NCiAgfQ0KICANCiAgIyBXaW5zb3JpemUgYmFnaWFuIGF0YXMNCiAgZm9yKGkgaW4gKG4taysxKTpuKXsNCiAgICB3aW5zb3JfZGF0YVtpXSA8LSBkYXRhX3NvcnRbbi1rXQ0KICB9DQogIA0KICByZXR1cm4od2luc29yX2RhdGEpDQp9DQoNCiMgRnVuZ3NpIHVudHVrIFdpbnNvcml6ZWQgbWVhbg0Kd2luc29yX21lYW4gPC0gZnVuY3Rpb24oZGF0YSwgYWxwaGEpew0KICANCiAgZGF0YV93aW5zb3IgPC0gd2luc29yX2Z1bmN0aW9uKGRhdGEsIGFscGhhKQ0KICBtZWFuX3cgPC0gbWVhbl9mdW5jdGlvbihkYXRhX3dpbnNvcikNCiAgDQogIHJldHVybihsaXN0KA0KICAgIHdpbnNvcl9kYXRhID0gZGF0YV93aW5zb3IsDQogICAgd2luc29yX21lYW4gPSBtZWFuX3cNCiAgKSkNCn0NCg0KIyBIYXNpbCBXaW5zb3JpemVkIGRlbmdhbiBhbHBoYSAwLjIgZGFuIDANCnJlc3VsdF8wMiA8LSB3aW5zb3JfbWVhbihkYXRhLCAwLjIpDQpyZXN1bHRfMDAgPC0gd2luc29yX21lYW4oZGF0YSwgMCkNCg0KIyBNZW5hbXBpbGthbiBoYXNpbA0KY2F0KCJXaW5zb3JpemVkIGRhdGEgYWxwaGEgMC4yIDoiLCByZXN1bHRfMDIkd2luc29yX2RhdGEsICJcbiIpDQoNCmNhdCgiV2luc29yaXplZCBtZWFuIGFscGhhIDAuMiA6IiwgcmVzdWx0XzAyJHdpbnNvcl9tZWFuLCAiXG5cbiIpDQoNCmNhdCgiV2luc29yaXplZCBkYXRhIGFscGhhIDAgOiIsIHJlc3VsdF8wMCR3aW5zb3JfZGF0YSwgIlxuIikNCg0KY2F0KCJXaW5zb3JpemVkIG1lYW4gYWxwaGEgMCA6IiwgcmVzdWx0XzAwJHdpbnNvcl9tZWFuLCAiXG4iKQ0KYGBgDQoNCkFkZCBhIG5ldyBjaHVuayBieSBjbGlja2luZyB0aGUgKkluc2VydCBDaHVuayogYnV0dG9uIG9uIHRoZSB0b29sYmFyIG9yIGJ5IHByZXNzaW5nICpDdHJsK0FsdCtJKi4NCg0KV2hlbiB5b3Ugc2F2ZSB0aGUgbm90ZWJvb2ssIGFuIEhUTUwgZmlsZSBjb250YWluaW5nIHRoZSBjb2RlIGFuZCBvdXRwdXQgd2lsbCBiZSBzYXZlZCBhbG9uZ3NpZGUgaXQgKGNsaWNrIHRoZSAqUHJldmlldyogYnV0dG9uIG9yIHByZXNzICpDdHJsK1NoaWZ0K0sqIHRvIHByZXZpZXcgdGhlIEhUTUwgZmlsZSkuDQoNClRoZSBwcmV2aWV3IHNob3dzIHlvdSBhIHJlbmRlcmVkIEhUTUwgY29weSBvZiB0aGUgY29udGVudHMgb2YgdGhlIGVkaXRvci4gQ29uc2VxdWVudGx5LCB1bmxpa2UgKktuaXQqLCAqUHJldmlldyogZG9lcyBub3QgcnVuIGFueSBSIGNvZGUgY2h1bmtzLiBJbnN0ZWFkLCB0aGUgb3V0cHV0IG9mIHRoZSBjaHVuayB3aGVuIGl0IHdhcyBsYXN0IHJ1biBpbiB0aGUgZWRpdG9yIGlzIGRpc3BsYXllZC4NCg0KDQo=