Visualise the likelihood function

ch08.r

@@ -0,0 +1,20 @@
+#############
+# Chapter 8.1 Maximum-likelihood Estimation
+
+## Visualizing the likelihood function
+
+# R: Visualize the likelooh function
+library(pracma)
+library(plot3D)
+
+N = 50
+S = 1:N
+theta = seq(0.1, 0.9, (0.9+0.1)/100)
+mesh = meshgrid(S, theta)
+S_grid = mesh$X
+theta_grid = mesh$Y
+L = S_grid * log(theta_grid) + (N-S_grid) * log(1-theta_grid)
+L = t(L)
+persp3D(S, theta, L, theta=65, phi=15)
+
+#############