Quero criar uma matriz de 38 linhas e duas colunas com coeficientes nas linhas ímpares e os intervalos de confiança para os parâmetros nas linhas pares.
Em seguida chamar o stargazer.
Meu código é:
library(quantreg)
y <- rnorm(200)
x <- matrix(rnorm(123),200)
fit <- rq(y~x,tau=seq(0.05, 0.95, by = 0.05),method="br")
summary(fit,se="nid")
Consegui avançar até aqui:
Repare que consigo colocar os coeficientes nas linhas ímpares sem problemas. O erro óbvio está quando tento construir e colocar o intervalo de confiança dos parâmetros nas linhas pares:
fit_Matrix<-matrix(0,38,2)
for(i in 1:19) for (j in seq(1, 38, by = 2)[i]) for (k in seq(2, 38, by = 2)[i]) {
fit_Matrix[j,1]<-summary(fit,se="nid")[[i]]$coefficients[1]
fit_Matrix[j,2]<-summary(fit,se="nid")[[i]]$coefficients[2]
fit_Matrix[k,1]<-summary(fit,se="nid")[[i]]$coefficients[1]-1.645*summary(fit,se="nid")[[i]]$coefficients[3];
summary(fit,se="nid")[[i]]$coefficients[1]+1.645*summary(fit,se="nid")[[i]]$coefficients[3]
fit_Matrix[k,2]<-summary(fit,se="nid")[[i]]$coefficients[2]-1.645*summary(fit,se="nid")[[1]]$coefficients[4];
summary(fit,se="nid")[[i]]$coefficients[2]-1.645*summary(fit,se="nid")[[i]]$coefficients[4]
}
A segunda parte do código está errada. Como faço para corrigir? Lembrando que gostaria de transporta a matriz para o LaTeX.
Obrigado, desde já.
lenght
é 14. Mas a informação que eu quero esta nosummary
que possui 19. Não é? Pelo que vi 14 é o tamanho da lista:str(fit)
.