Boa noite pessoal, estou tentando criar uma matriz com os intervalos de confiança que eu calculei para os um conjunto de dados. No entanto, os valores calculados para o intervalo de confiança ficaram ocultos na matriz.
ic<- function(x, sigma2, conf.level=0.95){
n <- length(x)
xbarra <- mean(x)
if (!missing(sigma2)){
z <- qnorm((1+conf.level)/2)
}else{
sigma2 <- var(x)
z <- qt((1+conf.level)/2, n-1)
}
return(mean(x)+c(-1, 1)*z*sqrt(sigma2)/sqrt(n))
}
tr<-gl(3,4)
resp<-rnorm(12,30,5)
d<-data.frame(tr,resp)
int.conf<-tapply(resp,tr,ic)
mat<-matrix(int.conf)
mat
# matriz obtida.
[,1]
Agradeço desde já pela ajuda.
Att.
Tiago.