Running examples in 'tcgsaseq-Ex.R' failed
The error most likely occurred in:

> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: vc_score_h
> ### Title: Computes variance component test statistic for homogeneous
> ###   trajectory
> ### Aliases: vc_score_h
> ### Keywords: internal
> 
> ### ** Examples
> 
> rm(list=ls())
> set.seed(123)
> 
> ##generate some fake data
> ########################
> ng <- 100
> nindiv <- 30
> nt <- 5
> nsample <- nindiv*nt
> tim <- matrix(rep(1:nt), nindiv, ncol=1, nrow=nsample)
> tim <- cbind(tim, tim^2)
> sigma <- 5
> b0 <- 10
> 
> #under the null:
> beta1 <- rnorm(n=ng, 0, sd=0)
> #under the (heterogen) alternative:
> beta1 <- rnorm(n=ng, 0, sd=0.1)
> #under the (homogen) alternative:
> beta1 <- rnorm(n=ng, 0.06, sd=0)
> 
> y.tilde <- b0 + rnorm(ng, sd = sigma)
> y <- t(matrix(rep(y.tilde, nsample), ncol=ng, nrow=nsample, byrow=TRUE) +
+       matrix(rep(beta1, each=nsample), ncol=ng, nrow=nsample, byrow=FALSE)*matrix(rep(tim, ng),
+                                                             ncol=ng, nrow=nsample, byrow=FALSE) +
+       matrix(rnorm(ng*nsample, sd = sigma), ncol=ng, nrow=nsample, byrow=FALSE)
+       )
Error in matrix(rep(tim, ng), ncol = ng, nrow = nsample, byrow = FALSE) : 
  data length differs from size of matrix: [30000 != 150 x 100]
Calls: t -> matrix
Execution halted
Running examples in ‘tcgsaseq-Ex.R’ failed
The error most likely occurred in:

> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: vc_score_h
> ### Title: Computes variance component test statistic for homogeneous
> ###   trajectory
> ### Aliases: vc_score_h
> ### Keywords: internal
> 
> ### ** Examples
> 
> rm(list=ls())
> set.seed(123)
> 
> ##generate some fake data
> ########################
> ng <- 100
> nindiv <- 30
> nt <- 5
> nsample <- nindiv*nt
> tim <- matrix(rep(1:nt), nindiv, ncol=1, nrow=nsample)
> tim <- cbind(tim, tim^2)
> sigma <- 5
> b0 <- 10
> 
> #under the null:
> beta1 <- rnorm(n=ng, 0, sd=0)
> #under the (heterogen) alternative:
> beta1 <- rnorm(n=ng, 0, sd=0.1)
> #under the (homogen) alternative:
> beta1 <- rnorm(n=ng, 0.06, sd=0)
> 
> y.tilde <- b0 + rnorm(ng, sd = sigma)
> y <- t(matrix(rep(y.tilde, nsample), ncol=ng, nrow=nsample, byrow=TRUE) +
+       matrix(rep(beta1, each=nsample), ncol=ng, nrow=nsample, byrow=FALSE)*matrix(rep(tim, ng),
+                                                             ncol=ng, nrow=nsample, byrow=FALSE) +
+       matrix(rnorm(ng*nsample, sd = sigma), ncol=ng, nrow=nsample, byrow=FALSE)
+       )
Error in matrix(rep(tim, ng), ncol = ng, nrow = nsample, byrow = FALSE) : 
  data length differs from size of matrix: [30000 != 150 x 100]
Calls: t -> matrix
Execution halted
Running examples in ‘tcgsaseq-Ex.R’ failed
The error most likely occurred in:

> ### Name: vc_score_h
> ### Title: Computes variance component test statistic for homogeneous
> ###   trajectory
> ### Aliases: vc_score_h
> ### Keywords: internal
> 
> ### ** Examples
> 
> rm(list=ls())
> set.seed(123)
> 
> ##generate some fake data
> ########################
> ng <- 100
> nindiv <- 30
> nt <- 5
> nsample <- nindiv*nt
> tim <- matrix(rep(1:nt), nindiv, ncol=1, nrow=nsample)
> tim <- cbind(tim, tim^2)
> sigma <- 5
> b0 <- 10
> 
> #under the null:
> beta1 <- rnorm(n=ng, 0, sd=0)
> #under the (heterogen) alternative:
> beta1 <- rnorm(n=ng, 0, sd=0.1)
> #under the (homogen) alternative:
> beta1 <- rnorm(n=ng, 0.06, sd=0)
> 
> y.tilde <- b0 + rnorm(ng, sd = sigma)
> y <- t(matrix(rep(y.tilde, nsample), ncol=ng, nrow=nsample, byrow=TRUE) +
+       matrix(rep(beta1, each=nsample), ncol=ng, nrow=nsample, byrow=FALSE)*matrix(rep(tim, ng),
+                                                             ncol=ng, nrow=nsample, byrow=FALSE) +
+       matrix(rnorm(ng*nsample, sd = sigma), ncol=ng, nrow=nsample, byrow=FALSE)
+       )
Error in matrix(rep(tim, ng), ncol = ng, nrow = nsample, byrow = FALSE) : 
  data length differs from size of matrix: [30000 != 150 x 100]
Calls: t -> matrix
Execution halted
Running examples in ‘tcgsaseq-Ex.R’ failed
The error most likely occurred in:

> ### Name: vc_score_h
> ### Title: Computes variance component test statistic for homogeneous
> ###   trajectory
> ### Aliases: vc_score_h
> ### Keywords: internal
> 
> ### ** Examples
> 
> rm(list=ls())
> set.seed(123)
> 
> ##generate some fake data
> ########################
> ng <- 100
> nindiv <- 30
> nt <- 5
> nsample <- nindiv*nt
> tim <- matrix(rep(1:nt), nindiv, ncol=1, nrow=nsample)
> tim <- cbind(tim, tim^2)
> sigma <- 5
> b0 <- 10
> 
> #under the null:
> beta1 <- rnorm(n=ng, 0, sd=0)
> #under the (heterogen) alternative:
> beta1 <- rnorm(n=ng, 0, sd=0.1)
> #under the (homogen) alternative:
> beta1 <- rnorm(n=ng, 0.06, sd=0)
> 
> y.tilde <- b0 + rnorm(ng, sd = sigma)
> y <- t(matrix(rep(y.tilde, nsample), ncol=ng, nrow=nsample, byrow=TRUE) +
+       matrix(rep(beta1, each=nsample), ncol=ng, nrow=nsample, byrow=FALSE)*matrix(rep(tim, ng),
+                                                             ncol=ng, nrow=nsample, byrow=FALSE) +
+       matrix(rnorm(ng*nsample, sd = sigma), ncol=ng, nrow=nsample, byrow=FALSE)
+       )
Error in matrix(rep(tim, ng), ncol = ng, nrow = nsample, byrow = FALSE) : 
  data length differs from size of matrix: [30000 != 150 x 100]
Calls: t -> matrix
Execution halted