Альтернативы для вывода регрессии с использованием miEconAids (кроме `stargazer`)

Question

Я ищу более симпатичный вывод из оценок с использованием библиотеки micEconAids.stargazer, кажется, не работает, если я не оцениваю уравнения "вручную", используя lm или другой пакет регрессии для, казалось бы, не связанных уравнений.Некоторые из методов, реализованных в пакете, я мог бы воспроизвести самостоятельно, но не все из них.

Итак, есть ли альтернативы, которые могли бы дать мне лучший результат, кроме выхода регрессии из summary?

Ниже я вставил код, чтобы проиллюстрировать проблему:

rm(list=ls())

library(micEconAids)

data(Blanciforti86)

priceNames = c( "pFood1", "pFood2", "pFood3", "pFood4" )
shareNames = c( "wFood1", "wFood2", "wFood3", "wFood4" )

laaidsResult1 = aidsEst( priceNames, shareNames, "xFood", data = Blanciforti86,
                         priceIndex = "S"  )

laaidsResult2 = aidsEst( priceNames, shareNames, "xFood", data = Blanciforti86,
                        priceIndex = "P"  )

laaidsResult3 = aidsEst( priceNames, shareNames, "xFood", data = Blanciforti86,
                        priceIndex = "SL"  )

summary(laaidsResult1)
summary(laaidsResult2)
summary(laaidsResult3)

#Stargazer does not work with `micEconAids`:

stargazer::stargazer(laaidsResult1, laaidsResult2, laaidsResult3)

#What are the alternatives for showing a prettier output?

#######################################################################################################################

#In the previous examples, an AIDS model was regressed.
#Below there is an example of how the models could/should look like -- in a three equation setup --
#using `lm` function, which works with `stargazer`.


reg1 = lm(Blanciforti86$wFood1 ~ Blanciforti86$pFood1)
reg2 = lm(Blanciforti86$wFood2 ~ Blanciforti86$pFood2)
reg3 = lm(Blanciforti86$wFood3 ~ Blanciforti86$pFood3)

stargazer::stargazer(reg1,reg2,reg3, type="text")

Answer 1

Вывод aidsEst содержит массу информации в структуре списка.Если есть определенный выходной набор, который вы хотите проверить, вы можете извлечь его и использовать broom::tidy() для просмотра в более хорошем формате.

Например, aidsEst() содержит $coef объект, который, в свою очередь, имеет $stat фрейм данных.
Использование tidy():

tidy(laaidsResult1$coef$stat)

   .rownames     Estimate  Std..Error    t.value     Pr...t..
1    alpha 1 -0.247298293 0.070079963 -3.5288017 6.798727e-04
2    alpha 2  0.109249097 0.057874966  1.8876745 6.252270e-02
3    alpha 3  0.268238436 0.034413442  7.7945832 1.560840e-11
4    alpha 4  0.869810761 0.090578360  9.6028539 3.635728e-15
5     beta 1  0.323989176 0.041041424  7.8941991 9.873279e-12
6     beta 2  0.055863164 0.033599356  1.6626261 1.001146e-01
7     beta 3 -0.078626097 0.020025657 -3.9262681 1.760820e-04
8     beta 4 -0.301226244 0.052987085 -5.6848993 1.848756e-07
9  gamma 1 1  0.104150200 0.020804442  5.0061521 3.018569e-06
10 gamma 1 2 -0.139880152 0.015303438 -9.1404395 3.094504e-14

Использование str(laaidsResult1) для поиска в выводе aidsEst для конкретных вещей, которые выИщите.