Блог

Я новичок в области дендрологии и программирования R. У меня есть данные о кольцах деревьев с 1884 по 2019 год и данные о месячных осадках с 1950 по 2019 год. Я хочу рассчитать корреляцию между ростом деревьев и месячными осадками с 1950 по 2019 год. Я также хочу провести корреляцию между ростом в 1-м году ежемесячных метеорологических переменных в течение этого года (январь-октябрь), а также с метеорологическими переменными годом ранее. Проблема в том, что я не знаю, как создать код, который я могу использовать для вычисления любого из вышеперечисленных, может кто-нибудь, пожалуйста, помогите.

Это то, что мне удалось сделать на данный момент:

 rm(list = ls())  # load in dplR package library(dplR)  # AF28OFINAL # read in measurements chaf28o6 lt;- read.rwl("F:/chaf28o6.rwl") library(readr) # rwl.report(chaf28o6) # plot(chaf28o6, plot.type="spag") # detrending all the series in the chaf2806 at once chaf28o6.rwl lt;- detrend(rwl = chaf28o6, method = "Friedman")  #descriptive statistics rwl.stats(chaf28o6[1:17]) #Building a Mean Value Chronology chaf28o6.crn lt;- chron(chaf28o6.rwl, prefix = "CAM") ## library(treeclim) Abies2800 lt;- chaf28o6.crn library(readr)    clim2800 lt;- read_csv("Prec.csv", col_types = cols(Year = col_number(),   Month = col_number(), Tp = col_number())) View(clim2800)  ## correlation function analysis## when I try to run the following code, I get the error message: Error in x[1, 1]:x[dim(x)[1], 1] : NA/NaN argument dcc(chrono = Abies2800,  climate = clim2800,  selection = 1:10,  method = "cor",  dynamic = "static",  win_size = 25,  win_offset = 1,  start_last = TRUE,  timespan = NULL,  var_names = NULL,  ci = 0.05,  boot = "stationary",  sb = TRUE  )