Я анализирую данные промыслового ряда для траулеров.Я использовал пакет strucchange для определения контрольных точек в атрибуте (увеличение мощности траулеров) с 1993 по 2017 год. Это ежегодная статистика, и я хочу найти переломные моменты в течение этого периода, соответствующие политикам управления рыболовством.Функция breakpoints () обнаружила три переломных момента.Но я столкнулся с проблемой при оценке доверительных интервалов для точек останова.Первый пункт - 1994, и КИ - это НС.Поэтому, когда я попытался использовать lines () для добавления CI в график, произошел сбой с отчетом об ошибке, который также отображался при печати объекта confint в r:
Error in if (any(obj$confint < 1) | any(obj$confint > obj$nobs)) warning(paste("Confidence intervals outside data time interval\n\t from ", :
missing value where TRUE/FALSE needed
Я пытался опустить NA изобъект confint (ci_rht_ts), но нашел ключ к его созданию.
Ниже приведены мои r коды:
идентифицируют точки останова в моих данных (r_ht.ts)
r_ht.ts<-na.omit(r_ht.ts) ## r_ht.ts is the original data with two columns, one is 'Date', the other is 'R_h_tr' (the increase rate of horsepower of trawlers). There are some NAs in some years of my original data (originally 1950 - 2017, no data from 1950 - 1993).
r_ht.ts$Date<-as.Date(r_ht.ts$Date,format="%d/%m/%Y")
r_ht.ts$Date<-as.POSIXct(r_ht.ts$Date)
str(r_ht.ts)
r_ht.ts<-ts(r_ht.ts$R_h_tr,start = c(1993,1),end = c(2017,1),frequency=1)
rht_ts <- breakpoints(r_ht.ts ~ 1, h = 0.1)# get the breakpoints
ci_rht_ts<-confint(rht_ts)
ci_rht_ts # the error pops up as mentioned above
plot(r_ht.ts,xlab="Year",ylab="Increase rate of horsepower of trawlers")
lines(rht_ts)
lines(ci_rht_ts)
ошибка появляется, когда я запускаю ci_rht_ts и строки (ci_rht_ts):
Доверительные интервалы для точек останова оптимального 4-сегментного раздела:
Вызов:
confint.breakpointsfull(object = rht_ts)
Точки останова на номере наблюдения:
2.5 % breakpoints 97.5 %
1 NA 2 NA
2 3 4 7
3 9 11 15
Соответствует перерывам:
Error in if (any(obj$confint < 1) | any(obj$confint > obj$nobs)) warning(paste("Confidence intervals outside data time interval\n\t from ", :
missing value where TRUE/FALSE needed
Я ожидаю игнорировать NA в объекте ограничения для первых точек останова, но построю CI для остальных двух.