/dports/finance/R-cran-plm/plm/R/ |
H A D | tool_ercomp.R | 308 CPZM <- solve(crossprod(Z)) 356 CPXM <- solve(crossprod(WX)) 393 …M["id", "eta"] <- O - sum(TG) - trace(crossprod(t(solve(crossprod(XBetaBlambda))), crossprod(XS… 396 … M["gp", "eta"] <- sum(TG) - trace(crossprod(t(CPZBlambda.solve), crossprod(ZBlambda, ZSeta))) 397 … M["gp", "lambda"] <- O - trace(crossprod(t(CPZBlambda.solve), crossprod(ZSlambda, Z))) 493 CPZM <- solve(crossprod(Z)) 500 CPZSmu <- crossprod(ZSmu, Z) 505 CPZW <- crossprod(ZW) 520 CPZBeta <- crossprod(ZBeta) 530 CPZBmu <- crossprod(ZBmu) [all …]
|
H A D | tool_vcovG.R | 554 Sl[ , , i-l] <- crossprod(X, E(u, ul)) %*% Xl 574 pane <- solve(crossprod(demX)) 576 mycov <- tcrossprod(crossprod(t(pane), salame), t(pane)) # == pane %*% salame %*% pane 1121 salame[ , , i] <- crossprod(xi, OmegaTi) %*% xi 1127 pane <- solve(crossprod(demX)) 1130 mycov <- tcrossprod(crossprod(t(pane), salame), t(pane)) # == pane %*% salame %*% pane 1200 res1s <- lapply(yX, function(x) x[ , 1L] - crossprod(t(x[ , -1L, drop = FALSE]), coef1s)) 1213 z <- crossprod(t(x[ , k, drop = FALSE]), t(y)) 1220 crossprod(x, crossprod(y, x)), 1225 vcovr <- B2 + crossprod(t(D), B2) + t(crossprod(t(D), B2)) + D %*% vcov1s %*% t(D) [all …]
|
/dports/math/R/R-4.1.2/tests/ |
H A D | reg-BLAS.R | 48 cprod(crossprod(m, m), crossprod(Cm, Cm), matrix(c(5,11,11,25), 2, 2) ) 49 cprod(crossprod(m), crossprod(Cm), matrix(c(5,11,11,25), 2, 2) ) 50 cprod(crossprod(m, cv), crossprod(Cm, Ccv), matrix(c(35,81), 2, 1) ) 51 cprod(crossprod(m, v), crossprod(Cm, Cv), matrix(c(35,81), 2, 1) ) 52 cprod(crossprod(cv, m), crossprod(Ccv, Cm), matrix(c(35,81), 1, 2) ) 53 cprod(crossprod(v, m), crossprod(Cv, Cm), matrix(c(35,81), 1, 2) ) 54 cprod(crossprod(cv, cv), crossprod(Ccv, Ccv), matrix(265,1,1) ) 55 cprod(crossprod(v, v), crossprod(Cv, Cv), matrix(265,1,1) ) 56 cprod(crossprod(rv, rv), crossprod(Crv, Crv), matrix(c(121,132,132,144), 2, 2) ) 96 cprod(crossprod(m1, cv1), crossprod(Cm1, Ccv1), matrix(c(Inf,182), 2, 1) ) [all …]
|
/dports/math/libRmath/R-4.1.1/tests/ |
H A D | reg-BLAS.R | 48 cprod(crossprod(m, m), crossprod(Cm, Cm), matrix(c(5,11,11,25), 2, 2) ) 49 cprod(crossprod(m), crossprod(Cm), matrix(c(5,11,11,25), 2, 2) ) 50 cprod(crossprod(m, cv), crossprod(Cm, Ccv), matrix(c(35,81), 2, 1) ) 51 cprod(crossprod(m, v), crossprod(Cm, Cv), matrix(c(35,81), 2, 1) ) 52 cprod(crossprod(cv, m), crossprod(Ccv, Cm), matrix(c(35,81), 1, 2) ) 53 cprod(crossprod(v, m), crossprod(Cv, Cm), matrix(c(35,81), 1, 2) ) 54 cprod(crossprod(cv, cv), crossprod(Ccv, Ccv), matrix(265,1,1) ) 55 cprod(crossprod(v, v), crossprod(Cv, Cv), matrix(265,1,1) ) 56 cprod(crossprod(rv, rv), crossprod(Crv, Crv), matrix(c(121,132,132,144), 2, 2) ) 96 cprod(crossprod(m1, cv1), crossprod(Cm1, Ccv1), matrix(c(Inf,182), 2, 1) ) [all …]
|
/dports/math/R-cran-car/car/R/ |
H A D | wcrossprod.R | 10 if (missing(y)) return(crossprod(x)) else return(crossprod(x, y)) 13 if (missing (y)) return(w[1]*crossprod(x)) else return(w[1]*crossprod(x, y)) 18 if (missing(y)) return(crossprod(x, w*x)) else return(crossprod(x, w*y)) 22 if (missing(y)) return(crossprod(x, w %*% x)) else return(crossprod(x, w %*% y))
|
/dports/math/R-cran-spdep/spdep/R/ |
H A D | SpatialFiltering.R | 121 M <- diag(1,nofreg) - X %*% solve(crossprod(X),t(X)) 137 cyMy <- crossprod(y, M) %*% y 138 cyMSMy <- crossprod(y, MSM) %*% y 178 res <- y - X %*% solve(crossprod(X), crossprod(X, y)) 180 if (((crossprod(res, S) %*% res) / crossprod(res)) < 0) acsign <- -1 188 Xcoeffs <- solve(crossprod(X), crossprod(X, y)) 196 crossprod(vec[,j], y)) 220 res <- y - xe %*% solve(crossprod(xe), crossprod(xe, y)) 222 mi <- (crossprod(res, S) %*% res) / crossprod(res) 245 M <- diag(1,nofreg) - X %*% solve(crossprod(X),t(X)) [all …]
|
H A D | kpgm_new.R | 159 SSE <- c(crossprod(et)) 185 J[2,1] <- crossprod(KP04de, KP04eo) + crossprod(KP04se) 186 J[1,2] <- - KP04c*(crossprod(KP04de) - KP04a*crossprod(KP04se)) 192 A2N <- crossprod(W) 208 JJI <- 1/crossprod(J1) 533 uu <- crossprod(u,u) 534 uwu <- crossprod(u,wu) 535 uwpuw <- crossprod(wu,wu) 536 uwwu <- crossprod(u,wwu) 537 wwupwu <- crossprod(wwu,wu) [all …]
|
H A D | aple.R | 18 if (is.complex(eig)) trWW <- Re(crossprod(eig)) 19 else trWW <- crossprod(eig) 27 W2 <- crossprod(W) + corterm 33 xwx <- crossprod(x, (pre$WU %*% x)) 34 xwwx <- crossprod(x, (pre$W2 %*% x))
|
/dports/finance/R-cran-urca/urca/R/ |
H A D | ablrtest.R | 28 M00 <- crossprod(z@Z0)/N 29 M11 <- crossprod(z@Z1)/N 30 MKK <- crossprod(z@ZK)/N 31 M01 <- crossprod(z@Z0, z@Z1)/N 32 M0K <- crossprod(z@Z0, z@ZK)/N 33 MK0 <- crossprod(z@ZK, z@Z0)/N 34 M10 <- crossprod(z@Z1, z@Z0)/N 48 Saa.b <- crossprod(RA, RA)/N 49 Sak.b <- crossprod(RA, RK)/N 50 Ska.b <- crossprod(RK, RA)/N [all …]
|
H A D | alrtest.R | 19 M00 <- crossprod(z@Z0)/N 20 M11 <- crossprod(z@Z1)/N 21 MKK <- crossprod(z@ZK)/N 22 M01 <- crossprod(z@Z0, z@Z1)/N 23 M0K <- crossprod(z@Z0, z@ZK)/N 24 MK0 <- crossprod(z@ZK, z@Z0)/N 25 M10 <- crossprod(z@Z1, z@Z0)/N 39 Saa.b <- crossprod(RA, RA)/N 40 Sak.b <- crossprod(RA, RK)/N 41 Ska.b <- crossprod(RK, RA)/N [all …]
|
H A D | blrtest.R | 23 M00 <- crossprod(z@Z0)/N 24 M11 <- crossprod(z@Z1)/N 25 MKK <- crossprod(z@ZK)/N 26 M01 <- crossprod(z@Z0, z@Z1)/N 27 M0K <- crossprod(z@Z0, z@ZK)/N 28 MK0 <- crossprod(z@ZK, z@Z0)/N 29 M10 <- crossprod(z@Z1, z@Z0)/N 30 M1K <- crossprod(z@Z1, z@ZK)/N 31 MK1 <- crossprod(z@ZK, z@Z1)/N
|
H A D | bh5lrtest.R | 29 M00 <- crossprod(z@Z0)/N 30 M11 <- crossprod(z@Z1)/N 31 MKK <- crossprod(z@ZK)/N 32 M01 <- crossprod(z@Z0, z@Z1)/N 33 M0K <- crossprod(z@Z0, z@ZK)/N 34 MK0 <- crossprod(z@ZK, z@Z0)/N 35 M10 <- crossprod(z@Z1, z@Z0)/N 36 M1K <- crossprod(z@Z1, z@ZK)/N 37 MK1 <- crossprod(z@ZK, z@Z1)/N
|
/dports/math/R-cran-pls/pls/R/ |
H A D | simpls.fit.R | 101 S <- crossprod(X, Y) 110 q.a <- eigen(crossprod(S), symmetric = TRUE)$vectors[,1] 112 q.a <- c(crossprod(S, eigen(S %*% t(S), 114 q.a <- q.a / sqrt(c(crossprod(q.a))) 122 tnorm <- sqrt(c(crossprod(t.a))) 125 p.a <- crossprod(X, t.a) # X block factor loadings 126 q.a <- crossprod(Y, t.a) # Y block factor loadings 129 v.a <- v.a - V %*% crossprod(V, p.a) # v.a orth to previous loadings 131 v.a <- v.a / sqrt(c(crossprod(v.a))) # normalize orthogonal loadings 132 S <- S - v.a %*% crossprod(v.a, S) # deflate S [all …]
|
H A D | kernelpls.fit.R | 115 XtY <- crossprod(X, Y) 120 w.a <- XtY / sqrt(c(crossprod(XtY))) 124 q <- eigen(crossprod(XtY), symmetric = TRUE)$vectors[,1] 126 w.a <- w.a / sqrt(c(crossprod(w.a))) 136 r.a <- r.a - colSums(crossprod(w.a, P[,1:(a-1), drop=FALSE]) %*% 145 tsq <- c(crossprod(t.a)) 146 p.a <- crossprod(X, t.a) / tsq 147 q.a <- crossprod(XtY, r.a) / tsq 160 u.a <- Y %*% q.a / c(crossprod(q.a)) # Ok for nresp == 1 ?? 162 if (a > 1) u.a <- u.a - TT %*% (crossprod(TT, u.a) / tsqs)
|
/dports/science/R-cran-kernlab/kernlab/R/ |
H A D | ipop.R | 101 smwinner <- chol(V + crossprod(H)) 104 smwa2 <- smwa1 - (H %*% solve(smwinner,solve(t(smwinner),crossprod(H,smwa1)))) 105 smwc2 <- smwc1 - (H %*% solve(smwinner,solve(t(smwinner),crossprod(H,smwc1)))) 118 mu <- as.vector(crossprod(z,g) + crossprod(v,w) + crossprod(s,t) + crossprod(p,q))/(2 * (m + n)) 133 H.dot.x <- H %*% crossprod(H,x) 138 sigma <- c - crossprod(A, y) - z + s + H.dot.x 145 x.dot.H.dot.x <- crossprod(x, H.dot.x) 148 primal.obj <- crossprod(c,x) + 0.5 * x.dot.H.dot.x 149 …dual.obj <- crossprod(b,y) - 0.5 * x.dot.H.dot.x + crossprod(l, z) - crossprod(u,s) - crossprod(r,… 201 newmu <- (crossprod(z,g) + crossprod(v,w) + crossprod(s,t) + crossprod(p,q))/(2 * (m + n)) [all …]
|
/dports/math/R-cran-irlba/irlba/R/ |
H A D | svdr.R | 104 B <- crossprod(x, Q) 108 Q <- qr.Q(qr(x %*% Q - rep(1, nrow(x)) %*% crossprod(center, Q))) 109 B <- crossprod(Q, x) - tcrossprod(crossprod(Q, rep(1, nrow(x))), center) 122 B <- crossprod(Q, x) 125 Q <- qr.Q(qr(x %*% Q - rep(1, nrow(x)) %*% crossprod(center, Q))) 126 B <- crossprod(Q, x) - tcrossprod(crossprod(Q, rep(1, nrow(x))), center)
|
/dports/science/xtb/xtb-6.4.1/src/ |
H A D | basic_geo.f90 | 122 call crossprod(ra,rb,na) 123 call crossprod(rb,rc,nb) 138 call crossprod(na,rb,rab) 139 call crossprod(nb,ra,rba) 140 call crossprod(na,rc,rac) 141 call crossprod(nb,rb,rbb) 142 call crossprod(nb,rc,rbc) 143 call crossprod(na,ra,raa) 204 call crossprod(re,rd,rn) 299 call crossprod(ra,rb,na) [all …]
|
/dports/math/R-cran-gmp/gmp/man/ |
H A D | matrix.Rd | 18 % \alias{crossprod,bigq,bigq-method} 19 % \alias{crossprod,bigq,ANY-method} 20 % \alias{crossprod,ANY,bigq-method} 22 % \alias{crossprod,bigz,ANY-method} 23 % \alias{crossprod,ANY,bigz-method} 35 \alias{crossprod} 36 \alias{crossprod.default} 37 \alias{crossprod.bigq} 38 \alias{crossprod.bigz} 132 crossprod(V)# scalar product [all …]
|
/dports/math/R-cran-lava/lava/R/ |
H A D | Objective.R | 15 val <- log(detC) + tcrossprod(z,crossprod(z,iC))[1] 19 W <- suppressMessages(crossprod(rbind(mu-xi))) 69 as.numeric(crossprod(u,iC%*%W)%*%D$dxi + 76 as.numeric(score0 + crossprod(u, iC)%*%D$dxi + 97 W <- crossprod(rbind(mu-xi)) 113 Grad <- Grad - (n/2*crossprod(D$dT,vec.iC)) 179 as.numeric(crossprod(u,iC%*%W)%*%D$dxi + 181 - iC %*% crossprod(rbind(u)) 191 as.numeric(crossprod(u[i,],iC%*%W)%*%dxi + 203 res <- crossprod(S) [all …]
|
/dports/math/R-cran-quadprog/quadprog/tests/ |
H A D | test2.R | 12 print(crv1 <- (crossprod(res$unc, Dmat)/2-dvec)%*%res$unc) 13 print(crv2 <- (crossprod(res$solution, Dmat)/2-dvec)%*%res$solution) 16 print(all(crossprod(Amat, res$solution) >= bvec)) 29 print(crv1 <- (crossprod(res$unc, Dmat)/2-dvec)%*%res$unc) 30 print(crv2 <- (crossprod(res$solution, Dmat)/2-dvec)%*%res$solution) 33 print(all(crossprod(Amat, res$solution) >= bvec))
|
H A D | test1.R | 26 print(crv1 <- (crossprod(res$unc, Dmat)/2-dvec)%*%res$unc) 27 print(crv2 <- (crossprod(res$solution, Dmat)/2-dvec)%*%res$solution) 31 print(all(crossprod(Amat, res$solution) >= bvec)) 46 print(crv1 <- (crossprod(res$unc, Dmat)/2-dvec)%*%res$unc) 47 print(crv2 <- (crossprod(res$solution, Dmat)/2-dvec)%*%res$solution) 51 print(all(crossprod(Amat, res$solution) >= bvec))
|
/dports/finance/R-cran-gmm/gmm/R/ |
H A D | gmm.R | 320 xx <- crossprod(dg, whx) 331 par <- solve(crossprod(hm,xm),crossprod(hm,ym)) 370 hmxm <- crossprod(hm,xm) 375 par <- par-solve(crossprod(hmxm, T1), crossprod(T1, mustar)) 379 value <- crossprod(gb, solve(w, gb)) 381 value <- crossprod(gb, w%*%gb) 393 obj <- crossprod(gbar,w)%*%gbar 423 obj <- crossprod(gbar,solve(w,gbar)) 425 obj <- crossprod(gbar,w%*%gbar) 534 w <- crossprod(gt)/n [all …]
|
/dports/math/R/R-4.1.2/src/library/base/man/ |
H A D | crossprod.Rd | 1 % File src/library/base/man/crossprod.Rd 6 \name{crossprod} 7 \alias{crossprod} 13 faster than) the call \code{t(x) \%*\% y} (\code{crossprod}) or 17 crossprod(x, y = NULL) 56 (z <- crossprod(1:4)) # = sum(1 + 2^2 + 3^2 + 4^2) 61 crossprod(t(x))) 66 identical(tcrossprod(v, m), crossprod(v, t(m))), 67 identical(crossprod(m, v2), t(m) \%*\% v2))
|
/dports/math/libRmath/R-4.1.1/src/library/base/man/ |
H A D | crossprod.Rd | 1 % File src/library/base/man/crossprod.Rd 6 \name{crossprod} 7 \alias{crossprod} 13 faster than) the call \code{t(x) \%*\% y} (\code{crossprod}) or 17 crossprod(x, y = NULL) 56 (z <- crossprod(1:4)) # = sum(1 + 2^2 + 3^2 + 4^2) 61 crossprod(t(x))) 66 identical(tcrossprod(v, m), crossprod(v, t(m))), 67 identical(crossprod(m, v2), t(m) \%*\% v2))
|
/dports/math/R-cran-SQUAREM/SQUAREM/R/ |
H A D | squarem.R | 115 sr2 <- crossprod(q1) 123 sq2 <- sqrt(crossprod(q2)) 125 sv2 <- crossprod(q2-q1) 126 srv <- crossprod(q1, q2-q1) 231 sr2 <- crossprod(q1) 239 sq2 <- sqrt(crossprod(q2)) 242 sv2 <- crossprod(q2-q1) 243 srv <- crossprod(q1, q2-q1) 335 res <- sqrt(c(crossprod(ptmp - par))) 474 res <- sqrt(c(crossprod(ptmp - par))) [all …]
|