- res = 0
- for (i in 1:d)
- {
- term = 0
- for (k in 1:K)
- term = term + p[k]*.G(link,1,λ[k],b[k])*β[i,k]
- res = res + (term - M1[i])^2
- for (j in 1:d)
- {
- term = 0
- for (k in 1:K)
- term = term + p[k]*.G(link,2,λ[k],b[k])*β2[i,j,k]
- res = res + (term - M2[i,j])^2
- for (l in 1:d)
- {
- term = 0
- for (k in 1:K)
- term = term + p[k]*.G(link,3,λ[k],b[k])*β3[i,j,l,k]
- res = res + (term - M3[i,j,l])^2
- }
- }
- }
- res
+ res = 0
+ for (i in 1:d)
+ {
+ term = 0
+ for (k in 1:K)
+ term = term + p[k]*.G(link,1,λ[k],b[k])*β[i,k]
+ res = res + (term - M1[i])^2
+ for (j in 1:d)
+ {
+ term = 0
+ for (k in 1:K)
+ term = term + p[k]*.G(link,2,λ[k],b[k])*β2[i,j,k]
+ res = res + (term - M2[i,j])^2
+ for (l in 1:d)
+ {
+ term = 0
+ for (k in 1:K)
+ term = term + p[k]*.G(link,3,λ[k],b[k])*β3[i,j,l,k]
+ res = res + (term - M3[i,j,l])^2
+ }
+ }
+ }
+ res