1

I am trying to use the stargazer package to display a Ridge Regression, but the error appears: unrecognized object type. Do any of you know how to fix it? Or know other packages that show a good table with the regression results in R? I appreciate any help you can provide.

#Here is the main code used for one of the regressions:

base_sem_NAs_padronizada <- base_sem_NAs %>% 
  mutate(
    HHI = as.numeric(scale(HHI)),
    log_pib = as.numeric(scale(log(pib_milhares))),
    pib_milhares = as.numeric(scale(pib_milhares)),
    log_pib_per_capita = as.numeric(scale(log(pib_per_capita))),
    pib_per_capita = as.numeric(scale(pib_per_capita)),
    incremento = as.numeric(scale(incremento)), 
    log_pib_agricola = as.numeric(scale(log(pib_agricola))),
    pib_agricola = as.numeric(scale(pib_agricola)),
    # log_pib_industria = as.numeric(scale(log(pib_industria))),
    pib_industria = as.numeric(scale(pib_industria)),
    populacao = as.numeric(scale(populacao)), 
    ano_de_eleicao = as.numeric(scale(ano_de_eleicao)),
    nuvem = as.numeric(scale(nuvem)),
    area_km2 = as.numeric(scale(area_km2)),
    homicidio = as.numeric(scale(homicidio)),
    idhm_2010 = as.numeric(scale(idhm_2010)),
    incremento = incremento - mean(incremento) #centralizar resposta
  )

reg1 <- lm(
  incremento ~ pib_per_capita + ano_de_eleicao + nuvem  + homicidio + idhm_2010 - 1, # + log_pib_industria,
  data = base_sem_NAs_padronizada
)

beta_original1 <- reg1$coefficients #para reg1
p1 <- length(beta_original1) #para reg1
sighat2_1 <- summary(reg1)$sigma^2 #para reg1

novo_k1 <- function(beta) p*sighat2_1/(t(beta)%*%beta) #para reg1
novo_beta1 <- function(k) { #para reg1
  cur_fit1 <- ridge::linearRidge(
    incremento ~ pib_per_capita + ano_de_eleicao + nuvem  + homicidio + idhm_2010 - 1, # + log_pib_industria,
    lambda = k,
    data = base_sem_NAs_padronizada,
    scaling = "none"
  )
  return(cur_fit1$coef)
}

kk1 <- novo_k1(beta_original1) #para reg1
for(j in 1:200) kk1 <- kk1 %>% novo_beta1() %>% novo_k1() #para reg1

ridge_fit1 <- ridge::linearRidge( #para reg1
  incremento ~ pib_per_capita + ano_de_eleicao + nuvem  + homicidio + idhm_2010 - 1, # + log_pib_industria,
  lambda = kk1,
  data = base_sem_NAs_padronizada,
  scaling = "none"
)

ridge_fit1 %>% 
  summary()

stargazer(ridge_fit1)

% Error: Unrecognized object type.

Phil
  • 7,287
  • 3
  • 36
  • 66

0 Answers0