[R-br] Função não reconhecida

Paulo Justiniano paulojus em leg.ufpr.br
Sábado Outubro 8 10:05:55 BRT 2011


é só olhar a mensagem de erro e o código...
o objeto z nao usando da funcao nao existe...

em
 w<-matrix(c(beta,sigma2),ncol=2)
 colnames(z)<-c("beta","sigma2")
 return(w)

deve ser
 w<-matrix(c(beta,sigma2),ncol=2)
 colnames(w)<-c("beta","sigma2")
 return(w)



On Sat, 8 Oct 2011, Mauro Sznelwar wrote:

> Alguém sabe dizer o que está errado por aqui?
>  
> set.seed(123)
>  b<-5
>  sigma2<-2
>  n<-20
>  x<-runif(n,10,20)
>  y<-b*x+rnorm(n,0,sqrt(sigma2))
>  plot(x,y)
>  ###############################################
>  #Función para estimar la distribución posterior
>  f1<-function(beta0,nu,a,b,m,x,y)
>  {
>  A<-0
>  for(i in 1:n)
>  {
>  for(j in 1:n)
>  {
>  A<-A+x[j]*y[i]*(x[j]*y[i]-x[i]*y[j])/(sum(x^2)+nu)
>  }
>  }
>  B0<-(nu*sum((y-beta0*x)^2))/(sum(x)+nu)
>  med<-(sum(x*y)+nu*beta0)/(sum(x^2)+nu)
>  dp<-sqrt((A+B0+b)/((n+a)*(sum(x^2)+nu)))
>  beta<-rt(m,(n+a)/2)*dp+med
>  gamma1<-(n+a)/2;
>  delta1<-(A+B0+b)/2
>  sigma2<-1/rgamma(m,gamma1,delta1)
>  w<-matrix(c(beta,sigma2),ncol=2)
>  colnames(z)<-c("beta","sigma2")
>  return(w)
>  }
>  
> ###############################################
>  #Asignamos los valores para las ditribuciones
>  # apriori y el número de simulaciones
>  
> beta0<-4
>  nu<-0.00001
>  a<-1
>  b<-1
>  m<-1000
>  
> > simula<-f1(beta0,nu,a,b,m,x,y)
> Erro em colnames(z) <- c("beta", "sigma2") : objeto 'z' não encontrado
> > simula
> Erro: objeto 'simula' não encontrado
> 
>


Mais detalhes sobre a lista de discussão R-br