diff --git a/GPy/likelihoods/Gaussian.py b/GPy/likelihoods/Gaussian.py
new file mode 100644
index 00000000..2397ce38
--- /dev/null
+++ b/GPy/likelihoods/Gaussian.py
@@ -0,0 +1,16 @@
+import numpy as np
+
+class Gaussian:
+    def __init__(self,data,variance=1.,normalise=False):
+        self.data = data
+        if normalise:
+            foo
+        self._variance = variance
+    def _get_params(self):
+        return np.asarray(self.variance)
+    def _set_params(self,x):
+        self._variance = x
+    def fit(self):
+        pass
+    def _gradients(self,foo):
+        return bar(foo)
diff --git a/GPy/likelihoods/__init__.py b/GPy/likelihoods/__init__.py
new file mode 100644
index 00000000..d1369c43
--- /dev/null
+++ b/GPy/likelihoods/__init__.py
@@ -0,0 +1,3 @@
+from EP import EP
+from Gaussian import Gaussian
+# TODO: from Laplace import Laplace