This function retrives the fixed effects from a femlm estimation. It is useful only when there are more than one cluster.
A femlm object.
A list containig the vectors of the fixed effects.
See also the main estimation function
summary.femlm to see the results with the appropriate standard-errors,
getFE to extract the cluster coefficients, and the functions
res2tex to visualize the results of multiple estimations.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
# Bilateral network nb = 20 n = nb**2 k = nb id1 = factor(rep(1:k, each=n/k)) id2 = factor(rep(1:(n/k), times=k)) d = rep(rnorm(k)**2, each=n/k) x = rnorm(n, 1, 5)**2 y = rnorm(n, -1, 5)**2 z = rpois(n, x*y+rnorm(n, sd = 3)**2) + rpois(n, 2) base = data.frame(x, y, z, id1, id2) # We want to use the ID's of each observation as a variable: we use the option cluster est_poisson = femlm(z~log(x)+log(y), base, family="poisson", cluster=c("id1", "id2")) # To get the FE: myFE = getFE(est_poisson)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.