我用winbugs编程迭代得到了参数的后验分布值，现在想要通过图形看拟合的效果。该怎么画图呢？是自己编程的吗？还是说winbugs里有相关的应用板块。下面是我的程序的主要部分，我就是想画出fit[j]和y10[j]的图。望大家踊跃帮忙。
model{
for(i in 1:22){
a2[i,1]<-0
a2[i,2]<-0
a2[i,3]<-0
a2[i,4]<-0
a3[i,1]<-0
a3[i,2]<-0
a3[i,3]<-0
b[i,1:4]~dmnorm(a2[i,1:4],Sigma2[,])   
a[i,1:3]~dmnorm(a3[i,1:3],Sigma3[,])
}
for(j in 1:103){
z.star[j]<-(alpha[1]+a[y4[j],1])+(alpha[2]+a[y4[j],2])*y13[j]+(alpha[3]+a[y4[j],3])*y13[j]*y13[j]+delta2*w2[j]   ##skew-normal
w2[j]~dnorm(0,1)I(0,)
y11[j]~dnorm(z.star[j],tau2)   
betai1[j]<-beta[1]+b[y4[j],1]     
betai2[j]<-beta[2]+b[y4[j],2]
betai3[j]<-beta[3]+b[y4[j],3]
betaij4[j]<-beta[4]+beta[5]*z.star[j]+b[y4[j],4]    #mu.not[1]*Z[j,2]+mu.not[2]*Z[j,3]
dm1[j]<-betai1[j]-step(betai2[j]-betaij4[j])*betai2[j]*y13[j]
dm2[j]<-betai3[j]-step(betai2[j]-betaij4[j])*betaij4[j]*y13[j]
dm3[j]<-exp(dm1[j])
dm4[j]<-exp(dm2[j])
dm5[j]<-dm3[j]+dm4[j]
mu[j]<-log(dm5[j])+delta*w[j]
w[j]~dnorm(0,1)I(0,)
y10[j]~dnorm(mu[j],tau)   
fit[j]<-mu[j]
res[j]<-y10[j]-fit[j]
}