獨斷論

Windows subsystems linux에 파일을 윈도우에서 복사하려고 할때 파일탐색기를 연다. 주소창에 역슬래시를 두번하고 wsl$를 치고 역슬래쉬하고 Ubuntu라고 입력한다. Ubuntu는 자신이 설치한 리눅스가 Ubuntu일 경우에 해당된다. \\wsl$\Ubuntu wsl의 하위디렉토리중 home을 선택하면 자신의 사용자명이 나올것이다

책 홈페이지: https://xcelab.net/rm/statistical-rethinking/ 2M1 풀이 # Define grid p_grid = seq(from=0, to=1, length.out=20) # Prior prior = rep(1, 20) # Likelihood # 2M1 (1) W = 3 N = 3 ## 2M1 (2) #W = 3 #N = 4 ## 2M1 (3) #W = 5 #N = 7 likelihood = dbinom(x=W, size=N, prob=p_grid) # Posterior posterior = likelihood * prior posterior = posterior / sum(posterior) # Plot plot(p_grid, posterior, type="b", ..

앞서 1과 2에서는 JAGS모델을 별도의 파일에 저장하여 R script에서 불러들였다. 이는 모델 전체를 이해하는데 도움이 되지 않는다. 여기서는 별도의 *.jags 파일을 작성하지 않고 R script에 직접 입력하는 방법을 알아보자. 입문1을 아래와같이 R script만으로 JAGS를 이용할수 있다. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 18 19 20 21 22 23 24 25 setwd("d:/tmp/rcode/") rm(list = ls()) x

간단한 회귀분석(linear regression)을 베이시언(Bayesian)으로 수행하여 보자 아래 데이터 파일을 R의 working directory에 내려받은 후에 R에서 읽는다. R에서 데이터를 아래와 같이 읽어들인다 1 2 3 4 5 rm(list = ls()) setwd("d:/tmp/rcode/") regdat

패키지 설치JAGS를 먼저 컴퓨터에 설치한다. 현재 버전은 4.3.0 http://mcmc-jags.sourceforge.net/ R에서 JAGS를 사용하기 위하여 rjags를 R에 설치한다. R console에서 아래 스크립트를 실행 > install.packages("rjags") 데이터 내려받기평균과 표준편차를 알지 못하는 가상의 데이터 example1.csv를 내려받아 R의 working directory에 저장한다. JAGS 모델 작성 이제 평균과 표준편차를 모르는 표준정규분포를 따르는 example1.csv 데이터의 평균과 표준편차를 베이시언Bayesian 추론하여보자. JAGS 모델은 다음과 같다. 12345678 model { for (i in 1:N) { x[i] ~ dnorm(mu, t..

Python qtconsole은 기본적으로 하얀 바탕에 검은 글씨로 되어 있어서 오래 코딩하다보면 눈이 아픈데 이걸 좀 바꿔보자.Anaconda prompt에서 jupyter qtconsole이라 입력하자(예전에는 ipython qtconsole이었는데 이제 바뀌었다).(base) C:\Users\your_name>jupyter qtconsole위와같은 희멀건한 창이 뜨고 글꼴도 매우 작다.이를 한번 원하는데로 바꿔보자. C:\Users\your_name\.jupyter\안에 jupyter_qtconsole_config.py가 있는지 보고 없으면 아래와 같이 anaconda prompt에서 실행한다.(base) C:\Users\your_name>jupyter qtconsole --generate-conf..

R은 버전이 바뀔때마다 라이브러리 디렉토리(폴더)도 바뀌어서 패키지 설치한걸 다시 또 설치해야하는 번거로움이 있는데 라이브러리 디렉토리를 영구적으로 지정해두면 R 버전이 업그레이드 되더라도 패키지를 다시 설치할 필요가 없다. R library 디렉토리를 영구적으로 지정하는 방법을 알아보자. notepad와 같은 편집기를 열고 아래와 같이 입력한다. 여기서 "yourname"은 자신의 윈도우 로그인 계정이름이다. .libPaths(c("C:/Users/yourname/Documents/R/Library", .libPaths())) 이 파일의 이름을 .Rprofile로 저장하는데, 주의해야할 것은 파일명이 점(.)으로 시작하고 확장자가 없다. 또한 디렉토리 표시가 역슬레쉬가 아니라 슬레쉬이다. 저장할 디렉토..

R 4.0이상은 더이상 지원하지 않는 것 같다 (2021. 2. 7) 통계소프트웨어 R은 주로 RStudio를 통합개발환경(IDE)로 사용하지만 비주얼스튜디오(visual studio)도 사용할수가 있다. 우선 Microsoft R open을 설처히야 한다. cran에서 제공하는 R과 달리 선형대수에 사용하는 수학관련된 라이브러리가 5~10배정도 빠르므로 GNU R을 사용하였을때 실행속도가 느리다면 Microsoft R Open을 사용해보는것도 좋다. Microsoft R Open 설치하기 왼쪽 링크 (https://mran.microsoft.com/download)에서 자신의 OS에 맞는걸 선택하여 클릭한다. 32bit는 제공하지 않는다. 현재 제공하는 OS는 윈도우, 우분투, 레드햇(페도라), 수세리..