[R-br] Função não reconhecida
Mauro Sznelwar
sznelwar em uol.com.br
Domingo Outubro 9 20:13:46 BRT 2011
Muito obrigado, era isto mesmo, desculpem minha ignorância!
Mauro,
na função f1 você quer dar nomes às colunas de 'z' mas em nenhum
momento você definiu 'z'.
Creio que o que você queria era, na verdade, 'colnames(w)' e não 'colnames(z)'.
Gustavo H. M. A. Rocha
Em 8 de outubro de 2011 00:43, Mauro Sznelwar <sznelwar em uol.com.br> escreveu:
> 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
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111009/fef850bd/attachment.html>
Mais detalhes sobre a lista de discussão R-br