ClojureとRの連携による統計的学習

>install.packages("Rserve", repos="http://cran.rstudio.com/")

>R CMD Rserve
R version 4.2.1 (2022-06-23) -- "Funny-Looking Kid"
Copyright (C) 2022 The R Foundation for Statistical Computing
Platform: aarch64-apple-darwin21.5.0 (64-bit)

R は、自由なソフトウェアであり、「完全に無保証」です。
一定の条件に従えば、自由にこれを再配布することができます。
配布条件の詳細に関しては、'license()' あるいは 'licence()' と入力してください。

R は多くの貢献者による共同プロジェクトです。
詳しくは 'contributors()' と入力してください。
また、R や R のパッケージを出版物で引用する際の形式については
'citation()' と入力してください。

'demo()' と入力すればデモをみることができます。
'help()' とすればオンラインヘルプが出ます。
'help.start()' で HTML ブラウザによるヘルプがみられます。
'q()' と入力すれば R を終了します。

Rserv started in daemon mode.
>

ps ax | grep Rserve　
 3021   ??  Ss     0:00.00 /opt/homebrew/Cellar/r/4.2.1_2/lib/R/bin/Rserve
 3059 s001  R+     0:00.00 grep Rserve

> kill 3021

(defproject cljisr-test01 "0.1.0-SNAPSHOT"
  :description "FIXME: write description"
  :url "http://example.com/FIXME"
  :license {:name "EPL-2.0 OR GPL-2.0-or-later WITH Classpath-exception-2.0"
            :url "https://www.eclipse.org/legal/epl-2.0/"}
  :dependencies [[org.clojure/clojure "1.11.1"]
                 [scicloj/clojisr "1.0.0-BETA20"]
                 [techascent/tech.ml.dataset "6.094"]]
  :repl-options {:init-ns cljisr-test01.core}
  :jvm-opts ["-Dclojure.tools.logging.factory=clojure.tools.logging.impl/jul-factory"])

(ns cljisr-test01.core)

(require '[clojisr.v1.r :as r :refer
            [r eval-r->java r->java java->r java->clj java->native-clj
             clj->java r->clj clj->r ->code r+ colon require-r]]
          '[clojisr.v1.robject :as robject]
          '[clojisr.v1.session :as session]
          '[tech.v3.dataset :as dataset])


(r "mean(rnorm(100000,mean=1.0,sd=3.0))").  ;; 0.9851342

(r
   "abc <- runif(1000);
          f <- function(x) {mean(log(x))};  ;; -1.001936
          f(abc)")

(def x (r "1+2"))

(->> x
      r->clj).  ;; [3.0]

(def f (r "function(x) x*10"))

(->> 5
      f
      r->clj).    ;; [50.0]

(->> "5*5"
      r
      f
      r->clj).    ;; [250.0]

(require-r '[graphics :refer [plot hist]])
 (require-r '[ggplot2 :refer [ggplot aes geom_point xlab ylab labs]])
 (require '[clojisr.v1.applications.plotting :refer
            [plot->svg plot->file plot->buffered-image]])

(spit "out/out-test01.svg" (plot->svg (fn []
              (->> rand
                   (repeatedly 30)
                   (reductions +)
                   (plot :xlab "t" :ylab "y" :type "l")))))


(spit "out/pout-test02.svg" (plot->svg
   (let [x (repeatedly 99 rand)
         y (map + x (repeatedly 99 rand))]
     (-> {:x x, :y y}
         dataset/->dataset
         (ggplot (aes :x x :y y :color '(+ x y) :size '(/ x y)))
         (r+ (geom_point) (xlab "x") (ylab "y"))))))

(defproject roj-test01 "0.1.0-SNAPSHOT"
  :description "FIXME: write description"
  :url "http://example.com/FIXME"
  :license {:name "EPL-2.0 OR GPL-2.0-or-later WITH Classpath-exception-2.0"
            :url "https://www.eclipse.org/legal/epl-2.0/"}
  :dependencies [[org.clojure/clojure "1.11.1"]
                 [rojure "0.2.0"]].      ;;<-add
  :repl-options {:init-ns roj-test01.core})

(ns roj-test01.core)

(use '(rojure core))

(def r (get-r))

(r-eval r "data(iris)")

(r-eval r "iris")

(r-eval r "x11()")

(r-eval r "plot(Sepal.Length ~ Sepal.Width, data = iris)")

(r-eval r "dev.off()")