Dúvida Inferência Bayesiana
Galera, qual comando que me da a distribuição gama inversa? Como que calculo intervalo de credibilidade? Abraços, Victor Eduardo
teoria de probabilidade: se X eh Gamma(a, b), entao Y=1/X eh Gamma(a, 1/b)... 2011/5/24 Victor Eduardo <victorduca08@gmail.com>:
Galera, qual comando que me da a distribuição gama inversa? Como que calculo intervalo de credibilidade?
Abraços,
Victor Eduardo
_______________________________________________ R-br mailing list R-br@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)
Oi Victor, Pode usar se inspirar no exemplo abaixo: library(TeachingDemos) x=c(62,58,61,60,51,61, 66, 63, 71, 59, 68, 50, 62, 62, 58, 54, 56, 58, 49, 65) p <- 0.95 s2=var(x) gl=length(x)-1 pinvchisq <- function(x, df=gl, scale=s2) { prob <- pchisq(scale/x, df) return(1-prob) } qinvchisq <- function(p, df=gl, scale) { qprob <- qchisq(1-p, df) return(scale/qprob) } hpd(qinvchisq,df=gl,scale=s2,conf=.95) Att, Greice. 2011/5/23 Benilton Carvalho <beniltoncarvalho@gmail.com>
teoria de probabilidade: se X eh Gamma(a, b), entao Y=1/X eh Gamma(a, 1/b)...
2011/5/24 Victor Eduardo <victorduca08@gmail.com>:
Galera, qual comando que me da a distribuição gama inversa? Como que calculo intervalo de credibilidade?
Abraços,
Victor Eduardo
_______________________________________________ R-br mailing list R-br@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) _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Vitor, vc pode gerar de uma Gamma comum e inverter, depois pode olhar os quantis. Att. Em 23 de maio de 2011 21:03, Greice Laureano <greice.laureano@gmail.com>escreveu:
Oi Victor,
Pode usar se inspirar no exemplo abaixo:
library(TeachingDemos) x=c(62,58,61,60,51,61, 66, 63, 71, 59, 68, 50, 62, 62, 58, 54, 56, 58, 49, 65) p <- 0.95 s2=var(x) gl=length(x)-1
pinvchisq <- function(x, df=gl, scale=s2) { prob <- pchisq(scale/x, df) return(1-prob) }
qinvchisq <- function(p, df=gl, scale) { qprob <- qchisq(1-p, df) return(scale/qprob) }
hpd(qinvchisq,df=gl,scale=s2,conf=.95)
Att, Greice.
2011/5/23 Benilton Carvalho <beniltoncarvalho@gmail.com>
teoria de probabilidade: se X eh Gamma(a, b), entao Y=1/X eh Gamma(a, 1/b)...
2011/5/24 Victor Eduardo <victorduca08@gmail.com>:
Galera, qual comando que me da a distribuição gama inversa? Como que calculo intervalo de credibilidade?
Abraços,
Victor Eduardo
_______________________________________________ R-br mailing list R-br@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) _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.
# Gamma inversa g <- 1/gamma(n,a,b) Fábio Mathias Corrêa Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e da Terra - DCET Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna CEP 45662-900. Ilhéus-Bahia Tel.: 73-3680-5076 ________________________________ De: Vinicius Brito Rocha <viniciusbritor@gmail.com> Para: r-br@listas.c3sl.ufpr.br Enviadas: Segunda-feira, 23 de Maio de 2011 21:09 Assunto: Re: [R-br] Dúvida Inferência Bayesiana Vitor, vc pode gerar de uma Gamma comum e inverter, depois pode olhar os quantis. Att. Em 23 de maio de 2011 21:03, Greice Laureano <greice.laureano@gmail.com> escreveu: Oi Victor,
Pode usar se inspirar no exemplo abaixo:
library(TeachingDemos) x=c(62,58,61,60,51,61, 66, 63, 71, 59, 68, 50, 62, 62, 58, 54, 56, 58, 49, 65) p <- 0.95 s2=var(x) gl=length(x)-1
pinvchisq <- function(x, df=gl, scale=s2) { prob <- pchisq(scale/x, df) return(1-prob) }
qinvchisq <- function(p, df=gl, scale) { qprob <- qchisq(1-p, df) return(scale/qprob) }
hpd(qinvchisq,df=gl,scale=s2,conf=.95)
Att, Greice.
2011/5/23 Benilton Carvalho <beniltoncarvalho@gmail.com>
teoria de probabilidade: se X eh Gamma(a, b), entao Y=1/X eh Gamma(a, 1/b)...
2011/5/24 Victor Eduardo <victorduca08@gmail.com>:
Galera, qual comando que me da a distribuição gama inversa? Como que calculo intervalo de credibilidade?
Abraços,
Victor Eduardo
_______________________________________________ R-br mailing list R-br@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) _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein. _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
# Gamma inversa g <- 1/rgamma(n,a,b) Fábio Mathias Corrêa Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e da Terra - DCET Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna CEP 45662-900. Ilhéus-Bahia Tel.: 73-3680-5076 ________________________________ De: Vinicius Brito Rocha <viniciusbritor@gmail.com> Para: r-br@listas.c3sl.ufpr.br Enviadas: Segunda-feira, 23 de Maio de 2011 21:09 Assunto: Re: [R-br] Dúvida Inferência Bayesiana Vitor, vc pode gerar de uma Gamma comum e inverter, depois pode olhar os quantis. Att. Em 23 de maio de 2011 21:03, Greice Laureano <greice.laureano@gmail.com> escreveu: Oi Victor,
Pode usar se inspirar no exemplo abaixo:
library(TeachingDemos) x=c(62,58,61,60,51,61, 66, 63, 71, 59, 68, 50, 62, 62, 58, 54, 56, 58, 49, 65) p <- 0.95 s2=var(x) gl=length(x)-1
pinvchisq <- function(x, df=gl, scale=s2) { prob <- pchisq(scale/x, df) return(1-prob) }
qinvchisq <- function(p, df=gl, scale) { qprob <- qchisq(1-p, df) return(scale/qprob) }
hpd(qinvchisq,df=gl,scale=s2,conf=.95)
Att, Greice.
2011/5/23 Benilton Carvalho <beniltoncarvalho@gmail.com>
teoria de probabilidade: se X eh Gamma(a, b), entao Y=1/X eh Gamma(a, 1/b)...
2011/5/24 Victor Eduardo <victorduca08@gmail.com>:
Galera, qual comando que me da a distribuição gama inversa? Como que calculo intervalo de credibilidade?
Abraços,
Victor Eduardo
_______________________________________________ R-br mailing list R-br@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) _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein. _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Valeu pela ajuda, galera!! Consegui aqui! Abraços! Em 24 de maio de 2011 06:27, Fabio Mathias Corrêa <fabio.ufla@yahoo.com.br>escreveu:
# Gamma inversa
g <- 1/gamma(n,a,b)
Fábio Mathias Corrêa
Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e da Terra - DCET
Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna CEP 45662-900. Ilhéus-Bahia
Tel.: 73-3680-5076 ------------------------------ *De:* Vinicius Brito Rocha <viniciusbritor@gmail.com> *Para:* r-br@listas.c3sl.ufpr.br *Enviadas:* Segunda-feira, 23 de Maio de 2011 21:09 *Assunto:* Re: [R-br] Dúvida Inferência Bayesiana
Vitor,
vc pode gerar de uma Gamma comum e inverter, depois pode olhar os quantis.
Att.
Em 23 de maio de 2011 21:03, Greice Laureano <greice.laureano@gmail.com>escreveu:
Oi Victor,
Pode usar se inspirar no exemplo abaixo:
library(TeachingDemos) x=c(62,58,61,60,51,61, 66, 63, 71, 59, 68, 50, 62, 62, 58, 54, 56, 58, 49, 65) p <- 0.95 s2=var(x) gl=length(x)-1
pinvchisq <- function(x, df=gl, scale=s2) { prob <- pchisq(scale/x, df) return(1-prob) }
qinvchisq <- function(p, df=gl, scale) { qprob <- qchisq(1-p, df) return(scale/qprob) }
hpd(qinvchisq,df=gl,scale=s2,conf=.95)
Att, Greice.
2011/5/23 Benilton Carvalho <beniltoncarvalho@gmail.com>
teoria de probabilidade: se X eh Gamma(a, b), entao Y=1/X eh Gamma(a, 1/b)...
2011/5/24 Victor Eduardo <victorduca08@gmail.com>:
Galera, qual comando que me da a distribuição gama inversa? Como que calculo intervalo de credibilidade?
Abraços,
Victor Eduardo
_______________________________________________ R-br mailing list R-br@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) _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional
www.aplicademic.blogspot.com http://twitter.com/viniciusbritor
"Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
participantes (5)
-
Benilton Carvalho -
Fabio Mathias Corrêa -
Greice Laureano -
Victor Eduardo -
Vinicius Brito Rocha