Description Usage Arguments Value Author(s) Examples
Compute the predicted response matrix for a new data set.
1 2 |
object |
an object of class lsgl, produced with |
x |
a data matrix of size N_\textrm{new} \times p. |
sparse.data |
if TRUE |
... |
ignored. |
Yhat |
the predicted response matrix (of size N_\textrm{new} \times K) |
Martin Vincent
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 | set.seed(100) # This may be removed, it ensures consistency of the daily tests
## Simulate from Y=XB+E, the dimension of Y is N x K, X is N x p, B is p x K
N <- 100 #number of samples
p <- 50 #number of features
K <- 25 #number of groups
B<-matrix(sample(c(rep(1,p*K*0.1),rep(0, p*K-as.integer(p*K*0.1)))),nrow=p,ncol=K)
X1<-matrix(rnorm(N*p,1,1),nrow=N,ncol=p)
Y1 <-X1%*%B+matrix(rnorm(N*K,0,1),N,K)
X2<-matrix(rnorm(N*p,1,1),nrow=N,ncol=p)
Y2 <-X2%*%B+matrix(rnorm(N*K,0,1),N,K)
#### Fit models using X1
lambda <- lsgl::lambda(X1, Y1, alpha = 1, d = 25L, lambda.min = 0.5, intercept = FALSE)
fit <- lsgl::fit(X1, Y1, alpha = 1, lambda = lambda, intercept = FALSE)
# Predict Y2 using the estimated models and X2
res <- predict(fit, X2)
# Frobenius norm \|Y2hat - Y2\|_F
sapply(res$Yhat, function(y) sqrt(sum((y - Y2)^2)))
# This is proportional to the errors compute with Err:
Err(fit, X2, Y2)*length(Y2)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.