File: predplot.R

package info (click to toggle)
r-cran-learnbayes 2.15-2
  • links: PTS, VCS
  • area: main
  • in suites: jessie, jessie-kfreebsd
  • size: 2,064 kB
  • sloc: sh: 16; makefile: 1
file content (10 lines) | stat: -rw-r--r-- 366 bytes parent folder | download | duplicates (4) 
1
2
3
4
5
6
7
8
9
10
 
predplot=function(prior,n,yobs)
{
	y=0:n; a=prior[1]; b=prior[2]
	probs=pbetap(prior,n,y)
      m=max(probs)*1.05
	plot(y,probs,type="h",ylab="Probability",ylim=c(0,m),
 main=paste("Predictive Dist., beta(",a,",",b,") prior, n=",n,
            ", yobs=",yobs),lwd=2,col="blue")
	points(yobs,0,pch=19,cex=2.5,col="red")
      text(yobs,m/8,"yobs",col="red")}