R: (RHRV) Ошибка в plot.window (…): нужны конечные значения xlim - PullRequest
0 голосов
/ 21 марта 2019

Я использую пакет RHRV в R, и у меня возникают проблемы с анализом временных интервалов короче 400 с, особенность в том, что это происходит только в определенных входных файлах, а не в других.

ДляНапример, я хочу проанализировать данные на 300 с, поэтому временное окно будет длиться 300 с, я могу сделать это для файла 1.txt, но я не могу сделать это для файла 2.txt, но если я увеличу временное окнодо 400 с (где это комментируется в фрагменте кода) это будет работать.Файлы данных имеют одинаковый формат.

Это файл сценария:

library(RHRV)
library(ggplot2)
library(expss)
library(knitr)
library(gridExtra)
library(tidyverse)

hrv.data <- CreateHRVData()
hrv.data <- SetVerbose(hrv.data, TRUE) 
hrv.data <- LoadBeatAscii(hrv.data,"C:/Users/User/Documents/2.txt")

hrv.data = BuildNIHR(hrv.data)
hrv.data = FilterNIHR(hrv.data)
hrv.data = InterpolateNIHR (hrv.data, freqhr = 4)
hrv.data_time <- ExtractTimeSegment(hrv.data, 0, 400) #It doesn't work for examples under 400 sec

spectral_analysis <- CreateFreqAnalysis(hrv.data_time)
spectral_analysis <- CalculatePowerBand(spectral_analysis , indexFreqAnalysis= 1,
                                        size = 300, shift = 30, type = "fourier",
                                        ULFmin = 0, ULFmax = 0.03, VLFmin = 0.03, VLFmax = 0.05,
                                        LFmin = 0.05, LFmax = 0.15, HFmin = 0.15, HFmax = 0.4 )

timedomain_analysis <- CreateTimeAnalysis(hrv.data_time, size = 300)

nonlinear_analysis <- CreateNonLinearAnalysis(hrv.data_time)
nonlinear_analysis <- CalculateCorrDim(nonlinear_analysis,
                                       indexNonLinearAnalysis = 1, 
                                       minEmbeddingDim = 2,
                                       maxEmbeddingDim = 8, timeLag = 1,
                                       minRadius = 1, maxRadius = 15,
                                       pointsRadius = 20, theilerWindow = 10,
                                       corrOrder = 2, doPlot = FALSE)

nonlinear_analysis <- CalculateSampleEntropy(nonlinear_analysis)
nonlinear_analysis <- EstimateSampleEntropy(nonlinear_analysis)
nonlinear_analysis <- PoincarePlot(nonlinear_analysis)


SDNN <- timedomain_analysis$TimeAnalysis[[1]]$SDNN
SDANN <-timedomain_analysis$TimeAnalysis[[1]]$SDANN
SDNNIDX <-timedomain_analysis$TimeAnalysis[[1]]$SDNNIDX
SDSD <- timedomain_analysis$TimeAnalysis[[1]]$SDSD
pNN50 <- timedomain_analysis$TimeAnalysis[[1]]$pNN50 
rMSSD <- timedomain_analysis$TimeAnalysis[[1]]$rMSSD 
IRRR <-timedomain_analysis$TimeAnalysis[[1]]$IRRR
MADRR <- timedomain_analysis$TimeAnalysis[[1]]$MADRR
TINN <- timedomain_analysis$TimeAnalysis[[1]]$TINN
HRVi <- timedomain_analysis$TimeAnalysis[[1]]$HRVi

HRV <- mean(spectral_analysis$FreqAnalysis[[1]]$HRV)
LF <- mean(spectral_analysis$FreqAnalysis[[1]]$LF)
HF <- mean(spectral_analysis$FreqAnalysis[[1]]$HF)
LFHF <- mean(spectral_analysis$FreqAnalysis[[1]]$LFHF)

entropy <- mean(nonlinear_analysis$NonLinearAnalysis[[1]]$sampleEntropy$statistic) 
SD1 <- nonlinear_analysis$NonLinearAnalysis[[1]]$PoincarePlot$SD1
SD2 <- nonlinear_analysis$NonLinearAnalysis[[1]]$PoincarePlot$SD2

hrv_results_table <- c(SDNN, SDANN, SDNNIDX, SDSD,pNN50, rMSSD, HRVi, MADRR, TINN, HRV, LF, HF, LFHF, entropy, SD1, SD2)

hrv_results_table_DF <- data.frame(hrv_results_table)
hrv_results_table_DF <- hrv_results_table_DF %>% `row.names<-`(c("SDNN","SDANN", "SDNNIDX", "SDSD", "pNN50","r-MSSD","HRVi", "MADRR", "TINN", "HRV", "LF", "HF", "LF/HF", "Sample entropy", "Poincare SD1", "Poincare SD2"))
hrv_results_table_DF <- round(hrv_results_table_DF,2)
name.width <- max(sapply(names(hrv_results_table_DF), nchar))
print(hrv_results_table_DF)

Это сообщение об ошибке:

nonlinear_analysis <- CalculateSampleEntropy(nonlinear_analysis) Computing the sample entropy of order 2 no non-missing arguments to min; returning Infno non-missing arguments to max; returning -Infno non-missing arguments to min; returning Infno non-missing arguments to max; returning -InfError in plot.window(...) : need finite 'xlim' values

Это входной файл счто у меня возникли проблемы: 1.txt 2.txt

Буду признателен за любую помощь, Петр

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...