Извлечь коэффициенты и t-значения из квантильной регрессии, rqpd

Вы можете проверить, работает ли (безымянный) пакет, который вы используете, с stargazer, но в противном случае получение звездочек не будет автоматическим. Обычно эту информацию можно извлечь с помощью summary().

# coefficients and t-values in matrix
cf_10 <- summary(x)$coefficients[1:8,c(1,3)]
# p-values
cf_10.pVals <- summary(x)$coefficients[1:8,4]

Если вы действительно хотите, вы можете реализовать алгоритм для вставки звездочек на коэффициенты.

Вероятно, есть способ получше, но его будет довольно просто прочитать:

# function to print stars
starPrinter <- function(pVal) {
  if(pVal < 0.01) return("***")
  if(pVal < 0.05) return("**")
  if(pVal < 0.1) return("*")

  return("")
}

# a matrix, with digits rounded to 3rd decimal
myTable <- round(cf_10, 3)
# get stars and put tvalues in parentheses
for(i in 1:nrow(myTable)) {
  myTable[i, 1] <- paste0(myTable[i, 1], starPrinter(cf_10.pVals[i]))
  myTable[i, 2] <- paste0("(", myTable[i, 2], ")")
}

Если вы хотите сделать это вектором, где у вас есть coef1, tval1, coef2, tval2 ...

myEstVector.cf_10 <- as.character(t(myTable))

Затем вы можете объединить несколько результатов в новую таблицу:

myNewTable <- cbind(myEstVector.cf_10, myEstVector.cf_25, myEstVector.cf_99)

Просто пытаюсь добавить в вышеупомянутую ветку. Я выполнил приведенный выше код и экспортировал результаты в формат файла csv.

# y is the dependent variable 
# x6 is a vector of 6 independent variables. 
#setting regression
r6.form <- y ~ x6 |as.factor(cty)
r6090 <- rqpd(r6.form, panel= panel(taus=c(0.1, 0.25, 0.5, 0.75, 0.9), tauw = rep(1/5,5)), 
              na.action = 'na.omit', data=mydata)

#Storing results 
c6090 <- summary(r6090, se = "boot", R=1000, covariance = TRUE)

#making a data frame for results
df6 = as.data.frame(c6090$coefficients)

#coefficients and t-values in a matrix
cf_6 <- df6[1:71, c(1,3)]

#p-values
cf_6.pval <- df6[1:71,4]

#function to print stars
starPrinter <- function(pVal) {
  if(pVal < 0.01) return("***")
  if(pVal < 0.05) return("**")
  if(pVal < 0.1) return("*")

  return("")
}

# a matrix, with digits rounded to 3rd decimal
myTable6 <- round(cf_6, 3)

# get stars and put tvalues in parentheses
for(i in 1:nrow(myTable6)) {
  myTable6[i, 1] <- paste0(myTable6[i, 1], starPrinter(cf_6.pval[i]))
  myTable6[i, 2] <- paste0("(", myTable6[i, 2], ")")
}

#creating a vector
table.cf_6 <- (t(myTable6))

#exporting it to CSV file
write.table(table.cf_6, "c6090.csv", sep =',', row.names = TRUE)

Это извлекает результаты в разные столбцы с именами переменных.