Symbolic Regression in R

최근에 검색을 하다가 기호적회귀분석이란 것을 알게 되었다.

 

이를 적용함에 있어 R패키지 보다 python 패키지가 좀 더 유용해 보였다.

 

따라서 python을 R에서 쓰는 식으로 예제코드를 생성하였다. 

 

세팅하는 부분은 아래 블로그를 참고하기 바란다.

 

https://rcj92.tistory.com/10

R tensorflow-gpu 설정하기

 

https://pycj92.tistory.com/23

기호적회귀분석(Symbolic Regression)

library(forecast)
library(reticulate)
data=decompose(AirPassengers)
data=na.omit(data.frame(x=data$x,seasonal=data$seasonal,trand=data$trend,random=data$random))
train=data[1:80,]
test=data[81:nrow(data),]
gpl=import('gplearn')
function_set=c('add','sub','mul','div')
model=gpl$genetic$SymbolicRegressor(population_size = 3000L, 
                              tournament_size = 5L,
                              generations = 25L,
                              stopping_criteria = 0.1,
                              function_set = function_set,
                              metric = 'rmse',
                              p_crossover = 0.65,
                              p_subtree_mutation = 0.15,
                              p_hoist_mutation = 0.05, 
                              p_point_mutation = 0.1,
                              verbose = 1,
                              random_state = 42L,
                              n_jobs = -1L)
model$fit(train,as.vector(apply(train,1,sum)))
plot(apply(data,1,sum))
points(c(apply(train,1,sum),model$predict(test)),col=c(rep(1,80),rep(2,52)))