# DEMONSTRATION OF SAMPLING FROM A PRIOR OVER FUNCTIONS DEFINED USING AN MLP source("mlp.r") X <- matrix (seq(-1,1,by=0.01), ncol=1) p <- 1 par(mfrow=c(2,2)) # ---------------------------------------------------------------------- set.seed(1) m <- 1 plot(0,0,type="n",xlim=c(-1,1),ylim=c(-10,10),xlab="x",ylab="f(x)") for (i in 1:10) { # w2_0 w2 w1_0 w1 w <- c (rnorm(1,0,1), rnorm(m,0,4/sqrt(m)), rnorm(m,0,3), rnorm(p*m,0,7)) f <- mlp.forward (X, m, w) $ o lines(X,f,type="l",col=i) } # ---------------------------------------------------------------------- set.seed(1) m <- 100 plot(0,0,type="n",xlim=c(-1,1),ylim=c(-10,10),xlab="x",ylab="f(x)") for (i in 1:10) { # w2_0 w2 w1_0 w1 w <- c (rnorm(1,0,1), rnorm(m,0,4/sqrt(m)), rnorm(m,0,3), rnorm(p*m,0,7)) f <- mlp.forward (X, m, w) $ o lines(X,f,type="l",col=i) } # ---------------------------------------------------------------------- set.seed(1) m <- 100 plot(0,0,type="n",xlim=c(-1,1),ylim=c(-10,10),xlab="x",ylab="f(x)") for (i in 1:10) { # w2_0 w2 w1_0 w1 w <- c (rnorm(1,0,1), rnorm(m,0,4/sqrt(m)), rnorm(m,0,3*5), rnorm(p*m,0,7*5)) f <- mlp.forward (X, m, w) $ o lines(X,f,type="l",col=i) } # ---------------------------------------------------------------------- set.seed(1) m <- 100 plot(0,0,type="n",xlim=c(-1,1),ylim=c(-10,10),xlab="x",ylab="f(x)") for (i in 1:10) { # w2_0 w2 w1_0 w1 w <- c (rnorm(1,0,4), rnorm(m,0,4/sqrt(m)), rnorm(m,0,3/5), rnorm(p*m,0,7/5)) f <- mlp.forward (X, m, w) $ o lines(X,f,type="l",col=i) }