search for: a_csmh

...e size in each stratum must be > 1") I <- dim(x)[1]; J <- dim(x)[2]; K <- dim(x)[3] df_GMH <- (I - 1) * (J - 1); df_SMH <- I - 1; df_CSMH <- 1 A_GMH <- cbind(diag(I-1), 0) %x% cbind(diag(J-1), 0) A_SMH <- cbind(diag(I-1), 0) %x% t(1:J) A_CSMH <- t(1:I) %x% t(1:J) Y_GMH <- matrix(0, nc = 1, nr = df_GMH) Y_SMH <- matrix(0, nc = 1, nr = df_SMH) Y_CSMH <- matrix(0, nc = 1, nr = df_CSMH) S_GMH <- matrix(0, nc = df_GMH, nr = df_GMH) S_SMH <- matrix(0, nc = df_SMH, nr = df_SMH) S_CSMH <- matrix(0, n...

....null(col_scores)) { col_scores =3D 1:J } df_GMH <- (I - 1) * (J - 1); df_SMH <- I - 1; df_CSMH <- 1 A_GMH <- cbind(diag(I-1), 0) %x% cbind(diag(J-1), 0) #A_SMH <- cbind(diag(I-1), 0) %x% t(1:J) A_SMH <- cbind(diag(I-1), 0) %x% t(col_scores) #A_CSMH <- t(1:I) %x% t(1:J) A_CSMH <- t(row_scores) %x% t(col_scores) Y_GMH <- matrix(0, nc =3D 1, nr =3D df_GMH) Y_SMH <- matrix(0, nc =3D 1, nr =3D df_SMH) Y_CSMH <- matrix(0, nc =3D 1, nr =3D df_CSMH) S_GMH <- matrix(0, nc =3D df_GMH, nr =3D df_GMH) S_SMH <-...