transpose

Jelena Markovic · Jelena Markovic · commit d2a3b9bf075f · 2017-11-07T23:04:22.000-08:00
diff --git a/selectiveInference/R/funs.randomized.R b/selectiveInference/R/funs.randomized.R
@@ -265,7 +265,7 @@ importance_weight = function(noise_scale,
 get_mean_cov = function(noise_scale, linear_term, offset_term){
     temp = solve(t(linear_term) %*% linear_term)
     cov = noise_scale^2*temp
-    mean = temp %*% linear_term %*% offset_term
+    mean = temp %*% t(linear_term) %*% offset_term
     return(list(mean=mean, cov=cov))
 }
 

Original file line number	Diff line number	Diff line change
`@@ -265,7 +265,7 @@ importance_weight = function(noise_scale,`
`265`	`265`	`get_mean_cov = function(noise_scale, linear_term, offset_term){`
`266`	`266`	`temp = solve(t(linear_term) %*% linear_term)`
`267`	`267`	`cov = noise_scale^2*temp`
`268`		`- mean = temp %% linear_term %% offset_term`
	`268`	`+ mean = temp %% t(linear_term) %% offset_term`
`269`	`269`	`return(list(mean=mean, cov=cov))`
`270`	`270`	`}`
`271`	`271`