We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent 5b53e9c commit 7d7c3aaCopy full SHA for 7d7c3aa
chapter10/linRegVb.m
@@ -46,7 +46,7 @@
46
KLw = -sum(log(diag(U)));
47
% q(alpha)
48
w2 = dot(Ew,Ew);
49
- invU = U\I;
+ invU = U'\I;
50
trS = dot(invU(:),invU(:));
51
b = b0+0.5*(w2+trS); % 10.95
52
Ealpha = a/b; % 10.102
chapter10/rvmRegVb.m
@@ -49,7 +49,7 @@
w2 = Ew.*Ew;
53
dgS = dot(invU,invU,2);
54
b = b0+0.5*(w2+dgS);
55
Ealpha = a./b;
0 commit comments