[R-br] Como calcular intervalos de confiança

[Benilton Carvalho]

uma extensao, baseada no codigo inicial do victor, que tambem funciona
para variacia desconhecida...

ic.media <- 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))
}

ex <- rnorm(100)
ic.media(ex)
ic.media(ex, 1)

b

2011/5/11 Victor Eduardo <victorduca08 em gmail.com>:
> Se quiser Intervalo de confianca para a media com variancia conhecida, você
> pode usar a função abaixo
>
> ic.media <- function(x,sigma2,conf.level=
> 0.95){
>
> n<-length(x)
>
> xbarra<-mean(x)
>
> z<-qnorm((1+conf.level)/2)
>
> LI <- xbarra - z*sqrt(sigma2)/sqrt(n)
>
> LS <- xbarra + z*sqrt(sigma2)/sqrt(n)
>
> return(c(LI,LS))
>
> }
>
>
> Se você não estiver interessado em calcular o intervalo de confiança, e sim
> apenas do resultado final, utilize o t.test()$conf.level
>
> Em 11 de maio de 2011 17:45, Victor Eduardo <victorduca08 em gmail.com>
> escreveu:
>>
>> Se quiser Intervalo de confianca para a media com variancia conhecida,
>> você pode usar a função abaixo
>>
>> ic.media <- function(x,sigma2,conf.level=0.95){
>>
>> n<-length(x)
>>
>> xbarra<-mean(x)
>>
>> z<-qnorm((1+conf.level)/2)
>>
>> LI <- xbarra - z*sqrt(sigma2)/sqrt(n)
>>
>> LS <- xbarra + z*sqrt(sigma2)/sqrt(n)
>>
>> return(c(LI,LS))
>>
>> }
>>
>>
>> Se você não estiver interessado em calcular o intervalo de confiança, e
>> sim apenas do resultado final, utilize o t.test()$conf.level
>>
>>
>>
>> Em 11 de maio de 2011 16:23, Cesar Rabak <cesar.rabak em gmail.com> escreveu:
>>>
>>> Em 11/5/2011 16:15, Claudio Carneiro escreveu:
>>>>
>>>> Olá a todos,
>>>>
>>>> Gostaria de saber como calcular intervalos de confiança para a média de
>>>> uma amostra.
>>>>
>>>
>>> http://lmgtfy.com/?q=calcular+intervalos+de+confian%C3%A7a+para+a+m%C3%A9dia+de+uma+amostra&l=1
>>>
>>> HTH
>>>
>>> --
>>> Cesar Rabak
>>> GNU/Linux User 52247.
>>> Get counted: http://counter.li.org/
>>> _______________________________________________
>>> R-br mailing list
>>> R-br em listas.c3sl.ufpr.br
>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>
>



-- 
Successful people ask better questions, and as a result, they get
better answers. (Tony Robbins)