diff --git a/python/likelihoods/likelihood_function.py b/python/likelihoods/likelihood_function.py
index 61b5c427..50f9b620 100644
--- a/python/likelihoods/likelihood_function.py
+++ b/python/likelihoods/likelihood_function.py
@@ -88,7 +88,6 @@ class student_t(likelihood_function):
         Need to find what the variance is at the latent points for a student t*normal p(y*|f*)p(f*)
         (((g((v+1)/2))/(g(v/2)*s*sqrt(v*pi)))*(1+(1/v)*((y-f)/s)^2)^(-(v+1)/2))
         *((1/(s*sqrt(2*pi)))*exp(-(1/(2*(s^2)))*((y-f)^2)))
-
         """
 
         #We want the variance around test points y which comes from int p(y*|f*)p(f*) df*
@@ -144,9 +143,6 @@ class student_t(likelihood_function):
         p_025 = stats.scoreatpercentile(student_t_samples, .025, axis=1)[:, None]
         p_975 = stats.scoreatpercentile(student_t_samples, .975, axis=1)[:, None]
 
-        p_025 = 1+p_025
-        p_975 = 1+p_975
-
         ##Alernenately we could sample from int p(y|f*)p(f*|x*) df*
         def t_gaussian(f, mu, var):
             return (((gamma((self.v+1)*0.5)) / (gamma(self.v*0.5)*self.sigma*np.sqrt(self.v*np.pi))) * ((1+(1/self.v)*(((mu-f)/self.sigma)**2))**(-(self.v+1)*0.5))