Rui Barradas já deu uma boa resposta sobre como calcular o intervalo de confiança, esta é para entender o que estava fazendo de errado. Porque não é apenas um erro de programação, mas de entendimento do IC.
A mensagem de erro
O erro que o R está retornando vem de usar um valor negativo para o desvio padrão na função rnorm
:
> rnorm(5, -.4, .1)
[1] -0.4070113 -0.3860859 -0.3772370 -0.4746763 -0.3817434
> rnorm(5, -.4, -.1)
[1] NaN NaN NaN NaN NaN
Warning message:
In rnorm(5, -0.4, -0.1) : NAs produced
Isso porque desvio padrão não assume valores negativos. Mas por que está usando rnorm
para calcular os limites? Por que usar o valor absoluto da média?
Cálculo do IC
Existe um motivo para os pacotes básicos do R não terem uma função para calcular IC. Ele é uma medida do quão confiável é a estimativa de um parâmetro, o modo com que é calculado depende do parâmetro e do modelo atribuído aos dados.
Sua fórmula está correta para mu estimado de uma amostra pequena (ou com grande variância) seguindo distribuição próxima da normal. Compare com o resultado de um teste t:
ic <- media + qt(c(0.025, 0.975), df = n-1) * sqrt(var/n)
> ic
[1] -0.8479307 -0.4110068
> t.test(dados$y)$conf.int
[1] -0.8479307 -0.4110068
attr(,"conf.level")
[1] 0.95
Mas esses são bons pressupostos para seus dados? Contraste com o resultado da função do @rui-barradas. E compare com o resultado usando diferentes métodos com a média estimada por bootstrap:
library(boot)
media.b <- boot(dados$y, function(x,i) mean(x[i]), 10000)
> boot.ci(media.b, conf = .95, type = c("norm","basic", "perc", "bca"))
BOOTSTRAP CONFIDENCE INTERVAL CALCULATIONS
Based on 10000 bootstrap replicates
CALL :
boot.ci(boot.out = media.b, conf = 0.95, type = c("norm", "basic",
"perc", "bca"))
Intervals :
Level Normal Basic
95% (-0.8366, -0.4206 ) (-0.8527, -0.4365 )
Level Percentile BCa
95% (-0.8224, -0.4063 ) (-0.8041, -0.3735 )
Calculations and Intervals on Original Scale
Adicionando limites aos dados
Uma vez que calculou o intervalo, bastaria adicionar/subtrair os valores de y:
df$lim_inf <- df$y - ic[1]
df$lim_sup <- df$y + ic[2]
Mas isso é errado. O intervalo de confiança aplica-se ao parâmetro estimado (a média, no caso), não aos valores individuais. Uma opção mais adequada é usar o erro padrão:
ep <- sqrt(var/n)
df$lim_inf <- df$y - ep
df$lim_sup <- df$y + ep
media<-abs(mean(df$y))
? Existe alguma restrição para a média ser não-negativa? Porque, sem calcular o valor absoluto dela, seu valor é negativo.rnorm(nrow(df), media/-0.832, media/-0.399)
: como é que o desvio padrão pode ser negativo?