<!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>
<DIV><FONT face=Arial color=#0000ff>Muito obrigado, era isto mesmo, desculpem
minha ignorância!</FONT></DIV>
<DIV><FONT face=Arial><BR><FONT size=2></FONT></FONT></DIV><FONT face=Arial
size=2>Mauro,<BR><BR>na função f1 você quer dar nomes às colunas de 'z' mas em
nenhum<BR>momento você definiu 'z'.<BR>Creio que o que você queria era, na
verdade, 'colnames(w)' e não 'colnames(z)'.<BR><BR>Gustavo H. M. A.
Rocha<BR><BR><BR><BR>Em 8 de outubro de 2011 00:43, Mauro Sznelwar <</FONT><A
href="mailto:sznelwar@uol.com.br"><FONT face=Arial
size=2>sznelwar@uol.com.br</FONT></A><FONT face=Arial size=2>>
escreveu:<BR>> Alguém sabe dizer o que está errado por aqui?<BR>><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>><BR></FONT></BODY></HTML>