Create a coefplot from matrix

Question

I am trying to create a coefplot (https://lrberge.github.io/fixest/reference/coefplot.html) from a matrix. (This is to show linear combinations of interaction terms, see below example code.) The documentation - see above link - indicates that this should be possible: object can be a matrix of coefficients table. Alternatively, I should be able to provide a vector of estimates, plus confidence interval information, to create the coefplot. I haven't been able to do so.

I have tried the following:

##packages and libraries required: 
##fixest
##biostat3

library(biostat3)
library(fixest)

##load example data
data(base_did)
base_inter = base_did

##Run simple regression with interactions
A <- lm(y ~factor(period)*treat, base_inter)

##Create linear combinations
C<-lincom(A,c("treat","treat+factor(period)2:treat","treat+factor(period)3:treat", "treat+factor(period)10:treat"))

##Transform to matrix
T<-as.matrix(C)

##Remove unnecessary information and dimension names
T<-T[,1:3]
dimnames(T) <- NULL

##run coefplot from fixest
fixest::coefplot(T[,1],ci_low=T[,2],ci_high=T[,3])

I get an error:

Error in data.frame(estimate = estimate, ci_low = ci_low, ci_high = ci_high,  : 
  arguments imply differing number of rows: 4, 0

I don't understand this, because the matrix T has three columns with four rows - so dimensions should be correct.

    [,1]       [,2]       [,3]     
[1,] -1.061001  -2.99068   0.8686775
[2,] -0.7100586 -2.639737  1.21962  
[3,] 1.458543   -0.4711357 3.388222 
[4,] 7.90124    5.971561   9.830918

Just a style note, by default T and F are shorthand for TRUE and FALSE. I'd strongly recommend not using T or F as object names--it could lead to bugs and also just looks weird to experienced users. — Gregor Thomas
– Gregor Thomas, Commented Jan 23, 2024 at 20:48
I've never used fixest::coefplot, but as a guess try naming your first argument, fixest::coefplot(estimate = T[,1], ci_low=T[,2], ci_high=T[,3]) — Gregor Thomas
– Gregor Thomas, Commented Jan 23, 2024 at 20:51

stefan · Accepted Answer · 2024-01-23 20:54:10Z

1

While your matrices C and T look like matrices of numbers they are actually matrices of lists. Hence, T[, 1], ... are lists too. To fix your issue you could use apply to unlist the columns before passing them to fixest::coefplot:

library(biostat3)
library(fixest)

T <- apply(T, 2, unlist)

## run coefplot from fixest
fixest::coefplot(
  T[, 1],
  ci_low = T[, 2],
  ci_high = T[, 3]
)

answered Jan 23, 2024 at 20:54

stefan

130k6 gold badges42 silver badges84 bronze badges

Sign up to request clarification or add additional context in comments.

Collectives™ on Stack Overflow

Create a coefplot from matrix

1 Answer 1

Comments

Your Answer

Hot Network Questions

Collectives™ on Stack Overflow

1 Answer 1

Comments

Your Answer

Sign up or log in

Post as a guest

Related