Voce precisa inspecionar o objeto que contem os resultados para ver como utilizá-lo aqui vai uma sugestão: str(dist2) names(dist2) dim(dist2$beta) par(mfrow=c(1,3)) with(dist2, { hist(beta[,1], prob=T); lines(density(beta[,1])) hist(beta[,2], prob=T); lines(density(beta[,2])) hist(sigma, prob=T); lines(density(sigma)) }) On Mon, 18 Jun 2012, Eliana Silva wrote:
a função regressao linear é y=c+betaX+erro
estou fazendo a seguinte função library(LearnBayes)
dist<-lm(y~x) summary(dist)
#Simular a distribuição posteriori beta e sigma2 #1.Matriz X x<-cbind(1,x)
dist2<-blinreg(y, x, 5000) apply(dist2$beta,2,mean) apply(dist2$beta,2,sd)
podem me dizer como encontro a distribuição de c,beta e erro?
Estava a pensar fazer um gráfico
No dia 18 de Junho de 2012 00:06, Eder David Borges da Silva <eder@leg.ufpr.br> escreveu: Bernardo,As notas de aulas neste site podem ajudar bastante. http://faculty.agecon.vt.edu/moeltner/AAEC5126.html Att
Em 17 de junho de 2012 19:26, Bernardo Rodrigues <rodriguesbernardo22@gmail.com> escreveu: Alguem me pode ajudar como é que se determina distribuições posteriori da regressão linear a partir de distribuição a priori nao informativa?
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.