Как получить имя встроенной функции в R

Я создал код, который подбирает разные типы копул к данным и выбирает лучший на основе наименьшего значения AIC. Код работает очень хорошо. Однако результат следующий:

Normal copula, dim. d = 3 
Dimension:  3 
Parameters:
  rho.1   = 0.8094
dispstr:  ex 

Затем я могу получить из результата размер и параметры, но не имя связки. То есть мне нужно извлечь имя копуала. Например,

    selected_copula@parameters
[1] 0.7833736

Я также хочу иметь возможность использовать selected_copula@copula, чтобы получить имя наиболее подходящей связки, но это не работает.

Функция fitCopula позволяет нам извлечь имя связки с помощью @copula, чего я не могу сделать с помощью своего кода.

Полный код выглядит следующим образом:

library(copula)
fit_select_copula <- function(copula_list, data) {
  fit_copula_models <- function(copula_list, data) {
    num_copulas <- length(copula_list)
    copula_results <- vector("list", length = num_copulas)

    for (i in seq_along(copula_list)) {
      copula <- copula_list[[i]]
      fitted_model <- fitCopula(copula,data)

      copula_results[[i]] <- list(
        copula = fitted_model@copula,
        loglik = fitted_model@loglik,
        aic = AIC(fitted_model)
      )
    }

    return(copula_results)
  }

  copula_results <- fit_copula_models(copula_list, data)

  aic_values <- sapply(copula_results, function(x) x$aic)
  lowest_aic_index <- which.min(aic_values)

  selected_copula <- copula_results[[lowest_aic_index]]$copula

  return(selected_copula)
}

# Example usage
data <- rCopula(1000, normalCopula(0.8, dim = 3))  # Replace with your actual data
copula_list <- list(
  frank = frankCopula(6, 3),
  clayton = claytonCopula(2, 3),
  normal = normalCopula(0.7, 3)
)

selected_copula <- fit_select_copula(copula_list, data)

поскольку у copula_list есть имена, просто сделайте names(copula_results) <- names(copula_list)

Onyambu 01.07.2024 21:44
Стоит ли изучать PHP в 2023-2024 годах?
Стоит ли изучать PHP в 2023-2024 годах?
Привет всем, сегодня я хочу высказать свои соображения по поводу вопроса, который я уже много раз получал в своем сообществе: "Стоит ли изучать PHP в...
Поведение ключевого слова "this" в стрелочной функции в сравнении с нормальной функцией
Поведение ключевого слова "this" в стрелочной функции в сравнении с нормальной функцией
В JavaScript одним из самых запутанных понятий является поведение ключевого слова "this" в стрелочной и обычной функциях.
Приемы CSS-макетирования - floats и Flexbox
Приемы CSS-макетирования - floats и Flexbox
Здравствуйте, друзья-студенты! Готовы совершенствовать свои навыки веб-дизайна? Сегодня в нашем путешествии мы рассмотрим приемы CSS-верстки - в...
Тестирование функциональных ngrx-эффектов в Angular 16 с помощью Jest
В системе управления состояниями ngrx, совместимой с Angular 16, появились функциональные эффекты. Это здорово и делает код определенно легче для...
Концепция локализации и ее применение в приложениях React ⚡️
Концепция локализации и ее применение в приложениях React ⚡️
Локализация - это процесс адаптации приложения к различным языкам и культурным требованиям. Это позволяет пользователям получить опыт, соответствующий...
Пользовательский скаляр GraphQL
Пользовательский скаляр GraphQL
Листовые узлы системы типов GraphQL называются скалярами. Достигнув скалярного типа, невозможно спуститься дальше по иерархии типов. Скалярный тип...
0
1
51
1
Перейти к ответу Данный вопрос помечен как решенный

Ответы 1

Ответ принят как подходящий

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

library(copula)
fit_select_copula <- function(copula_list, data) {
  fit_copula_models <- function(copula_list, data) {
    num_copulas <- length(copula_list)
    copula_names <- names(copula_list) #  <<< ADDED THIS  <<<<
    copula_results <- vector("list", length = num_copulas)
    
    for (i in seq_along(copula_list)) {
      copula <- copula_list[[i]]
      fitted_model <- fitCopula(copula,data)
      
      copula_results[[i]] <- list(
        copula_name = copula_names[[i]], #  <<< ADDED THIS  <<<<
        copula = fitted_model@copula,
        loglik = fitted_model@loglik,
        aic = AIC(fitted_model)
      )
    }
    
    return(copula_results)
  }
  
  copula_results <- fit_copula_models(copula_list, data)
  aic_values <- sapply(copula_results, function(x) x$aic)
  lowest_aic_index <- which.min(aic_values)
  
 ### ALTERED THIS ####
  selected_copula <- list(name=copula_results[[lowest_aic_index]]$copula_name,
                          content = copula_results[[lowest_aic_index]]$copula)
  
  return(selected_copula)
}

# Example usage
data <- rCopula(1000, normalCopula(0.8, dim = 3))  # Replace with your actual data
copula_list <- list(
  frank = frankCopula(6, 3),
  clayton = claytonCopula(2, 3),
  normal = normalCopula(0.7, 3)
)

selected_copula <- fit_select_copula(copula_list, data)

Другие вопросы по теме