<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content="text/html; charset=iso-8859-1">
<META content="MSHTML 6.00.2900.3429" name=GENERATOR>
<STYLE></STYLE>
</HEAD>
<BODY bgColor=#ffffff>
<DIV><FONT face=Arial color=#0000ff>Alguém sabe dizer o que está errado por
aqui?</FONT></DIV>
<DIV><FONT face=Arial size=2> <BR>
set.seed(123)<BR> b<-5<BR> sigma2<-2<BR> n<-20<BR> x<-runif(n,10,20)<BR> y<-b*x+rnorm(n,0,sqrt(sigma2))<BR> plot(x,y)<BR> ###############################################<BR> #Función
para estimar la distribución
posterior<BR> f1<-function(beta0,nu,a,b,m,x,y)<BR> {<BR> A<-0<BR> for(i
in 1:n)<BR> {<BR> for(j in
1:n)<BR> {<BR> A<-A+x[j]*y[i]*(x[j]*y[i]-x[i]*y[j])/(sum(x^2)+nu)<BR> }<BR> }<BR> B0<-(nu*sum((y-beta0*x)^2))/(sum(x)+nu)<BR> med<-(sum(x*y)+nu*beta0)/(sum(x^2)+nu)<BR> dp<-sqrt((A+B0+b)/((n+a)*(sum(x^2)+nu)))<BR> beta<-rt(m,(n+a)/2)*dp+med<BR> gamma1<-(n+a)/2;<BR> delta1<-(A+B0+b)/2<BR> sigma2<-1/rgamma(m,gamma1,delta1)<BR> w<-matrix(c(beta,sigma2),ncol=2)<BR> colnames(z)<-c("beta","sigma2")<BR> return(w)<BR> }<BR> <BR>
###############################################<BR> #Asignamos los valores
para las ditribuciones<BR> # apriori y el número de
simulaciones<BR> <BR>
beta0<-4<BR> nu<-0.00001<BR> a<-1<BR> b<-1<BR> m<-1000<BR> <BR>>
simula<-f1(beta0,nu,a,b,m,x,y)<BR>Erro em colnames(z) <- c("beta",
"sigma2") : objeto 'z' não encontrado<BR>> simula<BR>Erro: objeto 'simula'
não encontrado<BR></FONT></DIV></BODY></HTML>