Displaying 2 results from an estimated 2 matches for "ybar2".
Did you mean:
bar2
2010 Nov 18
0
On efficiency, Vectorize and loops
...)local({force(t);
function(u){
crossprod(xi[,t],sqrt(2)*cos(pi*outer(1:d,u))) +
crossprod(2^(1-(1:10))*Z[,t],sqrt(2)+sin(pi*outer(1:10,u)))
}
}))
######################
# remove later
Ybar<-function(u){
Ybar<-0
for (t in 1:n){
Ybar<-Ybar+Y[[t]](u)
}
Ybar<-1/n*Ybar
}
Ybar2<-function(u){
Ybar2<-0
for (t in 1:n){
Ybar2<-Ybar2+Y2[[t]](u)
}
Ybar2<-1/n*Ybar2
}
Ybar3<-function(u){
Ybar3<-0
for (t in 1:n){
Ybar3<-Ybar3+Y3[[t]](u)
}
Ybar3<-1/n*Ybar3
}
Ybar4<-function(u){
Ybar4<-0
for (t in 1:n){
Ybar4<-Ybar4+...
2010 Jul 18
2
loop troubles
...to have elements from
the previous posterior to be the new prior and now I cant figure out how to
do the code below in a loop. The data below is dummy data, I used a starting
mu of 1, and starting precision of 0.
bayes.analysis.treat<-function(mu0,p0){
n1 = 5
n2 = 10
n3 = 15
ybar1 = 12
ybar2 = 13
ybar3 = 14
sd1 = 2
sd2 = 3
sd3 = 4
#posterior 1
var1 = sd1^2 #sample variance
p1 = n1/var1 #sample precision
p1n = p0+p1
mu1 = ((p0)/(p1n)*mu0)+((p1)/(p1n))*ybar1
sigma1 = 1/sqrt(p1n)
#posterior 2
var2 = sd2^2 #sample variance
p2 = n2/var2 #sample precision
p2n = p1n+p2
mu2 = ((...