R_tilde の変更点 - アールメカブ

アールメカブ


R_tilde の変更点


[[Rの備忘録]]


[[Lattice Multivariate Data Visualization with R:http://www.springer.com/statistics/computational/book/978-0-387-75968-5]] 10.1  ### p.166 \section{非標準評価} に

 form <- ~ boxcox.trans(gcsescore, p) | gender
という自分では使ったことのないコードがある.ただし,ここでの事例は,動作しないということで紹介されている.

10.1  ### p.166 \section{非標準評価}




 boxcox.trans <- function(x, lambda) {
      if (lambda == 0) log(x) else (x^lambda - 1) / lambda
  }
 
 data(Chem97, package = "mlmRev")
 
 trellis.device(pdf, file = "Chem97BoxCox.pdf", 
                 width = 8, height = 6)
 
 for (p in seq(0, 3, by = 0.5)) {
      plot(qqmath(~boxcox.trans(gcsescore, p) | gender,
    data = Chem97, 
        groups = score, f.value = ppoints(100), 
        main = as.expression(substitute(lambda == v, 
                                        list(v = p)))))
  }
 dev.off()
 
# 以下は動作しない

 form <- ~ boxcox.trans(gcsescore, p) | gender
 qqboxcox <- function(lambda) {
      for (p in lambda)
          plot(qqmath(form, data = Chem97, 
             groups = score, f.value = ppoints(100), 
             main = as.expression(substitute(lambda == v, 
                                             list(v = p)))))
  }
 qqboxcox(lambda = seq(0, 3, by = 0.5))