오늘은 재미있는 코드를 들고왔습니다. 이번 포스팅은 요즘 제가 보고있는 Rcpp 책을 보다가 나온 코드임을 밝힙니다. 피보나치 수열 피보나치 수열은 수학에서 유명한 수열이죠. 이 수열은 다음과 같은 점화식을 써서 나타낼 수 있습니다. \[ F_1 = 1, F_2 = 1, F_{n+2} = F_{n+1} + F_n \] 즉, 임의의 항을 만들때 그 전과 그 전전 항을 더해서 만들어 나가는 꼴입니다. 이 수열은 컴공 알고리즘 시간에도 많이 나오는데, 다음과 같은 재귀 함수를 써써 짜는 것을 배웁니다. fibonacci
이번 포스팅에서 다룰 R함수는 바로 apply() 함수이다. 사실 R에는 여러가지 apply()로 끝나는 이름의 함수들이 많다. 예를들어, lappy(), sapply(), tapply(), 그리고 mapply() 등등 사실 이런 함수들의 사용법을 알고 있어도 필자도 매번 잊어버리고 도움말을 찾아본다. 일일이 이러한 함수의 사용법을 알고 있는 것 자체도 비효율적이라고 생각하는 측면도 있고, 따지고 보면 오늘 설명할 apply 함수에서 출발한 아류들이기 때문이다. 따라서 apply 함수를 확실하게 이해하고, 자신의 코드 상황에 맞춰 다른 함수들을 응용할 수 있게 된다면, 자신의 코드가 훨씬 간결해지는 것을 느낄 수 있을 것이다. 함수의 이름이 apply인 이유는 기존에 정의된 함수를 어떤 대상에 편하게 적..
Make your own R package! - 10. 설명서(Vignette)와 마크다운(markdown)에 대하여 Vignette이란 단어를 R 패키지를 만들며 처음 접했을 때, 생소한 단어라서 당황했던 기억이 있다. Vignette의 뜻은 책의 타이틀 페이지나 요약본이라는 뜻인데, 프랑스에서 전해진 단어라서 그런지 우리가 생각하는 비그넷이 아닌 “빈옛(트)”으로 g를 묵음처리해서 발음한다. 다음의 링크를 타고 한번 들어보도록 하자.(http://aha-dic.com/View.asp?Word=vignette) 오늘은 패키지의 설명서라고 할 수 있는 Vignette을 만드는 방법에 대하여 알아보도록 하자. 이제까지 포스팅을 하면서 사용한 함수들을 필자의 패키지에 넣어놓았는데, 이제 필자의 패키지에도 들..
Make your own R package! - 9. C언어와 연결된 함수 R 패키지에 넣기 이제까지 우리는 R와 C를 연결하는 방법에 대하여 알아보았다. 오늘은 이렇게 작성한 함수를 우리가 만든 패키지에 넣는 방법에 대하여 알아보도록 하자. 사용할 파일 준비하기 R 파일 준비하기 우리가 패키지에 넣을 함수는 파레토 분포의 확률밀도함수를 계산하는 dpareto 함수이다. 우리가 정의할 dpareto 함수는 R 사용자로부터 입력값을 넘겨 받아 .C 함수를 사용하여 전달받은 입력값들을 C로 짜여진 함수로 넘겨주게 된다. 필자의 경우 paretodens 라는 C함수로 넘겨줄 예정이다. 이렇게 R과 C사이의 교량 역할을 하는 함수를 보통 wrapper 함수라고 하는데, 이유는 뜯어보면 C에서 작동하는 함수이지만..
이전 포스팅에서 우리는 R함수와 C언어 짜여진 함수를 연결하는 방법에 대하여 알아보았다. 이전 포스팅에서는 R 함수를 이용하여 입력값의 길이가 다른 경우 lengthout 옵션을 이용하여 같은 길이의 벡터로 바꿔주었는데, 이 부분을 이번에는 C에서 구현해보도록 하자. C 함수에서의 recyling 구현 C 함수에서 R에서의 recycling을 구현하기 위해서는 % 연산자를 잘 사용하면 된다. 필자가 정의한 함수의 코드를 살펴보도록 하자. 아래의 코드에서 입력값들은 x, alpha, x_m 벡터와 각각의 길이를 담고있는 input_len 벡터, 그리고 가장 긴 벡터의 길이 max_len, log를 처리하기 위한 indicator, 결과를 전달받을 density 벡터를 의미한다. 가장 핵심이 되는 부분은 f..
Make your own R package! - 7. R과 C 언어 연결하기 오늘의 주제는 R와 C 언어를 연결하기 위한 준비단계로 이전 포스팅에서 짰던 dpareto함수를 C언어로 짜는 것을 알아보도록 하겠다. R을 주로 포스팅하는 필자가 C언어에 대하여 포스팅하는 이유는 C언어를 다룰수 있으면 R이 제일 취약한 부분인 루프(loop) 부분을 C로 넘겨서 계산을 할 수 있도록 연결해줄 수 있기 때문이다. C 언어로 함수 정의하기 이전 포스팅에서 다루었던 파레토 분포의 확률밀도함수 dpareto를 C언어를 사용하여 만들어 보도록 하자. 먼저 가장 기본이 되는 변수 하나씩 받을 수 있는 함수를 paretodens라는 이름으로 다음과 같이 정의할 수 있다. C의 경우 함수의 앞에 함수의 결과값이 어떠한 형식..
Make your own R package! - 6. 함수 벡터화(Vectorized function)에 대하여 오늘의 주제는 함수를 정의하는 방법에 대한 내용이다. 패키지에 함수를 정의해서 사용하는 방법은 사실 이전 포스팅에서 많이 다루었다. 하지만 만약 훗날 다른 사람들을 위하여 함수를 정의하고 배포할 경우 고려해야하는 사항들을 한번 알아보도록 하자. 이 포스팅은 필자가 수강한 학교 수업의 숙제에서 따온 내용이다. Pareto분포의 확률밀도함수 dpareto 정의하기 R은 여러 통계 분포에 관한 함수를 제공하고 있다. 예를 들어 정규분포의 경우, 확률밀도함수는 dnorm, 누적분포함수는 pnorm 등 앞의 분포의 이름 앞에 제공하는 함수의 성질을 간략한 약자를 사용하여 나타내고있다. 오늘 우리가 만들..
r4issactoast을 이용한 XY-plot 쉽게 그리기 r4issactoast 패키지 XY-plot 기능 추가 필자가 통계를 주로 설명하는 블로그를 운영하고 있다보니, 그래프를 그릴 일이 많다. 아무래도 시각적인 것이 글로 써져있는 것보다 훨씬 잘 눈에 들어오기 때문일 것이다. 그 중에서도 가장 많이 사용되는 것을 꼽으라면 단연 XY plot이다. 매번 XYplot을 그리는 입장에서 코드 쓰는게 여간 귀찮은게 아니었다. 따라서 좀 더 빨리 XY plot을 그리고 싶은 마음이 있었는데, 마침 예전에 필자가 기본 패키지를 사용하여 예쁜 XY 그래프를 그리는 것을 연재한 것이 생각나, 이번에 시간을 내어 그 때 사용했던 코드를 r4issactoast에 넣어놓았다. r4issactoast 패키지를 통한 XY..
Tensorflow + R - 2. placeholder와 노드(node) 개념 이해 본 텐서플로우 시리즈는 Sung Kim 교수님의 유튜브 강의(python 기반, 참고자료의 링크)를 기반으로 필자가 R 버전에 맞도록 코드를 수정 및 변형 시켰음을 알려드립니다. library(tensorflow) Hello, Tensorflow 저번 포스팅에서 우리는 Tensorflow상에 상수를 선언하는 방법에 대하여 알아보았다. 오늘은 텐서플로우에서 그래프의 개념과 Placeholder의 개념에 대하여 알아보도록 하겠다. 먼저, 아래의 코드와 같이 문자열을 텐서플로우 상에 정의해보도록 하자. # Say hello to your tensorflow hello
ggplot2 패키지 정리 - 1. stat_density, stat_function, geom_rug, geom_point, anotate, labs 이번 포스팅에서 사용될 R 패키지는 다음과 같다. library(RColorBrewer) library(magrittr) library(ggplot2) 이번 포스팅에서는 필자의 개인 참고를 위하여 최근 필자가 사용했던 ggplot2의 코드들을 정리해놓으려 한다. RColorBrewer 최근 필자가 RColorBrewer의 위력을 실감하였다. 패키지가 하는 일은 아주 간단하다. 색깔을 선택해준다. 사용법 # pick 5 colors with theme: "Accent" display.brewer.pal(5, "Accent") my_color
Make your own R package! - 5. Imports, Suggests, and README.rmd 저번 포스팅을 통해서 필자는 roxygen2 패키지를 사용하여 자동으로 man 폴더와 NAMESPACE 파일을 관리하는 방법에 대하여 알아보았다. 오늘 배울 내용은 우리의 패키지를 만드는 과정에서 다른 패키지들에 존재하는 함수를 빌려서 사용하는 방법을 알아보자. mle_gamma() 함수 현재 필자의 r4issactoast 패키지에는 mom_gamma() 함수가 들어있다. 이번 포스팅에서는 mle_gamma() 함수를 예로 설명해보도록 하겠다. 왜냐하면 mle_gamma() 함수는 아래의 코드와 같이 R의 stats 패키지의 uniroot() 함수를 사용해서 정의되었기 때문이다. 함수에 관한..
Make your own R package! - 4. 나만의 함수 추가하기 with roxygen2 이번 포스팅에서 사용될 R 패키지는 다음과 같다. library(devtools) library(roxygen2) 저번 포스팅을 통해서 필자는 DESCRIPTION 파일 설정에 대하여 알아보았다. 이번 시간에는 함수를 추가하는 방법에 대하여 알아보도록 하자. 추가할 함수 마련하기 필자가 r4issactoast패키지를 만들기로 결심한 것은 최근의 통계 포스팅 시리즈를 연재하면서 부터이다. 따라서 필자는 예전에 만들어놓은 감마분포의 MOME 함수를 패키지에 넣도록 하겠다. 함수는 다음과 같다.(함수와 관련된 내용은 이전 포스팅을 참고하자.) mom_gamma
