[통계, R] 연대기 만들고, 그래프 그리기

R을 활용한 연륜연대 분석

#1. 기본 패키지 불러오기 및 DATA 생성하기

##패키지 불러오기 ##
library(dplR)
library(utils)
library(graphics)
library(tcltk)
library(ggplot2)
library(usethis)
library(devtools)
library(dendroTools)
library(FactoMineR)
library(factoextra)
library(corrplot)
library(detrendeR)

detrender()

##작업공간 설정 ##

setwd("/R/")

## 데이터 불러오기 ##
#FH 파일 불러오기#
sapd.fh <- read.fh("sapd_raw_dated.fh")

#RWL 파일 불러오기#

sapd.rwl <- read.rwl("sapd_raw_dated.rwl")

## 데이터 불러오기 ##

sapd.rwl <- read.rwl('sapd.rwl', format = c("auto")) 

SCPDf 

## 옵션 값

format = c("auto", "tucson", "compact", "tridas", "heidelberg", "csv")

## 2번째 분석코어부터 24번째 분석코어까지만 추출 ##

sapd.s2 <- sapd.s[,c(2:24)] 

sapd.s2

##기본 통계 ##
rwl.report(sapd.fh)

rwl.stats(SCPDf)

## 기본 그래프 보기  ##

plot(sapd.fh, plot.type="seg")

plot(sapd.fh, plot.type="spag")

## 연대기 만들기, prefix=3글자 미만,  ##
sapd.fh.crn <- chron(sapd.fh, prefix = "설악산", prewhiten=TRUE)
sapd.fh.crn

## 색상 코드 ##
my.cols <- c("#3182BD","#9ECAE1","#DEEBF7","#31A354","#A1D99B","#E5F5E0")

## 연대기 그래프 ##
crn.plot(sapd.fh.crn,add.spline=TRUE,nyrs=32,
             crn.line.col=my.cols[5],
             spline.line.col=my.cols[4],
             samp.depth.col=my.cols[3],
             samp.depth.border.col=my.cols[2],
             abline.col=my.cols[1],
             crn.lwd=1.5,spline.lwd=3,
             abline.lwd=1,
             xlab="Year", ylab="RWI")




참고사이트

• https://rdrr.io/cran/dplR/