как использовать пакет верификации R для набора данных сезонного прогнозирования климата

Question

Я хочу использовать пакет verification в R. У меня есть ежедневные данные obs и прогноза в следующем формате:

• Данные Obs имеют размерность 120, 396, где 120 - лидерствовремя, а 396 - общее начальное время.
• Данные прогноза имеют измерение 120, 396, 10, где первые два измерения совпадают с измерением данных obs, а третье измерение является членом ансамбля.

Кроме того, я уже удалил суточное климатологическое значение из моих данных.Чтобы получить вероятностный навык, мне нужно преобразовать свои данные obs в набор двоичных данных, а данные прогнозирования - в вероятности.Для этого мне нужно использовать функцию проверки в пакете verification в R (я думаю, что мое понимание верно).Но когда я применяю эту функцию, она выдает некоторую ошибку как

If baseline is not included, baseline values  will be calculated from the sample obs.

Error in e1[[1L]] : subscript out of bounds

Дополнительно:

Предупреждающие сообщения:

1. In is.finite(obs) & is.finite(pred) : longer object length is not a multiple of shorter object length
2. In max(pred) : no non-missing arguments to max; returning -Inf
3. In min(pred) : no non-missing arguments to min; returning Inf

Я не знаю, как решить эту функцию.Кто-нибудь может мне помочь в этом отношении?

Answer 1

Я решил проблему чтения данных и получил несколько графиков из пакета подтверждения R.Теперь я хочу знать, как упорядочить данные в R, то есть я упорядочил данные наблюдений, например (120, 396), где 120 - это время выполнения, а 396 - общее начальное время.Данные прогноза имеют одно дополнительное измерение, например, член ансамбля (120, 396,10).Прежде чем поместить эти наборы данных в функцию проверки пакета верификации, я делаю так:

obs <- round(runif(100))
forecast <- runif(100)

Я не знаю, почему этот процесс необходим и что он делает, но это необходимо, прежде чем вводить в функцию проверки.Затем, чтобы использовать функцию проверки, которую я использую,

verify(obs, forecast, obs.type = "binary", frcst.type="prob")

Тогда все в порядке, и я могу использовать результат функции проверки для других вещей.Но я не знаю, как работает пакет верификации, и правильное ли расположение данных?Какое измерение берется в качестве измерения образца в приведенном выше наборе данных 120, 396?

Answer 2

Вы пробовали работать через виньетка ?