HW7

knitr::knit_hooks$set(webgl = hook_webgl) 
# 讀檔案
dta7 <- read.table("hs0.txt", header = TRUE)

# 設定ssq function
ssq <- function(mu, sigma, y) {sum(((y - mu) / sigma)^2)}
vssq <- Vectorize(ssq, c("mu", "sigma"))

# 設定最小化的vssq
functionVectorize <- function(data = NA){
  x <- seq(mean(data) - 5, mean(data) + 5, by = 0.1)
  y <- sd(data)
  z <- vssq(x, sd(data), data)
  rgl::plot3d(x, y, z, 
              xlab = "mu", ylab = "sigma(control)", zlab = "ssq")
}

functionVectorize(data = dta7$math)

You must enable Javascript to view this page properly.